diff --git a/test/integration/analytical_policy.hpp b/test/integration/analytical_policy.hpp
index 26edd267c..a0aad9a39 100644
--- a/test/integration/analytical_policy.hpp
+++ b/test/integration/analytical_policy.hpp
@@ -1780,12 +1780,6 @@ void test_analytical_oregonator(auto& solver, double tolerance = 1e-8)
 
   for (size_t i = 0; i < model_concentrations.size(); ++i)
   {
-    // EXPECT_NEAR(model_concentrations[i][_a], analytical_concentrations[i][0], tolerance)
-    //     << "Arrays differ at index (" << i << ", " << 0 << ")";
-    // EXPECT_NEAR(model_concentrations[i][_b], analytical_concentrations[i][1], tolerance)
-    //     << "Arrays differ at index (" << i << ", " << 1 << ")";
-    // EXPECT_NEAR(model_concentrations[i][_c], analytical_concentrations[i][2], tolerance)
-    //     << "Arrays differ at index (" << i << ", " << 2 << ")";
     double rel_diff = relative_difference(model_concentrations[i][_a], analytical_concentrations[i][0]);
     EXPECT_NEAR(0, rel_diff, tolerance) << "Arrays differ at index (" << i << ", " << 0 << ")";
     rel_diff = relative_difference(model_concentrations[i][_b], analytical_concentrations[i][1]);
diff --git a/test/integration/analytical_rosenbrock.cpp b/test/integration/analytical_rosenbrock.cpp
index e792b5bc9..ef4157fb5 100644
--- a/test/integration/analytical_rosenbrock.cpp
+++ b/test/integration/analytical_rosenbrock.cpp
@@ -113,7 +113,11 @@ TEST(AnalyticalExamples, Oregonator)
   using OregonatorTest = Oregonator<micm::Matrix<double>, SparseMatrixTest>;
 
   auto rosenbrock_solver = OregonatorTest::template CreateSolver<RosenbrockTest<OregonatorTest>, LinearSolverTest>(params, 1);
-  auto backward_euler_solver = OregonatorTest::template CreateSolver<BackwardEulerTest<OregonatorTest>, LinearSolverTest>(micm::BackwardEulerSolverParameters(), 1);
+  auto backward_euler_params = micm::BackwardEulerSolverParameters();
+  double abolute_tolerance = 1e-16;
+  backward_euler_params.absolute_tolerance_ = { abolute_tolerance, abolute_tolerance, abolute_tolerance };
+  backward_euler_params.relative_tolerance_ = 1e-10;
+  auto backward_euler_solver = OregonatorTest::template CreateSolver<BackwardEulerTest<OregonatorTest>, LinearSolverTest>(backward_euler_params, 1);
 
   test_analytical_oregonator(rosenbrock_solver, 1e-3);
   // test_analytical_oregonator(backward_euler_solver, 1e-1);