shares:time

shares:einstein

EXTENDS KEYWORD evolution_system ""
{
   "Conf_ADM" :: "Evolve using Conformal ADM equation in Conf_Hyp"
}

EXTENDS KEYWORD initial_data ""
{
"conf_adm_test1"  :: "a test metric, extrinsic curvature, and lapse "
}

EXTENDS KEYWORD slicing ""
{
   "ch_harmonic"   :: "harmonic slicing inside Conf_Hyp"
   "ch_harmonic2"  :: "harmonic2 slicing inside Conf_Hyp"
   "ch_oneplog"    :: "1+log slicing"
} ""

shares:gen_rel_perfect_fluid

#USES BOOLEAN mahc_evolution ""
#{
#}

shares: mahc_boundary

USES KEYWORD mahc_bound ""
{
}  ""


restricted:

KEYWORD conf_adm_method "Method for evolving Conf_ADM"
{
   "icn"      :: "Evolve Conf_ADM using Iterative Crank-Nicholson"
   "none"     :: "Do not use any method for conf_adm"
} "none"


private:

BOOLEAN Conf_setinit "should we set initial data for conf_adm?"
{
} "no"

BOOLEAN Conf_Storage "should we set storage for conf_adm?"
{
} "no"

BOOLEAN Conf_ADM_verbose "should Conf_ADM be loud?"
{
} "no"

BOOLEAN ch_evolveG "should we evolve the gammas?"
{
} "yes"

BOOLEAN ch_RfromH_A_eq  "should we use the hamiltonian constraint in the A_ij equations?"
{
} "yes"

BOOLEAN ch_RfromH_K_eq "should we use the hamiltonian constraint in the trK equation?"
{
} "yes"

BOOLEAN ch_subtract_trA "should we enforce A to be traceless?"
{
} "yes"

BOOLEAN ch_trAfromD "should we use the momentum constraint in the gamma evolution equation?"
{
} "yes"

BOOLEAN ch_forceKzero "should we force the trace of the extrinsic curvature to be 0?"
{
} "no"

STRING ch_maxsolve_method "method used by maximal equation solve"
{
   "bam" :: "use BAM to solve maximal equation"
} "bam"

BOOLEAN ch_RfromH_maximal  "should we use the hamiltonian constraint for scalar curv in maximal eq?"
{
} "yes"

KEYWORD ch_initlapse "how to set the initial lapse"
{
   "none"    :: "do not do anything to the lapse"
   "maximal" :: "solve the maximal equation for the lapse"
} "none"

REAL ch_maxsolve_tol "tolerance for maximal solve in Conf_ADM"
{
 : :: ""
} 1e-7

REAL ch_oneplog_const "constant in front of d_t alp slice cond."
{
 : :: ""
} 2.0

INT Conf_ICN_itnum "number of iterations for ICN"
{
 : :: ""
} 3

INT Conf_maxsolve_every "period (in timesteps) between maximal solves"
{
 : :: ""
} -1