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

fix(poly): increased default gamma and changed initial guess

Browse Source 
The default gamma value has been upped to 1e4 which is enough to strongly constrain the solution to have zero slope at the core region. Further, the initial guess has been changed from a series expansion of theta to a simple quadratic that is one at origin and zero at the polytrope radius. This is faster to evaluate and seems to work just as well.
This commit is contained in:
Emily Boudreaux
2025-03-26 12:34:30 -04:00
parent f56611b785
commit b2ddbc41e7
1 changed files with 6 additions and 6 deletions
Download Patch File Download Diff File Expand all files Collapse all files

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

@@ -104,7 +104,7 @@ void PolySolver::assembleNonlinearForm() {
compositeIntegrator->add_integrator(nonlinearIntegrator.release());
// Add the contraint term \gamma(\nabla \theta(0)\cdot\nabla v(0))^{2}
double gamma = config.get<double>("Poly:Solver:Constraint:Gamma", 1e2);
double gamma = config.get<double>("Poly:Solver:Constraint:Gamma", 1e4);
auto constraintIntegrator = std::make_unique<polyMFEMUtils::BilinearIntegratorWrapper>(
new polyMFEMUtils::ConstraintIntegrator(gamma, &mesh)
);
@@ -118,12 +118,12 @@ void PolySolver::solve(){
mfem::FunctionCoefficient initCoeff (
[this](const mfem::Vector &x) {
double r = x.Norml2();
double theta = laneEmden::thetaSerieseExpansion(r, n, 10);
return theta;
// double radius = Probe::getMeshRadius(mesh);
// double u = 1/radius;
// double theta = laneEmden::thetaSerieseExpansion(r, n, 10);
// return theta;
double radius = Probe::getMeshRadius(mesh);
double u = 1/radius;
// return -std::pow((u*r), 2)+1.0;
return -std::pow((u*r), 2)+1.0;
}
);
u->ProjectCoefficient(initCoeff);