diff --git a/tests/poly/polyTest.cpp b/tests/poly/polyTest.cpp
index 3d4a11d..92cf25f 100644
--- a/tests/poly/polyTest.cpp
+++ b/tests/poly/polyTest.cpp
@@ -1,40 +1,38 @@
 #include <gtest/gtest.h>
 
+#include<numbers>
+
 #include "quill/LogMacros.h"
 #include "mfem.hpp"
 
 #include "polySolver.h"
 #include "probe.h"
 #include "config.h"
+#include "meshIO.h"
 
 std::string CONFIG_FILENAME = std::string(getenv("MESON_SOURCE_ROOT")) + "/tests/testsConfig.yaml";
+std::string SPHERICAL_MESH = std::string(getenv("MESON_SOURCE_ROOT")) + "/src/resources/mesh/core.msh";
 
 class polyTest : public ::testing::Test {};
 
-TEST_F(polyTest, DefaultConstructor) {
-    Config& config = Config::getInstance();
-    config.loadConfig(CONFIG_FILENAME);
-    EXPECT_NO_THROW(PolySolver(1, 1));
-}
-
 TEST_F(polyTest, Solve) {
     Config& config = Config::getInstance();
     config.loadConfig(CONFIG_FILENAME);
-    PolySolver polytrope(1, 1);
-    EXPECT_NO_THROW(polytrope.solve());
     Probe::LogManager& logManager = Probe::LogManager::getInstance();
-    quill::Logger* logger = logManager.newFileLogger("polyTest.log", "polyTest");
-    LOG_INFO(logger, "Solving polytrope with n = 1.5");
-    mfem::Mesh& mesh = polytrope.getMesh();
-    mfem::GridFunction& solution = polytrope.getSolution();
-    Probe::glVisView(solution, mesh, "Polytrope solution");
+    quill::Logger* logger = logManager.getLogger("log");
+    LOG_INFO(logger, "Starting polytrope solve test 1...");
+    config.loadConfig(CONFIG_FILENAME);
 
-    // Get the mesh radius
+    MeshIO meshIO(SPHERICAL_MESH, std::numbers::pi);
+    mfem::Mesh& mesh = meshIO.GetMesh();
     double radius = Probe::getMeshRadius(mesh);
     LOG_INFO(logger, "Mesh radius: {}", radius);
 
-    // Get the ray solution
-    const std::vector<double> rayDirection = {0, 0};
-    int numSamples = 10;
-    // std::vector<double> samples = Probe::getRaySolution(solution, mesh, rayDirection, numSamples);
+    double polytropicIndex = config.get<double>("Tests:Poly:Index", 1);
+    LOG_INFO(logger, "Solving polytrope with n = {:0.2f}", polytropicIndex);
+
+    PolySolver polytrope(polytropicIndex, 1, mesh);
+    LOG_INFO(logger, "Solving polytrope...");
+    EXPECT_NO_THROW(polytrope.solve());
+    LOG_INFO(logger, "Polytrope solved.");
 }
\ No newline at end of file