#include <deal.II/base/mpi_remote_point_evaluation.h>

Classes
struct	CellData

Public Member Functions
	RemotePointEvaluation (const double tolerance=1e-6, const bool enforce_unique_mapping=false)

	~RemotePointEvaluation ()

void	reinit (const std::vector< Point< spacedim > > &points, const Triangulation< dim, spacedim > &tria, const Mapping< dim, spacedim > &mapping)

template<typename T >
void	evaluate_and_process (std::vector< T > &output, std::vector< T > &buffer, const std::function< void(const ArrayView< T > &, const CellData &)> &evaluation_function) const

template<typename T >
void	process_and_evaluate (const std::vector< T > &input, std::vector< T > &buffer, const std::function< void(const ArrayView< const T > &, const CellData &)> &evaluation_function) const

const std::vector< unsigned int > &	get_point_ptrs () const

bool	is_map_unique () const

const Triangulation< dim, spacedim > &	get_triangulation () const

const Mapping< dim, spacedim > &	get_mapping () const

bool	is_ready () const

Private Attributes
const double	tolerance

const bool	enforce_unique_mapping

boost::signals2::connection	tria_signal

bool	ready_flag

SmartPointer< const Triangulation< dim, spacedim > >	tria

SmartPointer< const Mapping< dim, spacedim > >	mapping

bool	unique_mapping

std::vector< unsigned int >	point_ptrs

std::vector< unsigned int >	recv_permutation

std::vector< unsigned int >	recv_ptrs

std::vector< unsigned int >	recv_ranks

CellData	cell_data

std::vector< unsigned int >	send_permutation

std::vector< unsigned int >	send_ranks

std::vector< unsigned int >	send_ptrs

Detailed Description

template<int dim, int spacedim = dim>
class Utilities::MPI::RemotePointEvaluation< dim, spacedim >

Helper class to access values on non-matching grids.

Note: The name of the fields are chosen with the method evaluate_and_process() in mind. Here, quantities are computed at specified arbitrary positioned points (and even on remote processes in the MPI universe) cell by cell and these values are sent to requesting processes, which receive the result and resort the result according to the points.

Definition at line 44 of file mpi_remote_point_evaluation.h.

Constructor & Destructor Documentation

◆ RemotePointEvaluation()

template<int dim, int spacedim>

Utilities::MPI::RemotePointEvaluation< dim, spacedim >::RemotePointEvaluation	(	const double	tolerance = `1e-6`,
		const bool	enforce_unique_mapping = `false`
	)

Constructor.

Parameters

tolerance	Tolerance in terms of unit cell coordinates for determining all cells around a point passed to the class during reinit(). Depending on the problem, it might be necessary to adjust the tolerance in order to be able to identify a cell. Floating point arithmetic implies that a point will, in general, not lie exactly on a vertex, edge, or face.
enforce_unique_mapping	Enforce unique mapping, i.e., (one-to-one) relation of points and cells.

Definition at line 39 of file mpi_remote_point_evaluation.cc.

◆ ~RemotePointEvaluation()

template<int dim, int spacedim>

Utilities::MPI::RemotePointEvaluation< dim, spacedim >::~RemotePointEvaluation

Destructor.

Definition at line 50 of file mpi_remote_point_evaluation.cc.

Member Function Documentation

◆ reinit()

template<int dim, int spacedim>

void Utilities::MPI::RemotePointEvaluation< dim, spacedim >::reinit	(	const std::vector< Point< spacedim > > &	points,
		const Triangulation< dim, spacedim > &	tria,
		const Mapping< dim, spacedim > &	mapping
	)

Set up internal data structures and communication pattern based on a list of points points and mesh description (tria and mapping).

Definition at line 60 of file mpi_remote_point_evaluation.cc.

◆ evaluate_and_process()

template<int dim, int spacedim>

template<typename T >

void Utilities::MPI::RemotePointEvaluation< dim, spacedim >::evaluate_and_process	(	std::vector< T > &	output,
		std::vector< T > &	buffer,
		const std::function< void(const ArrayView< T > &, const CellData &)> &	evaluation_function
	)		const

Evaluate function evaluation_function in the given points and triangulation. The result is stored in output.

Note: If the map of points to cells is not a one-to-one relation (is_map_unique()==false), the result needs to be processed with the help of get_point_ptrs(). This might be the case if a point coincides with a geometric entity (e.g., vertex) that is shared by multiple cells or a point is outside of the computational domain.

Definition at line 257 of file mpi_remote_point_evaluation.h.

◆ process_and_evaluate()

template<int dim, int spacedim>

template<typename T >

void Utilities::MPI::RemotePointEvaluation< dim, spacedim >::process_and_evaluate	(	const std::vector< T > &	input,
		std::vector< T > &	buffer,
		const std::function< void(const ArrayView< const T > &, const CellData &)> &	evaluation_function
	)		const

This method is the inverse of the method evaluate_and_process(). It makes the data at the points, provided by input, available in the function evaluation_function.

Definition at line 371 of file mpi_remote_point_evaluation.h.

◆ get_point_ptrs()

template<int dim, int spacedim>

const std::vector< unsigned int > & Utilities::MPI::RemotePointEvaluation< dim, spacedim >::get_point_ptrs

Return a CRS-like data structure to determine the position of the result corresponding a point and the amount.

Definition at line 162 of file mpi_remote_point_evaluation.cc.

◆ is_map_unique()

template<int dim, int spacedim>

bool Utilities::MPI::RemotePointEvaluation< dim, spacedim >::is_map_unique

Return if points and cells have a one-to-one relation. This is not the case if a point is not owned by any cell (the point is outside of the domain) or if multiple cells own the point (the point is positioned on a geometric entity shared by neighboring cells).

Definition at line 171 of file mpi_remote_point_evaluation.cc.

◆ get_triangulation()

template<int dim, int spacedim>

const Triangulation< dim, spacedim > & Utilities::MPI::RemotePointEvaluation< dim, spacedim >::get_triangulation

Return the Triangulation object used during reinit().

Definition at line 180 of file mpi_remote_point_evaluation.cc.

◆ get_mapping()

template<int dim, int spacedim>

const Mapping< dim, spacedim > & Utilities::MPI::RemotePointEvaluation< dim, spacedim >::get_mapping

Return the Mapping object used during reinit().

Definition at line 189 of file mpi_remote_point_evaluation.cc.

◆ is_ready()

template<int dim, int spacedim>

bool Utilities::MPI::RemotePointEvaluation< dim, spacedim >::is_ready

Return if the internal data structures have been set up and if yes whether they are still valid (and not invalidated due to changes of the Triangulation).

Definition at line 198 of file mpi_remote_point_evaluation.cc.