Cactus Code Thorn TmunuBase
Author(s)    : Erik Schnetter <schnetter@gmail.com>
Maintainer(s): Erik Schnetter <schnetter@gmail.com>
Licence      : LGPL
--------------------------------------------------------------------------
1. Purpose
Provide storage for the stress-energy tensor T_\mu\nu

# Interface definition for thorn TmunuBase
IMPLEMENTS: TmunuBase
USES INCLUDE HEADER: loop.hxx
PUBLIC:
CCTK_REAL eTtt TYPE=gf TAGS='index={0 0 0}' "T_00"
CCTK_REAL eTti TYPE=gf TAGS='index={0 0 0}' { eTtx eTty eTtz } "T_0i"
CCTK_REAL eTij TYPE=gf TAGS='index={0 0 0}' { eTxx eTxy eTxz eTyy eTyz eTzz } "T_ij"

# Parameter definitions for thorn TmunuBase

# Schedule definitions for thorn TmunuBase
SCHEDULE GROUP TmunuBase_SetTmunuVars AT initial
{
} "Schedule group for setting T_munu"
SCHEDULE GROUP TmunuBase_SetTmunuVars AT postregrid
{
} "Schedule group for setting T_munu"
SCHEDULE GROUP TmunuBase_SetTmunuVars AT poststep
{
} "Schedule group for setting T_munu"
SCHEDULE GROUP TmunuBase_SetTmunuVars IN ODESolvers_PostStep
{
} "Schedule group for setting T_munu"
SCHEDULE TmunuBase_ZeroTmunu IN TmunuBase_SetTmunuVars
{
  LANG: C
  WRITES: eTtt(everywhere) eTti(everywhere) eTij(everywhere)
} "Set T_munu to zero"
SCHEDULE GROUP TmunuBase_AddToTmunu IN TmunuBase_SetTmunuVars AFTER TmunuBase_ZeroTmunu
{
} "Add to T_munu here"

# Main make.code.defn file for thorn TmunuBase
# Source files in this directory
SRCS = tmunubase.cxx
# Subdirectories containing source files
SUBDIRS =

#include <loop.hxx>
#include <cctk.h>
#include <cctk_Arguments_Checked.h>
#include <cctk_Parameters.h>
namespace TmunuBase {
using namespace Loop;
extern "C" void TmunuBase_ZeroTmunu(CCTK_ARGUMENTS) {
  DECLARE_CCTK_ARGUMENTS_TmunuBase_ZeroTmunu;
  DECLARE_CCTK_PARAMETERS;
  const GF3D<CCTK_REAL, 0, 0, 0> eTtt_(cctkGH, eTtt);
  const GF3D<CCTK_REAL, 0, 0, 0> eTtx_(cctkGH, eTtx);
  const GF3D<CCTK_REAL, 0, 0, 0> eTty_(cctkGH, eTty);
  const GF3D<CCTK_REAL, 0, 0, 0> eTtz_(cctkGH, eTtz);
  const GF3D<CCTK_REAL, 0, 0, 0> eTxx_(cctkGH, eTxx);
  const GF3D<CCTK_REAL, 0, 0, 0> eTxy_(cctkGH, eTxy);
  const GF3D<CCTK_REAL, 0, 0, 0> eTxz_(cctkGH, eTxz);
  const GF3D<CCTK_REAL, 0, 0, 0> eTyy_(cctkGH, eTyy);
  const GF3D<CCTK_REAL, 0, 0, 0> eTyz_(cctkGH, eTyz);
  const GF3D<CCTK_REAL, 0, 0, 0> eTzz_(cctkGH, eTzz);
  loop_all<0, 0, 0>(cctkGH, [&](const PointDesc &p) { eTtt_(p.I) = 0; });
  loop_all<0, 0, 0>(cctkGH, [&](const PointDesc &p) { eTtx_(p.I) = 0; });
  loop_all<0, 0, 0>(cctkGH, [&](const PointDesc &p) { eTty_(p.I) = 0; });
  loop_all<0, 0, 0>(cctkGH, [&](const PointDesc &p) { eTtz_(p.I) = 0; });
  loop_all<0, 0, 0>(cctkGH, [&](const PointDesc &p) { eTxx_(p.I) = 0; });
  loop_all<0, 0, 0>(cctkGH, [&](const PointDesc &p) { eTxy_(p.I) = 0; });
  loop_all<0, 0, 0>(cctkGH, [&](const PointDesc &p) { eTxz_(p.I) = 0; });
  loop_all<0, 0, 0>(cctkGH, [&](const PointDesc &p) { eTyy_(p.I) = 0; });
  loop_all<0, 0, 0>(cctkGH, [&](const PointDesc &p) { eTyz_(p.I) = 0; });
  loop_all<0, 0, 0>(cctkGH, [&](const PointDesc &p) { eTzz_(p.I) = 0; });
}
} // namespace TmunuBase