NLCurlCurlIntegrator.h

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 MOOSE_MFEM_ENABLED

#pragma once
#include "MFEMVectorMagnitudeCoefficient.h"

namespace Moose::MFEM
{
class NLCurlCurlJacMatrixCoefficient : public mfem::MatrixCoefficient
{
public:
  NLCurlCurlJacMatrixCoefficient(mfem::Coefficient & k,
                                 mfem::Coefficient & curlu_dk_dcurlu,
                                 mfem::VectorCoefficient & curlu_vec,
                                 mfem::real_t curlu_zero_tol);

  void Eval(mfem::DenseMatrix & K,
            mfem::ElementTransformation & T,
            const mfem::IntegrationPoint & ip) override;
  void SetTime(mfem::real_t t) override;

protected:
  mfem::Coefficient & _k_coef;
  mfem::Coefficient & _curlu_dk_dcurlu_coef;
  const mfem::real_t _curlu_zero_tol;
  mfem::NormalizedVectorCoefficient _curlu_hat_coef;
};

class NLCurlCurlIntegrator : public mfem::NonlinearFormIntegrator
{
public:
  NLCurlCurlIntegrator(mfem::Coefficient & k,
                       mfem::Coefficient & curlu_dk_dcurlu,
                       mfem::VectorCoefficient & curlu_vec,
                       mfem::real_t curlu_zero_tol = 1e-32,
                       const mfem::IntegrationRule * ir = nullptr);

  virtual void AssembleElementVector(const mfem::FiniteElement & el,
                                     mfem::ElementTransformation & Tr,
                                     const mfem::Vector & elfun,
                                     mfem::Vector & elvect) override;
  virtual void AssembleElementGrad(const mfem::FiniteElement & el,
                                   mfem::ElementTransformation & Tr,
                                   const mfem::Vector & elfun,
                                   mfem::DenseMatrix & elmat) override;

protected:
  mfem::CurlCurlIntegrator _curlcurl_res_integ;      // (k(|curl u|) curl u, curl phi_j)
  NLCurlCurlJacMatrixCoefficient _curlcurl_jac_matrix_coef;
  mfem::CurlCurlIntegrator _curlcurl_jac_integ;
};
}

#endif