From e3719c2bc603c6e88fe41f9986980a62b5fb36d0 Mon Sep 17 00:00:00 2001 From: schnetter Date: Mon, 25 Apr 2005 10:43:02 +0000 Subject: Correct second derivative operator for first order dissipation git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/Dissipation/trunk@8 850bcc8b-0e4f-0410-8c26-8d28fbf1eda9 --- src/apply_dissipation.F77 | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/src/apply_dissipation.F77 b/src/apply_dissipation.F77 index b5d6939..74cf168 100644 --- a/src/apply_dissipation.F77 +++ b/src/apply_dissipation.F77 @@ -20,9 +20,9 @@ c $Header$ do i = 2, ni-1 rhs(i,j,k) = rhs(i,j,k) + epsdis - $ * (+ (var(i-1,j,k) - 2*var(i,j,k) - var(i+1,j,k)) / dx(1) - $ + (var(i,j-1,k) - 2*var(i,j,k) - var(i,j+1,k)) / dx(2) - $ + (var(i,j,k-1) - 2*var(i,j,k) - var(i,j,k+1)) / dx(3)) + $ * (+ (var(i-1,j,k) - 2*var(i,j,k) + var(i+1,j,k)) / dx(1) + $ + (var(i,j-1,k) - 2*var(i,j,k) + var(i,j+1,k)) / dx(2) + $ + (var(i,j,k-1) - 2*var(i,j,k) + var(i,j,k+1)) / dx(3)) end do end do -- cgit v1.2.3