CentroidMultiApp.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
 
// MOOSE includes
#include "CentroidMultiApp.h"
#include "MooseMesh.h"
#include "FEProblem.h"
 
// libMesh includes
#include "libmesh/parallel_algebra.h"
 
registerMooseObject("MooseApp", CentroidMultiApp);
 
// TODO: Deprecate and use Positions system
 
InputParameters
CentroidMultiApp::validParams()
{
  InputParameters params = TransientMultiApp::validParams();
  params += BlockRestrictable::validParams();
  params.addClassDescription(
      "Automatically generates Sub-App positions from centroids of elements in the parent app "
      " mesh.");
  params.suppressParameter<std::vector<Point>>("positions");
  params.suppressParameter<std::vector<FileName>>("positions_file");
  params.suppressParameter<std::vector<PositionsName>>("positions_objects");
  return params;
}
 
CentroidMultiApp::CentroidMultiApp(const InputParameters & parameters)
  : TransientMultiApp(parameters), BlockRestrictable(this)
{
}
 
void
CentroidMultiApp::fillPositions()
{
  MooseMesh & parent_app_mesh = _fe_problem.mesh();
 
  for (const auto & elem_ptr : parent_app_mesh.getMesh().active_local_element_ptr_range())
    if (hasBlocks(elem_ptr->subdomain_id()))
      _positions.push_back(elem_ptr->vertex_average());
 
  // Use the comm from the problem this MultiApp is part of
  libMesh::ParallelObject::comm().allgather(_positions);
 
  if (_positions.empty())
    mooseError("No positions found for CentroidMultiapp ", _name);
 
  // An attempt to try to make this parallel stable
  std::sort(_positions.begin(), _positions.end());
}