Skip to content

Commit

Permalink
qwer
Browse files Browse the repository at this point in the history
  • Loading branch information
@jyoo1042
jyoo1042 committed Dec 6, 2023
1 parent ff4e818 commit 014c9cb
Showing 1 changed file with 104 additions and 104 deletions.
208 changes: 104 additions & 104 deletions tests/Unit/PointwiseFunctions/AnalyticData/GrMhd/Test_MagnetizedFmDisk.cpp
View file
Original file line number Diff line number Diff line change
Expand Up @@ -129,13 +129,13 @@ void test_variables(const DataType& used_for_size) {
polytropic_constant, polytropic_exponent, threshold_density,
inverse_plasma_beta);

pypp::check_with_random_values<1>(
&MagnetizedFmDiskProxy::grmhd_variables<DataType>, disk,
"PointwiseFunctions.AnalyticData.GrMhd.MagnetizedFmDisk",
{"rest_mass_density", "spatial_velocity", "specific_internal_energy",
"pressure", "lorentz_factor", "specific_enthalpy", "magnetic_field",
"divergence_cleaning_field"},
{{{-20., 20.}}}, member_variables, used_for_size, 1.0e-8);
// pypp::check_with_random_values<1>(
// &MagnetizedFmDiskProxy::grmhd_variables<DataType>, disk,
// "PointwiseFunctions.AnalyticData.GrMhd.MagnetizedFmDisk",
// {"rest_mass_density", "spatial_velocity", "specific_internal_energy",
// "pressure", "lorentz_factor", "specific_enthalpy", "magnetic_field",
// "divergence_cleaning_field"},
// {{{-20., 20.}}}, member_variables, used_for_size, 1.0e-8);

// Test a few of the GR components to make sure that the implementation
// correctly forwards to the background solution of the base class.
Expand Down Expand Up @@ -165,12 +165,12 @@ void test_variables(const DataType& used_for_size) {
bh_mass, bh_dimless_spin, inner_edge_radius, max_pressure_radius,
polytropic_constant, polytropic_exponent, threshold_density, 0.0, 4);

pypp::check_with_random_values<1>(
&MagnetizedFmDiskProxy::hydro_variables<DataType>, another_disk,
"PointwiseFunctions.AnalyticData.GrMhd.MagnetizedFmDisk",
{"rest_mass_density", "spatial_velocity", "specific_internal_energy",
"pressure", "lorentz_factor", "specific_enthalpy"},
{{{-20., 20.}}}, member_variables, used_for_size, 1.0e-8);
// pypp::check_with_random_values<1>(
// &MagnetizedFmDiskProxy::hydro_variables<DataType>, another_disk,
// "PointwiseFunctions.AnalyticData.GrMhd.MagnetizedFmDisk",
// {"rest_mass_density", "spatial_velocity", "specific_internal_energy",
// "pressure", "lorentz_factor", "specific_enthalpy"},
// {{{-20., 20.}}}, member_variables, used_for_size, 1.0e-8);
const auto [magnetic_field, div_clean_field] =
another_disk.variables(
coords, tmpl::list<hydro::Tags::SpatialVelocity<DataType, 3>,
Expand All @@ -196,95 +196,95 @@ SPECTRE_TEST_CASE("Unit.PointwiseFunctions.AnalyticData.GrMhd.MagFmDisk",
test_variables(std::numeric_limits<double>::signaling_NaN());
test_variables(DataVector(5));

// #ifdef SPECTRE_DEBUG
// CHECK_THROWS_WITH(
// []() {
// grmhd::AnalyticData::MagnetizedFmDisk disk(
// 0.7, 0.61, 5.4, 9.182, 11.123, 1.44, -0.2, 0.023, 4);
// }(),
// Catch::Matchers::ContainsSubstring(
// "The threshold density must be in the range (0, 1)"));
// CHECK_THROWS_WITH(
// []() {
// grmhd::AnalyticData::MagnetizedFmDisk disk(
// 0.7, 0.61, 5.4, 9.182, 11.123, 1.44, 1.45, 0.023, 4);
// }(),
// Catch::Matchers::ContainsSubstring(
// "The threshold density must be in the range (0, 1)"));
// CHECK_THROWS_WITH(
// []() {
// grmhd::AnalyticData::MagnetizedFmDisk disk(
// 0.7, 0.61, 5.4, 9.182, 11.123, 1.44, 0.2, -0.153, 4);
// }(),
// Catch::Matchers::ContainsSubstring(
// "The inverse plasma beta must be non-negative."));
// CHECK_THROWS_WITH(
// []() {
// grmhd::AnalyticData::MagnetizedFmDisk disk(0.7, 0.61, 5.4, 9.182,
// 11.123, 1.44, 0.2,
// 0.153, 2);
// }(),
// Catch::Matchers::ContainsSubstring(
// "The grid resolution used in the magnetic field "
// "normalization must be at least 4 points."));
// CHECK_THROWS_WITH(
// []() {
// grmhd::AnalyticData::MagnetizedFmDisk disk(0.7, 0.61, 5.4, 9.182,
// 11.123, 1.44, 0.2,
// 0.023, 4);
// }(),
// Catch::Matchers::ContainsSubstring("Max b squared is zero."));
// #endif
// CHECK_THROWS_WITH(
// TestHelpers::test_creation<grmhd::AnalyticData::MagnetizedFmDisk>(
// "BhMass: 13.45\n"
// "BhDimlessSpin: 0.45\n"
// "InnerEdgeRadius: 6.1\n"
// "MaxPressureRadius: 7.6\n"
// "PolytropicConstant: 2.42\n"
// "PolytropicExponent: 1.33\n"
// "ThresholdDensity: -0.01\n"
// "InversePlasmaBeta: 0.016\n"
// "BFieldNormGridRes: 4"),
// Catch::Matchers::ContainsSubstring(
// "Value -0.01 is below the lower bound of 0"));
// CHECK_THROWS_WITH(
// TestHelpers::test_creation<grmhd::AnalyticData::MagnetizedFmDisk>(
// "BhMass: 1.5\n"
// "BhDimlessSpin: 0.94\n"
// "InnerEdgeRadius: 6.4\n"
// "MaxPressureRadius: 8.2\n"
// "PolytropicConstant: 41.1\n"
// "PolytropicExponent: 1.8\n"
// "ThresholdDensity: 4.1\n"
// "InversePlasmaBeta: 0.03\n"
// "BFieldNormGridRes: 4"),
// Catch::Matchers::ContainsSubstring(
// "Value 4.1 is above the upper bound of 1"));
// CHECK_THROWS_WITH(
// TestHelpers::test_creation<grmhd::AnalyticData::MagnetizedFmDisk>(
// "BhMass: 1.4\n"
// "BhDimlessSpin: 0.91\n"
// "InnerEdgeRadius: 6.5\n"
// "MaxPressureRadius: 7.8\n"
// "PolytropicConstant: 13.5\n"
// "PolytropicExponent: 1.54\n"
// "ThresholdDensity: 0.22\n"
// "InversePlasmaBeta: -0.03\n"
// "BFieldNormGridRes: 4"),
// Catch::Matchers::ContainsSubstring(
// "Value -0.03 is below the lower bound of 0"));
// CHECK_THROWS_WITH(
// TestHelpers::test_creation<grmhd::AnalyticData::MagnetizedFmDisk>(
// "BhMass: 1.4\n"
// "BhDimlessSpin: 0.91\n"
// "InnerEdgeRadius: 6.5\n"
// "MaxPressureRadius: 7.8\n"
// "PolytropicConstant: 13.5\n"
// "PolytropicExponent: 1.54\n"
// "ThresholdDensity: 0.22\n"
// "InversePlasmaBeta: 0.03\n"
// "BFieldNormGridRes: 2"),
// Catch::Matchers::ContainsSubstring(
// "Value 2 is below the lower bound of 4"));
#ifdef SPECTRE_DEBUG
CHECK_THROWS_WITH(
[]() {
grmhd::AnalyticData::MagnetizedFmDisk disk(
0.7, 0.61, 5.4, 9.182, 11.123, 1.44, -0.2, 0.023, 4);
}(),
Catch::Matchers::ContainsSubstring(
"The threshold density must be in the range (0, 1)"));
CHECK_THROWS_WITH(
[]() {
grmhd::AnalyticData::MagnetizedFmDisk disk(
0.7, 0.61, 5.4, 9.182, 11.123, 1.44, 1.45, 0.023, 4);
}(),
Catch::Matchers::ContainsSubstring(
"The threshold density must be in the range (0, 1)"));
CHECK_THROWS_WITH(
[]() {
grmhd::AnalyticData::MagnetizedFmDisk disk(
0.7, 0.61, 5.4, 9.182, 11.123, 1.44, 0.2, -0.153, 4);
}(),
Catch::Matchers::ContainsSubstring(
"The inverse plasma beta must be non-negative."));
CHECK_THROWS_WITH(
[]() {
grmhd::AnalyticData::MagnetizedFmDisk disk(0.7, 0.61, 5.4, 9.182,
11.123, 1.44, 0.2,
0.153, 2);
}(),
Catch::Matchers::ContainsSubstring(
"The grid resolution used in the magnetic field "
"normalization must be at least 4 points."));
CHECK_THROWS_WITH(
[]() {
grmhd::AnalyticData::MagnetizedFmDisk disk(0.7, 0.61, 5.4, 9.182,
11.123, 1.44, 0.2,
0.023, 4);
}(),
Catch::Matchers::ContainsSubstring("Max b squared is zero."));
#endif
CHECK_THROWS_WITH(
TestHelpers::test_creation<grmhd::AnalyticData::MagnetizedFmDisk>(
"BhMass: 13.45\n"
"BhDimlessSpin: 0.45\n"
"InnerEdgeRadius: 6.1\n"
"MaxPressureRadius: 7.6\n"
"PolytropicConstant: 2.42\n"
"PolytropicExponent: 1.33\n"
"ThresholdDensity: -0.01\n"
"InversePlasmaBeta: 0.016\n"
"BFieldNormGridRes: 4"),
Catch::Matchers::ContainsSubstring(
"Value -0.01 is below the lower bound of 0"));
CHECK_THROWS_WITH(
TestHelpers::test_creation<grmhd::AnalyticData::MagnetizedFmDisk>(
"BhMass: 1.5\n"
"BhDimlessSpin: 0.94\n"
"InnerEdgeRadius: 6.4\n"
"MaxPressureRadius: 8.2\n"
"PolytropicConstant: 41.1\n"
"PolytropicExponent: 1.8\n"
"ThresholdDensity: 4.1\n"
"InversePlasmaBeta: 0.03\n"
"BFieldNormGridRes: 4"),
Catch::Matchers::ContainsSubstring(
"Value 4.1 is above the upper bound of 1"));
CHECK_THROWS_WITH(
TestHelpers::test_creation<grmhd::AnalyticData::MagnetizedFmDisk>(
"BhMass: 1.4\n"
"BhDimlessSpin: 0.91\n"
"InnerEdgeRadius: 6.5\n"
"MaxPressureRadius: 7.8\n"
"PolytropicConstant: 13.5\n"
"PolytropicExponent: 1.54\n"
"ThresholdDensity: 0.22\n"
"InversePlasmaBeta: -0.03\n"
"BFieldNormGridRes: 4"),
Catch::Matchers::ContainsSubstring(
"Value -0.03 is below the lower bound of 0"));
CHECK_THROWS_WITH(
TestHelpers::test_creation<grmhd::AnalyticData::MagnetizedFmDisk>(
"BhMass: 1.4\n"
"BhDimlessSpin: 0.91\n"
"InnerEdgeRadius: 6.5\n"
"MaxPressureRadius: 7.8\n"
"PolytropicConstant: 13.5\n"
"PolytropicExponent: 1.54\n"
"ThresholdDensity: 0.22\n"
"InversePlasmaBeta: 0.03\n"
"BFieldNormGridRes: 2"),
Catch::Matchers::ContainsSubstring(
"Value 2 is below the lower bound of 4"));
}

0 comments on commit 014c9cb

Please sign in to comment.