OversetFindPointsGSLIB enables use of findpts for arbitrary number of overlapping grids. More...

#include <gslib.hpp>

Inheritance diagram for mfem::OversetFindPointsGSLIB:

Collaboration diagram for mfem::OversetFindPointsGSLIB:

Public Member Functions
	OversetFindPointsGSLIB ()

	OversetFindPointsGSLIB (MPI_Comm comm_)

void	Setup (Mesh &m, const int meshid, GridFunction *gfmax=NULL, const double bb_t=0.1, const double newt_tol=1.0e-12, const int npt_max=256)

void	FindPoints (const Vector &point_pos, Array< unsigned int > &point_id, int point_pos_ordering=Ordering::byNODES)

void	Interpolate (const Vector &point_pos, Array< unsigned int > &point_id, const GridFunction &field_in, Vector &field_out, int point_pos_ordering=Ordering::byNODES)

virtual void	Interpolate (const GridFunction &field_in, Vector &field_out)

void	Interpolate (const Vector &point_pos, const GridFunction &field_in, Vector &field_out, int point_pos_ordering=Ordering::byNODES)

void	Interpolate (Mesh &m, const Vector &point_pos, const GridFunction &field_in, Vector &field_out, int point_pos_ordering=Ordering::byNODES)

Public Member Functions inherited from mfem::FindPointsGSLIB
	FindPointsGSLIB ()

	FindPointsGSLIB (MPI_Comm comm_)

virtual	~FindPointsGSLIB ()

void	Setup (Mesh &m, const double bb_t=0.1, const double newt_tol=1.0e-12, const int npt_max=256)

void	FindPoints (const Vector &point_pos, int point_pos_ordering=Ordering::byNODES)

void	FindPoints (Mesh &m, const Vector &point_pos, int point_pos_ordering=Ordering::byNODES, const double bb_t=0.1, const double newt_tol=1.0e-12, const int npt_max=256)
	Setup FindPoints and search positions.

void	Interpolate (const Vector &point_pos, const GridFunction &field_in, Vector &field_out, int point_pos_ordering=Ordering::byNODES)

void	Interpolate (Mesh &m, const Vector &point_pos, const GridFunction &field_in, Vector &field_out, int point_pos_ordering=Ordering::byNODES)

virtual void	SetL2AvgType (AvgType avgtype_)

virtual void	SetDefaultInterpolationValue (double interp_value_)

virtual void	SetDistanceToleranceForPointsFoundOnBoundary (double bdr_tol_)

virtual void	SetGPUtoCPUFallback (bool mode)

virtual void	FreeData ()

virtual const Array< unsigned int > &	GetCode () const

virtual const Array< unsigned int > &	GetElem () const
	Return element number for each point found by FindPoints.

virtual const Array< unsigned int > &	GetProc () const
	Return MPI rank on which each point was found by FindPoints.

virtual const Vector &	GetReferencePosition () const
	Return reference coordinates for each point found by FindPoints.

virtual const Vector &	GetDist () const

virtual const Array< unsigned int > &	GetGSLIBElem () const

virtual const Vector &	GetGSLIBReferencePosition () const

void	GetAxisAlignedBoundingBoxes (Vector &aabb)

void	GetOrientedBoundingBoxes (DenseTensor &obbA, Vector &obbC, Vector &obbV)

virtual void	DistributePointInfoToOwningMPIRanks (Array< unsigned int > &recv_elem, Vector &recv_ref, Array< unsigned int > &recv_code)

virtual void	DistributeInterpolatedValues (const Vector &int_vals, const int vdim, const int ordering, Vector &field_out) const

Protected Attributes
bool	overset

unsigned int	u_meshid

Vector	distfint

Protected Attributes inherited from mfem::FindPointsGSLIB
Mesh *	mesh

Array< Mesh * >	mesh_split

Array< IntegrationRule * >	ir_split

Array< FiniteElementSpace * >	fes_rst_map

Array< GridFunction * >	gf_rst_map

FiniteElementCollection *	fec_map_lin

void *	fdataD

struct gslib::crystal *	cr

struct gslib::comm *	gsl_comm

int	dim

int	points_cnt

Array< unsigned int >	gsl_code

Array< unsigned int >	gsl_proc

Array< unsigned int >	gsl_elem

Array< unsigned int >	gsl_mfem_elem

Vector	gsl_mesh

Vector	gsl_ref

Vector	gsl_dist

Vector	gsl_mfem_ref

Array< unsigned int >	recv_proc

Array< unsigned int >	recv_index

bool	setupflag

double	default_interp_value

AvgType	avgtype

Array< int >	split_element_map

Array< int >	split_element_index

int	NE_split_total

int	mesh_points_cnt

double	bdr_tol

bool	gpu_to_cpu_fallback = false

struct {

bool setup_device = false

bool find_device = false

int local_hash_size

int dof1d

int dof1d_sol

int h_o_size

int h_nx

double newt_tol

struct gslib::crystal * cr

struct gslib::hash_data_3 * hash3

struct gslib::hash_data_2 * hash2

Vector bb

Vector wtend

Vector gll1d

Vector lagcoeff

Vector gll1d_sol

Vector lagcoeff_sol

Array< unsigned int > loc_hash_offset

Vector loc_hash_min

Vector loc_hash_fac

}	DEV

Additional Inherited Members
Public Types inherited from mfem::FindPointsGSLIB
enum	AvgType { NONE , ARITHMETIC , HARMONIC }

Protected Member Functions inherited from mfem::FindPointsGSLIB
virtual void	InterpolateH1 (const GridFunction &field_in, Vector &field_out)
	Use GSLIB for communication and interpolation.

virtual void	InterpolateGeneral (const GridFunction &field_in, Vector &field_out)

virtual void	SetupSplitMeshes ()

virtual void	SetupIntegrationRuleForSplitMesh (Mesh mesh, IntegrationRule irule, int order)

virtual void	GetNodalValues (const GridFunction *gf_in, Vector &node_vals)
	Get GridFunction value at the points expected by GSLIB.

virtual void	MapRefPosAndElemIndices ()

void	FindPointsLocal3 (const Vector &point_pos, int point_pos_ordering, Array< unsigned int > &gsl_code_dev_l, Array< unsigned int > &gsl_elem_dev_l, Vector &gsl_ref_l, Vector &gsl_dist_l, int npt)

void	FindPointsLocal2 (const Vector &point_pos, int point_pos_ordering, Array< unsigned int > &gsl_code_dev_l, Array< unsigned int > &gsl_elem_dev_l, Vector &gsl_ref_l, Vector &gsl_dist_l, int npt)

void	InterpolateLocal3 (const Vector &field_in, Array< int > &gsl_elem_dev_l, Vector &gsl_ref_l, Vector &field_out, int npt, int ncomp, int nel, int dof1dsol)

void	InterpolateLocal2 (const Vector &field_in, Array< int > &gsl_elem_dev_l, Vector &gsl_ref_l, Vector &field_out, int npt, int ncomp, int nel, int dof1dsol)

void	SetupDevice ()

void	FindPointsOnDevice (const Vector &point_pos, int point_pos_ordering=Ordering::byNODES)

void	InterpolateOnDevice (const Vector &field_in_evec, Vector &field_out, const int nel, const int ncomp, const int dof1dsol, const int ordering)

Detailed Description

OversetFindPointsGSLIB enables use of findpts for arbitrary number of overlapping grids.

The parameters in this class are the same as FindPointsGSLIB with the difference of additional inputs required to account for more than 1 mesh.

Definition at line 370 of file gslib.hpp.

Constructor & Destructor Documentation

◆ OversetFindPointsGSLIB() [1/2]

mfem::OversetFindPointsGSLIB::OversetFindPointsGSLIB ( )

inline

Definition at line 378 of file gslib.hpp.

◆ OversetFindPointsGSLIB() [2/2]

mfem::OversetFindPointsGSLIB::OversetFindPointsGSLIB ( MPI_Comm comm_ )

inline

Definition at line 382 of file gslib.hpp.

Member Function Documentation

◆ FindPoints()

void mfem::OversetFindPointsGSLIB::FindPoints	(	const Vector &	point_pos,
		Array< unsigned int > &	point_id,
		int	point_pos_ordering = Ordering::byNODES )

Searches positions given in physical space by point_pos. All output Arrays and Vectors are expected to have the correct size.

Parameters

[in]	point_pos	Positions to be found.
[in]	point_id	Index of the mesh that the point belongs to (corresponding to `meshid` in Setup).
[in]	point_pos_ordering	Ordering of the points: byNodes: (XXX...,YYY...,ZZZ) or byVDim: (XYZ,XYZ,....XYZ)

Definition at line 2477 of file gslib.cpp.

◆ Interpolate() [1/4]

void mfem::FindPointsGSLIB::Interpolate	(	const GridFunction &	field_in,
		Vector &	field_out )

virtual

Interpolation of field values at prescribed reference space positions.

Parameters

[in]	field_in	Function values that will be interpolated on the reference positions. Note: it is assumed that `field_in` is in H1 and in the same space as the mesh that was given to Setup().
[out]	field_out	Interpolated values. For points that are not found the value is set to default_interp_value.

Reimplemented from mfem::FindPointsGSLIB.

Definition at line 251 of file gslib.cpp.

◆ Interpolate() [2/4]

void mfem::OversetFindPointsGSLIB::Interpolate	(	const Vector &	point_pos,
		Array< unsigned int > &	point_id,
		const GridFunction &	field_in,
		Vector &	field_out,
		int	point_pos_ordering = Ordering::byNODES )

Search positions and interpolate

Definition at line 2557 of file gslib.cpp.

◆ Interpolate() [3/4]

void mfem::FindPointsGSLIB::Interpolate	(	const Vector &	point_pos,
		const GridFunction &	field_in,
		Vector &	field_out,
		int	point_pos_ordering = Ordering::byNODES )

Search positions and interpolate. The ordering (byNODES or byVDIM) of the output values in field_out corresponds to the ordering used in the input GridFunction field_in.

Definition at line 255 of file gslib.cpp.

◆ Interpolate() [4/4]

void mfem::FindPointsGSLIB::Interpolate	(	Mesh &	m,
		const Vector &	point_pos,
		const GridFunction &	field_in,
		Vector &	field_out,
		int	point_pos_ordering = Ordering::byNODES )

Setup FindPoints, search positions and interpolate. The ordering (byNODES or byVDIM) of the output values in field_out corresponds to the ordering used in the input GridFunction field_in.

Definition at line 261 of file gslib.cpp.

◆ Setup()

void mfem::OversetFindPointsGSLIB::Setup	(	Mesh &	m,
		const int	meshid,
		GridFunction *	gfmax = NULL,
		const double	bb_t = 0.1,
		const double	newt_tol = 1.0e-12,
		const int	npt_max = 256 )

Initializes the internal mesh in gslib, by sending the positions of the Gauss-Lobatto nodes of the input Mesh object m. Note: not tested with periodic meshes (L2). Note: the input mesh m must have Nodes set.

Parameters

[in]	m	Input mesh.
[in]	meshid	A unique # for each overlapping mesh. This id is used to make sure that points being searched are not looked for in the mesh that they belong to.
[in]	gfmax	(Optional) GridFunction in H1 that is used as a discriminator when one point is located in multiple meshes. The mesh that maximizes gfmax is chosen. For example, using the distance field based on the overlapping boundaries is helpful for convergence during Schwarz iterations.
[in]	bb_t	(Optional) Relative size of bounding box around each element.
[in]	newt_tol	(Optional) Newton tolerance for the gslib search methods.
[in]	npt_max	(Optional) Number of points for simultaneous iteration. This alters performance and memory footprint.

Definition at line 2382 of file gslib.cpp.

Member Data Documentation

◆ distfint

Vector mfem::OversetFindPointsGSLIB::distfint

protected

Definition at line 375 of file gslib.hpp.

◆ overset

bool mfem::OversetFindPointsGSLIB::overset

protected

Definition at line 373 of file gslib.hpp.

◆ u_meshid

unsigned int mfem::OversetFindPointsGSLIB::u_meshid

protected

Definition at line 374 of file gslib.hpp.

The documentation for this class was generated from the following files:

fem/gslib.hpp
fem/gslib.cpp

Public Member Functions

Protected Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ OversetFindPointsGSLIB() [1/2]

◆ OversetFindPointsGSLIB() [2/2]

Member Function Documentation

◆ FindPoints()

◆ Interpolate() [1/4]

◆ Interpolate() [2/4]

◆ Interpolate() [3/4]

◆ Interpolate() [4/4]

◆ Setup()

Member Data Documentation

◆ distfint

◆ overset

◆ u_meshid