Changeset 2027 for branches/newmole/source/mole_reactions.cpp

Timestamp:

05/06/08 14:47:46 (8 months ago)

Author:

rjrw

Message:

Add optional diagnostics for rate and chemical balance controlling
species update.

Remove logically dead code.

Files:

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

branches/newmole/source/mole_reactions.cpp

@@ -1893 +1893 @@
 {
         struct mole_reaction_s *rate;
+        enum { DEBUG_MOLE = false };
 
         DEBUG_ENTRY("mole_update_rks()");
@@ -1906 +1907 @@
         {
                 rate = p->second;
+                realnum oldrk = rate->rk;
                 ASSERT(rate->fun != NULL);
                 rate->rk = rate->a*rate->fun(rate);
+                if (DEBUG_MOLE)
+                {
+                        if (fabs(rate->rk-oldrk) > 0.1*rate->rk)
+                                fprintf(ioQQQ,"%s: %15.8g => %15.8g\n",
+                                                                rate->label,oldrk,rate->rk);
+                }
         }       
 }