(+) implemented poisoning: added three parameters

    - poison_derivatives
    - poison_dissipation
    - poison_value
    The first two set values of derivative GF and RHS respectively
    to poison_value in boundary_shiftout zone, which is specified by
    CoordBase shiftout parameters.


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

1 files changed, 13 insertions, 0 deletions

diff --git a/src/dissipation.c b/src/dissipation.c
index a3b2bb6..4dbacbf 100644
--- a/src/dissipation.c
+++ b/src/dissipation.c
@@ -147,6 +147,14 @@ void CCTK_FCALL CCTK_FNAME(dissipation_8_4_alt) (const CCTK_REAL *var,
                                              CCTK_REAL *rhs);
 void get_grid_offsets (CCTK_INT *offset);
 
+void CCTK_FCALL CCTK_FNAME(SBP_Poisoning) (  const CCTK_INT *ni, 
+                                             const CCTK_INT *nj, 
+                                             const CCTK_INT *nk, 
+                                             const CCTK_INT *bb, 
+                                             const CCTK_INT *offset, 
+                                             CCTK_REAL *dvar);
+
+
 static void
 apply (int const varindex, char const * const optstring, void * const arg);
 
@@ -330,4 +338,9 @@ apply (int const varindex, char const * const optstring, void * const arg)
       assert(0);
     }
   }
+
+  if (poison_dissipation) {
+      CCTK_FNAME(SBP_Poisoning)(&cctk_lsh[0], &cctk_lsh[1], &cctk_lsh[2],
+                                bbox, offset, rhsptr);
+  }
 }