diff options
author | ianhin <ianhin> | 2004-06-20 20:16:50 +0000 |
---|---|---|
committer | ianhin <ianhin> | 2004-06-20 20:16:50 +0000 |
commit | b5f570f18e2d12e5f33836b7a44c2bfdcf8d9f73 (patch) | |
tree | 977510294e413e6c7df2279114bf2f0445cdb139 /Auxiliary/Cactus/KrancNumericalTools | |
parent | 4050e520a67a6a937cee80922d81430b33fa4072 (diff) |
Extremely large hack to allow us to use D+ as a difference operator. See
accompanying email
----------------------------------------------------------------------
Diffstat (limited to 'Auxiliary/Cactus/KrancNumericalTools')
-rw-r--r-- | Auxiliary/Cactus/KrancNumericalTools/GenericFD/src/GenericFD.h | 20 |
1 files changed, 19 insertions, 1 deletions
diff --git a/Auxiliary/Cactus/KrancNumericalTools/GenericFD/src/GenericFD.h b/Auxiliary/Cactus/KrancNumericalTools/GenericFD/src/GenericFD.h index 9bc1002..d9061e1 100644 --- a/Auxiliary/Cactus/KrancNumericalTools/GenericFD/src/GenericFD.h +++ b/Auxiliary/Cactus/KrancNumericalTools/GenericFD/src/GenericFD.h @@ -424,4 +424,22 @@ #define D1_c2c4(gf,i,j,k) (fdweight_c2*D1_c2(gf, i,j,k) + fdweight_c4*D1_c4(gf,i,j,k)) #define D2_c2c4(gf,i,j,k) (fdweight_c2*D2_c2(gf, i,j,k) + fdweight_c4*D2_c4(gf,i,j,k)) #define D3_c2c4(gf,i,j,k) (fdweight_c2*D3_c2(gf, i,j,k) + fdweight_c4*D3_c4(gf,i,j,k)) - + + +/*****************************************************/ +/* */ +/* One-sided derivatives */ +/* (only for C atm) */ +/******************************************************/ + +#define Dplus1(gf,i,j,k) \ + ((gf[CCTK_GFINDEX3D(cctkGH,i+1,j,k)] \ + - gf[CCTK_GFINDEX3D(cctkGH,i,j,k)]) * dxi) + +#define Dplus2(gf,i,j,k) \ + ((gf[CCTK_GFINDEX3D(cctkGH,i,j+1,k)] \ + - gf[CCTK_GFINDEX3D(cctkGH,i,j,k)]) * dyi) + +#define Dplus3(gf,i,j,k) \ + ((gf[CCTK_GFINDEX3D(cctkGH,i,j,k+1)] \ + - gf[CCTK_GFINDEX3D(cctkGH,i,j,k)]) * dzi) |