diff options
Diffstat (limited to 'src/petsc_flat_solver.c')
-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++) { |