diff options
Diffstat (limited to 'Examples/SimpleWaveCaKernel/src/CaKernel__calc_rhs.code')
-rw-r--r-- | Examples/SimpleWaveCaKernel/src/CaKernel__calc_rhs.code | 13 |
1 files changed, 9 insertions, 4 deletions
diff --git a/Examples/SimpleWaveCaKernel/src/CaKernel__calc_rhs.code b/Examples/SimpleWaveCaKernel/src/CaKernel__calc_rhs.code index 35ba87e..40d5b56 100644 --- a/Examples/SimpleWaveCaKernel/src/CaKernel__calc_rhs.code +++ b/Examples/SimpleWaveCaKernel/src/CaKernel__calc_rhs.code @@ -1,8 +1,11 @@ +#undef KRANC_DIFF_FUNCTIONS #define KRANC_C - #include "Differencing.h" #include "GenericFD.h" +#define KRANC_GFOFFSET3D(u,i,j,k) I3D(u,i,j,k) + + /* Define macros used in calculations */ #define INITVALUE (42) #define QAD(x) (SQR(SQR(x))) @@ -10,8 +13,6 @@ #define SQR(x) ((x) * (x)) #define CUB(x) ((x) * (x) * (x)) -#define CCTK_GFINDEX3D(u,i,j,k) I3D(u,i,j,k) - CAKERNEL_calc_rhs_Begin /* Include user-supplied include files */ @@ -60,11 +61,15 @@ CAKERNEL_calc_rhs_Begin /* Include user supplied include files */ /* Precompute derivatives */ + CCTK_REAL const PDstandard2nd11phi = PDstandard2nd11(phi); + CCTK_REAL const PDstandard2nd22phi = PDstandard2nd22(phi); + CCTK_REAL const PDstandard2nd33phi = PDstandard2nd33(phi); /* Calculate temporaries and grid functions */ CCTK_REAL phirhsL = piL; - CCTK_REAL pirhsL = phiL; + CCTK_REAL pirhsL = PDstandard2nd11phi + PDstandard2nd22phi + + PDstandard2nd33phi; /* Copy local copies back to grid functions */ I3D(phirhs,0,0,0) = phirhsL; |