From 3cb4af668f8aae3585fad040ee5afcb7254f5fb5 Mon Sep 17 00:00:00 2001 From: Erik Schnetter Date: Mon, 5 Nov 2007 11:21:08 -0600 Subject: Internal changes, add debug statements Reverse order of canonical indices of symmetric tensors. Add debug statements for BSSN constraints. --- par/ks_bssn_1.par | 75 +++++++++++++++++++++++++++++++++++++++++++++++-------- par/ks_bssn_2.par | 65 +++++++++++++++++++++++++++++++++++++++++------ 2 files changed, 123 insertions(+), 17 deletions(-) (limited to 'par') diff --git a/par/ks_bssn_1.par b/par/ks_bssn_1.par index 38c2850..8a13b48 100644 --- a/par/ks_bssn_1.par +++ b/par/ks_bssn_1.par @@ -1,12 +1,41 @@ ActiveThorns = " - Boundary CartGrid3D CoordBase IOUtil InitBase MoL SymBase - ADMBase ADMCoupling CoordGauge StaticConformal + Boundary + CartGrid3D + CoordBase + IOUtil + InitBase + MoL + SymBase + + ADMAnalysis + ADMBase + ADMConstraints + ADMCoupling + ADMMacros + CoordGauge + SpaceMask + StaticConformal + NaNChecker - Carpet CarpetIOASCII CarpetIOBasic CarpetIOScalar CarpetLib CarpetReduce LoopControl + + Carpet + CarpetIOASCII + CarpetIOBasic + CarpetIOScalar + CarpetLib + CarpetReduce + LoopControl + ReflectionSymmetry - Exact Formaline GSL + + Exact + Formaline + GSL + Slab + GenericFD + ML_BSSN " @@ -25,9 +54,12 @@ CoordBase::dx = 0.2 CoordBase::dy = 0.2 CoordBase::dz = 0.2 -ReflectionSymmetry::reflection_x = yes -ReflectionSymmetry::reflection_y = yes -ReflectionSymmetry::reflection_z = yes +ReflectionSymmetry::reflection_x = yes +ReflectionSymmetry::reflection_y = yes +ReflectionSymmetry::reflection_z = yes +ReflectionSymmetry::avoid_origin_x = no +ReflectionSymmetry::avoid_origin_y = no +ReflectionSymmetry::avoid_origin_z = no CoordBase::boundary_shiftout_x_lower = 1 CoordBase::boundary_shiftout_y_lower = 1 @@ -52,15 +84,20 @@ ADMBase::initial_shift = "exact" ADMBase::initial_dtlapse = "zero" ADMBase::initial_dtshift = "zero" -Exact::exact_model = "Kerr/Kerr-Schild" -Exact::Kerr_KerrSchild__mass = 1.0 -Exact::Kerr_KerrSchild__spin = 0.0 +#Exact::exact_model = "Kerr/Kerr-Schild" +#Exact::Kerr_KerrSchild__mass = 1.0 +#Exact::Kerr_KerrSchild__spin = 0.0 + +Exact::exact_model = "Schwarzschild/BL" +Exact::Schwarzschild_BL__mass = 1.0 ML_BSSN::harmonicN = 0 # 1+log ML_BSSN::harmonicF = 2.0 # 1+log ML_BSSN::ShiftGammaCoeff = 0.75 ML_BSSN::BetaDriver = 0.5 +ADMMacros::spatial_order = 4 + IO::out_dir = $parfile @@ -78,6 +115,15 @@ IOScalar::outScalar_vars = " ADMBase::curv ADMBase::lapse ADMBase::shift + ADMAnalysis::trace_of_K + ADMAnalysis::ricci_scalar + ADMConstraints::Hamiltonian + ADMConstraints::momentum + ML_BSSN::log_confac + ML_BSSN::metric + ML_BSSN::trace_curv + ML_BSSN::curv + ML_BSSN::Gamma ML_BSSN::Ham ML_BSSN::mom ML_BSSN::cons_detg @@ -97,6 +143,15 @@ IOASCII::out1D_vars = " ADMBase::curv ADMBase::lapse ADMBase::shift + ADMAnalysis::trace_of_K + ADMAnalysis::ricci_scalar + ADMConstraints::Hamiltonian + ADMConstraints::momentum + ML_BSSN::log_confac + ML_BSSN::metric + ML_BSSN::trace_curv + ML_BSSN::curv + ML_BSSN::Gamma ML_BSSN::Ham ML_BSSN::mom ML_BSSN::cons_detg diff --git a/par/ks_bssn_2.par b/par/ks_bssn_2.par index 957c1bd..ac12b53 100644 --- a/par/ks_bssn_2.par +++ b/par/ks_bssn_2.par @@ -1,12 +1,40 @@ ActiveThorns = " - Boundary CartGrid3D CoordBase IOUtil InitBase MoL SymBase - ADMBase ADMCoupling CoordGauge StaticConformal + Boundary + CartGrid3D + CoordBase + IOUtil + InitBase + MoL + SymBase + + ADMBase + ADMConstraints + ADMCoupling + ADMMacros + CoordGauge + SpaceMask + StaticConformal + NaNChecker - Carpet CarpetIOASCII CarpetIOBasic CarpetIOScalar CarpetLib CarpetReduce LoopControl + + Carpet + CarpetIOASCII + CarpetIOBasic + CarpetIOScalar + CarpetLib + CarpetReduce + LoopControl + ReflectionSymmetry - Exact Formaline GSL + + Exact + Formaline + GSL + Slab + GenericFD + ML_BSSN " @@ -25,9 +53,12 @@ CoordBase::dx = 0.1 CoordBase::dy = 0.1 CoordBase::dz = 0.1 -ReflectionSymmetry::reflection_x = yes -ReflectionSymmetry::reflection_y = yes -ReflectionSymmetry::reflection_z = yes +ReflectionSymmetry::reflection_x = yes +ReflectionSymmetry::reflection_y = yes +ReflectionSymmetry::reflection_z = yes +ReflectionSymmetry::avoid_origin_x = no +ReflectionSymmetry::avoid_origin_y = no +ReflectionSymmetry::avoid_origin_z = no CoordBase::boundary_shiftout_x_lower = 1 CoordBase::boundary_shiftout_y_lower = 1 @@ -61,6 +92,8 @@ ML_BSSN::harmonicF = 2.0 # 1+log ML_BSSN::ShiftGammaCoeff = 0.75 ML_BSSN::BetaDriver = 0.5 +ADMMacros::spatial_order = 4 + IO::out_dir = $parfile @@ -78,6 +111,15 @@ IOScalar::outScalar_vars = " ADMBase::curv ADMBase::lapse ADMBase::shift + ADMAnalysis::trace_of_K + ADMAnalysis::ricci_scalar + ADMConstraints::Hamiltonian + ADMConstraints::momentum + ML_BSSN::log_confac + ML_BSSN::metric + ML_BSSN::trace_curv + ML_BSSN::curv + ML_BSSN::Gamma ML_BSSN::Ham ML_BSSN::mom ML_BSSN::cons_detg @@ -97,6 +139,15 @@ IOASCII::out1D_vars = " ADMBase::curv ADMBase::lapse ADMBase::shift + ADMAnalysis::trace_of_K + ADMAnalysis::ricci_scalar + ADMConstraints::Hamiltonian + ADMConstraints::momentum + ML_BSSN::log_confac + ML_BSSN::metric + ML_BSSN::trace_curv + ML_BSSN::curv + ML_BSSN::Gamma ML_BSSN::Ham ML_BSSN::mom ML_BSSN::cons_detg -- cgit v1.2.3