[cig-commits] r6354 - mc/3D/CitcomS/trunk/lib

tan2 at geodynamics.org tan2 at geodynamics.org
Thu Mar 22 12:58:11 PDT 2007


Author: tan2
Date: 2007-03-22 12:58:11 -0700 (Thu, 22 Mar 2007)
New Revision: 6354

Modified:
   mc/3D/CitcomS/trunk/lib/Stokes_flow_Incomp.c
Log:
Make the output of 1st iteration consistent with later iterations

Modified: mc/3D/CitcomS/trunk/lib/Stokes_flow_Incomp.c
===================================================================
--- mc/3D/CitcomS/trunk/lib/Stokes_flow_Incomp.c	2007-03-22 19:57:57 UTC (rev 6353)
+++ mc/3D/CitcomS/trunk/lib/Stokes_flow_Incomp.c	2007-03-22 19:58:11 UTC (rev 6354)
@@ -234,9 +234,15 @@
   convergent=0;
 
   if (E->control.print_convergence && E->parallel.me==0)  {
-    fprintf(E->fp,"AhatP (%03d) after %g seconds with div/v=%.3e for step %d\n",count,CPU_time0()-time0,E->monitor.incompressibility,E->monitor.solution_cycles); /**/
+    fprintf(E->fp,"AhatP (%03d) after %g seconds with div/v=%.3e dv/v=%.3e"
+            " and dp/p=%.3e for step %d\n",
+            count, CPU_time0()-time0, E->monitor.incompressibility,
+            0.0, 0.0, E->monitor.solution_cycles);
     fflush(E->fp);
-    fprintf(stderr,"AhatP (%03d) after %g seconds with div/v=%.3e for step %d\n",count,CPU_time0()-time0,E->monitor.incompressibility,E->monitor.solution_cycles); /**/
+    fprintf(stderr,"AhatP (%03d) after %g seconds with div/v=%.3e dv/v=%.3e"
+            " and dp/p=%.3e for step %d\n",
+            count, CPU_time0()-time0, E->monitor.incompressibility,
+            0.0, 0.0, E->monitor.solution_cycles);
   }
 
   dpressure = 1.0;



More information about the cig-commits mailing list