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        2168 : MFEMNewtonNonlinearSolver::validParams()
      19             : {
      20        2168 :   InputParameters params = Moose::MFEM::NonlinearSolverBase::validParams();
      21        2168 :   params.addClassDescription("MFEM native nonlinear solver using Newton's method.");
      22        2168 :   return params;
      23           0 : }
      24             : 
      25          35 : MFEMNewtonNonlinearSolver::MFEMNewtonNonlinearSolver(const InputParameters & parameters)
      26          35 :   : Moose::MFEM::NonlinearSolverBase(parameters)
      27             : {
      28          35 :   ConstructSolver();
      29          35 : }
      30             : 
      31             : void
      32          35 : MFEMNewtonNonlinearSolver::ConstructSolver()
      33             : {
      34          35 :   auto solver = std::make_unique<mfem::NewtonSolver>(getMFEMProblem().getComm());
      35          70 :   solver->iterative_mode = getParam<bool>("use_initial_guess");
      36          70 :   solver->SetRelTol(getParam<mfem::real_t>("rel_tol"));
      37          70 :   solver->SetAbsTol(getParam<mfem::real_t>("abs_tol"));
      38          70 :   solver->SetMaxIter(getParam<unsigned int>("max_its"));
      39          70 :   solver->SetPrintLevel(getParam<unsigned int>("print_level"));
      40          35 :   _solver = std::move(solver);
      41          35 : }
      42             : 
      43             : void
      44         371 : MFEMNewtonNonlinearSolver::SetOperator(const mfem::Operator & op)
      45             : {
      46         371 :   GetSolver().SetOperator(op);
      47         371 : }
      48             : 
      49             : void
      50         371 : MFEMNewtonNonlinearSolver::SetLinearSolver(mfem::Solver & solver)
      51             : {
      52         371 :   static_cast<mfem::NewtonSolver &>(GetSolver()).SetSolver(solver);
      53         371 : }
      54             : 
      55             : void
      56         371 : MFEMNewtonNonlinearSolver::Mult(const mfem::Vector & rhs, mfem::Vector & x)
      57             : {
      58         371 :   GetSolver().Mult(rhs, x);
      59         371 : }
      60             : #endif