Line data    Source code

       1             : //* This file is part of the MOOSE framework
       2             : //* https://mooseframework.inl.gov
       3             : //*
       4             : //* All rights reserved, see COPYRIGHT for full restrictions
       5             : //* https://github.com/idaholab/moose/blob/master/COPYRIGHT
       6             : //*
       7             : //* Licensed under LGPL 2.1, please see LICENSE for details
       8             : //* https://www.gnu.org/licenses/lgpl-2.1.html
       9             : 
      10             : #ifdef MOOSE_MFEM_ENABLED
      11             : 
      12             : #include "MFEMNewtonNonlinearSolver.h"
      13             : #include "MFEMProblem.h"
      14             : 
      15             : registerMooseObject("MooseApp", MFEMNewtonNonlinearSolver);
      16             : 
      17             : InputParameters
      18        2212 : MFEMNewtonNonlinearSolver::validParams()
      19             : {
      20        2212 :   InputParameters params = Moose::MFEM::NonlinearSolverBase::validParams();
      21        2212 :   params.addClassDescription("MFEM native nonlinear solver using Newton's method.");
      22        2212 :   return params;
      23           0 : }
      24             : 
      25          43 : MFEMNewtonNonlinearSolver::MFEMNewtonNonlinearSolver(const InputParameters & parameters)
      26          43 :   : Moose::MFEM::NonlinearSolverBase(parameters)
      27             : {
      28          43 :   ConstructSolver();
      29          43 : }
      30             : 
      31             : void
      32          43 : MFEMNewtonNonlinearSolver::ConstructSolver()
      33             : {
      34          43 :   auto solver = std::make_unique<mfem::NewtonSolver>(getMFEMProblem().getComm());
      35          86 :   solver->iterative_mode = getParam<bool>("use_initial_guess");
      36          86 :   solver->SetRelTol(getParam<mfem::real_t>("rel_tol"));
      37          86 :   solver->SetAbsTol(getParam<mfem::real_t>("abs_tol"));
      38          86 :   solver->SetMaxIter(getParam<unsigned int>("max_its"));
      39          86 :   solver->SetPrintLevel(getParam<unsigned int>("print_level"));
      40          43 :   _solver = std::move(solver);
      41          43 : }
      42             : 
      43             : void
      44         671 : MFEMNewtonNonlinearSolver::SetLinearSolver(mfem::Solver & solver)
      45             : {
      46         671 :   static_cast<mfem::NewtonSolver &>(GetSolver()).SetSolver(solver);
      47         671 : }
      48             : #endif