From 534a9cdf4be4159ca9fcae078e7762ab3f1686c2 Mon Sep 17 00:00:00 2001 From: Emily Boudreaux Date: Mon, 12 May 2025 14:27:41 -0400 Subject: [PATCH] feat(polySolver): added polytropic index as a constructor argument for polytrope operator --- src/poly/solver/private/polySolver.cpp | 5 +++-- 1 file changed, 3 insertions(+), 2 deletions(-) diff --git a/src/poly/solver/private/polySolver.cpp b/src/poly/solver/private/polySolver.cpp index 407748a..21dfc6b 100644 --- a/src/poly/solver/private/polySolver.cpp +++ b/src/poly/solver/private/polySolver.cpp @@ -119,7 +119,8 @@ void PolySolver::assembleBlockSystem() { std::move(forms->Q), std::move(forms->D), std::move(forms->f), - blockOffsets); + blockOffsets, + m_polytropicIndex); } mfem::Array PolySolver::computeBlockOffsets() const { @@ -270,7 +271,7 @@ void PolySolver::setInitialGuess() const { // 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 laneEmden::thetaSeriesExpansion(r, m_polytropicIndex, 10) + 3.0; + return laneEmden::thetaSeriesExpansion(r, m_polytropicIndex, 10); } );