FeatureVolumeVectorPostprocessor.h

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://www.mooseframework.org
//*
//* 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 "GeneralVectorPostprocessor.h"
#include "MooseVariableDependencyInterface.h"
#include "BoundaryRestrictable.h"
 
// Forward Declarations
class FeatureVolumeVectorPostprocessor;
class FeatureFloodCount;
 
template <>
InputParameters validParams<FeatureVolumeVectorPostprocessor>();
 
class FeatureVolumeVectorPostprocessor : public GeneralVectorPostprocessor,
                                         public MooseVariableDependencyInterface,
                                         public BoundaryRestrictable
{
public:
  FeatureVolumeVectorPostprocessor(const InputParameters & parameters);
 
  virtual void initialize() override;
  virtual void execute() override;
  virtual void finalize() override;
 
  Real getFeatureVolume(unsigned int feature_id) const;
 
protected:
  const bool _single_feature_per_elem;
  const bool _output_centroids;
 
  const FeatureFloodCount & _feature_counter;
 
  VectorPostprocessorValue & _var_num;
  VectorPostprocessorValue & _feature_volumes;
  VectorPostprocessorValue & _intersects_bounds;
  VectorPostprocessorValue & _intersects_specified_bounds;
  VectorPostprocessorValue & _percolated;
 
  bool _is_boundary_restricted;
 
private:
  void accumulateVolumes(const Elem * elem,
                         const std::vector<unsigned int> & var_to_features,
                         std::size_t num_features);
 
  void accumulateBoundaryFaces(const Elem * elem,
                               const std::vector<unsigned int> & var_to_features,
                               std::size_t num_features,
                               unsigned int side);
 
  Real computeIntegral(std::size_t var_index) const;
 
  Real computeFaceIntegral(std::size_t var_index) const;
 
  const std::vector<MooseVariableFEBase *> & _vars;
  std::vector<const VariableValue *> _coupled_sln;
 
  MooseMesh & _mesh;
  Assembly & _assembly;
  const MooseArray<Point> & _q_point;
  const QBase * const & _qrule;
  const MooseArray<Real> & _JxW;
  const MooseArray<Real> & _coord;
  const QBase * const & _qrule_face;
  const MooseArray<Real> & _JxW_face;
};