Only use one sided finite differences near multipatch boundaries. Should fix

the problems we have with periodic boundary conditions. Also added the missing
declarations of fortran routines with Kreiss-Oliger type dissipation for
certain cases.


git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@58 f69c4107-0314-4c4f-9ad4-17e986b73f4a

2 files changed, 95 insertions, 4 deletions

diff --git a/src/call_derivs.c b/src/call_derivs.c
index c4a0fb8..a49930c 100644
--- a/src/call_derivs.c
+++ b/src/call_derivs.c
@@ -17,6 +17,9 @@ void DiffGv ( const CCTK_POINTER_TO_CONST cctkGH_, const CCTK_INT dir,
   CCTK_REAL delta;
   CCTK_INT ierr;
   CCTK_INT lsh[3], bbox[6], bb[2], nghostzones[3];
+  CCTK_INT symtable, pen_sym_handle;
+  CCTK_INT symbnd[6];
+
   int nelements;
   void CCTK_FCALL CCTK_FNAME(deriv_gf_2_1)(const CCTK_REAL *var,
                                        const CCTK_INT *ni,
@@ -123,22 +126,61 @@ void DiffGv ( const CCTK_POINTER_TO_CONST cctkGH_, const CCTK_INT dir,
     loc_order = order;
   }
 
+  symtable = SymmetryTableHandleForGrid (cctkGH);
+  if (symtable<0) {
+    CCTK_WARN(0,"symtable is out of bounds");
+  }
+
+  ierr = Util_TableGetIntArray (symtable, 6, symbnd, "symmetry_handle");
+  if (ierr!=6) {
+    CCTK_WARN(0,"Util_TableGetIntArray returned error");
+  }
+
+  pen_sym_handle = SymmetryHandleOfName ( "multipatch" );
+
   switch(dir) {
   case 0: {
     delta = CCTK_DELTA_SPACE(0);
-    bb[0] = cctk_bbox[0]; bb[1] = cctk_bbox[1];
+    if ( symbnd[0] == pen_sym_handle && cctk_bbox[0] == 1 ) {
+      bb[0] = 1;
+    } else {
+      bb[0] = 0;
+    }
+    if ( symbnd[1] == pen_sym_handle && cctk_bbox[1] == 1 ) {
+      bb[1] = 1;
+    } else {
+      bb[1] = 0;
+    }
     gsize = cctk_nghostzones[0];
     break;
   }
   case 1: {
     delta = CCTK_DELTA_SPACE(1);
-    bb[0] = cctk_bbox[2]; bb[1] = cctk_bbox[3];
+    if ( symbnd[2] == pen_sym_handle && cctk_bbox[2] == 1) {
+      bb[0] = 1;
+    } else {
+      bb[0] = 0;
+    }
+    if ( symbnd[3] == pen_sym_handle && cctk_bbox[3] == 1) {
+      bb[1] = 1;
+    } else {
+      bb[1] = 0;
+    }
     gsize = cctk_nghostzones[1];
     break;
   }
   case 2: {
     delta = CCTK_DELTA_SPACE(2);
-    bb[0] = cctk_bbox[4]; bb[1] = cctk_bbox[5];
+    if ( symbnd[4] == pen_sym_handle && cctk_bbox[4] == 1) {
+      bb[0] = 1;
+    } else {
+      bb[0] = 0;
+    }
+    if ( symbnd[4] == pen_sym_handle && cctk_bbox[5] == 1) {
+      bb[1] = 1;
+    } else {
+      bb[1] = 0;
+    }
     gsize = cctk_nghostzones[2];
     break;
   }
diff --git a/src/dissipation.c b/src/dissipation.c
index cc93b33..48f7ba4 100644
--- a/src/dissipation.c
+++ b/src/dissipation.c
@@ -7,6 +7,9 @@
 #include "cctk_Arguments.h"
 #include "cctk_Parameters.h"
 
+#include "util_ErrorCodes.h"
+#include "util_Table.h"
+
 void CCTK_FCALL CCTK_FNAME(dissipation_4_3_opt) (const CCTK_REAL *var,
                                              const CCTK_INT *lsh,
                                              const CCTK_INT *gsh,
@@ -17,6 +20,16 @@ void CCTK_FCALL CCTK_FNAME(dissipation_4_3_opt) (const CCTK_REAL *var,
                                              const CCTK_REAL *epsdis,
                                              const CCTK_REAL *diss_fraction,
                                              CCTK_REAL *rhs);
+void CCTK_FCALL CCTK_FNAME(dissipation_4_3_alt) (const CCTK_REAL *var,
+                                             const CCTK_INT *lsh,
+                                             const CCTK_INT *gsh,
+                                             const CCTK_INT *lbnd,
+                                             const CCTK_INT *bbox,
+                                             const CCTK_INT *gsize,
+                                             const CCTK_REAL *dx,
+                                             const CCTK_REAL *epsdis,
+                                             const CCTK_REAL *diss_fraction,
+                                             CCTK_REAL *rhs);
 void CCTK_FCALL CCTK_FNAME(dissipation_6_5_opt) (const CCTK_REAL *var,
                                              const CCTK_INT *ni,
                                              const CCTK_INT *nj,
@@ -64,6 +77,15 @@ void CCTK_FCALL CCTK_FNAME(dissipation_4_2) (const CCTK_REAL *var,
                                              const CCTK_REAL *dx,
                                              const CCTK_REAL *epsdis,
                                              CCTK_REAL *rhs);
+void CCTK_FCALL CCTK_FNAME(dissipation_4_2_alt) (const CCTK_REAL *var,
+                                             const CCTK_INT *ni,
+                                             const CCTK_INT *nj,
+                                             const CCTK_INT *nk,
+                                             const CCTK_INT *bbox,
+                                             const CCTK_INT *gsize,
+                                             const CCTK_REAL *dx,
+                                             const CCTK_REAL *epsdis,
+                                             CCTK_REAL *rhs);
 void CCTK_FCALL CCTK_FNAME(dissipation_6_3) (const CCTK_REAL *var,
                                              const CCTK_INT *ni,
                                              const CCTK_INT *nj,
@@ -73,6 +95,15 @@ void CCTK_FCALL CCTK_FNAME(dissipation_6_3) (const CCTK_REAL *var,
                                              const CCTK_REAL *dx,
                                              const CCTK_REAL *epsdis,
                                              CCTK_REAL *rhs);
+void CCTK_FCALL CCTK_FNAME(dissipation_6_3_alt) (const CCTK_REAL *var,
+                                             const CCTK_INT *ni,
+                                             const CCTK_INT *nj,
+                                             const CCTK_INT *nk,
+                                             const CCTK_INT *bbox,
+                                             const CCTK_INT *gsize,
+                                             const CCTK_REAL *dx,
+                                             const CCTK_REAL *epsdis,
+                                             CCTK_REAL *rhs);
 void CCTK_FCALL CCTK_FNAME(dissipation_8_4) (const CCTK_REAL *var,
                                              const CCTK_INT *ni,
                                              const CCTK_INT *nj,
@@ -120,6 +151,8 @@ apply (int const varindex, char const * const optstring, void * const arg)
   int n;
   int d;
   int ierr;
+  CCTK_INT symtable, pen_sym_handle;
+  CCTK_INT symbnd[6];
   
   assert (varindex >= 0);
   
@@ -128,8 +161,24 @@ apply (int const varindex, char const * const optstring, void * const arg)
     gsize[d] = cctk_nghostzones[d];
   }
   
+  symtable = SymmetryTableHandleForGrid (cctkGH);
+  if (symtable<0) {
+    CCTK_WARN(0,"symtable is out of bounds");
+  }
+
+  ierr = Util_TableGetIntArray (symtable, 6, symbnd, "symmetry_handle");
+  if (ierr!=6) {
+    CCTK_WARN(0,"Util_TableGetIntArray returned error");
+  }
+
+  pen_sym_handle = SymmetryHandleOfName ( "multipatch" );
+
   for (d=0; d<6; ++d) {
-    bbox[d] = cctk_bbox[d];
+    if ( symbnd[d] == pen_sym_handle && cctk_bbox[d] == 1 ) {
+      bbox[d] = 1;
+    } else {
+      bbox[d] = 0;
+    }
   }
 
   rhsindex = MoLQueryEvolvedRHS (varindex);