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 "Executioner.h"
      13             : #include "MFEMProblemSolve.h"
      14             : #include "MFEMProblem.h"
      15             : 
      16             : InputParameters
      17       17958 : MFEMProblemSolve::validParams()
      18             : {
      19       17958 :   InputParameters params = emptyInputParameters();
      20       35916 :   params.addClassDescription("Solve object for MFEM problems.");
      21       71832 :   params.addParam<std::string>("device", "Run app on the chosen device.");
      22       71832 :   MooseEnum assembly_levels("legacy full element partial none", "legacy", true);
      23       53874 :   params.addParam<MooseEnum>("assembly_level", assembly_levels, "Matrix assembly level.");
      24       35916 :   return params;
      25       17958 : }
      26             : 
      27         349 : MFEMProblemSolve::MFEMProblemSolve(
      28             :     Executioner & ex,
      29         349 :     std::vector<std::shared_ptr<Moose::MFEM::ProblemOperatorBase>> & problem_operators)
      30             :   : SolveObject(ex),
      31         349 :     _mfem_problem(dynamic_cast<MFEMProblem &>(_problem)),
      32         349 :     _problem_operators(problem_operators)
      33             : {
      34         349 :   if (const auto compute_device = _app.getComputeDevice())
      35         276 :     _app.setMFEMDevice(*compute_device, Moose::PassKey<MFEMProblemSolve>());
      36             :   else
      37         292 :     _app.setMFEMDevice(isParamValid("device")    ? getParam<std::string>("device")
      38          68 :                        : _app.isUltimateMaster() ? "cpu"
      39             :                                                  : "",
      40         422 :                        Moose::PassKey<MFEMProblemSolve>());
      41         349 : }
      42             : 
      43             : bool
      44         877 : MFEMProblemSolve::solve()
      45             : {
      46             :   // FixedPointSolve::solve() is libMesh specific, so we need
      47             :   // to include all steps therein relevant to the MFEM backend here.
      48             : 
      49         877 :   bool converged = true;
      50             : 
      51             :   // need to back up multi-apps even when not doing fixed point iteration for recovering from failed
      52             :   // multiapp solve
      53         877 :   _mfem_problem.backupMultiApps(EXEC_MULTIAPP_FIXED_POINT_BEGIN);
      54         877 :   _mfem_problem.backupMultiApps(EXEC_TIMESTEP_BEGIN);
      55         877 :   _mfem_problem.backupMultiApps(EXEC_TIMESTEP_END);
      56         877 :   _mfem_problem.backupMultiApps(EXEC_MULTIAPP_FIXED_POINT_END);
      57             : 
      58             :   // Solve step begins
      59         877 :   _executioner.preSolve();
      60         877 :   _mfem_problem.execTransfers(EXEC_TIMESTEP_BEGIN);
      61             : 
      62         877 :   _mfem_problem.execute(EXEC_MULTIAPP_FIXED_POINT_BEGIN);
      63         877 :   _mfem_problem.execTransfers(EXEC_MULTIAPP_FIXED_POINT_BEGIN);
      64         877 :   _mfem_problem.execMultiApps(EXEC_MULTIAPP_FIXED_POINT_BEGIN, true);
      65         877 :   _mfem_problem.outputStep(EXEC_MULTIAPP_FIXED_POINT_BEGIN);
      66             : 
      67         877 :   _mfem_problem.execMultiApps(EXEC_TIMESTEP_BEGIN, true);
      68         877 :   _mfem_problem.execute(EXEC_TIMESTEP_BEGIN);
      69         877 :   _mfem_problem.outputStep(EXEC_TIMESTEP_BEGIN);
      70             : 
      71             :   // Update warehouse active objects
      72         877 :   _mfem_problem.updateActiveObjects();
      73             : 
      74         877 :   if (_mfem_problem.shouldSolve())
      75        1658 :     for (const auto & problem_operator : _problem_operators)
      76         829 :       problem_operator->Solve();
      77         877 :   _mfem_problem.displaceMesh();
      78             : 
      79             :   // Execute user objects, transfers, and multiapps at timestep end
      80         877 :   _mfem_problem.onTimestepEnd();
      81         877 :   _mfem_problem.execute(EXEC_TIMESTEP_END);
      82         877 :   _mfem_problem.execTransfers(EXEC_TIMESTEP_END);
      83         877 :   _mfem_problem.execMultiApps(EXEC_TIMESTEP_END, true);
      84         877 :   _executioner.postSolve();
      85             :   // Solve step ends
      86             : 
      87         877 :   if (converged)
      88             :   {
      89             :     // Fixed point iteration loop ends right above
      90         877 :     _mfem_problem.execute(EXEC_MULTIAPP_FIXED_POINT_END);
      91         877 :     _mfem_problem.execTransfers(EXEC_MULTIAPP_FIXED_POINT_END);
      92         877 :     _mfem_problem.execMultiApps(EXEC_MULTIAPP_FIXED_POINT_END, true);
      93         877 :     _mfem_problem.outputStep(EXEC_MULTIAPP_FIXED_POINT_END);
      94             :   }
      95             : 
      96         877 :   return converged;
      97             : }
      98             : #endif