diff options
Diffstat (limited to 'src/sor_confmetric.c')
-rw-r--r-- | src/sor_confmetric.c | 11 |
1 files changed, 7 insertions, 4 deletions
diff --git a/src/sor_confmetric.c b/src/sor_confmetric.c index b258c91..4f5a2e4 100644 --- a/src/sor_confmetric.c +++ b/src/sor_confmetric.c @@ -20,9 +20,12 @@ #include "CactusBase/CartGrid3D/src/Symmetry.h" #include "CactusElliptic/EllBase/src/Ell_DBstructure.h" - #define SQR(a) ((a)*(a)) +void sor_confmetric_3d(cGH *GH, int *MetricPsiI, int conformal, + int FieldIndex, int MIndex, int NIndex, + CCTK_REAL *AbsTol, CCTK_REAL *RelTol); + /*@@ @routine sor_confmetric @date Tue Sep 26 11:28:08 2000 @@ -95,7 +98,7 @@ void sor_confmetric_3d(cGH *GH, int *MetricPsiI, int conformal, CCTK_REAL ac; CCTK_REAL ae,aw,an,as,at,ab; CCTK_REAL ane, anw, ase, asw, ate, atw, abe, abw; - CCTK_REAL atn, ats, abn, abs; + CCTK_REAL atn, ats, abn, absol; /* Miscellaneous */ int sum_handle=-1; @@ -338,7 +341,7 @@ void sor_confmetric_3d(cGH *GH, int *MetricPsiI, int conformal, atn = uyz[ijpk]+uyz[ijkp]; ats =-uyz[ijkp]-uyz[ijmk]; abn =-uyz[ijkm]-uyz[ijpk]; - abs = uyz[ijkm]+uyz[ijmk]; + absol = uyz[ijkm]+uyz[ijmk]; residual = ac * var[ijk] + ae *var[ipjk] + aw*var[imjk] @@ -355,7 +358,7 @@ void sor_confmetric_3d(cGH *GH, int *MetricPsiI, int conformal, + ate*var[ipjkp] + atw*var[imjkp] + abe*var[ipjkm] + abw*var[imjkm] + atn*var[ijpkp] + ats*var[ijmkp] - + abn*var[ijpkm] + abs*var[ijmkm]; + + abn*var[ijpkm] + absol*var[ijmkm]; if (Nstorage) residual +=Nlin[ijk]; |