diff options
Diffstat (limited to 'src/sor_flat.F')
-rw-r--r-- | src/sor_flat.F | 4 |
1 files changed, 2 insertions, 2 deletions
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 |