Changeset 1851 for branches/newmole/source/dense.cpp

Timestamp:

03/16/08 03:11:06 (10 months ago)

Author:

rjrw

Message:

Deal with FPE in dense.cpp

Use iteratively improved matrix solver in ion_solver

Change criteria for effectively homogeneous matrix in ion_solver

Files:

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

branches/newmole/source/dense.cpp

@@ -56 +56 @@
                                                                 sum,
                                                                 dense.gas_phase[nelem],
-                                                                (sum-dense.gas_phase[nelem])/sum);
+                                                                (sum-dense.gas_phase[nelem])/dense.gas_phase[nelem]);
                                 fprintf(ioQQQ,"Conservation should be better than this\n");
                                 lgOK=false;