1 files changed, 67 insertions, 0 deletions
diff --git a/src/macro/CDZCDK_guts.h b/src/macro/CDZCDK_guts.h
new file mode 100644
index 0000000..20a99db
--- /dev/null
+++ b/src/macro/CDZCDK_guts.h
@@ -0,0 +1,67 @@
+/*@@
+  @header   CDZCDK_guts.h
+  @date     August 98
+  @author   Gabrielle Allen
+  @desc
+  Macro to calculate the first covariant derivatives of the 
+  extrinsic curvature with respect to y
+  @enddesc
+@@*/
+
+#ifndef CDZCDK_GUTS
+#define CDZCDK_GUTS
+
+#include "../../packages/CactusEinstein/Einstein/src/macro/DZDK_guts.h"
+#include "../../packages/CactusEinstein/Einstein/src/macro/CHR2_guts.h"
+
+#ifdef FCODE 
+
+      CDZCDK_CDZCDKXX = DZDK_DZDKXX - 2D0*(CHR2_XXZ*CDZCDK_KXX + 
+     &  CHR2_YXZ*CDZCDK_KXY + CHR2_ZXZ*CDZCDK_KXZ)
+
+      CDZCDK_CDZCDKYY = DZDK_DZDKYY - 2D0*(CHR2_XYZ*CDZCDK_KXY + 
+     &  CHR2_YYZ*CDZCDK_KYY + CHR2_ZYZ*CDZCDK_KYZ)
+
+      CDZCDK_CDZCDKZZ = DZDK_DZDKZZ - 2D0*(CHR2_XZZ*CDZCDK_KXZ + 
+     &  CHR2_YZZ*CDZCDK_KYZ + CHR2_ZZZ*CDZCDK_KZZ) 
+
+      CDZCDK_CDZCDKXY = DZDK_DZDKXY-(CHR2_XYZ*CDZCDK_KXX) - 
+     &  (CHR2_XXZ + CHR2_YYZ)*CDZCDK_KXY  - CHR2_ZYZ*CDZCDK_KXZ - 
+     &  CHR2_YXZ*CDZCDK_KYY - CHR2_ZXZ*CDZCDK_KYZ
+
+      CDZCDK_CDZCDKXZ = DZDK_DZDKXZ-(CHR2_XZZ*CDZCDK_KXX) -
+     &  CHR2_YZZ*CDZCDK_KXY - (CHR2_XXZ + CHR2_ZZZ)*CDZCDK_KXZ 
+     &  - CHR2_YXZ*CDZCDK_KYZ - CHR2_ZXZ*CDZCDK_KZZ
+
+      CDZCDK_CDZCDKYZ = DZDK_DZDKYZ-(CHR2_XZZ*CDZCDK_KXY) - 
+     &  CHR2_XYZ*CDZCDK_KXZ - CHR2_YZZ*CDZCDK_KYY - 
+     & (CHR2_YYZ + CHR2_ZZZ)*CDZCDK_KYZ - CHR2_ZYZ*CDZCDK_KZZ
+
+#endif
+
+#ifdef CCODE
+
+      CDZCDK_CDZCDKXX = DZDK_DZDKXX - 2D0*(CHR2_XXZ*CDZCDK_KXX + 
+       CHR2_YXZ*CDZCDK_KXY + CHR2_ZXZ*CDZCDK_KXZ);
+
+      CDZCDK_CDZCDKYY = DZDK_DZDKYY - 2D0*(CHR2_XYZ*CDZCDK_KXY + 
+       CHR2_YYZ*CDZCDK_KYY + CHR2_ZYZ*CDZCDK_KYZ);
+
+      CDZCDK_CDZCDKZZ = DZDK_DZDKZZ - 2D0*(CHR2_XZZ*CDZCDK_KXZ + 
+       CHR2_YZZ*CDZCDK_KYZ + CHR2_ZZZ*CDZCDK_KZZ); 
+
+      CDZCDK_CDZCDKXY = DZDK_DZDKXY-(CHR2_XYZ*CDZCDK_KXX) - 
+       (CHR2_XXZ + CHR2_YYZ)*CDZCDK_KXY  - CHR2_ZYZ*CDZCDK_KXZ - 
+       CHR2_YXZ*CDZCDK_KYY - CHR2_ZXZ*CDZCDK_KYZ;
+
+      CDZCDK_CDZCDKXZ = DZDK_DZDKXZ-(CHR2_XZZ*CDZCDK_KXX) -
+       CHR2_YZZ*CDZCDK_KXY - (CHR2_XXZ + CHR2_ZZZ)*CDZCDK_KXZ 
+       - CHR2_YXZ*CDZCDK_KYZ - CHR2_ZXZ*CDZCDK_KZZ;
+
+      CDZCDK_CDZCDKYZ = DZDK_DZDKYZ-(CHR2_XZZ*CDZCDK_KXY) - 
+       CHR2_XYZ*CDZCDK_KXZ - CHR2_YZZ*CDZCDK_KYY - 
+      (CHR2_YYZ + CHR2_ZZZ)*CDZCDK_KYZ - CHR2_ZYZ*CDZCDK_KZZ;
+
+#endif
+
+#endif