eschnett/cactusamrex - Change O4RDTKKLC3OK5MSXNMTMOXOONLV3XSMH4OPZMANRBVYRFVJJ5KGQC

ADMBase: New thorn

Created by Erik Schnetter on December 16, 2019

O4RDTKKLC3OK5MSXNMTMOXOONLV3XSMH4OPZMANRBVYRFVJJ5KGQC

Dependencies

In channels

main

master

Change contents

File addition: ADMBase (drwxr-xr-x)
[1.0]

File addition: README (-rw-r--r--)

[0.10]

Cactus Code Thorn ADMBase
Author(s)    : Erik Schnetter <schnetter@gmail.com>
Maintainer(s): Erik Schnetter <schnetter@gmail.com>
Licence      : LGPL
--------------------------------------------------------------------------
1. Purpose
Provide storage for the ADM variables

File addition: interface.ccl (-rw-r--r--)

[0.10]

# Interface definition for thorn ADMBase
IMPLEMENTS: ADMBase
USES INCLUDE HEADER: loop.hxx
CCTK_REAL metric TYPE=gf TAGS='index={0 0 0}' { gxx gxy gxz gyy gyz gzz } "ADM 3-metric g_ij"
CCTK_REAL curv TYPE=gf TAGS='index={0 0 0}' { kxx kxy kxz kyy kyz kzz } "ADM extrinsic curvature K_ij"
CCTK_REAL lapse TYPE=gf TAGS='index={0 0 0}' { alp } "ADM lapse function alpha"
CCTK_REAL dtlapse TYPE=gf TAGS='index={0 0 0}' { dtalp } "Time derivative of ADM lapse function"
CCTK_REAL shift TYPE=gf TAGS='index={0 0 0}' { betax betay betaz } "ADM shift vector beta^i"
CCTK_REAL dtshift TYPE=gf TAGS='index={0 0 0}' { dtbetax dtbetay dtbetaz } "Time derivative of ADM shift vector"

File addition: param.ccl (-rw-r--r--)

[0.10]

# Parameter definitions for thorn ADMBase
RESTRICTED:
KEYWORD initial_data "Initial metric and extrinsic curvature datasets"
{
  "Cartesian Minkowski" :: "Minkowski geometry in Cartesian coordinates"
} "Cartesian Minkowski"
KEYWORD initial_lapse "Initial lapse value"
{
  "one" :: "Uniform lapse"
} "one"
KEYWORD initial_shift "Initial shift value"
{
  "zero" :: "Shift is zero"
} "zero"
KEYWORD initial_dtlapse "Initial dtlapse value"
{
  "zero" :: "Dtlapse is zero"
} "zero"
KEYWORD initial_dtshift "Initial dtshift value"
{
  "zero" :: "Dtshift is zero"
} "zero"

File addition: schedule.ccl (-rw-r--r--)

[0.10]

# Schedule definitions for thorn ADMBase
if (CCTK_EQUALS(initial_data, "Cartesian Minkowski")) {
  SCHEDULE ADMBase_initial_data IN initial
  {
    LANG: C
    WRITES: metric(everywhere) curv(everywhere)
  } "Set up Cartesian Minkowski initial data"
}
if (CCTK_EQUALS(initial_lapse, "one")) {
  SCHEDULE ADMBase_initial_lapse IN initial
  {
    LANG: C
    WRITES: lapse(everywhere)
  } "Set lapse to one"
}
if (CCTK_EQUALS(initial_dtlapse, "zero")) {
  SCHEDULE ADMBase_initial_dtlapse IN initial
  {
    LANG: C
    WRITES: dtlapse(everywhere)
  } "Set dtlapse to zero"
}
if (CCTK_EQUALS(initial_shift, "zero")) {
  SCHEDULE ADMBase_initial_shift IN initial
  {
    LANG: C
    WRITES: shift(everywhere)
  } "Set shift to zero"
}
if (CCTK_EQUALS(initial_dtshift, "zero")) {
  SCHEDULE ADMBase_initial_dtshift IN initial
  {
    LANG: C
    WRITES: dtshift(everywhere)
  } "Set dtshift to zero"
}

File addition: src (drwxr-xr-x)
[0.10]

File addition: initial.cxx (-rw-r--r--)

[0.2509]

#include <loop.hxx>
#include <cctk.h>
#include <cctk_Arguments_Checked.h>
#include <cctk_Parameters.h>
namespace ADMBase {
extern "C" void ADMBase_initial_data(CCTK_ARGUMENTS) {
  DECLARE_CCTK_ARGUMENTS_ADMBase_initial_data;
  DECLARE_CCTK_PARAMETERS;
  const Loop::GF3D<CCTK_REAL, 0, 0, 0> gxx_(cctkGH, gxx);
  const Loop::GF3D<CCTK_REAL, 0, 0, 0> gxy_(cctkGH, gxy);
  const Loop::GF3D<CCTK_REAL, 0, 0, 0> gxz_(cctkGH, gxz);
  const Loop::GF3D<CCTK_REAL, 0, 0, 0> gyy_(cctkGH, gyy);
  const Loop::GF3D<CCTK_REAL, 0, 0, 0> gyz_(cctkGH, gyz);
  const Loop::GF3D<CCTK_REAL, 0, 0, 0> gzz_(cctkGH, gzz);
  const Loop::GF3D<CCTK_REAL, 0, 0, 0> kxx_(cctkGH, kxx);
  const Loop::GF3D<CCTK_REAL, 0, 0, 0> kxy_(cctkGH, kxy);
  const Loop::GF3D<CCTK_REAL, 0, 0, 0> kxz_(cctkGH, kxz);
  const Loop::GF3D<CCTK_REAL, 0, 0, 0> kyy_(cctkGH, kyy);
  const Loop::GF3D<CCTK_REAL, 0, 0, 0> kyz_(cctkGH, kyz);
  const Loop::GF3D<CCTK_REAL, 0, 0, 0> kzz_(cctkGH, kzz);
  Loop::loop_all<0, 0, 0>(cctkGH,
                          [&](const Loop::PointDesc &p) { gxx_(p.I) = 1; });
  Loop::loop_all<0, 0, 0>(cctkGH,
                          [&](const Loop::PointDesc &p) { gxy_(p.I) = 0; });
  Loop::loop_all<0, 0, 0>(cctkGH,
                          [&](const Loop::PointDesc &p) { gxz_(p.I) = 0; });
  Loop::loop_all<0, 0, 0>(cctkGH,
                          [&](const Loop::PointDesc &p) { gyy_(p.I) = 1; });
  Loop::loop_all<0, 0, 0>(cctkGH,
                          [&](const Loop::PointDesc &p) { gyz_(p.I) = 0; });
  Loop::loop_all<0, 0, 0>(cctkGH,
                          [&](const Loop::PointDesc &p) { gzz_(p.I) = 1; });
  Loop::loop_all<0, 0, 0>(cctkGH,
                          [&](const Loop::PointDesc &p) { kxx_(p.I) = 0; });
  Loop::loop_all<0, 0, 0>(cctkGH,
                          [&](const Loop::PointDesc &p) { kxy_(p.I) = 0; });
  Loop::loop_all<0, 0, 0>(cctkGH,
                          [&](const Loop::PointDesc &p) { kxz_(p.I) = 0; });
  Loop::loop_all<0, 0, 0>(cctkGH,
                          [&](const Loop::PointDesc &p) { kyy_(p.I) = 0; });
  Loop::loop_all<0, 0, 0>(cctkGH,
                          [&](const Loop::PointDesc &p) { kyz_(p.I) = 0; });
  Loop::loop_all<0, 0, 0>(cctkGH,
                          [&](const Loop::PointDesc &p) { kzz_(p.I) = 0; });
}
extern "C" void ADMBase_initial_lapse(CCTK_ARGUMENTS) {
  DECLARE_CCTK_ARGUMENTS_ADMBase_initial_lapse;
  DECLARE_CCTK_PARAMETERS;
  const Loop::GF3D<CCTK_REAL, 0, 0, 0> alp_(cctkGH, alp);
  Loop::loop_all<0, 0, 0>(cctkGH,
                          [&](const Loop::PointDesc &p) { alp_(p.I) = 1; });
}
extern "C" void ADMBase_initial_dtlapse(CCTK_ARGUMENTS) {
  DECLARE_CCTK_ARGUMENTS_ADMBase_initial_dtlapse;
  DECLARE_CCTK_PARAMETERS;
  const Loop::GF3D<CCTK_REAL, 0, 0, 0> dtalp_(cctkGH, dtalp);
  Loop::loop_all<0, 0, 0>(cctkGH,
                          [&](const Loop::PointDesc &p) { dtalp_(p.I) = 0; });
}
extern "C" void ADMBase_initial_shift(CCTK_ARGUMENTS) {
  DECLARE_CCTK_ARGUMENTS_ADMBase_initial_shift;
  DECLARE_CCTK_PARAMETERS;
  const Loop::GF3D<CCTK_REAL, 0, 0, 0> betax_(cctkGH, betax);
  const Loop::GF3D<CCTK_REAL, 0, 0, 0> betay_(cctkGH, betay);
  const Loop::GF3D<CCTK_REAL, 0, 0, 0> betaz_(cctkGH, betaz);
  Loop::loop_all<0, 0, 0>(cctkGH,
                          [&](const Loop::PointDesc &p) { betax_(p.I) = 0; });
  Loop::loop_all<0, 0, 0>(cctkGH,
                          [&](const Loop::PointDesc &p) { betay_(p.I) = 0; });
  Loop::loop_all<0, 0, 0>(cctkGH,
                          [&](const Loop::PointDesc &p) { betaz_(p.I) = 0; });
}
extern "C" void ADMBase_initial_dtshift(CCTK_ARGUMENTS) {
  DECLARE_CCTK_ARGUMENTS_ADMBase_initial_dtshift;
  DECLARE_CCTK_PARAMETERS;
  const Loop::GF3D<CCTK_REAL, 0, 0, 0> dtbetax_(cctkGH, dtbetax);
  const Loop::GF3D<CCTK_REAL, 0, 0, 0> dtbetay_(cctkGH, dtbetay);
  const Loop::GF3D<CCTK_REAL, 0, 0, 0> dtbetaz_(cctkGH, dtbetaz);
  Loop::loop_all<0, 0, 0>(cctkGH,
                          [&](const Loop::PointDesc &p) { dtbetax_(p.I) = 0; });
  Loop::loop_all<0, 0, 0>(cctkGH,
                          [&](const Loop::PointDesc &p) { dtbetay_(p.I) = 0; });
  Loop::loop_all<0, 0, 0>(cctkGH,
                          [&](const Loop::PointDesc &p) { dtbetaz_(p.I) = 0; });
}
} // namespace ADMBase

File addition: make.code.defn (-rw-r--r--)

[0.2509]

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

Insertion in azure-pipelines/carpetx.th at line 127 [4.11]
[2.41]
[3.878]
```
CarpetX/ADMBase
```

ADMBase: New thorn

Dependencies

In channels

Change contents

File addition: ADMBase (drwxr-xr-x)

File addition: README (-rw-r--r--)

File addition: interface.ccl (-rw-r--r--)

File addition: param.ccl (-rw-r--r--)

File addition: schedule.ccl (-rw-r--r--)

File addition: src (drwxr-xr-x)

File addition: initial.cxx (-rw-r--r--)

File addition: make.code.defn (-rw-r--r--)

Insertion in azure-pipelines/carpetx.th at line 127 [4.11]