ElementAndTraceScalarHDGAssemblyHelper.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
 
#pragma once
 
#include "MooseArray.h"
#include "MooseFunctorForward.h"
#include "MooseTypes.h"
#include "ThreeMaterialPropertyInterface.h"
 
#include <array>
#include <set>
#include <string>
#include <vector>
 
template <typename>
class MooseVariableFE;
class Assembly;
class MooseObject;
class MooseVariableDependencyInterface;
class SystemBase;
class TransientInterface;
struct ADResidualsPacket;
 
class ElementAndTraceScalarHDGAssemblyHelper : public ThreeMaterialPropertyInterface
{
public:
  static InputParameters validParams();
 
  ElementAndTraceScalarHDGAssemblyHelper(const MooseObject * moose_obj,
                                         MooseVariableDependencyInterface * mvdi,
                                         const TransientInterface * ti,
                                         SystemBase & sys,
                                         const Assembly & assembly,
                                         THREAD_ID tid,
                                         const std::set<SubdomainID> & blocks_ids,
                                         const std::set<BoundaryID> & boundary_ids);
 
  void resizeResiduals();
 
  virtual void scalarVolume() = 0;
 
  virtual void scalarFace() = 0;
 
  virtual void lmFace() = 0;
 
  virtual void scalarDirichlet(const Moose::Functor<Real> & dirichlet_value) = 0;
 
  void lmDirichlet(const Moose::Functor<Real> & dirichlet_value);
 
  void lmPrescribedFlux(const Moose::Functor<Real> & flux_value);
 
  std::array<ADResidualsPacket, 2> taggingData() const;
 
  std::set<std::string> additionalROVariables();
 
  virtual ~ElementAndTraceScalarHDGAssemblyHelper() = default;
 
protected:
  const TransientInterface & _ti;
 
  const MooseVariableFE<Real> & _u_var;
 
  const MooseVariableFE<Real> & _u_face_var;
 
  const std::vector<dof_id_type> & _u_dof_indices;
 
  const std::vector<dof_id_type> & _lm_u_dof_indices;
 
  const MooseArray<ADReal> & _u_sol;
 
  const MooseArray<ADReal> & _lm_u_sol;
 
  const MooseArray<std::vector<Real>> & _scalar_phi;
 
  const MooseArray<std::vector<RealVectorValue>> & _grad_scalar_phi;
 
  const MooseArray<std::vector<Real>> & _scalar_phi_face;
 
  const MooseArray<std::vector<Real>> & _lm_phi_face;
 
  const Real & _elem_volume;
 
  const Real & _side_area;
 
  const Elem * const & _current_elem;
 
  const unsigned int & _current_side;
 
  const Elem * const & _current_side_elem;
 
  const MooseArray<Real> & _JxW;
 
  const libMesh::QBase * const & _qrule;
 
  const MooseArray<Point> & _q_point;
 
  const MooseArray<Real> & _JxW_face;
 
  const libMesh::QBase * const & _qrule_face;
 
  const MooseArray<Point> & _q_point_face;
 
  const MooseArray<Point> & _normals;
 
  DenseVector<ADReal> _scalar_re;
 
  DenseVector<ADReal> _lm_re;
};
 
inline void
ElementAndTraceScalarHDGAssemblyHelper::resizeResiduals()
{
  _scalar_re.resize(_u_dof_indices.size());
  _lm_re.resize(_lm_u_dof_indices.size());
}