diff options
| author | lanfer <lanfer@1d96b42b-98df-4a6a-9d84-1b24288d4588> | 1999-09-15 14:12:09 +0000 |
|---|---|---|
| committer | lanfer <lanfer@1d96b42b-98df-4a6a-9d84-1b24288d4588> | 1999-09-15 14:12:09 +0000 |
| commit | 24f2fc0261923ba90b81944790a4ccd0637e4735 (patch) | |
| tree | 95788841ad168c9741bb42638f56f29635354ad3 | |
| parent | 20a12762a605dd0e6284510d7c30f92034cb5bae (diff) | |
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
| -rw-r--r-- | src/petsc_flat_solver.c | 14 |
1 files 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<kmax;k++) { |