Add BSSN RHS.

Remove Mathematica output trees; keep only Cactus version of trees.

1 files changed, 95 insertions, 25 deletions

diff --git a/ML_BSSN/interface.ccl b/ML_BSSN/interface.ccl
index d68366c..5dfcf8c 100644
--- a/ML_BSSN/interface.ccl
+++ b/ML_BSSN/interface.ccl
@@ -31,13 +31,13 @@ CCTK_INT FUNCTION Boundary_SelectVarForBC(CCTK_POINTER_TO_CONST IN GH, CCTK_INT
 USES FUNCTION Boundary_SelectVarForBC
 
 public:
-CCTK_REAL cons_detg type=GF timelevels=1 tags='tensortypealias="Scalar" tensorweight=1.0000000000000000000'
+CCTK_REAL cons_detg type=GF timelevels=1 tags='tensortypealias="Scalar" tensorweight=2.0000000000000000000'
 {
   cS
 } "cons_detg"
 
 public:
-CCTK_REAL cons_Gamma type=GF timelevels=1 tags='tensortypealias="U" tensorweight=1.0000000000000000000'
+CCTK_REAL cons_Gamma type=GF timelevels=1 tags='tensortypealias="U" tensorweight=0.66666666666666666667'
 {
   cXt1,
   cXt2,
@@ -51,7 +51,91 @@ CCTK_REAL cons_traceA type=GF timelevels=1 tags='tensortypealias="Scalar" tensor
 } "cons_traceA"
 
 public:
-CCTK_REAL curv type=GF timelevels=1 tags='tensortypealias="DD_sym" tensorweight=1.0000000000000000000'
+CCTK_REAL Ham type=GF timelevels=1 tags='tensortypealias="Scalar" tensorweight=1.0000000000000000000'
+{
+  H
+} "Ham"
+
+public:
+CCTK_REAL mom type=GF timelevels=1 tags='tensortypealias="D" tensorweight=1.0000000000000000000'
+{
+  M1,
+  M2,
+  M3
+} "mom"
+
+public:
+CCTK_REAL curvrhs type=GF timelevels=1 tags='tensortypealias="DD_sym" tensorweight=-0.66666666666666666667'
+{
+  At11rhs,
+  At21rhs,
+  At31rhs,
+  At22rhs,
+  At32rhs,
+  At33rhs
+} "curvrhs"
+
+public:
+CCTK_REAL dtlapserhs type=GF timelevels=1 tags='tensortypealias="Scalar" tensorweight=1.0000000000000000000'
+{
+  dtalpharhs
+} "dtlapserhs"
+
+public:
+CCTK_REAL dtshiftrhs type=GF timelevels=1 tags='tensortypealias="U" tensorweight=1.0000000000000000000'
+{
+  dtbeta1rhs,
+  dtbeta2rhs,
+  dtbeta3rhs
+} "dtshiftrhs"
+
+public:
+CCTK_REAL Gammarhs type=GF timelevels=1 tags='tensortypealias="U" tensorweight=0.66666666666666666667'
+{
+  Xt1rhs,
+  Xt2rhs,
+  Xt3rhs
+} "Gammarhs"
+
+public:
+CCTK_REAL lapserhs type=GF timelevels=1 tags='tensortypealias="Scalar" tensorweight=1.0000000000000000000'
+{
+  alpharhs
+} "lapserhs"
+
+public:
+CCTK_REAL log_confacrhs type=GF timelevels=1 tags='tensortypealias="Scalar" tensorweight=0.16666666666666666667'
+{
+  phirhs
+} "log_confacrhs"
+
+public:
+CCTK_REAL metricrhs type=GF timelevels=1 tags='tensortypealias="DD_sym" tensorweight=-0.66666666666666666667'
+{
+  gt11rhs,
+  gt21rhs,
+  gt31rhs,
+  gt22rhs,
+  gt32rhs,
+  gt33rhs
+} "metricrhs"
+
+public:
+CCTK_REAL shiftrhs type=GF timelevels=1 tags='tensortypealias="U" tensorweight=1.0000000000000000000'
+{
+  beta1rhs,
+  beta2rhs,
+  beta3rhs
+} "shiftrhs"
+
+public:
+CCTK_REAL trace_curvrhs type=GF timelevels=1 tags='tensortypealias="Scalar" tensorweight=1.0000000000000000000'
+{
+  trKrhs
+} "trace_curvrhs"
+
+public:
+CCTK_REAL curv type=GF timelevels=3 tags='tensortypealias="DD_sym" tensorweight=-0.66666666666666666667'
 {
   At11,
   At21,
@@ -62,13 +146,13 @@ CCTK_REAL curv type=GF timelevels=1 tags='tensortypealias="DD_sym" tensorweight=
 } "curv"
 
 public:
-CCTK_REAL dtlapse type=GF timelevels=1 tags='tensortypealias="Scalar" tensorweight=1.0000000000000000000'
+CCTK_REAL dtlapse type=GF timelevels=3 tags='tensortypealias="Scalar" tensorweight=1.0000000000000000000'
 {
   dtalpha
 } "dtlapse"
 
 public:
-CCTK_REAL dtshift type=GF timelevels=1 tags='tensortypealias="U" tensorweight=1.0000000000000000000'
+CCTK_REAL dtshift type=GF timelevels=3 tags='tensortypealias="U" tensorweight=1.0000000000000000000'
 {
   dtbeta1,
   dtbeta2,
@@ -76,7 +160,7 @@ CCTK_REAL dtshift type=GF timelevels=1 tags='tensortypealias="U" tensorweight=1.
 } "dtshift"
 
 public:
-CCTK_REAL Gamma type=GF timelevels=1 tags='tensortypealias="U" tensorweight=1.0000000000000000000'
+CCTK_REAL Gamma type=GF timelevels=3 tags='tensortypealias="U" tensorweight=0.66666666666666666667'
 {
   Xt1,
   Xt2,
@@ -84,25 +168,19 @@ CCTK_REAL Gamma type=GF timelevels=1 tags='tensortypealias="U" tensorweight=1.00
 } "Gamma"
 
 public:
-CCTK_REAL Ham type=GF timelevels=1 tags='tensortypealias="Scalar" tensorweight=1.0000000000000000000'
-{
-  H
-} "Ham"
-
-public:
-CCTK_REAL lapse type=GF timelevels=1 tags='tensortypealias="Scalar" tensorweight=1.0000000000000000000'
+CCTK_REAL lapse type=GF timelevels=3 tags='tensortypealias="Scalar" tensorweight=1.0000000000000000000'
 {
   alpha
 } "lapse"
 
 public:
-CCTK_REAL log_confac type=GF timelevels=1 tags='tensortypealias="Scalar" tensorweight=1.0000000000000000000'
+CCTK_REAL log_confac type=GF timelevels=3 tags='tensortypealias="Scalar" tensorweight=0.16666666666666666667'
 {
   phi
 } "log_confac"
 
 public:
-CCTK_REAL metric type=GF timelevels=1 tags='tensortypealias="DD_sym" tensorweight=1.0000000000000000000'
+CCTK_REAL metric type=GF timelevels=3 tags='tensortypealias="DD_sym" tensorweight=-0.66666666666666666667'
 {
   gt11,
   gt21,
@@ -113,15 +191,7 @@ CCTK_REAL metric type=GF timelevels=1 tags='tensortypealias="DD_sym" tensorweigh
 } "metric"
 
 public:
-CCTK_REAL mom type=GF timelevels=1 tags='tensortypealias="D" tensorweight=1.0000000000000000000'
-{
-  M1,
-  M2,
-  M3
-} "mom"
-
-public:
-CCTK_REAL shift type=GF timelevels=1 tags='tensortypealias="U" tensorweight=1.0000000000000000000'
+CCTK_REAL shift type=GF timelevels=3 tags='tensortypealias="U" tensorweight=1.0000000000000000000'
 {
   beta1,
   beta2,
@@ -129,7 +199,7 @@ CCTK_REAL shift type=GF timelevels=1 tags='tensortypealias="U" tensorweight=1.00
 } "shift"
 
 public:
-CCTK_REAL trace_curv type=GF timelevels=1 tags='tensortypealias="Scalar" tensorweight=1.0000000000000000000'
+CCTK_REAL trace_curv type=GF timelevels=3 tags='tensortypealias="Scalar" tensorweight=1.0000000000000000000'
 {
   trK
 } "trace_curv"