diff options
author | schnetter <schnetter@850bcc8b-0e4f-0410-8c26-8d28fbf1eda9> | 2005-04-25 10:43:02 +0000 |
---|---|---|
committer | schnetter <schnetter@850bcc8b-0e4f-0410-8c26-8d28fbf1eda9> | 2005-04-25 10:43:02 +0000 |
commit | e3719c2bc603c6e88fe41f9986980a62b5fb36d0 (patch) | |
tree | 9b61f9fa3677a000dd4709691513fb27687c827b | |
parent | b3267c3188ef8be05d13a33e215b42277ccaf56b (diff) |
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
-rw-r--r-- | src/apply_dissipation.F77 | 6 |
1 files 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 |