ProjectedStatefulMaterialNodalPatchRecovery.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
 
// MOOSE includes
#include "ElementUserObject.h"
#include "SerialAccess.h"
 
#include <unordered_map>
 
class ProjectedStatefulMaterialNodalPatchRecoveryBase : public ElementUserObject
{
public:
  static InputParameters validParams();
 
  ProjectedStatefulMaterialNodalPatchRecoveryBase(const InputParameters & parameters)
    : ElementUserObject(parameters)
  {
  }
 
  virtual Real nodalPatchRecovery(const Point & p,
                                  const std::vector<dof_id_type> & elem_ids,
                                  std::size_t component) const = 0;
};
 
template <typename T, bool is_ad>
class ProjectedStatefulMaterialNodalPatchRecoveryTempl
  : public ProjectedStatefulMaterialNodalPatchRecoveryBase
{
public:
  static InputParameters validParams();
 
  ProjectedStatefulMaterialNodalPatchRecoveryTempl(const InputParameters & parameters);
 
  virtual Real nodalPatchRecovery(const Point & p,
                                  const std::vector<dof_id_type> & elem_ids,
                                  std::size_t component) const override;
 
  virtual void initialSetup() override;
 
  virtual void initialize() override;
  virtual void execute() override;
  virtual void threadJoin(const UserObject &) override;
  virtual void finalize() override;
 
private:
  RealEigenVector evaluateBasisFunctions(const Point & q_point) const;
 
  typedef std::pair<RealEigenMatrix, std::vector<RealEigenVector>> ElementData;
 
  unsigned int _qp;
 
  std::size_t _n_components;
 
  const unsigned int _patch_polynomial_order;
 
  const std::vector<std::vector<unsigned int>> _multi_index;
 
  const unsigned int _q;
 
  const GenericMaterialProperty<T, is_ad> & _prop;
 
  std::unordered_map<dof_id_type, ElementData> _abs;
 
  const SubdomainID & _current_subdomain_id;
 
  std::unordered_map<SubdomainID, std::vector<MaterialBase *>> _required_materials;
};
 
typedef ProjectedStatefulMaterialNodalPatchRecoveryTempl<Real, false>
    ProjectedStatefulMaterialNodalPatchRecoveryReal;
typedef ProjectedStatefulMaterialNodalPatchRecoveryTempl<Real, true>
    ADProjectedStatefulMaterialNodalPatchRecoveryReal;
typedef ProjectedStatefulMaterialNodalPatchRecoveryTempl<RealVectorValue, false>
    ProjectedStatefulMaterialNodalPatchRecoveryRealVectorValue;
typedef ProjectedStatefulMaterialNodalPatchRecoveryTempl<RealVectorValue, true>
    ADProjectedStatefulMaterialNodalPatchRecoveryRealVectorValue;
typedef ProjectedStatefulMaterialNodalPatchRecoveryTempl<RankTwoTensor, false>
    ProjectedStatefulMaterialNodalPatchRecoveryRankTwoTensor;
typedef ProjectedStatefulMaterialNodalPatchRecoveryTempl<RankTwoTensor, true>
    ADProjectedStatefulMaterialNodalPatchRecoveryRankTwoTensor;
typedef ProjectedStatefulMaterialNodalPatchRecoveryTempl<RankFourTensor, false>
    ProjectedStatefulMaterialNodalPatchRecoveryRankFourTensor;
typedef ProjectedStatefulMaterialNodalPatchRecoveryTempl<RankFourTensor, true>
    ADProjectedStatefulMaterialNodalPatchRecoveryRankFourTensor;
 
// Prevent implicit instantiation in other translation units where these classes are used
extern template class ProjectedStatefulMaterialNodalPatchRecoveryTempl<Real, false>;
extern template class ProjectedStatefulMaterialNodalPatchRecoveryTempl<Real, true>;
extern template class ProjectedStatefulMaterialNodalPatchRecoveryTempl<RealVectorValue, false>;
extern template class ProjectedStatefulMaterialNodalPatchRecoveryTempl<RealVectorValue, true>;
extern template class ProjectedStatefulMaterialNodalPatchRecoveryTempl<RankTwoTensor, false>;
extern template class ProjectedStatefulMaterialNodalPatchRecoveryTempl<RankTwoTensor, true>;
extern template class ProjectedStatefulMaterialNodalPatchRecoveryTempl<RankFourTensor, false>;
extern template class ProjectedStatefulMaterialNodalPatchRecoveryTempl<RankFourTensor, true>;