Changeset 2156

Timestamp:

07/04/08 16:18:11 (5 months ago)

Author:

peter

Message:

source/grains_qheat.cpp:

Bug-fix: fix crash in grains_qheat.in due to very steep gradient in Phi.

Files:

• trunk/source/grains_qheat.cpp (modified) (2 diffs)

trunk/source/grains_qheat.cpp

@@ -1760 +1760 @@
                 p[k] = (sum + trap1*delu[k])/(2.*Lambda[k] - Phi[0]*delu[k]);
 
-                ASSERT( p[k] > 0. );
+                // total failure -> force next step to be smaller
+                if( p[k] <= 0. )
+                        return 3.*QHEAT_TOL;
         }
 
@@ -1766 +1768 @@
         p2k = (sum2 + trap12*step)/(2.*Lambda[k] - Phi[0]*step);
 
-        ASSERT( p2k > 0. );
+        // total failure -> force next step to be smaller
+        if( p2k <= 0. )
+                return 3.*QHEAT_TOL;
 
         RelErrPk = fabs(p2k-p[k])/p[k];