diff options
Diffstat (limited to 'src/macro/WAVEG_guts.h')
| -rw-r--r-- | src/macro/WAVEG_guts.h | 113 |
1 files changed, 113 insertions, 0 deletions
diff --git a/src/macro/WAVEG_guts.h b/src/macro/WAVEG_guts.h new file mode 100644 index 0000000..512daba --- /dev/null +++ b/src/macro/WAVEG_guts.h @@ -0,0 +1,113 @@ +/*@@ + @header WAVEG_guts.h + @date Jul 98 + @author Gabrielle Allen + @desc + Macro to calculate wave operator acting on the metric + + That is g^lm g_ij,lm + + @enddesc +@@*/ + +#ifndef WAVEG_GUTS +#define WAVEG_GUTS + +#include "../../packages/CactusEinstein/Einstein/src/macro/UPPERMET_guts.h" +#include "../../packages/CactusEinstein/Einstein/src/macro/DDG_guts.h" + +#ifdef FCODE + + WAVEG_DDGXX = UPPERMET_UXX*DDG_DXXDGXX + + & 2D0*UPPERMET_UXY*DDG_DXYDGXX + + & 2D0*UPPERMET_UXZ*DDG_DXZDGXX + + & UPPERMET_UYY*DDG_DYYDGXX + + & 2D0*UPPERMET_UYZ*DDG_DYZDGXX + + & UPPERMET_UZZ*DDG_DZZDGXX + + WAVEG_DDGXY = UPPERMET_UXX*DDG_DXXDGXY + + & 2D0*UPPERMET_UXY*DDG_DXYDGXY + + & 2D0*UPPERMET_UXZ*DDG_DXZDGXY + + & UPPERMET_UYY*DDG_DYYDGXY + + & 2D0*UPPERMET_UYZ*DDG_DYZDGXY + + & UPPERMET_UZZ*DDG_DZZDGXY + + WAVEG_DDGXZ = UPPERMET_UXX*DDG_DXXDGXZ + + & 2D0*UPPERMET_UXY*DDG_DXYDGXZ + + & 2D0*UPPERMET_UXZ*DDG_DXZDGXZ + + & UPPERMET_UYY*DDG_DYYDGXZ + + & 2D0*UPPERMET_UYZ*DDG_DYZDGXZ + + & UPPERMET_UZZ*DDG_DZZDGXZ + + WAVEG_DDGYY = UPPERMET_UXX*DDG_DXXDGYY + + & 2D0*UPPERMET_UXY*DDG_DXYDGYY + + & 2D0*UPPERMET_UXZ*DDG_DXZDGYY + + & UPPERMET_UYY*DDG_DYYDGYY + + & 2D0*UPPERMET_UYZ*DDG_DYZDGYY + + & UPPERMET_UZZ*DDG_DZZDGYY + + WAVEG_DDGYZ = UPPERMET_UXX*DDG_DXXDGYZ + + & 2D0*UPPERMET_UXY*DDG_DXYDGYZ + + & 2D0*UPPERMET_UXZ*DDG_DXZDGYZ + + & UPPERMET_UYY*DDG_DYYDGYZ + + & 2D0*UPPERMET_UYZ*DDG_DYZDGYZ + + & UPPERMET_UZZ*DDG_DZZDGYZ + + WAVEG_DDGZZ = UPPERMET_UXX*DDG_DXXDGZZ + + & 2D0*UPPERMET_UXY*DDG_DXYDGZZ + + & 2D0*UPPERMET_UXZ*DDG_DXZDGZZ + + & UPPERMET_UYY*DDG_DYYDGZZ + + & 2D0*UPPERMET_UYZ*DDG_DYZDGZZ + + & UPPERMET_UZZ*DDG_DZZDGZZ + +#endif + +#ifdef CCODE + + WAVEG_DDGXX = UPPERMET_UXX*DDG_DXXDGXX + + 2*UPPERMET_UXY*DDG_DXYDGXX + + 2*UPPERMET_UXZ*DDG_DXZDGXX + + UPPERMET_UYY*DDG_DYYDGXX + + 2*UPPERMET_UYZ*DDG_DYZDGXX + + UPPERMET_UZZ*DDG_DYYDGXX; + + WAVEG_DDGXY = UPPERMET_UXX*DDG_DXXDGXY + + 2*UPPERMET_UXY*DDG_DXYDGXY + + 2*UPPERMET_UXZ*DDG_DXZDGXY + + UPPERMET_UYY*DDG_DYYDGXY + + 2*UPPERMET_UYZ*DDG_DYZDGXY + + UPPERMET_UZZ*DDG_DYYDGXY; + + WAVEG_DDGXZ = UPPERMET_UXX*DDG_DXXDGXZ + + 2*UPPERMET_UXY*DDG_DXYDGXZ + + 2*UPPERMET_UXZ*DDG_DXZDGXZ + + UPPERMET_UYY*DDG_DYYDGXZ + + 2*UPPERMET_UYZ*DDG_DYZDGXZ + + UPPERMET_UZZ*DDG_DYYDGXZ; + + WAVEG_DDGYY = UPPERMET_UXX*DDG_DXXDGYY + + 2*UPPERMET_UXY*DDG_DXYDGYY + + 2*UPPERMET_UXZ*DDG_DXZDGYY + + UPPERMET_UYY*DDG_DYYDGYY + + 2*UPPERMET_UYZ*DDG_DYZDGYY + + UPPERMET_UZZ*DDG_DYYDGYY; + + WAVEG_DDGYZ = UPPERMET_UXX*DDG_DXXDGYZ + + 2*UPPERMET_UXY*DDG_DXYDGYZ + + 2*UPPERMET_UXZ*DDG_DXZDGYZ + + UPPERMET_UYY*DDG_DYYDGYZ + + 2*UPPERMET_UYZ*DDG_DYZDGYZ + + UPPERMET_UZZ*DDG_DYYDGYZ; + + WAVEG_DDGZZ = UPPERMET_UXX*DDG_DXXDGZZ + + 2*UPPERMET_UXY*DDG_DXYDGZZ + + 2*UPPERMET_UXZ*DDG_DXZDGZZ + + UPPERMET_UYY*DDG_DYYDGZZ + + 2*UPPERMET_UYZ*DDG_DYZDGZZ + + UPPERMET_UZZ*DDG_DYYDGZZ; + + +#endif + +#endif + |