Added test for J_1 and increases accuracy required.

phbasler · Apr 25, 2024 · 2cb99e9 · 2cb99e9
1 parent dce589a
commit 2cb99e9
Showing 1 changed file with 15 additions and 3 deletions.
diff --git a/tests/unittests/Test-gw.cpp b/tests/unittests/Test-gw.cpp
@@ -26,9 +26,21 @@ TEST_CASE("Test I_\alpha", "[gw]")
 
   BounceActionInt BACalc;
 
-  REQUIRE(BACalc.BesselI(3, 1) == Approx(0.0221684249).epsilon(5e-2));
-  REQUIRE(BACalc.BesselI(1, 3) == Approx(3.953370217).epsilon(5e-2));
-  REQUIRE(BACalc.BesselI(1, 1.5) == Approx(0.9816664).epsilon(5e-2));
+  REQUIRE(BACalc.BesselI(3, 1) == Approx(0.0221684249).epsilon(1e-8));
+  REQUIRE(BACalc.BesselI(1, 3) == Approx(3.953370217).epsilon(1e-8));
+  REQUIRE(BACalc.BesselI(1, 1.5) == Approx(0.9816664285779).epsilon(1e-8));
+}
+
+TEST_CASE("Test J_1", "[gw]")
+{
+  // Tests bounce solver with analytical derivative
+  using namespace BSMPT;
+
+  BounceActionInt BACalc;
+
+  REQUIRE(BACalc.BesselJ(1) == Approx(0.440050585744).epsilon(1e-8));
+  REQUIRE(BACalc.BesselJ(3) == Approx(0.3390589585259).epsilon(1e-8));
+  REQUIRE(BACalc.BesselJ(1.5) == Approx(0.5579365079).epsilon(1e-8));
 }
 
 TEST_CASE("Solve bounce equation with analytical derivative", "[gw]")