From e1a9bf40216e6e1542849b6595c43b5e882cd50e Mon Sep 17 00:00:00 2001 From: pollney Date: Thu, 2 Mar 2000 10:28:41 +0000 Subject: 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 --- src/sor_confmetric.F | 4 ++-- 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 -- cgit v1.2.3