Implements aabb_pe_finder given a dofmap.

Data Types
type	aabb_pe_finder_t
	Implements global interpolation for arbitrary points in the domain. More...

Functions/Subroutines
subroutine	aabb_pe_finder_init (this, x, y, z, nelv, xh, comm, padding)
	Initialize the global interpolation object on a set of coordinates.

subroutine	aabb_pe_finder_free (this)
	Destructor.

subroutine	aabb_pe_finder_find_candidates (this, my_point, pe_candidates)
	Find pe candidates.

subroutine	aabb_pe_finder_find_candidates_batch (this, points, n_points, points_at_pe, n_points_pe)

Variables
integer, parameter, public	glob_map_size = 4096
	Minimum number of total boxes in the aabb tree.

Function/Subroutine Documentation

◆ aabb_pe_finder_find_candidates()

subroutine aabb_pe_finder::aabb_pe_finder_find_candidates	(	class(aabb_pe_finder_t), intent(inout)	this,
		type(point_t), intent(in)	my_point,
		type(stack_i4_t), intent(inout)	pe_candidates
	)

private

Parameters

my_point	Point to find candidates for.
pe_candidates	Candidates for the point.

Definition at line 238 of file aabb_pe_finder.f90.

subroutine aabb_pe_finder::aabb_pe_finder_find_candidates_batch	(	class(aabb_pe_finder_t), intent(inout)	this,
		real(kind=rp), dimension(3,n_points), intent(in)	points,
		integer, intent(in)	n_points,
		type(stack_i4_t), dimension(0:(this%pe_size-1)), intent(inout)	points_at_pe,
		integer, dimension(0:(this%pe_size-1)), intent(inout)	n_points_pe
	)

private

Check which ranks might have this point

Definition at line 253 of file aabb_pe_finder.f90.

subroutine aabb_pe_finder::aabb_pe_finder_free ( class(aabb_pe_finder_t), intent(inout) this )

private

Definition at line 226 of file aabb_pe_finder.f90.

subroutine aabb_pe_finder::aabb_pe_finder_init	(	class(aabb_pe_finder_t), intent(inout)	this,
		real(kind=rp), dimension(:), intent(in), target	x,
		real(kind=rp), dimension(:), intent(in), target	y,
		real(kind=rp), dimension(:), intent(in), target	z,
		integer, intent(in)	nelv,
		type(space_t), intent(in), target	xh,
		type(mpi_comm), intent(in), optional	comm,
		real(kind=dp), intent(in)	padding
	)

Parameters

x	x-coordinates.
y	y-coordinates.
z	z-coordinates.
gdim	Geometric dimension.
nelv	Number of elements of the mesh in which to search for the points.
Xh	Space on which to interpolate.
tol	Tolerance for Newton iterations.

At least 2 boxes

Traverse the local tree and find ther top boxes

Set the boxes to be empty

Get boxes from all ranks

Create global tree for each rank

Definition at line 92 of file aabb_pe_finder.f90.

integer, parameter, public aabb_pe_finder::glob_map_size = 4096

Definition at line 59 of file aabb_pe_finder.f90.