From f1b684b628b47565e995c9fc05d7e4d7fb2d38c1 Mon Sep 17 00:00:00 2001
From: lanfer <lanfer@1d96b42b-98df-4a6a-9d84-1b24288d4588>
Date: Wed, 15 Sep 1999 09:25:47 +0000
Subject: misc.

git-svn-id: http://svn.cactuscode.org/arrangements/CactusElliptic/EllPETSc/trunk@9 1d96b42b-98df-4a6a-9d84-1b24288d4588
---
 src/petsc_confmetric_solver.c | 15 ++++++++++-----
 1 file 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++)
-- 
cgit v1.2.3