Changeset 2118 for branches/newmole/source/mole_solve.cpp

Timestamp:

06/22/08 10:51:55 (5 months ago)

Author:

rjrw

Message:

Use internal timestep limiting (...in an attempt...) to improve
robustness of molecular solver.

Files:

• branches/newmole/source/mole_solve.cpp (modified) (2 diffs)

branches/newmole/source/mole_solve.cpp

@@ -149 +149 @@
                 MoleMap.updateMolecules( newmols );
                 
-                // fprintf(stdout,"Mole zone %ld -- %7.2f error %15.8g nBad %d\n",nzone,fnzone,error,nBad);
+                // fprintf(stdout,"Mole zone %ld -- %7.2f error %15.8g nBad %d loop %d\n",nzone,fnzone,error,nBad,i);
         
                 if (error < BIGERROR && nBad == 0 && nPrevBad == 0) /* Stop early if good enough */
@@ -361 +361 @@
         else
         {
-                lgConserve = true; 
+                // Without pseudo-timestepping, need to apply conservation
+                // constraint to maintain non-singular matrix
+                lgConserve = false; 
         }