From 1e90b806a47b7314602c5cecb19ad868f789e556 Mon Sep 17 00:00:00 2001
From: Remy Prechelt <prechelt@hawaii.edu>
Date: Mon, 29 Jun 2020 15:18:27 -1000
Subject: [PATCH] Add tests for SetMagneticField.
---
Environment/testEnvironment.cc | 26 ++++++++++++++++++++++++--
1 file changed, 24 insertions(+), 2 deletions(-)
diff --git a/Environment/testEnvironment.cc b/Environment/testEnvironment.cc
index d2434aac2..ae993ad13 100644
--- a/Environment/testEnvironment.cc
+++ b/Environment/testEnvironment.cc
@@ -246,7 +246,18 @@ TEST_CASE("UniformMagneticField w/ Homogeneous Medium") {
QuantityVector B0(0_T, 0_T, 0_T);
// create our atmospheric model
- AtmModel const medium(B0, 19.2_g / cube(1_cm), protonComposition);
+ AtmModel medium(B0, 19.2_g / cube(1_cm), protonComposition);
+
+ // create a new magnetic field vector
+ QuantityVector B1(23_T, 57_T, -4_T);
+
+ // and update this atmospheric model
+ medium.SetMagneticField(B1);
+
+ // and test at several locations
+ REQUIRE(B1 == medium.GetMagneticField(Point(gCS, -10_m, 4_m, 35_km)));
+ REQUIRE(B1 == medium.GetMagneticField(Point(gCS, 1000_km, -1000_km, 1000_km)));
+ REQUIRE(B1 == medium.GetMagneticField(Point(gCS, 0_m, 0_m, 0_m)));
}
TEST_CASE("UniformMagneticField w/ FlatExponential") {
@@ -270,8 +281,19 @@ TEST_CASE("UniformMagneticField w/ FlatExponential") {
QuantityVector B0(23_T, 57_T, -4_T);
// create our atmospheric model
- AtmModel const medium(B0, gOrigin, axis, rho0, lambda, protonComposition);
+ AtmModel medium(B0, gOrigin, axis, rho0, lambda, protonComposition);
// check that the returned magnetic field is correct
REQUIRE(B0 == medium.GetMagneticField(Point(gCS, -10_m, 4_m, 35_km)));
+
+ // create a new magnetic field vector
+ QuantityVector B1(23_T, 57_T, -4_T);
+
+ // and update this atmospheric model
+ medium.SetMagneticField(B1);
+
+ // and test at several locations
+ REQUIRE(B1 == medium.GetMagneticField(Point(gCS, -10_m, 4_m, 35_km)));
+ REQUIRE(B1 == medium.GetMagneticField(Point(gCS, 1000_km, -1000_km, 1000_km)));
+ REQUIRE(B1 == medium.GetMagneticField(Point(gCS, 0_m, 0_m, 0_m)));
}
--
GitLab