diff options
author | jese <jese@b2a53a04-0f4f-0410-87ed-f9f25ced00cf> | 2008-05-07 15:10:05 +0000 |
---|---|---|
committer | jese <jese@b2a53a04-0f4f-0410-87ed-f9f25ced00cf> | 2008-05-07 15:10:05 +0000 |
commit | b8badf6c13e0cc80e21e34f3c076c7d726a995b3 (patch) | |
tree | adecf4fc512c7da917da403ec54c946ddf008711 /src/Newton.c | |
parent | ecba7a774c9b7c77528199d1418b5a8157389862 (diff) |
This is much slower, but much more accurate and does not give diffrent values after the last post-mass-solving call to the Newton solver.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/TwoPunctures/trunk@86 b2a53a04-0f4f-0410-87ed-f9f25ced00cf
Diffstat (limited to 'src/Newton.c')
-rw-r--r-- | src/Newton.c | 2 |
1 files changed, 2 insertions, 0 deletions
diff --git a/src/Newton.c b/src/Newton.c index 058a443..643404e 100644 --- a/src/Newton.c +++ b/src/Newton.c @@ -412,6 +412,7 @@ Newton (CCTK_POINTER_TO_CONST cctkGH, derivs v, CCTK_REAL tol, int itmax) { DECLARE_CCTK_PARAMETERS; + int ntotal = n1 * n2 * n3 * nvar, ii, j, it; CCTK_REAL *F, dmax, normres; derivs u, dv; @@ -436,6 +437,7 @@ Newton (CCTK_POINTER_TO_CONST cctkGH, for (j = 0; j < ntotal; j++) dv.d0[j] = 0; printf ("Newton: it=%d \t |F|=%e\n", it, (double)dmax); + printf ("bare mass: mp=%g \t mm=%g\n", (double) par_m_plus, (double) par_m_minus); fflush(stdout); ii = bicgstab (cctkGH, |