diff options
author | lanfer <lanfer@1d96b42b-98df-4a6a-9d84-1b24288d4588> | 1999-09-15 09:25:47 +0000 |
---|---|---|
committer | lanfer <lanfer@1d96b42b-98df-4a6a-9d84-1b24288d4588> | 1999-09-15 09:25:47 +0000 |
commit | f1b684b628b47565e995c9fc05d7e4d7fb2d38c1 (patch) | |
tree | d0657c33d768a12210da4a1ff9de2928f5bd9c66 | |
parent | f931c69d368ae51b5d3b8f483cd73f5a4231ad6a (diff) |
misc.
git-svn-id: http://svn.cactuscode.org/arrangements/CactusElliptic/EllPETSc/trunk@9 1d96b42b-98df-4a6a-9d84-1b24288d4588
-rw-r--r-- | src/petsc_confmetric_solver.c | 15 |
1 files changed, 10 insertions, 5 deletions
diff --git a/src/petsc_confmetric_solver.c b/src/petsc_confmetric_solver.c index 49fe8a2..3cf0476 100644 --- a/src/petsc_confmetric_solver.c +++ b/src/petsc_confmetric_solver.c @@ -95,7 +95,7 @@ void petsc_confmetric_solver(cGH *GH, int *MetricPsiI, int MetricPsiISize, int its; /* Number of iterations */ /* Loopers */ - int i,j,k,l,m,n; + int i,j,k,l,m,n,ll; /* loop limits put in for stencil_w !=1 Ed Evans 1/19/98 */ int imin,imax,jmin,jmax,kmin,kmax,interior; @@ -169,9 +169,6 @@ void petsc_confmetric_solver(cGH *GH, int *MetricPsiI, int MetricPsiISize, PetscTolStyle=0; tolerance =AbsTol[0]; - printf("\n**** PETSC !!!!****\n\n"); - - /* Get the link to pugh Extension */ pughGH = (pGH*)GH->extensions[CCTK_GHExtensionHandle("PUGH")]; if (!pughGH) CCTK_WARN(0,"ETERNAL ERROR: Cannot find PUGH Extension Handle\n"); @@ -733,7 +730,7 @@ void petsc_confmetric_solver(cGH *GH, int *MetricPsiI, int MetricPsiISize, /*$asb = uyz(i,j,k-1)+uyz(i,j-1,k)$*/ a(0,-1,-1) = Uyz(i,j,k-1)+Uyz(i,j-1,k); } /* end nabla_form=3 */ - + /* M phi */ if (Mstorage) a(0,0,0) = a(0,0,0) - Mlin[ijk]; @@ -767,6 +764,14 @@ void petsc_confmetric_solver(cGH *GH, int *MetricPsiI, int MetricPsiISize, direction, but not impossible. */ +#ifdef DEBUG + if (ijk==1918) { + for (l=-1;l<=1;l++) + for (m=-1;m<=1;m++) + for (n=-1;n<=1;n++) printf(" (%d %d %d): %f \n",l,m,n,a(l,m,n)); + printf("\n"); + } +#endif if (octant) for (l=-1;l<=0;l++) |