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