Changeset 684 for branches/newmole/source/mole_co_solve.cpp

Timestamp:

12/05/06 13:55:14 (2 years ago)

Author:

rjrw

Message:

Patch back to deal with SEGVs on previous step.

Still need to a) reorganize startup b) improve treatment of
decoupling within the network.

Files:

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

branches/newmole/source/mole_co_solve.cpp

@@ -50 +50 @@
         float abundan;
         struct chem_element_s *element;
+        double sum;
 
         CO_step();                  /* Calculate the matrix elements */
+
+        /* Ugly hack to deal with species which have become uncoupled */
+        for (i=0;i<mole.num_comole_calc;i++)
+        {
+                sum = 0.;
+                for (j=0;j<mole.num_comole_calc;j++)
+                {
+                        sum = sum+fabs(mole.c[i][j]);
+                }
+                if (sum < SMALLFLOAT && fabs(mole.b[i]) < SMALLFLOAT)
+                {
+                        mole.c[i][i] = 1.;
+                }
+        }
 
         cartot_mol = dense.xMolecules[ipCARBON] + findspecies("C")->hevmol + findspecies("C+")->hevmol;