CurvTmunu: add thorn

Author: lwJi

CurvTmunu/README

@@ -0,0 +1,9 @@
+Cactus Code Thorn CurvTmunu
+Author(s)    : Liwei Ji <jiliwei.phys@gmail.com>
+Maintainer(s): Liwei Ji <jiliwei.phys@gmail.com>
+Licence      : LGPL
+--------------------------------------------------------------------------
+
+1. Purpose
+
+Provide storage for the stress-energy tensor T_\mu\nu

CurvTmunu/configuration.ccl

@@ -0,0 +1,3 @@
+# Configuration definition for thorn CurvTmunu
+
+REQUIRES Loop

CurvTmunu/interface.ccl

@@ -0,0 +1,14 @@
+# Interface definition for thorn CurvTmunu
+
+IMPLEMENTS: CurvTmunu
+
+USES INCLUDE HEADER: loop_device.hxx
+
+
+PUBLIC:
+
+CCTK_REAL eTtt TYPE=gf CENTERING={ccc} TAGS='checkpoint="no"' "T_00"
+
+CCTK_REAL eTti TYPE=gf CENTERING={ccc} TAGS='checkpoint="no"' { eTtx eTty eTtz } "T_0i"
+
+CCTK_REAL eTij TYPE=gf CENTERING={ccc} TAGS='checkpoint="no"' { eTxx eTxy eTxz eTyy eTyz eTzz } "T_ij"

CurvTmunu/param.ccl

@@ -0,0 +1 @@
+# Parameter definitions for thorn CurvTmunu

CurvTmunu/schedule.ccl

@@ -0,0 +1,30 @@
+# Schedule definitions for thorn CurvTmunu
+
+SCHEDULE GROUP CurvTmunu_SetTmunuVars AT initial AFTER CurvADM_SetADMVars
+{
+} "Schedule group for setting T_munu"
+
+if (CCTK_IsThornActive("ODESolvers")) {
+  SCHEDULE GROUP CurvTmunu_SetTmunuVars IN ODESolvers_PostStep AFTER CurvADM_SetADMVars
+  {
+  } "Schedule group for setting T_munu"
+} else {
+  SCHEDULE GROUP CurvTmunu_SetTmunuVars AT postregrid AFTER CurvADM_SetADMVars
+  {
+  } "Schedule group for setting T_munu"
+
+  SCHEDULE GROUP CurvTmunu_SetTmunuVars AT poststep AFTER CurvADM_SetADMVars
+  {
+  } "Schedule group for setting T_munu"
+}
+
+
+SCHEDULE CurvTmunu_ZeroTmunu IN CurvTmunu_SetTmunuVars
+{
+  LANG: C
+  WRITES: eTtt(everywhere) eTti(everywhere) eTij(everywhere)
+} "Set T_munu to zero"
+
+SCHEDULE GROUP CurvTmunu_AddToTmunu IN CurvTmunu_SetTmunuVars AFTER CurvTmunu_ZeroTmunu
+{
+} "Add to T_munu here"

CurvTmunu/src/make.code.defn

@@ -0,0 +1,7 @@
+# Main make.code.defn file for thorn CurvTmunu
+
+# Source files in this directory
+SRCS = tmunu.cxx
+
+# Subdirectories containing source files
+SUBDIRS =

CurvTmunu/src/tmunu.cxx

@@ -0,0 +1,33 @@
+#include <loop_device.hxx>
+
+#include <cctk.h>
+#include <cctk_Arguments.h>
+#include <cctk_Parameters.h>
+
+#include <array>
+
+namespace CurvTmunu {
+using namespace std;
+using namespace Loop;
+
+extern "C" void CurvTmunu_ZeroTmunu(CCTK_ARGUMENTS) {
+  DECLARE_CCTK_ARGUMENTSX_CurvTmunu_ZeroTmunu;
+  DECLARE_CCTK_PARAMETERS;
+
+  grid.loop_all_device<1, 1, 1>(grid.nghostzones,
+                                [=] CCTK_DEVICE(const PointDesc &p)
+                                    CCTK_ATTRIBUTE_ALWAYS_INLINE {
+                                      eTtt(p.I) = 0;
+                                      eTtx(p.I) = 0;
+                                      eTty(p.I) = 0;
+                                      eTtz(p.I) = 0;
+                                      eTxx(p.I) = 0;
+                                      eTxy(p.I) = 0;
+                                      eTxz(p.I) = 0;
+                                      eTyy(p.I) = 0;
+                                      eTyz(p.I) = 0;
+                                      eTzz(p.I) = 0;
+                                    });
+}
+
+} // namespace CurvTmunu