ElementsAlongPlane.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 "ElementsAlongPlane.h"
#include "ElementsIntersectedByPlane.h"
#include "MooseMesh.h"

registerMooseObject("MooseApp", ElementsAlongPlane);

InputParameters
ElementsAlongPlane::validParams()
{
  InputParameters params = GeneralVectorPostprocessor::validParams();

  params.addRequiredParam<Point>("point", "Point in the plane");
  params.addRequiredParam<Point>("normal", "Normal vector to the plane");
  params.addClassDescription(
      "Outputs the IDs of every element intersected by a user-defined plane");
  return params;
}

ElementsAlongPlane::ElementsAlongPlane(const InputParameters & parameters)
  : GeneralVectorPostprocessor(parameters),
    _p0(getParam<Point>("point")),
    _normal(getParam<Point>("normal")),
    _elem_ids(declareVector("elem_ids"))
{
  _fe_problem.mesh().errorIfDistributedMesh("ElementsAlongPlane");
}

void
ElementsAlongPlane::initialize()
{
  _elem_ids.clear();
}

void
ElementsAlongPlane::execute()
{
  std::vector<const Elem *> intersected_elems;

  Moose::elementsIntersectedByPlane(_p0, _normal, _fe_problem.mesh(), intersected_elems);

  unsigned int num_elems = intersected_elems.size();

  _elem_ids.resize(num_elems);

  for (unsigned int i = 0; i < num_elems; i++)
    _elem_ids[i] = intersected_elems[i]->id();
}