NodalPatchRecoveryBase.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"
 
class ElementSubdomainModifierBase;
 
class NodalPatchRecoveryBase : public ElementUserObject
{
public:
  static InputParameters validParams();
 
  NodalPatchRecoveryBase(const InputParameters & parameters);
 
  virtual Real nodalPatchRecovery(const Point & p, const std::vector<dof_id_type> & elem_ids) const;
 
  const RealEigenVector getCoefficients(const std::vector<dof_id_type> & elem_ids) const;
 
  const RealEigenVector getCachedCoefficients(const std::vector<dof_id_type> & elem_ids);
 
  void initialize() override;
  void execute() override;
  void threadJoin(const UserObject &) override;
  void finalize() override;
 
  const std::vector<std::vector<unsigned int>> & multiIndex() const { return _multi_index; }
 
protected:
  virtual Real computeValue() = 0;
 
  unsigned int _qp;
 
private:
  std::unordered_map<processor_id_type, std::vector<dof_id_type>>
  gatherRequestList(const std::vector<dof_id_type> & specific_elems);
 
  std::unordered_map<processor_id_type, std::vector<dof_id_type>> gatherRequestList();
 
  RealEigenVector evaluateBasisFunctions(const Point & q_point) const;
 
  void sync(const std::vector<dof_id_type> & specific_elems);
 
  void sync();
 
  void
  syncHelper(const std::unordered_map<processor_id_type, std::vector<dof_id_type>> & query_ids);
 
  void
  addToQuery(const libMesh::Elem * elem,
             std::unordered_map<processor_id_type, std::vector<dof_id_type>> & query_ids) const;
 
  const unsigned int _patch_polynomial_order;
 
  const std::vector<std::vector<unsigned int>> _multi_index;
 
  const unsigned int _q;
 
  std::map<dof_id_type, RealEigenMatrix> _Ae;
 
  std::map<dof_id_type, RealEigenVector> _be;
 
  mutable std::vector<dof_id_type> _cached_elem_ids;
  mutable RealEigenVector _cached_coef;
 
  bool _distributed_mesh;
 
  std::vector<int> _proc_ids;
 
  friend class ElementSubdomainModifierBase;
};