Add new parameter use_dissipation_near_excision.

Add OpenMP parallelisation.


git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/Dissipation/trunk@30 850bcc8b-0e4f-0410-8c26-8d28fbf1eda9

1 files changed, 6 insertions, 2 deletions

diff --git a/src/apply_dissipation.F77 b/src/apply_dissipation.F77
index d67a02e..b4e3545 100644
--- a/src/apply_dissipation.F77
+++ b/src/apply_dissipation.F77
@@ -2,12 +2,12 @@ c     $Header$
 
 #include "cctk.h"
 
-      subroutine apply_dissipation (var, rhs, ni, nj, nk, dx, dt, order, epsdis)
+      subroutine apply_dissipation (var, rhs, ni, nj, nk, dx, order, epsdis)
       implicit none
       
       integer   ni, nj, nk
       CCTK_REAL var(ni,nj,nk), rhs(ni,nj,nk)
-      CCTK_REAL dx(3), dt
+      CCTK_REAL dx(3)
       CCTK_INT  order
       CCTK_REAL epsdis(ni,nj,nk)
       
@@ -15,6 +15,7 @@ c     $Header$
       
       if (order .eq. 1) then
          
+c$omp parallel do
          do k = 2, nk-1
             do j = 2, nj-1
                do i = 2, ni-1
@@ -29,6 +30,7 @@ c     $Header$
          
       else if (order .eq. 3) then
          
+c$omp parallel do
          do k = 3, nk-2
             do j = 3, nj-2
                do i = 3, ni-2
@@ -44,6 +46,7 @@ c     $Header$
          
       else if (order .eq. 5) then
          
+c$omp parallel do
         do k = 4, nk-3
           do j = 4, nj-3
             do i = 4, ni-3
@@ -59,6 +62,7 @@ c     $Header$
         
       else if (order .eq. 7) then
          
+c$omp parallel do
         do k = 5, nk-4
           do j = 5, nj-4
             do i = 5, ni-4