From 0cc8e5bdfa48b3fb476d61fb4056d4a0e185ec84 Mon Sep 17 00:00:00 2001 From: Emily Boudreaux Date: Fri, 18 Apr 2025 11:18:55 -0400 Subject: [PATCH] feat(polyInitialGuess): use numeric derivitive for phi initial guess --- src/poly/solver/private/polySolver.cpp | 10 +--------- 1 file changed, 1 insertion(+), 9 deletions(-) diff --git a/src/poly/solver/private/polySolver.cpp b/src/poly/solver/private/polySolver.cpp index 99ea033..d3e70cd 100644 --- a/src/poly/solver/private/polySolver.cpp +++ b/src/poly/solver/private/polySolver.cpp @@ -270,16 +270,8 @@ void PolySolver::setInitialGuess() const { return -std::pow((u*r), 2)+1.0; } ); - mfem::VectorFunctionCoefficient phiInitGuess (m_mesh->SpaceDimension(), - [this](const mfem::Vector &x, mfem::Vector &v) { - double radius = Probe::getMeshRadius(*m_mesh); - double u = -1/std::pow(radius,2); - v(0) = 2*x(0)*u; - v(1) = 2*x(1)*u; - v(2) = 2*x(2)*u; - } - ); m_theta->ProjectCoefficient(thetaInitGuess); + mfem::GradientGridFunctionCoefficient phiInitGuess (m_theta.get()); m_phi->ProjectCoefficient(phiInitGuess); if (m_config.get("Poly:Solver:ViewInitialGuess", false)) { Probe::glVisView(*m_theta, *m_mesh, "θ init");