Changeset 2055 for branches/newmole/source/mole_eval_balance.cpp

Timestamp:

05/13/08 17:01:00 (8 months ago)

Author:

rjrw

Message:

Move further towards Householder projection method for of conservation
constraints

Some other tidyings (move from explicit pointers to multi_arr)

Files:

• branches/newmole/source/mole_eval_balance.cpp (modified) (5 diffs)

branches/newmole/source/mole_eval_balance.cpp

@@ -29 +29 @@
 /*=================================================================*/
 
-void mole_eval_balance(long int num_total, double *b, double **c)
+void mole_eval_balance(long int num_total, double *b, bool lgJac, multi_arr<double,2> &c)
 {
         long int i, j, nelem, ion, ion2;
@@ -44 +44 @@
         {
                 mole.zone[i].src = mole.zone[i].snk = 0.;
-                if (c)
+                if (lgJac)
                         for( j=0; j < num_total; j++ )
                         {
@@ -129 +129 @@
                 }
                 
-                if (c)
+                if (lgJac)
                 {
                         for(j=0;j<rate->nreactants;j++)
@@ -149 +149 @@
         }
 #if !defined(NDEBUG)
-        if (c)
+        if (lgJac)
         {
                 for (i=0;i<mole.num_calc;i++) 
@@ -206 +206 @@
                         fprintf(stdout,"]\n");
                 }
-                if (c)
+                if (lgJac)
                         fprintf(stdout,"Diag %g err %g rat %g\n",c[debug_species->index][debug_species->index],
                                                         b[debug_species->index],b[debug_species->index]/SDIV(-c[debug_species->index][debug_species->index]));