From 24f2fc0261923ba90b81944790a4ccd0637e4735 Mon Sep 17 00:00:00 2001 From: lanfer Date: Wed, 15 Sep 1999 14:12:09 +0000 Subject: a flat petsc solver, needs to be tested git-svn-id: http://svn.cactuscode.org/arrangements/CactusElliptic/EllPETSc/trunk@13 1d96b42b-98df-4a6a-9d84-1b24288d4588 --- src/petsc_flat_solver.c | 14 +++++++------- 1 file changed, 7 insertions(+), 7 deletions(-) diff --git a/src/petsc_flat_solver.c b/src/petsc_flat_solver.c index 1bbefc0..ccd1399 100644 --- a/src/petsc_flat_solver.c +++ b/src/petsc_flat_solver.c @@ -306,13 +306,13 @@ void petsc_flat(cGH *GH, int *FieldIndex, int *MIndex, int *NIndex, /* Initialize the stencil array */ for (I=0;I<27;I++) a[I] = 0.0; - a( 0, 0, 0) = -10.38; - a( 1, 0, 0) = 1.56; /* ae */ - a(-1, 0, 0) = 1.56; /* aw */ - a( 0, 1, 0) = 1.56; /* an */ - a( 0,-1, 0) = 1.56; /* as */ - a( 0, 0, 1) = 1.56; /* at */ - a( 0, 0,-1) = 1.56; /* ab */ + a( 0, 0, 0) =-6.0; + a( 1, 0, 0) = 1.0; /* ae */ + a(-1, 0, 0) = 1.0; /* aw */ + a( 0, 1, 0) = 1.0; /* an */ + a( 0,-1, 0) = 1.0; /* as */ + a( 0, 0, 1) = 1.0; /* at */ + a( 0, 0,-1) = 1.0; /* ab */ for (k=kmin;k