Changeset 2078 for branches/newmole/source/ion_solver.cpp

Timestamp:

05/17/08 07:00:31 (6 months ago)

Author:

rjrw

Message:

Deal with FPE in grains_lte model

Files:

• branches/newmole/source/ion_solver.cpp (modified) (1 diff)

branches/newmole/source/ion_solver.cpp

@@ -417 +417 @@
                 d = SDIV(d);
                 double ratio1 = a/b , ratio2 = c/d , fratio1=fabs(a/b),fratio2=fabs(c/d);
-                if( fabs(ratio1-ratio2)/MAX2(fratio1,fratio2) <DBL_EPSILON*1e4 )
+                if( fabs(ratio1-ratio2) <= DBL_EPSILON*1e4*MAX2(fratio1,fratio2) )
                 {
                         lgHomogeneous = true;