diff options
author | pollney <pollney@fa3da13c-9f13-4301-a575-cf5b8c5e1907> | 2000-03-02 10:28:41 +0000 |
---|---|---|
committer | pollney <pollney@fa3da13c-9f13-4301-a575-cf5b8c5e1907> | 2000-03-02 10:28:41 +0000 |
commit | e1a9bf40216e6e1542849b6595c43b5e882cd50e (patch) | |
tree | 7c2459148ee94fe10a861935f27ed339f78d0448 | |
parent | af75009bc85763049d5b6b2609ba743bc3f476f2 (diff) |
Switch signs for M and N in the EllBase interface so that it now solves
L\phi + M\phi + N\phi = 0
(in line with BAM_Elliptic).
git-svn-id: http://svn.cactuscode.org/arrangements/CactusElliptic/EllSOR/trunk@53 fa3da13c-9f13-4301-a575-cf5b8c5e1907
-rw-r--r-- | src/sor_confmetric.F | 4 | ||||
-rw-r--r-- | src/sor_flat.F | 4 |
2 files changed, 4 insertions, 4 deletions
diff --git a/src/sor_confmetric.F b/src/sor_confmetric.F index 37dbf77..8c7867a 100644 --- a/src/sor_confmetric.F +++ b/src/sor_confmetric.F @@ -235,7 +235,7 @@ c Start loop with Red Black & -uyy(i,j+1,k) -2.*uyy(i,j,k) -uyy(i,j-1,k) & -uzz(i,j,k+1) -2.*uzz(i,j,k) -uzz(i,j,k-1) if (Mlinear_storage.eq.1) then - ac = ac - Mlinear(i,j,k) + ac = ac + Mlinear(i,j,k) endif ae = uxx(i+1,j,k)+uxx(i,j,k) @@ -268,7 +268,7 @@ c Start loop with Red Black & + atn*var(i,j+1,k+1) + ats*var(i,j-1,k+1) & + abn*var(i,j+1,k-1) + asb*var(i,j-1,k-1) if (Nsource_storage.eq.1) then - residual = residual - Nsource(i,j,k) + residual = residual + Nsource(i,j,k) endif c Accumulate to the total Norm of the residual diff --git a/src/sor_flat.F b/src/sor_flat.F index 803e0ed..fad105b 100644 --- a/src/sor_flat.F +++ b/src/sor_flat.F @@ -200,7 +200,7 @@ c Start loop with Red Black ac = ac_orig if (Mlinear_storage.eq.1) then - ac = ac - Mlinear(i,j,k) + ac = ac + Mlinear(i,j,k) endif residual = ac*var(i,j,k) @@ -209,7 +209,7 @@ c Start loop with Red Black & + at*var(i,j,k+1) + ab*var(i,j,k-1) if (Nsource_storage.eq.1) then - residual = residual - Nsource(i,j,k) + residual = residual + Nsource(i,j,k) endif c Accumulate to the total Norm of the residual |