NodalPatchRecovery.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 "AuxKernel.h"
#include "FEProblem.h"
 
class NodalPatchRecovery : public AuxKernel
{
public:
  static InputParameters validParams();
 
  NodalPatchRecovery(const InputParameters & parameters);
  virtual ~NodalPatchRecovery() {}
  template <typename T>
  const MaterialProperty<T> & getMaterialProperty(const std::string & name);
  template <typename T, bool is_ad>
  const GenericMaterialProperty<T, is_ad> & getGenericMaterialProperty(const std::string & name);
 
  template <typename T>
  const MaterialProperty<T> & getMaterialPropertyOld(const std::string & name);
  template <typename T>
  const MaterialProperty<T> & getMaterialPropertyOlder(const std::string & name);
 
protected:
  virtual std::vector<std::vector<unsigned int>> nChooseK(unsigned int N, unsigned int K);
 
  virtual std::vector<std::vector<unsigned int>> multiIndex(unsigned int dim, unsigned int order);
 
  virtual void computePVector(Point q_point);
 
  virtual void accumulateAMatrix();
 
  virtual void accumulateBVector(Real val);
 
  virtual void reinitPatch();
 
  virtual void compute() override;
 
private:
  const unsigned int _patch_polynomial_order;
 
  FEProblemBase & _fe_problem;
  std::vector<std::vector<unsigned int>> _multi_index;
  DenseVector<Number> _P;
  DenseMatrix<Number> _A;
  DenseVector<Number> _B;
  DenseVector<Number> _coef;
};
 
template <typename T>
const MaterialProperty<T> &
NodalPatchRecovery::getMaterialProperty(const std::string & name)
{
  return MaterialPropertyInterface::getMaterialProperty<T>(name);
}
 
template <typename T, bool is_ad>
const GenericMaterialProperty<T, is_ad> &
NodalPatchRecovery::getGenericMaterialProperty(const std::string & name)
{
  return MaterialPropertyInterface::getGenericMaterialProperty<T, is_ad>(name);
}
 
template <typename T>
const MaterialProperty<T> &
NodalPatchRecovery::getMaterialPropertyOld(const std::string & name)
{
  return MaterialPropertyInterface::getMaterialPropertyOld<T>(name);
}
 
template <typename T>
const MaterialProperty<T> &
NodalPatchRecovery::getMaterialPropertyOlder(const std::string & name)
{
  return MaterialPropertyInterface::getMaterialPropertyOlder<T>(name);
}