From 695441e597abcf7de5ccb20813e07121fc5fb243 Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Andr=C3=A9=20Schmidt?= <upwli@student.kit.edu>
Date: Wed, 17 Nov 2021 10:38:08 +0000
Subject: [PATCH] Update testMagneticField.cpp
---
tests/media/testMagneticField.cpp | 13 ++++++++++---
1 file changed, 10 insertions(+), 3 deletions(-)
diff --git a/tests/media/testMagneticField.cpp b/tests/media/testMagneticField.cpp
index 56742bca5..4ca395af8 100644
--- a/tests/media/testMagneticField.cpp
+++ b/tests/media/testMagneticField.cpp
@@ -14,6 +14,7 @@
#include <corsika/media/UniformMagneticField.hpp>
#include <corsika/media/IMagneticFieldModel.hpp>
#include <corsika/media/VolumeTreeNode.hpp>
+#include <corsika/media/WMM.hpp>
#include <SetupTestTrajectory.hpp>
#include <corsika/setup/SetupTrajectory.hpp>
@@ -85,7 +86,13 @@ TEST_CASE("UniformMagneticField w/ Homogeneous Medium") {
setup::testing::make_track<setup::Trajectory>(line, tEnd);
// create earth magnetic field vector
- MagneticFieldVector Earth_B = get_wmm(gCS, 2022.5, 100_km, -80, -120);
- CHECK(Earth_B.getComponents(gCS) == MagneticFieldVector{gCS, 5814.9658886215_nT,
- 14802.9663839328_nT, -49755.3119939183_nT}.getComponents(gCS));
+ MagneticFieldVector Earth_B_1 = get_wmm(gCS, 2022.5, 100_km, -80, -120);
+ MagneticFieldVector Earth_B_2 = get_wmm(gCS, 2022.5, 0_km, 0, 120);
+ MagneticFieldVector Earth_B_3 = get_wmm(gCS, 2020, 100_km, 80, 0);
+ CHECK(Earth_B_1.getComponents(gCS) ==
+ MagneticFieldVector{gCS, 5815_nT, 14803_nT, -49755.3_nT}.getComponents(gCS));
+ CHECK(Earth_B_2.getComponents(gCS) ==
+ MagneticFieldVector{gCS, 39684.7_nT, -42.2_nT, 10809.5_nT}.getComponents(gCS));
+ CHECK(Earth_B_3.getComponents(gCS) ==
+ MagneticFieldVector{gCS, 6570.4_nT, -146.3_nT, -54606_nT}.getComponents(gCS));
}
--
GitLab