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