[wildfire230]

Hi All,


I am having some interesting issues where I believe Paraview is losing precision with some of my data. A simple way to reproduce my case is to use funkySetFields to initialize two data sets as follows:

setC1
{
field C1; // Field to initialize
expression "exp(pos().x) + 1.0";
keepPatches true;
variables ();
}
setC2
{
field C2; // Field to initialize
expression "exp(pos().x)";
keepPatches true;
variables ();
}


You can see that I am initializing C1=exp(x)+1, and C2=exp(x).


My domain goes from x = -25 to +25, so that the exponential does become quite large, like 10^10. I am using a write precision of 16. Then, in Paraview, I take C1 - C2 using a calculator filter, and the result contains errors on the order of O(10), whereas if we really have 16 digits of precision, I believe the errors should be O(1e-5) in this case. Am I missing something? Does Paraview's calculator only do single-point precision??


EDIT: Ok, I have narrowed down the issue a bit. It seems that whenever you run "paraFoam", all the fields are loaded as type float. Is there any way to load as type double?? Alternately, is there a way to run foamToVTK and write as type double? It seems like this also defaults to float.




I would greatly appreciate any thoughts.


Best regards.