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

style(polySolver): removed redundant lines

Browse Source
This commit is contained in:
Emily Boudreaux
2025-05-11 14:37:26 -04:00
parent d78f09623c
commit 454d49c3d3
1 changed files with 7 additions and 8 deletions
Download Patch File Download Diff File Expand all files Collapse all files

15
src/poly/solver/private/polySolver.cpp
View File

@@ -193,8 +193,7 @@ void PolySolver::solve() const {
state_vector.GetBlock(0) = thetaVec; // NOLINT(*-slicing) state_vector.GetBlock(0) = thetaVec; // NOLINT(*-slicing)
state_vector.GetBlock(1) = phiVec; // NOLINT(*-slicing) state_vector.GetBlock(1) = phiVec; // NOLINT(*-slicing)
mfem::Vector zero_rhs(block_offsets.Last());
mfem::Vector zero_rhs(block_offsets.Last()); // TODO: Is this right?
zero_rhs = 0.0; zero_rhs = 0.0;
const solverBundle sb = setupNewtonSolver(); const solverBundle sb = setupNewtonSolver();
@@ -271,7 +270,7 @@ void PolySolver::setInitialGuess() const {
// return (-1.0/radius) * r + 1; // return (-1.0/radius) * r + 1;
// return -std::pow((u*r), 2)+1.0; // The series expansion is a better guess; however, this is cheaper and ensures that the value at the surface is very close to zero in a way that the series expansion does not // return -std::pow((u*r), 2)+1.0; // The series expansion is a better guess; however, this is cheaper and ensures that the value at the surface is very close to zero in a way that the series expansion does not
return laneEmden::thetaSeriesExpansion(r, m_polytropicIndex, 10); return laneEmden::thetaSeriesExpansion(r, m_polytropicIndex, 10) + 3.0;
} }
); );
@@ -326,14 +325,14 @@ void PolySolver::saveAndViewSolution(const mfem::BlockVector& state_vector) cons
// --- Extract the Solution --- // --- Extract the Solution ---
if (m_config.get<bool>("Poly:Output:1D:Save", true)) { if (m_config.get<bool>("Poly:Output:1D:Save", true)) {
auto solutionPath = m_config.get<std::string>("Poly:Output:1D:Path", "polytropeSolution_1D.csv"); const auto solutionPath = m_config.get<std::string>("Poly:Output:1D:Path", "polytropeSolution_1D.csv");
auto derivSolPath = "d" + solutionPath; auto derivSolPath = "d" + solutionPath;
auto rayCoLatitude = m_config.get<double>("Poly:Output:1D:RayCoLatitude", 0.0); const auto rayCoLatitude = m_config.get<double>("Poly:Output:1D:RayCoLatitude", 0.0);
auto rayLongitude = m_config.get<double>("Poly:Output:1D:RayLongitude", 0.0); const auto rayLongitude = m_config.get<double>("Poly:Output:1D:RayLongitude", 0.0);
auto raySamples = m_config.get<int>("Poly:Output:1D:RaySamples", 100); const auto raySamples = m_config.get<int>("Poly:Output:1D:RaySamples", 100);
std::vector rayDirection = {rayCoLatitude, rayLongitude}; const std::vector rayDirection = {rayCoLatitude, rayLongitude};
Probe::getRaySolution(x_theta, *m_feTheta, rayDirection, raySamples, solutionPath); Probe::getRaySolution(x_theta, *m_feTheta, rayDirection, raySamples, solutionPath);
// Probe::getRaySolution(x_phi, *m_fePhi, rayDirection, raySamples, derivSolPath); // Probe::getRaySolution(x_phi, *m_fePhi, rayDirection, raySamples, derivSolPath);