Merge pull request #60 from EinsteinToolkit/eschnett/kerrschild

KerrSchildX: New thorn

Author: eschnett

KerrSchildX/README

@@ -0,0 +1,9 @@
+Cactus Code Thorn KerrSchildX
+Author(s)    : Erik Schnetter <[email protected]>
+Maintainer(s): Erik Schnetter <[email protected]>
+Licence      : LGPL
+--------------------------------------------------------------------------
+
+1. Purpose
+
+Set up Kerr-Schild initial conditions

KerrSchildX/configuration.ccl

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

KerrSchildX/interface.ccl

@@ -0,0 +1,5 @@
+# Interface definition for thorn KerrSchildX
+
+IMPLEMENTS: KerrSchildX
+
+INHERITS: ADMBaseX CoordinatesX

KerrSchildX/par/kerrschild.par

@@ -0,0 +1,43 @@
+ActiveThorns = "
+    ADMBaseX
+    CarpetX
+    CoordinatesX
+    Formaline
+    IOUtil
+    KerrSchildX
+"
+
+Cactus::cctk_show_schedule = yes   # no
+CarpetX::verbose = no
+
+Cactus::presync_mode = "mixed-error"
+
+CarpetX::xmin = -4.0
+CarpetX::ymin = -4.0
+CarpetX::zmin = -4.0
+
+CarpetX::xmax = +4.0
+CarpetX::ymax = +4.0
+CarpetX::zmax = +4.0
+
+CarpetX::ncells_x = 16
+CarpetX::ncells_y = 16
+CarpetX::ncells_z = 16
+
+CarpetX::blocking_factor_x = 16
+CarpetX::blocking_factor_y = 16
+CarpetX::blocking_factor_z = 16
+
+ADMBaseX::initial_data = "kerrschildx"
+ADMBaseX::initial_lapse = "kerrschildx"
+ADMBaseX::initial_shift = "kerrschildx"
+ADMBaseX::initial_dtlapse = "kerrschildx"
+ADMBaseX::initial_dtshift = "kerrschildx"
+
+Cactus::cctk_itlast = 0
+
+IO::out_dir = $parfile
+CarpetX::out_tsv_every = 1
+CarpetX::out_tsv_vars = "ADMBaseX::metric ADMBaseX::curv ADMBaseX::lapse ADMBaseX::shift ADMBaseX::dtlapse ADMBaseX::dtshift"
+CarpetX::out_silo_every = 1
+CarpetX::out_silo_vars = "ADMBaseX::metric ADMBaseX::curv ADMBaseX::lapse ADMBaseX::shift ADMBaseX::dtlapse ADMBaseX::dtshift"

KerrSchildX/param.ccl

@@ -0,0 +1,47 @@
+# Parameter definitions for thorn KerrSchild
+
+SHARES: ADMBaseX
+
+EXTENDS KEYWORD initial_data 
+{
+  "kerrschildx" :: "Kerr-Schild black hole"
+} 
+
+EXTENDS KEYWORD initial_lapse
+{
+  "kerrschildx" :: "Kerr-Schild black hole"
+}
+
+EXTENDS KEYWORD initial_shift
+{
+  "kerrschildx" :: "Kerr-Schild black hole"
+}
+
+EXTENDS KEYWORD initial_dtlapse
+{
+  "kerrschildx" :: "Kerr-Schild black hole"
+}
+
+EXTENDS KEYWORD initial_dtshift
+{
+  "kerrschildx" :: "Kerr-Schild black hole"
+}
+
+
+
+PRIVATE:
+
+CCTK_REAL mass "Mass"
+{
+  0.0:* :: ""
+} 1.0
+
+CCTK_REAL spin "Spin parameter (with dimension `mass`)"
+{
+  *:* :: "Must be -mass < spin < mass"
+} 0.0
+
+CCTK_REAL epsilon "Avoid division by zero"
+{
+  0:* :: ""
+} 1.0e-8

KerrSchildX/schedule.ccl

@@ -0,0 +1,26 @@
+# Schedule definitions for thorn KerrSchild
+
+if (CCTK_Equals(initial_data, "kerrschildx") ||
+    CCTK_Equals(initial_lapse, "kerrschildx") ||
+    CCTK_Equals(initial_shift, "kerrschildx") ||
+    CCTK_Equals(initial_dtlapse, "kerrschildx") ||
+    CCTK_Equals(initial_dtshift, "kerrschildx"))
+{
+  SCHEDULE KerrSchildX_ParamCheck AT paramcheck
+  {
+    LANG: C
+    OPTIONS: global
+  } "Check Kerr-Schild parameters"
+
+  SCHEDULE KerrSchildX_InitialData IN ADMBaseX_SetADMVars
+  {
+    LANG: C
+    READS: CoordinatesX::vertex_coords(everywhere)
+    WRITES: ADMBaseX::metric(everywhere)
+    WRITES: ADMBaseX::curv(everywhere)
+    WRITES: ADMBaseX::lapse(everywhere)
+    WRITES: ADMBaseX::shift(everywhere)
+    WRITES: ADMBaseX::dtlapse(everywhere)
+    WRITES: ADMBaseX::dtshift(everywhere)
+  } "Create Kerr-Schild initial data"
+}