From 0c2a5e3d9a4c021b099dd586e67567c390b50755 Mon Sep 17 00:00:00 2001 From: Emily Boudreaux Date: Fri, 14 Mar 2025 08:02:30 -0400 Subject: [PATCH] refactor(polyMFEMUtils): removed unused methods --- src/poly/utils/private/polyMFEMUtils.cpp | 58 +----------------------- 1 file changed, 1 insertion(+), 57 deletions(-) diff --git a/src/poly/utils/private/polyMFEMUtils.cpp b/src/poly/utils/private/polyMFEMUtils.cpp index 40231f7..a97ed6a 100644 --- a/src/poly/utils/private/polyMFEMUtils.cpp +++ b/src/poly/utils/private/polyMFEMUtils.cpp @@ -428,62 +428,6 @@ namespace polyMFEMUtils { LOG_INFO(logger, "Getting element dofs for zero slope constraint...done"); LOG_INFO(logger, "Building location of zero slope constraint...done"); } - - // void ZeroSlopeNewtonSolver::ProcessNewState(const mfem::Vector &x) const { - // LOG_INFO(logger, "Processing new state for zero slope constraint..."); - // if (zeroSlopeElemID < 0) { - // LOG_ERROR(logger, "ZeroSlopeNewtonSolver::ProcessNewState: zero slope element ID is not set"); - // MFEM_ABORT("ZeroSlopeNewtonSolver::ProcessNewState: zero slope element ID is not set"); - // } - - // mfem::NonlinearForm *nlf = dynamic_cast(const_cast(oper)); - // if (!nlf) { - // LOG_ERROR(logger, "ZeroSlopeNewtonSolver::ProcessNewState: input operator is not a NonlinearForm"); - // MFEM_ABORT("ZeroSlopeNewtonSolver::ProcessNewState: input operator is not a NonlinearForm"); - // } - - // mfem::FiniteElementSpace *fes = nlf->FESpace(); - // if (!fes) { - // LOG_ERROR(logger, "ZeroSlopeNewtonSolver::ProcessNewState: input operator does not have a finite element space"); - // MFEM_ABORT("ZeroSlopeNewtonSolver::ProcessNewState: input operator does not have a finite element space"); - // } - - // mfem::Mesh *mesh = fes->GetMesh(); - // if (!mesh) { - // LOG_ERROR(logger, "ZeroSlopeNewtonSolver::ProcessNewState: input operator does not have a mesh"); - // MFEM_ABORT("ZeroSlopeNewtonSolver::ProcessNewState: input operator does not have a mesh"); - // } - - // mfem::ElementTransformation *T = mesh->GetElementTransformation(zeroSlopeElemID); - // if (!T) { - // LOG_ERROR(logger, "ZeroSlopeNewtonSolver::ProcessNewState: element transformation is not found"); - // MFEM_ABORT("ZeroSlopeNewtonSolver::ProcessNewState: element transformation is not found"); - // } - - // mfem::Vector grad_u(3); - - // mfem::GridFunction u_gf(fes); - // DEPRECATION_WARNING_OFF - // u_gf.SetData(x.GetData()); - // DEPRECATION_WARNING_ON - - // T->SetIntPoint(&zeroSlopeIP); - // u_gf.GetGradient(*T, grad_u); - - // int dof; - // LOG_DEBUG(logger, "Adjusting the residual to enforce the zero slope constraint by {:0.4E}...", -alpha*grad_u[0]); - // double rNorm = r.Norml2(); - // LOG_INFO(logger, "||r_B|| = {:0.4E}", rNorm); - // for (int i = 0; i < zeroSlopeDofs.Size(); i++) { - // dof = zeroSlopeDofs[i]; - // r[dof] -= alpha * grad_u[0]; - // r[dof] -= alpha * grad_u[1]; - // r[dof] -= alpha * grad_u[2]; - // } - // rNorm = r.Norml2(); - // LOG_INFO(logger, "||r_A|| = {:0.4E}", rNorm); - // // This still is not working; however, I think I am close. I also need to modify the jacobain. - // } void ZeroSlopeNewtonSolver::Mult(const mfem::Vector &b, mfem::Vector &x) const { using namespace mfem; @@ -508,7 +452,7 @@ namespace polyMFEMUtils { { r -= b; } - // ComputeConstrainedResidual(x, r); + ComputeConstrainedResidual(x, r); norm0 = norm = initial_norm = Norm(r); if (print_options.first_and_last && !print_options.iterations)