diff options
-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 |