MFEMNewtonNonlinearSolver.h

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://mooseframework.inl.gov
//*
//* All rights reserved, see COPYRIGHT for full restrictions
//* https://github.com/idaholab/moose/blob/master/COPYRIGHT
//*
//* Licensed under LGPL 2.1, please see LICENSE for details
//* https://www.gnu.org/licenses/lgpl-2.1.html
 
#ifdef MOOSE_MFEM_ENABLED
 
#pragma once
 
#include "MFEMNonlinearSolverBase.h"
 
namespace Moose::MFEM
{
class DampedNewtonSolver : public mfem::NewtonSolver
{
public:
  DampedNewtonSolver(MPI_Comm comm) : mfem::NewtonSolver(comm) {}
 
  void SetOperator(const mfem::Operator & op) override;
 
  void SetDampingFactor(mfem::real_t damping_factor) { _damping_factor = damping_factor; }
 
  void SetBacktracking(unsigned int max_its,
                       mfem::real_t contraction_factor,
                       mfem::real_t sufficient_decrease);
 
  mfem::real_t ComputeScalingFactor(const mfem::Vector & x, const mfem::Vector & b) const override;
 
protected:
  mfem::real_t _damping_factor = 1.0;
  unsigned int _ls_max_its = 0;
  mfem::real_t _ls_contraction_factor = 0.5;
  mfem::real_t _ls_sufficient_decrease = 1.0e-4;
  mutable mfem::Vector _ls_x_trial;
  mutable mfem::Vector _ls_r_trial;
};
} // namespace Moose::MFEM
 
class MFEMNewtonNonlinearSolver : public Moose::MFEM::NonlinearSolverBase
{
public:
  static InputParameters validParams();
 
  MFEMNewtonNonlinearSolver(const InputParameters & parameters);
 
  void ConstructSolver() override;
 
  void SetLinearSolver(mfem::Solver & solver) override;
  bool RequiresGradient() const override { return true; }
  bool RequiresExternalLinearSolver() const override { return true; }
};
 
#endif