[cig-commits] r11905 - mc/3D/CitcomCU/trunk/src

[tan2 at geodynamics.org]

Author: tan2
Date: 2008-05-02 16:41:00 -0700 (Fri, 02 May 2008)
New Revision: 11905

Modified:
   mc/3D/CitcomCU/trunk/src/General_matrix_functions.c
Log:
Fix for conjugate gradient solver, contributed by Magali Billen.


Modified: mc/3D/CitcomCU/trunk/src/General_matrix_functions.c
===================================================================
--- mc/3D/CitcomCU/trunk/src/General_matrix_functions.c	2008-05-02 23:39:49 UTC (rev 11904)
+++ mc/3D/CitcomCU/trunk/src/General_matrix_functions.c	2008-05-02 23:41:00 UTC (rev 11905)
@@ -413,6 +413,7 @@
 
 	if(!(E->control.NMULTIGRID || E->control.EMULTIGRID))
 	{
+		valid = (residual < acc) ? 0 : 1;
 		cycles = E->control.v_steps_low;
 		time = CPU_time0();
 		residual = conj_grad(E, D1, r, Au, acc, &cycles, high_lev);