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 MFEM_ENABLED
      11             : 
      12             : #include "MFEMDivDivKernel.h"
      13             : #include "MFEMProblem.h"
      14             : 
      15             : registerMooseObject("MooseApp", MFEMDivDivKernel);
      16             : 
      17             : InputParameters
      18        8648 : MFEMDivDivKernel::validParams()
      19             : {
      20        8648 :   InputParameters params = MFEMKernel::validParams();
      21        8648 :   params.addClassDescription(
      22             :       "Adds the domain integrator to an MFEM problem for the bilinear form "
      23             :       "$(k\\vec\\nabla \\cdot \\vec u, \\vec\\nabla \\cdot \\vec v)_\\Omega$ "
      24             :       "arising from the weak form of the grad-div operator "
      25             :       "$-\\vec\\nabla \\left( k \\vec\\nabla \\cdot \\vec u \\right)$.");
      26             : 
      27        8648 :   params.addParam<MFEMScalarCoefficientName>(
      28             :       "coefficient", "1.", "Name of property k to multiply the integrator by");
      29             : 
      30        8648 :   return params;
      31           0 : }
      32             : 
      33           9 : MFEMDivDivKernel::MFEMDivDivKernel(const InputParameters & parameters)
      34             :   : MFEMKernel(parameters),
      35             :     // FIXME: The MFEM bilinear form can also handle vector and matrix
      36             :     // coefficients, so ideally we'd handle all three too.
      37           9 :     _coef(getScalarCoefficient(getParam<MFEMScalarCoefficientName>("coefficient")))
      38             : {
      39           9 : }
      40             : 
      41             : mfem::BilinearFormIntegrator *
      42           9 : MFEMDivDivKernel::createBFIntegrator()
      43             : {
      44           9 :   return new mfem::DivDivIntegrator(_coef);
      45             : }
      46             : 
      47             : #endif