BondStatusConvergedPostprocessorPD.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 "BondStatusConvergedPostprocessorPD.h"

registerMooseObject("PeridynamicsApp", BondStatusConvergedPostprocessorPD);

InputParameters
BondStatusConvergedPostprocessorPD::validParams()
{
  InputParameters params = ElementPostprocessor::validParams();
  params.addClassDescription(
      "Postprocessor to check whether the bond status is converged within a time step");

  return params;
}

BondStatusConvergedPostprocessorPD::BondStatusConvergedPostprocessorPD(
    const InputParameters & parameters)
  : ElementPostprocessor(parameters),
    _bond_status_var(&_subproblem.getStandardVariable(_tid, "bond_status")),
    _pdmesh(dynamic_cast<PeridynamicsMesh &>(_mesh)),
    _dim(_pdmesh.dimension())
{
}

void
BondStatusConvergedPostprocessorPD::initialize()
{
  _num_bond_status_updated = 0;
}

void
BondStatusConvergedPostprocessorPD::execute()
{
  if ((_bond_status_var->getElementalValueOld(_current_elem) -
       _bond_status_var->getElementalValue(_current_elem)) > 0.5)
    _num_bond_status_updated += 1;
}

Real
BondStatusConvergedPostprocessorPD::getValue() const
{
  return _num_bond_status_updated;
}

void
BondStatusConvergedPostprocessorPD::finalize()
{
  gatherSum(_num_bond_status_updated);
}

void
BondStatusConvergedPostprocessorPD::threadJoin(const UserObject & uo)
{
  const auto & pps = static_cast<const BondStatusConvergedPostprocessorPD &>(uo);
  _num_bond_status_updated += pps._num_bond_status_updated;
}