
[Sponsors] 
March 17, 2013, 09:17 
is double precision default of OF?

#1 
Senior Member
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,209
Rep Power: 18 
can anyone answer me?
Thanks. 

March 17, 2013, 18:23 

#2 
Senior Member
Håkon Strandenes
Join Date: Dec 2011
Location: Norway
Posts: 111
Rep Power: 11 
Yes, by default it is double precision. But you can change this at compile time by setting WM_PRECISION_OPTION=SP instead of WM_PRECISION_OPTION=DP in etc/bashrc before you compile. Anyways, that is not something I would recommend.


March 18, 2013, 02:13 

#3 
Senior Member
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,907
Rep Power: 27 
One comment: the calculation is done in double precision by default, but results are saved to disk typically with a precision of 6 digits. Use the binary format to preserve precision.
Best,
__________________
Alberto Passalacqua GeekoCFD  A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as in both physical and virtual formats. OpenQBMM  An opensource implementation of quadraturebased moment methods. To obtain more accurate answers, please specify the version of OpenFOAM you are using. 

March 18, 2013, 05:44 

#4 
Senior Member
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,209
Rep Power: 18 
thank you dear Hacon and Alberto.
double precision is 14 digit after point,isn't? When tolerances are set for example e20 does it have any effect on accuracy of calculations since it's lower than DP criteria? Thanks again. 

March 18, 2013, 18:56 

#5  
Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 9,556
Blog Entries: 39
Rep Power: 97 
Greetings to all!
Quote:
In essence:
If I'm not mistaken, this is where the continuity indicators help figure out how good or bad the solution is... but that's all I can remember Best regards, Bruno
__________________


March 21, 2013, 13:30 

#6 
Senior Member
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,209
Rep Power: 18 
thanks but I wonder then whats its relation to zero of machine thats about 1e17 that must reached by continuity residuals?Is it a limit while DP can work with numbers very smaller than 1e17?


March 21, 2013, 13:36 

#7 
Senior Member
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,209
Rep Power: 18 
I was curious about DP but didn't find it in etc/bashrc file.


March 24, 2013, 13:22 

#8  
Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 9,556
Blog Entries: 39
Rep Power: 97 
Quick answers:
Quote:
Quote:
Quote:
edit: "eps" is actually the name of the command on MATLAB that gives us this machine precision... it looks like it's not an acronym but it actually stands for "epsilon"
__________________
Last edited by wyldckat; March 24, 2013 at 13:49. Reason: see "edit:" 

March 24, 2013, 13:49 

#9  
Senior Member
Håkon Strandenes
Join Date: Dec 2011
Location: Norway
Posts: 111
Rep Power: 11 
Quote:
See this page for more information: http://en.wikipedia.org/wiki/Double...gpoint_format 

March 29, 2013, 02:04 

#10 
Senior Member
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,209
Rep Power: 18 
thank you both so much for explanations.
But still a doubt is remaind about the relation between machine eps(around 1e16) and DP(that calculates until e308 in low values) or SP. 

March 29, 2013, 17:39 

#11 
Senior Member
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,209
Rep Power: 18 
hakon and bruno,please complete your useful explanations with giving a comparison between DP bound and machine epsilon.
Why machine epsilon is restricted to e17 when calculations could be done untile e308 according to DP accuracy? 

March 29, 2013, 19:58 

#12  
Super Moderator
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 9,556
Blog Entries: 39
Rep Power: 97 
Quote:
Anyway, to try to answer your question, I'll go through a somewhat of summary bullet point presentation:
__________________


March 30, 2013, 05:20 

#13 
Senior Member
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,209
Rep Power: 18 
thank you Bruno
I grasped it. There isn't any comparison description in wiki.no,i don't remember such topics in the courses.professor assumed such subject too simple that he want to describe.so my friends i asked wasn't certain like me.although we had courses like numerical analysis in bsc but they're too far to remember! Have a nice day and good asleep. 

March 30, 2013, 08:49 

#14 
Senior Member
Håkon Strandenes
Join Date: Dec 2011
Location: Norway
Posts: 111
Rep Power: 11 
I started on writing a short explanation yesterday, so I thought I might post it even tough this thread should be laid dead:
The magic lies in the definition of epsilon, which according to the Wikipedia article "gives an upper bound on the relative error due to rounding in floating point arithmetic" (my emphasis). To use my previous example, if you add two numbers of different magnitude, you might get roundoff errors. This also goes for other operations, such as multiplication/division. For the sake of simplicity assume that one of the numbers is one (exactly). First, let the other number also be one (exactly). In double precision that is: Code:
1+1=2.00000000000000000000 Then let us assume that we are to do a division, 1/10, witch we know (from numerics) is exactly 0.1. In double precision that is actually not the case: Code:
1/10 = 0.10000000000000000555 We might as well do this with some other numbers: Code:
1e67/1e66 = 0.09999999999999999167 Lets do some additions: Code:
1e100+2e100 = 3.00000000000000005998e100 Do you get the concept? The crux is that we are dealing with relative errors due to roundoffs. The last example is an good one, because we are dealing with numbers that are small (order of magnitude of 10^100), and thus get small absolute errors. The key in the definition of epsilon is that it (i repeat myself) "gives an upper bound on the relative error due to rounding in floating point arithmetic". In all of our examples the relative errors are smaller than epsilon, and that is what it is all about. Disclaimer: I have not considered that numbers like 1e100 are not possible to represent exactly in double precision, and hence the error calculated here is a result of first converting 1e100 (in the exact sense) to a binary number (witch is 1.00000000000000001999e100) and then doing the same thing with the other numbers before actually performing the operation. Nevertheless, it illustrate the concept fairly well without adding to much complexity. I'm just a silly hydrodynamics engineer, so if some computer geeks wants to blow my head off and come with some exact calculations and considerations, feel free to do that. 

March 30, 2013, 09:20 

#15 
Senior Member
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,209
Rep Power: 18 
thank you so much dear hakon for precise writing.that illuminated the concept more.this topic is alive for refers and add thoughts.I hope it is useful to anyone who has some questions related machine epsilon and double precision.


Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Continuing User Defined Real Gas Model issues  aeroman  FLUENT  6  April 8, 2016 03:34 
Parallel User Defined Real Gas Model  aeroman  FLUENT  4  July 1, 2015 06:09 
Missing math.h header  Travis  FLUENT  4  January 15, 2009 12:48 
what's wrong about my code for 2d burgers equation  morxio  Main CFD Forum  3  April 27, 2007 10:38 
REAL GAS UDF  brian  FLUENT  6  September 11, 2006 08:23 