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