diff options
author | knarf <knarf@b2a53a04-0f4f-0410-87ed-f9f25ced00cf> | 2004-11-17 08:54:32 +0000 |
---|---|---|
committer | knarf <knarf@b2a53a04-0f4f-0410-87ed-f9f25ced00cf> | 2004-11-17 08:54:32 +0000 |
commit | 05f3af0337f3a4abfeb32f004ef3d9be5a7bf2dd (patch) | |
tree | 82ca84952a9ef78bcf95c2c54d26a990a40f17e4 /src/Newton.c | |
parent | 9624179def3e1e7b1d9a771b0d236d8bf9e4a53f (diff) |
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
Diffstat (limited to 'src/Newton.c')
-rw-r--r-- | src/Newton.c | 11 |
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); |