Changeset 1693 for branches/newmole/source

Timestamp:

12/18/07 02:38:17 (11 months ago)

Author:

rjrw

Message:

Previous fix for

limit_supra

breaks

grains_lte, h2_cr_grains, h2_solomon, h2_t2000, h2_t500, helike_n,
pdr_leiden_hack_v3

-- back out the fix for the moment, will try something else...

Files:

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

branches/newmole/source/mole_eval_balance.cpp

@@ -29 +29 @@
 /*=================================================================*/
 
+#define FABSSWITCH(a) (a)
 
 void mole_eval_balance(long int num_total, double *b, double **c)
@@ -90 +91 @@
                                 if(i!=j)
                                 {
-                                        rate_deriv[i] *= fabs(rate->reactants[j]->hevmol);
-                                }
-                        }
-                }
-                
-                rate_tot = rate_deriv[0]*fabs(rate->reactants[0]->hevmol);
+                                        rate_deriv[i] *= FABSSWITCH(rate->reactants[j]->hevmol);
+                                }
+                        }
+                }
+                
+                rate_tot = rate_deriv[0]*FABSSWITCH(rate->reactants[0]->hevmol);
                 
                 for(i=0;i<rate->nreactants;i++)