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"

#include <array>

// Forward Declarations
class FeatureFloodCount;

class FeatureVolumeVectorPostprocessor : public GeneralVectorPostprocessor,
                                         public MooseVariableDependencyInterface,
                                         public BoundaryRestrictable
{
public:
  static InputParameters validParams();

  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;

  std::array<VectorPostprocessorValue *, 3> _centroid;
};