MFEMCurlAux.C

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 MFEM_ENABLED

#include "MFEMCurlAux.h"
#include "MFEMProblem.h"

registerMooseObject("MooseApp", MFEMCurlAux);

/*
Class to set an H(div) auxvariable to be the curl of a H(curl) vector variable.
*/
InputParameters
MFEMCurlAux::validParams()
{
  InputParameters params = MFEMAuxKernel::validParams();
  params.addClassDescription(
      "Calculates the curl of an H(curl) conforming ND source variable and stores the result"
      " on an H(div) conforming RT result auxvariable");
  params.addRequiredParam<VariableName>("source",
                                        "Vector H(curl) MFEMVariable to take the curl of.");
  params.addParam<mfem::real_t>("scale_factor", 1.0, "Factor to scale result auxvariable by.");
  return params;
}

MFEMCurlAux::MFEMCurlAux(const InputParameters & parameters)
  : MFEMAuxKernel(parameters),
    _source_var_name(getParam<VariableName>("source")),
    _source_var(*getMFEMProblem().getProblemData().gridfunctions.Get(_source_var_name)),
    _scale_factor(getParam<mfem::real_t>("scale_factor")),
    _curl(_source_var.ParFESpace(), _result_var.ParFESpace())
{
  _curl.Assemble();
  _curl.Finalize();
}

// Computes the auxvariable.
void
MFEMCurlAux::execute()
{
  _result_var = 0.0;
  _curl.AddMult(_source_var, _result_var, _scale_factor);
}

#endif