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>;