From 776e8a14c5294578e10fff98333262ac5494297d Mon Sep 17 00:00:00 2001 From: tiglio Date: Thu, 18 May 2006 17:35:58 +0000 Subject: Added a small ifdef to make it work with PUGH (as it was the derivs were wrong near real boundaries if using PUGH). One can do something better but I did not want to mess up the MP case. It will not work with reflection symmetries and PUGH as it is, though, so this is a temporary solution only. git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@69 f69c4107-0314-4c4f-9ad4-17e986b73f4a --- src/call_derivs.c | 12 +++++++++++- 1 file changed, 11 insertions(+), 1 deletion(-) diff --git a/src/call_derivs.c b/src/call_derivs.c index a49930c..a6da543 100644 --- a/src/call_derivs.c +++ b/src/call_derivs.c @@ -4,6 +4,7 @@ #include "util_ErrorCodes.h" #include "util_Table.h" +#include "cctk_DefineThorn.h" void DiffGv ( const CCTK_POINTER_TO_CONST cctkGH_, const CCTK_INT dir, const CCTK_REAL *var, CCTK_REAL *dvar, @@ -136,7 +137,15 @@ void DiffGv ( const CCTK_POINTER_TO_CONST cctkGH_, const CCTK_INT dir, CCTK_WARN(0,"Util_TableGetIntArray returned error"); } - pen_sym_handle = SymmetryHandleOfName ( "multipatch" ); + + pen_sym_handle = SymmetryHandleOfName ( "multipatch" ); + +#ifdef CACTUSPUGH_PUGH // a little hack to make this work with pugh (without reflection sym, though) + for (int i=0; i<6; i++){ + symbnd[i] = 1; + pen_sym_handle=1; + } +#endif switch(dir) { case 0: { @@ -188,6 +197,7 @@ void DiffGv ( const CCTK_POINTER_TO_CONST cctkGH_, const CCTK_INT dir, CCTK_WARN (0, "Wrong direction specified"); } + if ( CCTK_Equals(norm_type,"Diagonal") ) { if ( CCTK_Equals(operator_type,"Minimal Bandwidth") ) { switch(loc_order) { -- cgit v1.2.3