aboutsummaryrefslogtreecommitdiff
path: root/ML_ADM
diff options
context:
space:
mode:
authorIan Hinder <ian.hinder@aei.mpg.de>2013-09-12 22:09:32 +0200
committerIan Hinder <ian.hinder@aei.mpg.de>2013-09-12 22:09:32 +0200
commitcaadfb1214b142be4512d86994fdfdac9a45206b (patch)
treeb20c527044ad74f0cb953296f456c8c429cb9cb7 /ML_ADM
parent8c4ac493e00d4cd3978d4504372b3c37b85c64ea (diff)
Regenerate code (Kranc commit 3237af2566109f099ee78a7effb511a65dc18ae0)
Diffstat (limited to 'ML_ADM')
-rw-r--r--ML_ADM/src/Boundaries.cc728
-rw-r--r--ML_ADM/src/RegisterSymmetries.cc34
2 files changed, 381 insertions, 381 deletions
diff --git a/ML_ADM/src/Boundaries.cc b/ML_ADM/src/Boundaries.cc
index f150f48..bc3f0e5 100644
--- a/ML_ADM/src/Boundaries.cc
+++ b/ML_ADM/src/Boundaries.cc
@@ -32,6 +32,17 @@ extern "C" void ML_ADM_SelectBoundConds(CCTK_ARGUMENTS)
CCTK_INT ierr CCTK_ATTRIBUTE_UNUSED = 0;
+ if (CCTK_EQUALS(ML_metric_bound, "none" ) ||
+ CCTK_EQUALS(ML_metric_bound, "static") ||
+ CCTK_EQUALS(ML_metric_bound, "flat" ) ||
+ CCTK_EQUALS(ML_metric_bound, "zero" ) )
+ {
+ ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, -1,
+ "ML_ADM::ML_metric", ML_metric_bound);
+ if (ierr < 0)
+ CCTK_WARN(0, "Failed to register ML_metric_bound BC for ML_ADM::ML_metric!");
+ }
+
if (CCTK_EQUALS(ML_curv_bound, "none" ) ||
CCTK_EQUALS(ML_curv_bound, "static") ||
CCTK_EQUALS(ML_curv_bound, "flat" ) ||
@@ -54,17 +65,6 @@ extern "C" void ML_ADM_SelectBoundConds(CCTK_ARGUMENTS)
CCTK_WARN(0, "Failed to register ML_lapse_bound BC for ML_ADM::ML_lapse!");
}
- if (CCTK_EQUALS(ML_metric_bound, "none" ) ||
- CCTK_EQUALS(ML_metric_bound, "static") ||
- CCTK_EQUALS(ML_metric_bound, "flat" ) ||
- CCTK_EQUALS(ML_metric_bound, "zero" ) )
- {
- ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, -1,
- "ML_ADM::ML_metric", ML_metric_bound);
- if (ierr < 0)
- CCTK_WARN(0, "Failed to register ML_metric_bound BC for ML_ADM::ML_metric!");
- }
-
if (CCTK_EQUALS(ML_shift_bound, "none" ) ||
CCTK_EQUALS(ML_shift_bound, "static") ||
CCTK_EQUALS(ML_shift_bound, "flat" ) ||
@@ -76,6 +76,72 @@ extern "C" void ML_ADM_SelectBoundConds(CCTK_ARGUMENTS)
CCTK_WARN(0, "Failed to register ML_shift_bound BC for ML_ADM::ML_shift!");
}
+ if (CCTK_EQUALS(g11_bound, "none" ) ||
+ CCTK_EQUALS(g11_bound, "static") ||
+ CCTK_EQUALS(g11_bound, "flat" ) ||
+ CCTK_EQUALS(g11_bound, "zero" ) )
+ {
+ ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, -1,
+ "ML_ADM::g11", g11_bound);
+ if (ierr < 0)
+ CCTK_WARN(0, "Failed to register g11_bound BC for ML_ADM::g11!");
+ }
+
+ if (CCTK_EQUALS(g12_bound, "none" ) ||
+ CCTK_EQUALS(g12_bound, "static") ||
+ CCTK_EQUALS(g12_bound, "flat" ) ||
+ CCTK_EQUALS(g12_bound, "zero" ) )
+ {
+ ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, -1,
+ "ML_ADM::g12", g12_bound);
+ if (ierr < 0)
+ CCTK_WARN(0, "Failed to register g12_bound BC for ML_ADM::g12!");
+ }
+
+ if (CCTK_EQUALS(g13_bound, "none" ) ||
+ CCTK_EQUALS(g13_bound, "static") ||
+ CCTK_EQUALS(g13_bound, "flat" ) ||
+ CCTK_EQUALS(g13_bound, "zero" ) )
+ {
+ ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, -1,
+ "ML_ADM::g13", g13_bound);
+ if (ierr < 0)
+ CCTK_WARN(0, "Failed to register g13_bound BC for ML_ADM::g13!");
+ }
+
+ if (CCTK_EQUALS(g22_bound, "none" ) ||
+ CCTK_EQUALS(g22_bound, "static") ||
+ CCTK_EQUALS(g22_bound, "flat" ) ||
+ CCTK_EQUALS(g22_bound, "zero" ) )
+ {
+ ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, -1,
+ "ML_ADM::g22", g22_bound);
+ if (ierr < 0)
+ CCTK_WARN(0, "Failed to register g22_bound BC for ML_ADM::g22!");
+ }
+
+ if (CCTK_EQUALS(g23_bound, "none" ) ||
+ CCTK_EQUALS(g23_bound, "static") ||
+ CCTK_EQUALS(g23_bound, "flat" ) ||
+ CCTK_EQUALS(g23_bound, "zero" ) )
+ {
+ ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, -1,
+ "ML_ADM::g23", g23_bound);
+ if (ierr < 0)
+ CCTK_WARN(0, "Failed to register g23_bound BC for ML_ADM::g23!");
+ }
+
+ if (CCTK_EQUALS(g33_bound, "none" ) ||
+ CCTK_EQUALS(g33_bound, "static") ||
+ CCTK_EQUALS(g33_bound, "flat" ) ||
+ CCTK_EQUALS(g33_bound, "zero" ) )
+ {
+ ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, -1,
+ "ML_ADM::g33", g33_bound);
+ if (ierr < 0)
+ CCTK_WARN(0, "Failed to register g33_bound BC for ML_ADM::g33!");
+ }
+
if (CCTK_EQUALS(K11_bound, "none" ) ||
CCTK_EQUALS(K11_bound, "static") ||
CCTK_EQUALS(K11_bound, "flat" ) ||
@@ -153,72 +219,6 @@ extern "C" void ML_ADM_SelectBoundConds(CCTK_ARGUMENTS)
CCTK_WARN(0, "Failed to register alpha_bound BC for ML_ADM::alpha!");
}
- if (CCTK_EQUALS(g11_bound, "none" ) ||
- CCTK_EQUALS(g11_bound, "static") ||
- CCTK_EQUALS(g11_bound, "flat" ) ||
- CCTK_EQUALS(g11_bound, "zero" ) )
- {
- ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, -1,
- "ML_ADM::g11", g11_bound);
- if (ierr < 0)
- CCTK_WARN(0, "Failed to register g11_bound BC for ML_ADM::g11!");
- }
-
- if (CCTK_EQUALS(g12_bound, "none" ) ||
- CCTK_EQUALS(g12_bound, "static") ||
- CCTK_EQUALS(g12_bound, "flat" ) ||
- CCTK_EQUALS(g12_bound, "zero" ) )
- {
- ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, -1,
- "ML_ADM::g12", g12_bound);
- if (ierr < 0)
- CCTK_WARN(0, "Failed to register g12_bound BC for ML_ADM::g12!");
- }
-
- if (CCTK_EQUALS(g13_bound, "none" ) ||
- CCTK_EQUALS(g13_bound, "static") ||
- CCTK_EQUALS(g13_bound, "flat" ) ||
- CCTK_EQUALS(g13_bound, "zero" ) )
- {
- ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, -1,
- "ML_ADM::g13", g13_bound);
- if (ierr < 0)
- CCTK_WARN(0, "Failed to register g13_bound BC for ML_ADM::g13!");
- }
-
- if (CCTK_EQUALS(g22_bound, "none" ) ||
- CCTK_EQUALS(g22_bound, "static") ||
- CCTK_EQUALS(g22_bound, "flat" ) ||
- CCTK_EQUALS(g22_bound, "zero" ) )
- {
- ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, -1,
- "ML_ADM::g22", g22_bound);
- if (ierr < 0)
- CCTK_WARN(0, "Failed to register g22_bound BC for ML_ADM::g22!");
- }
-
- if (CCTK_EQUALS(g23_bound, "none" ) ||
- CCTK_EQUALS(g23_bound, "static") ||
- CCTK_EQUALS(g23_bound, "flat" ) ||
- CCTK_EQUALS(g23_bound, "zero" ) )
- {
- ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, -1,
- "ML_ADM::g23", g23_bound);
- if (ierr < 0)
- CCTK_WARN(0, "Failed to register g23_bound BC for ML_ADM::g23!");
- }
-
- if (CCTK_EQUALS(g33_bound, "none" ) ||
- CCTK_EQUALS(g33_bound, "static") ||
- CCTK_EQUALS(g33_bound, "flat" ) ||
- CCTK_EQUALS(g33_bound, "zero" ) )
- {
- ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, -1,
- "ML_ADM::g33", g33_bound);
- if (ierr < 0)
- CCTK_WARN(0, "Failed to register g33_bound BC for ML_ADM::g33!");
- }
-
if (CCTK_EQUALS(beta1_bound, "none" ) ||
CCTK_EQUALS(beta1_bound, "static") ||
CCTK_EQUALS(beta1_bound, "flat" ) ||
@@ -252,6 +252,25 @@ extern "C" void ML_ADM_SelectBoundConds(CCTK_ARGUMENTS)
CCTK_WARN(0, "Failed to register beta3_bound BC for ML_ADM::beta3!");
}
+ if (CCTK_EQUALS(ML_metric_bound, "radiative"))
+ {
+ /* select radiation boundary condition */
+ static CCTK_INT handle_ML_metric_bound CCTK_ATTRIBUTE_UNUSED = -1;
+ if (handle_ML_metric_bound < 0) handle_ML_metric_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE);
+ if (handle_ML_metric_bound < 0) CCTK_WARN(0, "could not create table!");
+ if (Util_TableSetReal(handle_ML_metric_bound , ML_metric_bound_limit, "LIMIT") < 0)
+ CCTK_WARN(0, "could not set LIMIT value in table!");
+ if (Util_TableSetReal(handle_ML_metric_bound ,ML_metric_bound_speed, "SPEED") < 0)
+ CCTK_WARN(0, "could not set SPEED value in table!");
+
+ ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, handle_ML_metric_bound,
+ "ML_ADM::ML_metric", "Radiation");
+
+ if (ierr < 0)
+ CCTK_WARN(0, "Failed to register Radiation BC for ML_ADM::ML_metric!");
+
+ }
+
if (CCTK_EQUALS(ML_curv_bound, "radiative"))
{
/* select radiation boundary condition */
@@ -290,25 +309,6 @@ extern "C" void ML_ADM_SelectBoundConds(CCTK_ARGUMENTS)
}
- if (CCTK_EQUALS(ML_metric_bound, "radiative"))
- {
- /* select radiation boundary condition */
- static CCTK_INT handle_ML_metric_bound CCTK_ATTRIBUTE_UNUSED = -1;
- if (handle_ML_metric_bound < 0) handle_ML_metric_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE);
- if (handle_ML_metric_bound < 0) CCTK_WARN(0, "could not create table!");
- if (Util_TableSetReal(handle_ML_metric_bound , ML_metric_bound_limit, "LIMIT") < 0)
- CCTK_WARN(0, "could not set LIMIT value in table!");
- if (Util_TableSetReal(handle_ML_metric_bound ,ML_metric_bound_speed, "SPEED") < 0)
- CCTK_WARN(0, "could not set SPEED value in table!");
-
- ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, handle_ML_metric_bound,
- "ML_ADM::ML_metric", "Radiation");
-
- if (ierr < 0)
- CCTK_WARN(0, "Failed to register Radiation BC for ML_ADM::ML_metric!");
-
- }
-
if (CCTK_EQUALS(ML_shift_bound, "radiative"))
{
/* select radiation boundary condition */
@@ -328,6 +328,120 @@ extern "C" void ML_ADM_SelectBoundConds(CCTK_ARGUMENTS)
}
+ if (CCTK_EQUALS(g11_bound, "radiative"))
+ {
+ /* select radiation boundary condition */
+ static CCTK_INT handle_g11_bound CCTK_ATTRIBUTE_UNUSED = -1;
+ if (handle_g11_bound < 0) handle_g11_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE);
+ if (handle_g11_bound < 0) CCTK_WARN(0, "could not create table!");
+ if (Util_TableSetReal(handle_g11_bound , g11_bound_limit, "LIMIT") < 0)
+ CCTK_WARN(0, "could not set LIMIT value in table!");
+ if (Util_TableSetReal(handle_g11_bound ,g11_bound_speed, "SPEED") < 0)
+ CCTK_WARN(0, "could not set SPEED value in table!");
+
+ ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_g11_bound,
+ "ML_ADM::g11", "Radiation");
+
+ if (ierr < 0)
+ CCTK_WARN(0, "Failed to register Radiation BC for ML_ADM::g11!");
+
+ }
+
+ if (CCTK_EQUALS(g12_bound, "radiative"))
+ {
+ /* select radiation boundary condition */
+ static CCTK_INT handle_g12_bound CCTK_ATTRIBUTE_UNUSED = -1;
+ if (handle_g12_bound < 0) handle_g12_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE);
+ if (handle_g12_bound < 0) CCTK_WARN(0, "could not create table!");
+ if (Util_TableSetReal(handle_g12_bound , g12_bound_limit, "LIMIT") < 0)
+ CCTK_WARN(0, "could not set LIMIT value in table!");
+ if (Util_TableSetReal(handle_g12_bound ,g12_bound_speed, "SPEED") < 0)
+ CCTK_WARN(0, "could not set SPEED value in table!");
+
+ ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_g12_bound,
+ "ML_ADM::g12", "Radiation");
+
+ if (ierr < 0)
+ CCTK_WARN(0, "Failed to register Radiation BC for ML_ADM::g12!");
+
+ }
+
+ if (CCTK_EQUALS(g13_bound, "radiative"))
+ {
+ /* select radiation boundary condition */
+ static CCTK_INT handle_g13_bound CCTK_ATTRIBUTE_UNUSED = -1;
+ if (handle_g13_bound < 0) handle_g13_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE);
+ if (handle_g13_bound < 0) CCTK_WARN(0, "could not create table!");
+ if (Util_TableSetReal(handle_g13_bound , g13_bound_limit, "LIMIT") < 0)
+ CCTK_WARN(0, "could not set LIMIT value in table!");
+ if (Util_TableSetReal(handle_g13_bound ,g13_bound_speed, "SPEED") < 0)
+ CCTK_WARN(0, "could not set SPEED value in table!");
+
+ ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_g13_bound,
+ "ML_ADM::g13", "Radiation");
+
+ if (ierr < 0)
+ CCTK_WARN(0, "Failed to register Radiation BC for ML_ADM::g13!");
+
+ }
+
+ if (CCTK_EQUALS(g22_bound, "radiative"))
+ {
+ /* select radiation boundary condition */
+ static CCTK_INT handle_g22_bound CCTK_ATTRIBUTE_UNUSED = -1;
+ if (handle_g22_bound < 0) handle_g22_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE);
+ if (handle_g22_bound < 0) CCTK_WARN(0, "could not create table!");
+ if (Util_TableSetReal(handle_g22_bound , g22_bound_limit, "LIMIT") < 0)
+ CCTK_WARN(0, "could not set LIMIT value in table!");
+ if (Util_TableSetReal(handle_g22_bound ,g22_bound_speed, "SPEED") < 0)
+ CCTK_WARN(0, "could not set SPEED value in table!");
+
+ ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_g22_bound,
+ "ML_ADM::g22", "Radiation");
+
+ if (ierr < 0)
+ CCTK_WARN(0, "Failed to register Radiation BC for ML_ADM::g22!");
+
+ }
+
+ if (CCTK_EQUALS(g23_bound, "radiative"))
+ {
+ /* select radiation boundary condition */
+ static CCTK_INT handle_g23_bound CCTK_ATTRIBUTE_UNUSED = -1;
+ if (handle_g23_bound < 0) handle_g23_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE);
+ if (handle_g23_bound < 0) CCTK_WARN(0, "could not create table!");
+ if (Util_TableSetReal(handle_g23_bound , g23_bound_limit, "LIMIT") < 0)
+ CCTK_WARN(0, "could not set LIMIT value in table!");
+ if (Util_TableSetReal(handle_g23_bound ,g23_bound_speed, "SPEED") < 0)
+ CCTK_WARN(0, "could not set SPEED value in table!");
+
+ ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_g23_bound,
+ "ML_ADM::g23", "Radiation");
+
+ if (ierr < 0)
+ CCTK_WARN(0, "Failed to register Radiation BC for ML_ADM::g23!");
+
+ }
+
+ if (CCTK_EQUALS(g33_bound, "radiative"))
+ {
+ /* select radiation boundary condition */
+ static CCTK_INT handle_g33_bound CCTK_ATTRIBUTE_UNUSED = -1;
+ if (handle_g33_bound < 0) handle_g33_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE);
+ if (handle_g33_bound < 0) CCTK_WARN(0, "could not create table!");
+ if (Util_TableSetReal(handle_g33_bound , g33_bound_limit, "LIMIT") < 0)
+ CCTK_WARN(0, "could not set LIMIT value in table!");
+ if (Util_TableSetReal(handle_g33_bound ,g33_bound_speed, "SPEED") < 0)
+ CCTK_WARN(0, "could not set SPEED value in table!");
+
+ ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_g33_bound,
+ "ML_ADM::g33", "Radiation");
+
+ if (ierr < 0)
+ CCTK_WARN(0, "Failed to register Radiation BC for ML_ADM::g33!");
+
+ }
+
if (CCTK_EQUALS(K11_bound, "radiative"))
{
/* select radiation boundary condition */
@@ -461,120 +575,6 @@ extern "C" void ML_ADM_SelectBoundConds(CCTK_ARGUMENTS)
}
- if (CCTK_EQUALS(g11_bound, "radiative"))
- {
- /* select radiation boundary condition */
- static CCTK_INT handle_g11_bound CCTK_ATTRIBUTE_UNUSED = -1;
- if (handle_g11_bound < 0) handle_g11_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE);
- if (handle_g11_bound < 0) CCTK_WARN(0, "could not create table!");
- if (Util_TableSetReal(handle_g11_bound , g11_bound_limit, "LIMIT") < 0)
- CCTK_WARN(0, "could not set LIMIT value in table!");
- if (Util_TableSetReal(handle_g11_bound ,g11_bound_speed, "SPEED") < 0)
- CCTK_WARN(0, "could not set SPEED value in table!");
-
- ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_g11_bound,
- "ML_ADM::g11", "Radiation");
-
- if (ierr < 0)
- CCTK_WARN(0, "Failed to register Radiation BC for ML_ADM::g11!");
-
- }
-
- if (CCTK_EQUALS(g12_bound, "radiative"))
- {
- /* select radiation boundary condition */
- static CCTK_INT handle_g12_bound CCTK_ATTRIBUTE_UNUSED = -1;
- if (handle_g12_bound < 0) handle_g12_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE);
- if (handle_g12_bound < 0) CCTK_WARN(0, "could not create table!");
- if (Util_TableSetReal(handle_g12_bound , g12_bound_limit, "LIMIT") < 0)
- CCTK_WARN(0, "could not set LIMIT value in table!");
- if (Util_TableSetReal(handle_g12_bound ,g12_bound_speed, "SPEED") < 0)
- CCTK_WARN(0, "could not set SPEED value in table!");
-
- ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_g12_bound,
- "ML_ADM::g12", "Radiation");
-
- if (ierr < 0)
- CCTK_WARN(0, "Failed to register Radiation BC for ML_ADM::g12!");
-
- }
-
- if (CCTK_EQUALS(g13_bound, "radiative"))
- {
- /* select radiation boundary condition */
- static CCTK_INT handle_g13_bound CCTK_ATTRIBUTE_UNUSED = -1;
- if (handle_g13_bound < 0) handle_g13_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE);
- if (handle_g13_bound < 0) CCTK_WARN(0, "could not create table!");
- if (Util_TableSetReal(handle_g13_bound , g13_bound_limit, "LIMIT") < 0)
- CCTK_WARN(0, "could not set LIMIT value in table!");
- if (Util_TableSetReal(handle_g13_bound ,g13_bound_speed, "SPEED") < 0)
- CCTK_WARN(0, "could not set SPEED value in table!");
-
- ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_g13_bound,
- "ML_ADM::g13", "Radiation");
-
- if (ierr < 0)
- CCTK_WARN(0, "Failed to register Radiation BC for ML_ADM::g13!");
-
- }
-
- if (CCTK_EQUALS(g22_bound, "radiative"))
- {
- /* select radiation boundary condition */
- static CCTK_INT handle_g22_bound CCTK_ATTRIBUTE_UNUSED = -1;
- if (handle_g22_bound < 0) handle_g22_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE);
- if (handle_g22_bound < 0) CCTK_WARN(0, "could not create table!");
- if (Util_TableSetReal(handle_g22_bound , g22_bound_limit, "LIMIT") < 0)
- CCTK_WARN(0, "could not set LIMIT value in table!");
- if (Util_TableSetReal(handle_g22_bound ,g22_bound_speed, "SPEED") < 0)
- CCTK_WARN(0, "could not set SPEED value in table!");
-
- ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_g22_bound,
- "ML_ADM::g22", "Radiation");
-
- if (ierr < 0)
- CCTK_WARN(0, "Failed to register Radiation BC for ML_ADM::g22!");
-
- }
-
- if (CCTK_EQUALS(g23_bound, "radiative"))
- {
- /* select radiation boundary condition */
- static CCTK_INT handle_g23_bound CCTK_ATTRIBUTE_UNUSED = -1;
- if (handle_g23_bound < 0) handle_g23_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE);
- if (handle_g23_bound < 0) CCTK_WARN(0, "could not create table!");
- if (Util_TableSetReal(handle_g23_bound , g23_bound_limit, "LIMIT") < 0)
- CCTK_WARN(0, "could not set LIMIT value in table!");
- if (Util_TableSetReal(handle_g23_bound ,g23_bound_speed, "SPEED") < 0)
- CCTK_WARN(0, "could not set SPEED value in table!");
-
- ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_g23_bound,
- "ML_ADM::g23", "Radiation");
-
- if (ierr < 0)
- CCTK_WARN(0, "Failed to register Radiation BC for ML_ADM::g23!");
-
- }
-
- if (CCTK_EQUALS(g33_bound, "radiative"))
- {
- /* select radiation boundary condition */
- static CCTK_INT handle_g33_bound CCTK_ATTRIBUTE_UNUSED = -1;
- if (handle_g33_bound < 0) handle_g33_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE);
- if (handle_g33_bound < 0) CCTK_WARN(0, "could not create table!");
- if (Util_TableSetReal(handle_g33_bound , g33_bound_limit, "LIMIT") < 0)
- CCTK_WARN(0, "could not set LIMIT value in table!");
- if (Util_TableSetReal(handle_g33_bound ,g33_bound_speed, "SPEED") < 0)
- CCTK_WARN(0, "could not set SPEED value in table!");
-
- ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_g33_bound,
- "ML_ADM::g33", "Radiation");
-
- if (ierr < 0)
- CCTK_WARN(0, "Failed to register Radiation BC for ML_ADM::g33!");
-
- }
-
if (CCTK_EQUALS(beta1_bound, "radiative"))
{
/* select radiation boundary condition */
@@ -632,6 +632,23 @@ extern "C" void ML_ADM_SelectBoundConds(CCTK_ARGUMENTS)
}
+ if (CCTK_EQUALS(ML_metric_bound, "scalar"))
+ {
+ /* select scalar boundary condition */
+ static CCTK_INT handle_ML_metric_bound CCTK_ATTRIBUTE_UNUSED = -1;
+ if (handle_ML_metric_bound < 0) handle_ML_metric_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE);
+ if (handle_ML_metric_bound < 0) CCTK_WARN(0, "could not create table!");
+ if (Util_TableSetReal(handle_ML_metric_bound ,ML_metric_bound_scalar, "SCALAR") < 0)
+ CCTK_WARN(0, "could not set SCALAR value in table!");
+
+ ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, handle_ML_metric_bound,
+ "ML_ADM::ML_metric", "scalar");
+
+ if (ierr < 0)
+ CCTK_WARN(0, "Failed to register Scalar BC for ML_ADM::ML_metric!");
+
+ }
+
if (CCTK_EQUALS(ML_curv_bound, "scalar"))
{
/* select scalar boundary condition */
@@ -666,23 +683,6 @@ extern "C" void ML_ADM_SelectBoundConds(CCTK_ARGUMENTS)
}
- if (CCTK_EQUALS(ML_metric_bound, "scalar"))
- {
- /* select scalar boundary condition */
- static CCTK_INT handle_ML_metric_bound CCTK_ATTRIBUTE_UNUSED = -1;
- if (handle_ML_metric_bound < 0) handle_ML_metric_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE);
- if (handle_ML_metric_bound < 0) CCTK_WARN(0, "could not create table!");
- if (Util_TableSetReal(handle_ML_metric_bound ,ML_metric_bound_scalar, "SCALAR") < 0)
- CCTK_WARN(0, "could not set SCALAR value in table!");
-
- ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, 1, handle_ML_metric_bound,
- "ML_ADM::ML_metric", "scalar");
-
- if (ierr < 0)
- CCTK_WARN(0, "Failed to register Scalar BC for ML_ADM::ML_metric!");
-
- }
-
if (CCTK_EQUALS(ML_shift_bound, "scalar"))
{
/* select scalar boundary condition */
@@ -700,6 +700,108 @@ extern "C" void ML_ADM_SelectBoundConds(CCTK_ARGUMENTS)
}
+ if (CCTK_EQUALS(g11_bound, "scalar"))
+ {
+ /* select scalar boundary condition */
+ static CCTK_INT handle_g11_bound CCTK_ATTRIBUTE_UNUSED = -1;
+ if (handle_g11_bound < 0) handle_g11_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE);
+ if (handle_g11_bound < 0) CCTK_WARN(0, "could not create table!");
+ if (Util_TableSetReal(handle_g11_bound ,g11_bound_scalar, "SCALAR") < 0)
+ CCTK_WARN(0, "could not set SCALAR value in table!");
+
+ ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_g11_bound,
+ "ML_ADM::g11", "scalar");
+
+ if (ierr < 0)
+ CCTK_WARN(0, "Error in registering Scalar BC for ML_ADM::g11!");
+
+ }
+
+ if (CCTK_EQUALS(g12_bound, "scalar"))
+ {
+ /* select scalar boundary condition */
+ static CCTK_INT handle_g12_bound CCTK_ATTRIBUTE_UNUSED = -1;
+ if (handle_g12_bound < 0) handle_g12_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE);
+ if (handle_g12_bound < 0) CCTK_WARN(0, "could not create table!");
+ if (Util_TableSetReal(handle_g12_bound ,g12_bound_scalar, "SCALAR") < 0)
+ CCTK_WARN(0, "could not set SCALAR value in table!");
+
+ ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_g12_bound,
+ "ML_ADM::g12", "scalar");
+
+ if (ierr < 0)
+ CCTK_WARN(0, "Error in registering Scalar BC for ML_ADM::g12!");
+
+ }
+
+ if (CCTK_EQUALS(g13_bound, "scalar"))
+ {
+ /* select scalar boundary condition */
+ static CCTK_INT handle_g13_bound CCTK_ATTRIBUTE_UNUSED = -1;
+ if (handle_g13_bound < 0) handle_g13_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE);
+ if (handle_g13_bound < 0) CCTK_WARN(0, "could not create table!");
+ if (Util_TableSetReal(handle_g13_bound ,g13_bound_scalar, "SCALAR") < 0)
+ CCTK_WARN(0, "could not set SCALAR value in table!");
+
+ ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_g13_bound,
+ "ML_ADM::g13", "scalar");
+
+ if (ierr < 0)
+ CCTK_WARN(0, "Error in registering Scalar BC for ML_ADM::g13!");
+
+ }
+
+ if (CCTK_EQUALS(g22_bound, "scalar"))
+ {
+ /* select scalar boundary condition */
+ static CCTK_INT handle_g22_bound CCTK_ATTRIBUTE_UNUSED = -1;
+ if (handle_g22_bound < 0) handle_g22_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE);
+ if (handle_g22_bound < 0) CCTK_WARN(0, "could not create table!");
+ if (Util_TableSetReal(handle_g22_bound ,g22_bound_scalar, "SCALAR") < 0)
+ CCTK_WARN(0, "could not set SCALAR value in table!");
+
+ ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_g22_bound,
+ "ML_ADM::g22", "scalar");
+
+ if (ierr < 0)
+ CCTK_WARN(0, "Error in registering Scalar BC for ML_ADM::g22!");
+
+ }
+
+ if (CCTK_EQUALS(g23_bound, "scalar"))
+ {
+ /* select scalar boundary condition */
+ static CCTK_INT handle_g23_bound CCTK_ATTRIBUTE_UNUSED = -1;
+ if (handle_g23_bound < 0) handle_g23_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE);
+ if (handle_g23_bound < 0) CCTK_WARN(0, "could not create table!");
+ if (Util_TableSetReal(handle_g23_bound ,g23_bound_scalar, "SCALAR") < 0)
+ CCTK_WARN(0, "could not set SCALAR value in table!");
+
+ ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_g23_bound,
+ "ML_ADM::g23", "scalar");
+
+ if (ierr < 0)
+ CCTK_WARN(0, "Error in registering Scalar BC for ML_ADM::g23!");
+
+ }
+
+ if (CCTK_EQUALS(g33_bound, "scalar"))
+ {
+ /* select scalar boundary condition */
+ static CCTK_INT handle_g33_bound CCTK_ATTRIBUTE_UNUSED = -1;
+ if (handle_g33_bound < 0) handle_g33_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE);
+ if (handle_g33_bound < 0) CCTK_WARN(0, "could not create table!");
+ if (Util_TableSetReal(handle_g33_bound ,g33_bound_scalar, "SCALAR") < 0)
+ CCTK_WARN(0, "could not set SCALAR value in table!");
+
+ ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_g33_bound,
+ "ML_ADM::g33", "scalar");
+
+ if (ierr < 0)
+ CCTK_WARN(0, "Error in registering Scalar BC for ML_ADM::g33!");
+
+ }
+
if (CCTK_EQUALS(K11_bound, "scalar"))
{
/* select scalar boundary condition */
@@ -819,108 +921,6 @@ extern "C" void ML_ADM_SelectBoundConds(CCTK_ARGUMENTS)
}
- if (CCTK_EQUALS(g11_bound, "scalar"))
- {
- /* select scalar boundary condition */
- static CCTK_INT handle_g11_bound CCTK_ATTRIBUTE_UNUSED = -1;
- if (handle_g11_bound < 0) handle_g11_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE);
- if (handle_g11_bound < 0) CCTK_WARN(0, "could not create table!");
- if (Util_TableSetReal(handle_g11_bound ,g11_bound_scalar, "SCALAR") < 0)
- CCTK_WARN(0, "could not set SCALAR value in table!");
-
- ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_g11_bound,
- "ML_ADM::g11", "scalar");
-
- if (ierr < 0)
- CCTK_WARN(0, "Error in registering Scalar BC for ML_ADM::g11!");
-
- }
-
- if (CCTK_EQUALS(g12_bound, "scalar"))
- {
- /* select scalar boundary condition */
- static CCTK_INT handle_g12_bound CCTK_ATTRIBUTE_UNUSED = -1;
- if (handle_g12_bound < 0) handle_g12_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE);
- if (handle_g12_bound < 0) CCTK_WARN(0, "could not create table!");
- if (Util_TableSetReal(handle_g12_bound ,g12_bound_scalar, "SCALAR") < 0)
- CCTK_WARN(0, "could not set SCALAR value in table!");
-
- ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_g12_bound,
- "ML_ADM::g12", "scalar");
-
- if (ierr < 0)
- CCTK_WARN(0, "Error in registering Scalar BC for ML_ADM::g12!");
-
- }
-
- if (CCTK_EQUALS(g13_bound, "scalar"))
- {
- /* select scalar boundary condition */
- static CCTK_INT handle_g13_bound CCTK_ATTRIBUTE_UNUSED = -1;
- if (handle_g13_bound < 0) handle_g13_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE);
- if (handle_g13_bound < 0) CCTK_WARN(0, "could not create table!");
- if (Util_TableSetReal(handle_g13_bound ,g13_bound_scalar, "SCALAR") < 0)
- CCTK_WARN(0, "could not set SCALAR value in table!");
-
- ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_g13_bound,
- "ML_ADM::g13", "scalar");
-
- if (ierr < 0)
- CCTK_WARN(0, "Error in registering Scalar BC for ML_ADM::g13!");
-
- }
-
- if (CCTK_EQUALS(g22_bound, "scalar"))
- {
- /* select scalar boundary condition */
- static CCTK_INT handle_g22_bound CCTK_ATTRIBUTE_UNUSED = -1;
- if (handle_g22_bound < 0) handle_g22_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE);
- if (handle_g22_bound < 0) CCTK_WARN(0, "could not create table!");
- if (Util_TableSetReal(handle_g22_bound ,g22_bound_scalar, "SCALAR") < 0)
- CCTK_WARN(0, "could not set SCALAR value in table!");
-
- ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_g22_bound,
- "ML_ADM::g22", "scalar");
-
- if (ierr < 0)
- CCTK_WARN(0, "Error in registering Scalar BC for ML_ADM::g22!");
-
- }
-
- if (CCTK_EQUALS(g23_bound, "scalar"))
- {
- /* select scalar boundary condition */
- static CCTK_INT handle_g23_bound CCTK_ATTRIBUTE_UNUSED = -1;
- if (handle_g23_bound < 0) handle_g23_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE);
- if (handle_g23_bound < 0) CCTK_WARN(0, "could not create table!");
- if (Util_TableSetReal(handle_g23_bound ,g23_bound_scalar, "SCALAR") < 0)
- CCTK_WARN(0, "could not set SCALAR value in table!");
-
- ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_g23_bound,
- "ML_ADM::g23", "scalar");
-
- if (ierr < 0)
- CCTK_WARN(0, "Error in registering Scalar BC for ML_ADM::g23!");
-
- }
-
- if (CCTK_EQUALS(g33_bound, "scalar"))
- {
- /* select scalar boundary condition */
- static CCTK_INT handle_g33_bound CCTK_ATTRIBUTE_UNUSED = -1;
- if (handle_g33_bound < 0) handle_g33_bound = Util_TableCreate(UTIL_TABLE_FLAGS_CASE_INSENSITIVE);
- if (handle_g33_bound < 0) CCTK_WARN(0, "could not create table!");
- if (Util_TableSetReal(handle_g33_bound ,g33_bound_scalar, "SCALAR") < 0)
- CCTK_WARN(0, "could not set SCALAR value in table!");
-
- ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, 1, handle_g33_bound,
- "ML_ADM::g33", "scalar");
-
- if (ierr < 0)
- CCTK_WARN(0, "Error in registering Scalar BC for ML_ADM::g33!");
-
- }
-
if (CCTK_EQUALS(beta1_bound, "scalar"))
{
/* select scalar boundary condition */
@@ -977,6 +977,11 @@ extern "C" void ML_ADM_SelectBoundConds(CCTK_ARGUMENTS)
/* template for entries in parameter file:
+#$bound$#ML_ADM::ML_metric_bound = "skip"
+#$bound$#ML_ADM::ML_metric_bound_speed = 1.0
+#$bound$#ML_ADM::ML_metric_bound_limit = 0.0
+#$bound$#ML_ADM::ML_metric_bound_scalar = 0.0
+
#$bound$#ML_ADM::ML_curv_bound = "skip"
#$bound$#ML_ADM::ML_curv_bound_speed = 1.0
#$bound$#ML_ADM::ML_curv_bound_limit = 0.0
@@ -987,16 +992,41 @@ extern "C" void ML_ADM_SelectBoundConds(CCTK_ARGUMENTS)
#$bound$#ML_ADM::ML_lapse_bound_limit = 0.0
#$bound$#ML_ADM::ML_lapse_bound_scalar = 0.0
-#$bound$#ML_ADM::ML_metric_bound = "skip"
-#$bound$#ML_ADM::ML_metric_bound_speed = 1.0
-#$bound$#ML_ADM::ML_metric_bound_limit = 0.0
-#$bound$#ML_ADM::ML_metric_bound_scalar = 0.0
-
#$bound$#ML_ADM::ML_shift_bound = "skip"
#$bound$#ML_ADM::ML_shift_bound_speed = 1.0
#$bound$#ML_ADM::ML_shift_bound_limit = 0.0
#$bound$#ML_ADM::ML_shift_bound_scalar = 0.0
+#$bound$#ML_ADM::g11_bound = "skip"
+#$bound$#ML_ADM::g11_bound_speed = 1.0
+#$bound$#ML_ADM::g11_bound_limit = 0.0
+#$bound$#ML_ADM::g11_bound_scalar = 0.0
+
+#$bound$#ML_ADM::g12_bound = "skip"
+#$bound$#ML_ADM::g12_bound_speed = 1.0
+#$bound$#ML_ADM::g12_bound_limit = 0.0
+#$bound$#ML_ADM::g12_bound_scalar = 0.0
+
+#$bound$#ML_ADM::g13_bound = "skip"
+#$bound$#ML_ADM::g13_bound_speed = 1.0
+#$bound$#ML_ADM::g13_bound_limit = 0.0
+#$bound$#ML_ADM::g13_bound_scalar = 0.0
+
+#$bound$#ML_ADM::g22_bound = "skip"
+#$bound$#ML_ADM::g22_bound_speed = 1.0
+#$bound$#ML_ADM::g22_bound_limit = 0.0
+#$bound$#ML_ADM::g22_bound_scalar = 0.0
+
+#$bound$#ML_ADM::g23_bound = "skip"
+#$bound$#ML_ADM::g23_bound_speed = 1.0
+#$bound$#ML_ADM::g23_bound_limit = 0.0
+#$bound$#ML_ADM::g23_bound_scalar = 0.0
+
+#$bound$#ML_ADM::g33_bound = "skip"
+#$bound$#ML_ADM::g33_bound_speed = 1.0
+#$bound$#ML_ADM::g33_bound_limit = 0.0
+#$bound$#ML_ADM::g33_bound_scalar = 0.0
+
#$bound$#ML_ADM::K11_bound = "skip"
#$bound$#ML_ADM::K11_bound_speed = 1.0
#$bound$#ML_ADM::K11_bound_limit = 0.0
@@ -1032,36 +1062,6 @@ extern "C" void ML_ADM_SelectBoundConds(CCTK_ARGUMENTS)
#$bound$#ML_ADM::alpha_bound_limit = 0.0
#$bound$#ML_ADM::alpha_bound_scalar = 0.0
-#$bound$#ML_ADM::g11_bound = "skip"
-#$bound$#ML_ADM::g11_bound_speed = 1.0
-#$bound$#ML_ADM::g11_bound_limit = 0.0
-#$bound$#ML_ADM::g11_bound_scalar = 0.0
-
-#$bound$#ML_ADM::g12_bound = "skip"
-#$bound$#ML_ADM::g12_bound_speed = 1.0
-#$bound$#ML_ADM::g12_bound_limit = 0.0
-#$bound$#ML_ADM::g12_bound_scalar = 0.0
-
-#$bound$#ML_ADM::g13_bound = "skip"
-#$bound$#ML_ADM::g13_bound_speed = 1.0
-#$bound$#ML_ADM::g13_bound_limit = 0.0
-#$bound$#ML_ADM::g13_bound_scalar = 0.0
-
-#$bound$#ML_ADM::g22_bound = "skip"
-#$bound$#ML_ADM::g22_bound_speed = 1.0
-#$bound$#ML_ADM::g22_bound_limit = 0.0
-#$bound$#ML_ADM::g22_bound_scalar = 0.0
-
-#$bound$#ML_ADM::g23_bound = "skip"
-#$bound$#ML_ADM::g23_bound_speed = 1.0
-#$bound$#ML_ADM::g23_bound_limit = 0.0
-#$bound$#ML_ADM::g23_bound_scalar = 0.0
-
-#$bound$#ML_ADM::g33_bound = "skip"
-#$bound$#ML_ADM::g33_bound_speed = 1.0
-#$bound$#ML_ADM::g33_bound_limit = 0.0
-#$bound$#ML_ADM::g33_bound_scalar = 0.0
-
#$bound$#ML_ADM::beta1_bound = "skip"
#$bound$#ML_ADM::beta1_bound_speed = 1.0
#$bound$#ML_ADM::beta1_bound_limit = 0.0
diff --git a/ML_ADM/src/RegisterSymmetries.cc b/ML_ADM/src/RegisterSymmetries.cc
index 0d05f97..0f972fb 100644
--- a/ML_ADM/src/RegisterSymmetries.cc
+++ b/ML_ADM/src/RegisterSymmetries.cc
@@ -19,67 +19,67 @@ extern "C" void ML_ADM_RegisterSymmetries(CCTK_ARGUMENTS)
sym[0] = 1;
sym[1] = 1;
sym[2] = 1;
- SetCartSymVN(cctkGH, sym, "ML_ADM::K11");
+ SetCartSymVN(cctkGH, sym, "ML_ADM::g11");
sym[0] = -1;
sym[1] = -1;
sym[2] = 1;
- SetCartSymVN(cctkGH, sym, "ML_ADM::K12");
+ SetCartSymVN(cctkGH, sym, "ML_ADM::g12");
sym[0] = -1;
sym[1] = 1;
sym[2] = -1;
- SetCartSymVN(cctkGH, sym, "ML_ADM::K13");
+ SetCartSymVN(cctkGH, sym, "ML_ADM::g13");
sym[0] = 1;
sym[1] = 1;
sym[2] = 1;
- SetCartSymVN(cctkGH, sym, "ML_ADM::K22");
+ SetCartSymVN(cctkGH, sym, "ML_ADM::g22");
sym[0] = 1;
sym[1] = -1;
sym[2] = -1;
- SetCartSymVN(cctkGH, sym, "ML_ADM::K23");
-
- sym[0] = 1;
- sym[1] = 1;
- sym[2] = 1;
- SetCartSymVN(cctkGH, sym, "ML_ADM::K33");
+ SetCartSymVN(cctkGH, sym, "ML_ADM::g23");
sym[0] = 1;
sym[1] = 1;
sym[2] = 1;
- SetCartSymVN(cctkGH, sym, "ML_ADM::alpha");
+ SetCartSymVN(cctkGH, sym, "ML_ADM::g33");
sym[0] = 1;
sym[1] = 1;
sym[2] = 1;
- SetCartSymVN(cctkGH, sym, "ML_ADM::g11");
+ SetCartSymVN(cctkGH, sym, "ML_ADM::K11");
sym[0] = -1;
sym[1] = -1;
sym[2] = 1;
- SetCartSymVN(cctkGH, sym, "ML_ADM::g12");
+ SetCartSymVN(cctkGH, sym, "ML_ADM::K12");
sym[0] = -1;
sym[1] = 1;
sym[2] = -1;
- SetCartSymVN(cctkGH, sym, "ML_ADM::g13");
+ SetCartSymVN(cctkGH, sym, "ML_ADM::K13");
sym[0] = 1;
sym[1] = 1;
sym[2] = 1;
- SetCartSymVN(cctkGH, sym, "ML_ADM::g22");
+ SetCartSymVN(cctkGH, sym, "ML_ADM::K22");
sym[0] = 1;
sym[1] = -1;
sym[2] = -1;
- SetCartSymVN(cctkGH, sym, "ML_ADM::g23");
+ SetCartSymVN(cctkGH, sym, "ML_ADM::K23");
sym[0] = 1;
sym[1] = 1;
sym[2] = 1;
- SetCartSymVN(cctkGH, sym, "ML_ADM::g33");
+ SetCartSymVN(cctkGH, sym, "ML_ADM::K33");
+
+ sym[0] = 1;
+ sym[1] = 1;
+ sym[2] = 1;
+ SetCartSymVN(cctkGH, sym, "ML_ADM::alpha");
sym[0] = -1;
sym[1] = 1;