[cig-commits] commit 1944 by heister to /var/svn/dealii/aspect

[dealii.demon at gmail.com]

Revision 1944

do not use == 0.0 for doubles but a tolerance instead

U   trunk/aspect/source/simulator/assembly.cc


http://www.dealii.org/websvn/revision.php?repname=Aspect+Repository&path=%2F&rev=1944&peg=1944

Diff:
Modified: trunk/aspect/source/simulator/assembly.cc
===================================================================
--- trunk/aspect/source/simulator/assembly.cc	2013-10-08 15:13:56 UTC (rev 1943)
+++ trunk/aspect/source/simulator/assembly.cc	2013-10-08 15:18:50 UTC (rev 1944)
@@ -595,7 +595,9 @@
                      const double                        cell_diameter,
                      const TemperatureOrComposition     &temperature_or_composition) const
   {
-    if (global_u_infty == 0 || global_entropy_variation == 0 || global_field_variation == 0)
+    if (std::abs(global_u_infty) < 1e-50
+	|| std::abs(global_entropy_variation) < 1e-50
+	|| std::abs(global_field_variation) < 1e-50)
       return 5e-3 * cell_diameter;
 
     double max_residual = 0;