SideIntegralPostprocessor.C

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
 
#include "SideIntegralPostprocessor.h"
 
#include "libmesh/quadrature.h"
 
InputParameters
SideIntegralPostprocessor::validParams()
{
  InputParameters params = SidePostprocessor::validParams();
  return params;
}
 
SideIntegralPostprocessor::SideIntegralPostprocessor(const InputParameters & parameters)
  : SidePostprocessor(parameters), _qp(0), _integral_value(0), _qp_integration(true)
{
}
 
void
SideIntegralPostprocessor::initialSetup()
{
  SidePostprocessor::initialSetup();
 
  if (!_qp_integration && _mesh.isFiniteVolumeInfoDirty())
    errorNoFaceInfo();
}
 
void
SideIntegralPostprocessor::initialize()
{
  _integral_value = 0;
}
 
void
SideIntegralPostprocessor::execute()
{
  _integral_value += computeIntegral();
}
 
Real
SideIntegralPostprocessor::getValue() const
{
  return _integral_value;
}
 
void
SideIntegralPostprocessor::threadJoin(const UserObject & y)
{
  const auto & pps = static_cast<const SideIntegralPostprocessor &>(y);
  _integral_value += pps._integral_value;
}
 
Real
SideIntegralPostprocessor::computeIntegral()
{
  Real sum = 0;
  if (_qp_integration)
    for (_qp = 0; _qp < _qrule->n_points(); _qp++)
      sum += _JxW[_qp] * _coord[_qp] * computeQpIntegral();
  else
  {
    // Finite volume functors integration is over FaceInfo, not quadrature points
    getFaceInfos();
 
    for (auto & fi : _face_infos)
      sum += fi->faceArea() * fi->faceCoord() * computeFaceInfoIntegral(fi);
  }
  return sum;
}
 
void
SideIntegralPostprocessor::finalize()
{
  gatherSum(_integral_value);
}