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 "MFEMGradAux.h"
      13             : #include "MFEMProblem.h"
      14             : 
      15             : registerMooseObject("MooseApp", MFEMGradAux);
      16             : 
      17             : InputParameters
      18        2534 : MFEMGradAux::validParams()
      19             : {
      20        2534 :   InputParameters params = MFEMAuxKernel::validParams();
      21        5068 :   params.addClassDescription(
      22             :       "Calculates the gradient of an H1 conforming source variable and stores the result"
      23             :       " on an H(curl) conforming ND result auxvariable");
      24       10136 :   MFEMExecutedObject::addRequiredDependencyParam<VariableName>(
      25             :       params, "source", "Scalar H1 MFEMVariable to take the gradient of.");
      26        7602 :   params.addParam<mfem::real_t>("scale_factor", 1.0, "Factor to scale result auxvariable by.");
      27        2534 :   return params;
      28           0 : }
      29             : 
      30         218 : MFEMGradAux::MFEMGradAux(const InputParameters & parameters)
      31             :   : MFEMAuxKernel(parameters),
      32         218 :     _source_var_name(getParam<VariableName>("source")),
      33         218 :     _source_var(*getMFEMProblem().getGridFunction(_source_var_name)),
      34         436 :     _scale_factor(getParam<mfem::real_t>("scale_factor")),
      35         436 :     _grad(_source_var.ParFESpace(), _result_var.ParFESpace())
      36             : {
      37         218 :   _sequence = _source_var.GetSequence() + _result_var.GetSequence();
      38         218 :   _grad.Assemble();
      39         218 :   _grad.Finalize();
      40         218 : }
      41             : 
      42             : // Computes the auxvariable.
      43             : void
      44         218 : MFEMGradAux::execute()
      45             : {
      46         218 :   update();
      47         218 :   _grad.AddMult(_source_var, _result_var = 0, _scale_factor);
      48         218 : }
      49             : 
      50             : void
      51         218 : MFEMGradAux::update()
      52             : {
      53         218 :   if (long sequence = _source_var.GetSequence() + _result_var.GetSequence() > _sequence)
      54             :   {
      55          26 :     _sequence = sequence;
      56          26 :     _grad.Update();
      57          26 :     _grad.Assemble();
      58          26 :     _grad.Finalize();
      59             :   }
      60         218 : }
      61             : 
      62             : #endif