SlopeReconstructionBase.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 "BCUserObject.h"
#include "ElementLoopUserObject.h"
 
using libMesh::RealGradient;
 
class SlopeReconstructionBase : public ElementLoopUserObject
{
public:
  static InputParameters validParams();
 
  SlopeReconstructionBase(const InputParameters & parameters);
 
  virtual void initialize();
  virtual void finalize();
 
  virtual void computeElement();
 
  virtual const std::vector<RealGradient> & getElementSlope(dof_id_type elementid) const;
 
  virtual const std::vector<Real> & getElementAverageValue(dof_id_type elementid) const;
 
  virtual const std::vector<Real> & getBoundaryAverageValue(dof_id_type elementid,
                                                            unsigned int side) const;
 
  virtual const Point & getSideCentroid(dof_id_type elementid, dof_id_type neighborid) const;
 
  virtual const Point & getBoundarySideCentroid(dof_id_type elementid, unsigned int side) const;
 
  virtual const Point & getSideNormal(dof_id_type elementid, dof_id_type neighborid) const;
 
  virtual const Point & getBoundarySideNormal(dof_id_type elementid, unsigned int side) const;
 
  virtual const Real & getSideArea(dof_id_type elementid, dof_id_type neighborid) const;
 
  virtual const Real & getBoundarySideArea(dof_id_type elementid, unsigned int side) const;
 
  virtual void reconstructElementSlope() = 0;
 
  virtual void meshChanged();
 
protected:
  virtual void serialize(std::string & serialized_buffer);
  virtual void deserialize(std::vector<std::string> & serialized_buffers);
 
  std::map<dof_id_type, std::vector<RealGradient>> & _rslope;
 
  std::map<dof_id_type, std::vector<Real>> & _avars;
 
  std::map<std::pair<dof_id_type, unsigned int>, std::vector<Real>> & _bnd_avars;
 
  std::map<std::pair<dof_id_type, dof_id_type>, Point> & _side_centroid;
 
  std::map<std::pair<dof_id_type, unsigned int>, Point> & _bnd_side_centroid;
 
  std::map<std::pair<dof_id_type, dof_id_type>, Real> & _side_area;
 
  std::map<std::pair<dof_id_type, unsigned int>, Real> & _bnd_side_area;
 
  std::map<std::pair<dof_id_type, dof_id_type>, Point> & _side_normal;
 
  std::map<std::pair<dof_id_type, unsigned int>, Point> & _bnd_side_normal;
 
  const MooseArray<Point> & _q_point_face;
  const QBase * const & _qrule_face;
  const MooseArray<Real> & _JxW_face;
  const MooseArray<Point> & _normals_face;
 
  const unsigned int & _side;
 
  const Elem * const & _side_elem;
  const Real & _side_volume;
 
  const Elem * const & _neighbor_elem;
 
  bool _side_geoinfo_cached;
 
private:
  static Threads::spin_mutex _mutex;
};