tboudreaux/SERiF
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

test(tests/poly): removed default constructor test

Browse Source
This commit is contained in:
Emily Boudreaux
2025-03-05 12:57:35 -05:00
parent ea2661dbbb
commit b0252927e0
1 changed files with 16 additions and 18 deletions
Download Patch File Download Diff File Expand all files Collapse all files

34
tests/poly/polyTest.cpp
View File

@@ -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.");
}