implement initial guess, matter rescaling, debug output files

git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/TwoPunctures/trunk@31 b2a53a04-0f4f-0410-87ed-f9f25ced00cf

1 files changed, 8 insertions, 3 deletions

diff --git a/src/Newton.c b/src/Newton.c
index 2eefea1..076b926 100644
--- a/src/Newton.c
+++ b/src/Newton.c
@@ -421,9 +421,6 @@ Newton (CCTK_POINTER_TO_CONST cctkGH,
   allocate_derivs (&u, ntotal);
 
 /*         TestRelax(nvar, n1, n2, n3, v, dv.d0); */
-  for (j = 0; j < ntotal; j++)
-    v.d0[j] = 0;
-
   it = 0;
   dmax = 1;
   while (dmax > tol && it < itmax)
@@ -456,6 +453,14 @@ Newton (CCTK_POINTER_TO_CONST cctkGH,
     }
     it += 1;
   }
+  if (itmax==0)
+  {
+      F_of_v (cctkGH, nvar, n1, n2, n3, v, F, u);
+      dmax = -1;
+      for (j = 0; j < ntotal; j++)
+        if (fabs (F[j]) > dmax)
+          dmax = fabs (F[j]);
+  }
   printf ("Newton: it=%d \t |F|=%e \n", it, dmax);
   fflush(stdout);