axom::primal Namespace Reference

Classes
class	BoundingBox
struct	NonChar
	Type trait to avoid outputting chars when a value is expected This avoids unintentionally outputting system beeps. More...
struct	NonChar< char >
struct	NonChar< unsigned char >
class	NumericArray
	A simple statically sized array of data with component-wise operators. More...
class	OrientedBoundingBox
	OrientedBoundingBox represents an oriented bounding box defined by its centroid, axes, and extents. More...
class	Plane
	Defines an oriented plane within a 2-D or 3-D Euclidean space and provides associated operators, such as, projection, signed distance, orientation, etc. More...
class	Point
	The point class represents a point, \( p \in \mathcal{R}^d \) . It provides access methods to set and query the point coordinates. More...
class	Polygon
	Represents a polygon defined by an array of points. More...
class	Ray
class	Segment
class	Sphere
	Defines an oriented Sphere in 2-D (i.e., a circle) or 3-D given by its center, \( \mathcal{X} \) and radius \( \mathcal{R} \). The Sphere object provides associated operations on a sphere, such as, signed distance and orientation. More...
class	Tetrahedron
	Represents a tetrahedral geometric shape defined by four points. More...
class	Triangle
	Represents a triangular geometric shape defined by three points. More...
class	Vector
	Represents a vector, \( v \in \mathcal{R}^d \). It provides access methods for setting and querying the vector components as well as vector math operators, e.g., adding, subtracting, dot_product and cross_product. More...

Typedefs
Pre-defined point types
typedef Point< double, 2 >	Point2D
typedef Point< double, 3 >	Point3D
Pre-defined Vector types
typedef Vector< double, 2 >	Vector2D
typedef Vector< double, 3 >	Vector3D

Enumerations
enum	OrientationResult { ON_BOUNDARY, ON_POSITIVE_SIDE, ON_NEGATIVE_SIDE }
	Enumerates possible return values for orientation tests. More...

Functions
template<typename T , int SIZE>
NumericArray< T, SIZE >	operator- (const NumericArray< T, SIZE > &lhs, const NumericArray< T, SIZE > &rhs)
	Performs component-wise subtraction of two numeric arrays. More...
template<typename T , int NDIMS>
std::ostream &	operator<< (std::ostream &os, const Polygon< T, NDIMS > &poly)
	Overloaded output operator for polygons. More...
template<typename T , int NDIMS>
std::ostream &	operator<< (std::ostream &os, const Ray< T, NDIMS > &ray)
	Overloaded output operator for rays. More...
template<typename T , int NDIMS>
std::ostream &	operator<< (std::ostream &os, const Segment< T, NDIMS > &seg)
	Overloaded output operator for Segment. More...
template<typename T , int NDIMS>
std::ostream &	operator<< (std::ostream &os, const Tetrahedron< T, NDIMS > &tet)
	Free functions implementing Tetrahedron's operators. More...
template<typename T , int NDIMS>
std::ostream &	operator<< (std::ostream &os, const Triangle< T, NDIMS > &tri)
	Overloaded output operator for triangles. More...
template<typename T >
Polygon< T, 3 >	clip (const Triangle< T, 3 > &tri, const BoundingBox< T, 3 > &bbox)
	Clips a 3D triangle against an axis-aligned bounding box in 3D. More...
template<typename T , int NDIMS>
Point< T, NDIMS >	closest_point (const Point< T, NDIMS > &P, const Triangle< T, NDIMS > &tri, int *loc=nullptr)
	Computes the closest point from a point, P, to a given triangle. More...
template<typename T , int NDIMS>
Point< T, NDIMS >	closest_point (const Point< T, NDIMS > &pt, const OrientedBoundingBox< T, NDIMS > &obb)
	Computes the closest point from a point to a given OBB. More...
template<typename T , int NDIMS>
OrientedBoundingBox< T, NDIMS >	merge_boxes (const OrientedBoundingBox< T, NDIMS > &l, const OrientedBoundingBox< T, NDIMS > &r)
	Creates an oriented bounding box which contains the passed in OBBs. More...
template<typename T , int NDIMS>
BoundingBox< T, NDIMS >	merge_boxes (const BoundingBox< T, NDIMS > &l, const BoundingBox< T, NDIMS > &r)
	Constructor. Creates a bounding box which contains the passed in bounding boxes. More...
template<typename T , int NDIMS>
AXOM_HOST_DEVICE BoundingBox< T, NDIMS >	compute_bounding_box (const Triangle< T, NDIMS > &tri)
	Creates a bounding box around a Triangle. More...
template<typename T >
bool	in_sphere (const Point< T, 2 > &q, const Point< T, 2 > &p0, const Point< T, 2 > &p1, const Point< T, 2 > &p2)
	Tests whether a query point lies inside a 2D triangle's circumcircle. More...
template<typename T >
bool	in_sphere (const Point< T, 3 > &q, const Point< T, 3 > &p0, const Point< T, 3 > &p1, const Point< T, 3 > &p2, const Point< T, 3 > &p3)
	Tests whether a query point lies inside a 3D tetrahedron's circumsphere. More...
template<typename T >
int	orientation (const Point< T, 3 > &p, const Triangle< T, 3 > &tri)
	Computes the orientation of the given point, p, with respect to a supplied oriented triangle. More...
template<typename T >
int	orientation (const Point< T, 2 > &p, const Segment< T, 2 > &seg)
	Computes the orientation of the given point, p, with respect to a supplied oriented segment. More...
double	squared_distance (const double *A, const double *B, int N)
	Computes the squared distance from point A to point B, represented by arrays of length N. More...
template<typename T , int NDIMS>
double	squared_distance (const Point< T, NDIMS > &A, const Point< T, NDIMS > &B)
	Computes the squared distance from point A to point B. More...
template<typename T , int NDIMS>
double	squared_distance (const Point< T, NDIMS > &P, const BoundingBox< T, NDIMS > &B)
	Computes the minimum squared distance from a query point, P, to a given axis-aligned bounding box B. More...
template<typename T , int NDIMS>
double	squared_distance (const Point< T, NDIMS > &P, const Segment< T, NDIMS > &S)
	Computes the minimum squared distance from a query point, P, to the given segment, S. More...
template<typename T , int NDIMS>
double	squared_distance (const Point< T, NDIMS > &P, const Triangle< T, NDIMS > &tri)
	Computes the minimum squared distance from a query point, P, to the closest point on the given triangle. More...
Forward Declared Overloaded Operators
template<typename T , int NDIMS>
bool	operator== (const BoundingBox< T, NDIMS > &lhs, const BoundingBox< T, NDIMS > &rhs)
	Equality comparison operator for bounding boxes. Two bounding boxes are equal when they have the same bounds. More...
template<typename T , int NDIMS>
bool	operator!= (const BoundingBox< T, NDIMS > &lhs, const BoundingBox< T, NDIMS > &rhs)
	Inequality comparison operator for bounding boxes. Two bounding boxes are unequal when they have different bounds. More...
template<typename T , int NDIMS>
std::ostream &	operator<< (std::ostream &os, const BoundingBox< T, NDIMS > &bb)
	Overloaded output operator for bounding boxes. More...
template<typename T , int SIZE>
bool	operator== (const NumericArray< T, SIZE > &lhs, const NumericArray< T, SIZE > &rhs)
	Checks if two numeric arrays are component-wise equal. More...
template<typename T , int SIZE>
bool	operator!= (const NumericArray< T, SIZE > &lhs, const NumericArray< T, SIZE > &rhs)
	Checks if two numeric arrays are not component-wise equal. More...
template<typename T , int SIZE>
NumericArray< T, SIZE >	operator+ (const NumericArray< T, SIZE > &lhs, const NumericArray< T, SIZE > &rhs)
	Performs component-wise addition of two numeric arrays. More...
template<typename T , int SIZE>
AXOM_HOST_DEVICE NumericArray< T, SIZE >	operator- (const NumericArray< T, SIZE > &lhs, const NumericArray< T, SIZE > &rhs)
	Performs component-wise subtraction of two numeric arrays. More...
template<typename T , int SIZE>
NumericArray< T, SIZE >	operator- (const NumericArray< T, SIZE > &arr)
	Unary negation of a numeric array instance. More...
template<typename T , int SIZE>
NumericArray< T, SIZE >	operator* (const NumericArray< T, SIZE > &arr, double scalar)
	Scalar multiplication a numeric array; Scalar on rhs. More...
template<typename T , int SIZE>
NumericArray< T, SIZE >	operator* (double scalar, const NumericArray< T, SIZE > &arr)
	Scalar multiplication a numeric array; Scalar on lhs. More...
template<typename T , int SIZE>
NumericArray< T, SIZE >	operator* (const NumericArray< T, SIZE > &lhs, const NumericArray< T, SIZE > &rhs)
	Component-wise multiplication of NumericArrays. More...
template<typename T , int SIZE>
NumericArray< T, SIZE >	operator/ (const NumericArray< T, SIZE > &lhs, const NumericArray< T, SIZE > &rhs)
	Component-wise division of NumericArrays. More...
template<typename T , int SIZE>
NumericArray< T, SIZE >	operator/ (const NumericArray< T, SIZE > &arr, double scalar)
	Scalar division of NumericArray; Scalar on rhs. More...
template<typename T , int SIZE>
NumericArray< T, SIZE >	abs (const NumericArray< T, SIZE > &arr)
	Coordinate-wise absolute value on the NumericArray. More...
template<typename T , int SIZE>
std::ostream &	operator<< (std::ostream &os, const NumericArray< T, SIZE > &arr)
	Overloaded output operator for numeric arrays. More...
template<typename T , int NDIMS>
bool	operator== (const OrientedBoundingBox< T, NDIMS > &lhs, const OrientedBoundingBox< T, NDIMS > &rhs)
	Equality comparison operator for oriented bounding boxes. More...
template<typename T , int NDIMS>
bool	operator!= (const OrientedBoundingBox< T, NDIMS > &lhs, const OrientedBoundingBox< T, NDIMS > &rhs)
	Inequality comparison operator for oriented bounding boxes. More...
template<typename T , int NDIMS>
std::ostream &	operator<< (std::ostream &os, const OrientedBoundingBox< T, NDIMS > &box)
	Overloaded output operator for oriented bounding boxes. More...
template<typename T , int NDIMS>
OrientedBoundingBox< T, NDIMS >	compute_oriented_bounding_box (const Point< T, NDIMS > *pts, int n)
	Creates a bounding box which contains the collection of passed in points. More...
template<typename T , int NDIMS>
std::ostream &	operator<< (std::ostream &os, const Plane< T, NDIMS > &p)
template<typename T , int NDIMS>
bool	operator== (const Point< T, NDIMS > &lhs, const Point< T, NDIMS > &rhs)
	Equality comparison operator for points. More...
template<typename T , int NDIMS>
bool	operator!= (const Point< T, NDIMS > &lhs, const Point< T, NDIMS > &rhs)
	Inequality comparison operator for points. More...
template<typename T , int NDIMS>
std::ostream &	operator<< (std::ostream &os, const Point< T, NDIMS > &pt)
	Overloaded output operator for points. More...
template<typename T , int NDIMS>
std::ostream &	operator<< (std::ostream &os, const Sphere< T, NDIMS > &s)
template<typename T , int NDIMS>
Vector< T, NDIMS >	operator+ (const Vector< T, NDIMS > &A, const Vector< T, NDIMS > &B)
	Adds vectors A, B and stores the result into a new vector C. More...
template<typename T , int NDIMS>
Vector< T, NDIMS >	operator- (const Vector< T, NDIMS > &A, const Vector< T, NDIMS > &B)
	Subtracts vectors A, B and stores the result into a new vector C. More...
template<typename T , int NDIMS>
Vector< T, NDIMS >	operator- (const Vector< T, NDIMS > &vec1)
	Unary negation of a vector instance. More...
template<typename T , int NDIMS>
Vector< T, NDIMS >	operator* (const Vector< T, NDIMS > &vec, const T scalar)
	Scalar multiplication of vector; Scalar on rhs. More...
template<typename T , int NDIMS>
Vector< T, NDIMS >	operator* (const T scalar, const Vector< T, NDIMS > &vec)
	Scalar multiplication of vector; Scalar on lhs. More...
template<typename T , int NDIMS>
Vector< T, NDIMS >	operator/ (const Vector< T, NDIMS > &vec, const T scalar)
	Scalar division of vector; Scalar on rhs. More...
template<typename T , int NDIMS>
std::ostream &	operator<< (std::ostream &os, const Vector< T, NDIMS > &vec)
	Overloaded output operator for vectors. More...
Triangle Intersection Routines
template<typename T >
AXOM_HOST_DEVICE bool	intersect (const Triangle< T, 3 > &t1, const Triangle< T, 3 > &t2, bool includeBoundary=false, double EPS=1E-08)
	Tests if 3D Triangles t1 and t2 intersect. More...
template<typename T >
bool	intersect (const Triangle< T, 2 > &t1, const Triangle< T, 2 > &t2, bool includeBoundary=false, double EPS=1E-08)
	Tests if 2D Triangles t1 and t2 intersect. More...
template<typename T >
bool	intersect (const Triangle< T, 3 > &tri, const BoundingBox< T, 3 > &bb)
	Determines if a triangle and a bounding box intersect. More...
template<typename T >
bool	intersect (const Triangle< T, 3 > &tri, const Ray< T, 3 > &ray)
	Determines if a 3D triangle intersects a 3D ray. More...
template<typename T >
bool	intersect (const Triangle< T, 3 > &tri, const Ray< T, 3 > &ray, T &t)
	Determines if a 3D triangle intersects a 3D ray. More...
template<typename T >
bool	intersect (const Triangle< T, 3 > &tri, const Ray< T, 3 > &ray, T &t, Point< double, 3 > &p)
	Determines if a 3D triangle intersects a 3D ray. More...
template<typename T >
bool	intersect (const Triangle< T, 3 > &tri, const Segment< T, 3 > &seg)
	Determines if a 3D triangle intersects a 3D segment. More...
template<typename T >
bool	intersect (const Triangle< T, 3 > &tri, const Segment< T, 3 > &seg, T &t)
	Determines if a 3D triangle intersects a 3D segment. More...
template<typename T >
bool	intersect (const Triangle< T, 3 > &tri, const Segment< T, 3 > &seg, T &t, Point< double, 3 > &p)
	Determines if a 3D triangle intersects a 3D segment. More...
Ray Intersection Routines
template<typename T >
bool	intersect (const Ray< T, 2 > &R, const Segment< T, 2 > &S, Point< T, 2 > &ip)
	Computes the intersection of the given ray, R, with the segment, S. More...
template<typename T , int DIM>
bool	intersect (const Ray< T, DIM > &R, const BoundingBox< T, DIM > &bb, Point< T, DIM > &ip)
	Computes the intersection of the given ray, R, with the Box, bb. More...
Segment Intersection Routines
template<typename T , int DIM>
bool	intersect (const Segment< T, DIM > &S, const BoundingBox< T, DIM > &bb, Point< T, DIM > &ip)
	Computes the intersection of the given segment, S, with the Box, bb. If an intersection is found, output parameter ip contains an intersection point. More...
Axis-Aligned Bounding Box Intersection Routines
template<typename T , int DIM>
bool	intersect (const BoundingBox< T, DIM > &bb1, const BoundingBox< T, DIM > &bb2)
	Determines if two axis aligned bounding boxes intersect. More...
Sphere Intersection Routines
template<typename T , int DIM>
bool	intersect (const Sphere< T, DIM > &s1, const Sphere< T, DIM > &s2, double TOL=1.e-9)
	Determines if two spheres intersect. More...
Oriented Bounding Box Intersection Routines
template<typename T >
bool	intersect (const OrientedBoundingBox< T, 1 > &b1, const OrientedBoundingBox< T, 1 > &b2)
template<typename T >
bool	intersect (const OrientedBoundingBox< T, 2 > &b1, const OrientedBoundingBox< T, 2 > &b2)
	Determines if a 2D OBB intersects a 2D OBB. More...
template<typename T >
bool	intersect (const OrientedBoundingBox< T, 3 > &b1, const OrientedBoundingBox< T, 3 > &b2, double EPS=1E-4)
	Determines if a 3D OBB intersects a 3D OBB. More...
Bezier Curve Intersection Routines
template<typename T , int NDIMS>
bool	intersect (const BezierCurve< T, NDIMS > &c1, const BezierCurve< T, NDIMS > &c2, std::vector< T > &sp, std::vector< T > &tp, double tol=1E-8)
	Tests if two Bezier Curves c1 and c2 intersect. More...

Typedef Documentation

Point2D

typedef Point<double, 2> axom::primal::Point2D

Point3D

typedef Point<double, 3> axom::primal::Point3D

Vector2D

typedef Vector<double, 2> axom::primal::Vector2D

Vector3D

typedef Vector<double, 3> axom::primal::Vector3D

Enumeration Type Documentation

OrientationResult

enum axom::primal::OrientationResult

Enumerates possible return values for orientation tests.

Enumerator
ON_BOUNDARY	primitive is on the boundary of a primitive
ON_POSITIVE_SIDE	primitive is on the positive side of a primitive
ON_NEGATIVE_SIDE	primitive is on the negative side of a primitive

Function Documentation

operator==() [1/4]

template<typename T , int NDIMS>

bool axom::primal::operator==	(	const BoundingBox< T, NDIMS > &	lhs,
		const BoundingBox< T, NDIMS > &	rhs
	)

Equality comparison operator for bounding boxes. Two bounding boxes are equal when they have the same bounds.

Free functions implementing comparison and arithmetic operators.

References axom::primal::BoundingBox< T, NDIMS >::getMax(), and axom::primal::BoundingBox< T, NDIMS >::getMin().

operator!=() [1/4]

template<typename T , int NDIMS>

bool axom::primal::operator!=	(	const BoundingBox< T, NDIMS > &	lhs,
		const BoundingBox< T, NDIMS > &	rhs
	)

Inequality comparison operator for bounding boxes. Two bounding boxes are unequal when they have different bounds.

operator<<() [1/12]

template<typename T , int NDIMS>

std::ostream & axom::primal::operator<<	(	std::ostream &	os,
		const BoundingBox< T, NDIMS > &	bb
	)

Overloaded output operator for bounding boxes.

operator==() [2/4]

template<typename T , int SIZE>

bool axom::primal::operator==	(	const NumericArray< T, SIZE > &	lhs,
		const NumericArray< T, SIZE > &	rhs
	)

Checks if two numeric arrays are component-wise equal.

Free functions implementing comparison and arithmetic operators.

Parameters

[in]	lhs	numeric array instance on the left-hand side.
[in]	rhs	numeric array instance on the right-hand side.

Returns

status true if lhs==rhs, otherwise, false.

operator!=() [2/4]

template<typename T , int SIZE>

bool axom::primal::operator!=	(	const NumericArray< T, SIZE > &	lhs,
		const NumericArray< T, SIZE > &	rhs
	)

Checks if two numeric arrays are not component-wise equal.

Parameters

[in]	lhs	numeric array instance on the left-hand side.
[in]	rhs	numeric array instance on the right-hand side.

Returns

status true if lhs!=rhs, otherwise, false.

operator+() [1/2]

template<typename T , int SIZE>

NumericArray< T, SIZE > axom::primal::operator+ ( const NumericArray< T, SIZE > &  lhs, const NumericArray< T, SIZE > &  rhs  )

inline

Performs component-wise addition of two numeric arrays.

Parameters

[in]	lhs	numeric array instance on the left-hand side.
[in]	rhs	numeric array instance on the right-hand side.

Returns

C resulting numeric array from the component-wise addition.

operator-() [1/5]

template<typename T , int SIZE>

AXOM_HOST_DEVICE NumericArray<T, SIZE> axom::primal::operator- ( const NumericArray< T, SIZE > &  lhs, const NumericArray< T, SIZE > &  rhs  )

inline

Performs component-wise subtraction of two numeric arrays.

Parameters

[in]	lhs	numeric array instance on the left-hand side.
[in]	rhs	numeric array instance on the right-hand side.

Returns

C resulting numeric array from component-wise subtraction.

operator-() [2/5]

template<typename T , int SIZE>

NumericArray< T, SIZE > axom::primal::operator- ( const NumericArray< T, SIZE > &  arr )

inline

Unary negation of a numeric array instance.

Parameters

[in]

arr

numeric array instance on the left-hand side.

Returns

C resulting numeric array from unary negation.

operator*() [1/5]

template<typename T , int SIZE>

NumericArray< T, SIZE > axom::primal::operator* ( const NumericArray< T, SIZE > &  arr, double  scalar  )

inline

Scalar multiplication a numeric array; Scalar on rhs.

Parameters

[in]	arr	numeric array instance.
[in]	scalar	user-supplied scalar.

Returns

C resutling numeric array, \( \ni: C_i = scalar*arr_i, \forall i\)

operator*() [2/5]

template<typename T , int SIZE>

NumericArray< T, SIZE > axom::primal::operator* ( double  scalar, const NumericArray< T, SIZE > &  arr  )

inline

Scalar multiplication a numeric array; Scalar on lhs.

Parameters

[in]	scalar	user-supplied scalar.
[in]	arr	numeric array instance.

Returns

C resulting numeric array, \( \ni: C_i = scalar*arr_i, \forall i\)

operator*() [3/5]

template<typename T , int SIZE>

NumericArray< T, SIZE > axom::primal::operator* ( const NumericArray< T, SIZE > &  lhs, const NumericArray< T, SIZE > &  rhs  )

inline

Component-wise multiplication of NumericArrays.

Parameters

[in]	lhs	numeric array instance on the left-hand side.
[in]	rhs	numeric array instance on the right-hand side.

Returns

C resulting numeric array, \( \ni: C_i = lhs_i * rhs_i, \forall i\)

operator/() [1/3]

template<typename T , int SIZE>

NumericArray< T, SIZE > axom::primal::operator/ ( const NumericArray< T, SIZE > &  lhs, const NumericArray< T, SIZE > &  rhs  )

inline

Component-wise division of NumericArrays.

Parameters

[in]	lhs	numeric array instance on the left-hand side.
[in]	rhs	numeric array instance on the right-hand side.

Returns

C resulting numeric array, \( \ni: C_i = lhs_i / rhs_i, \forall i\)

Precondition

\( rhs_i != 0.0, \forall i \)

operator/() [2/3]

template<typename T , int SIZE>

NumericArray< T, SIZE > axom::primal::operator/ ( const NumericArray< T, SIZE > &  arr, double  scalar  )

inline

Scalar division of NumericArray; Scalar on rhs.

Parameters

[in]	arr	numeric array instance
[in]	scalar	user-supplied scalar

Returns

C resulting numeric array, \( \ni: C_i = arr_i/scalar, \forall i\)

Precondition

scalar != 0.0

abs()

template<typename T , int SIZE>

NumericArray< T, SIZE > axom::primal::abs ( const NumericArray< T, SIZE > &  arr )

inline

Coordinate-wise absolute value on the NumericArray.

Parameters

[in]

arr

numeric array instance

Precondition

std::abs is defined for template type T

Returns

A NumericArray whose coordinates are the absolute value of arr

References axom::utilities::abs().

operator<<() [2/12]

template<typename T , int SIZE>

std::ostream & axom::primal::operator<<	(	std::ostream &	os,
		const NumericArray< T, SIZE > &	arr
	)

Overloaded output operator for numeric arrays.

Parameters

[in]	os	C++ output stream
[in]	arr	numeric array instance.

operator-() [3/5]

template<typename T , int SIZE>

NumericArray<T, SIZE> axom::primal::operator- ( const NumericArray< T, SIZE > &  lhs, const NumericArray< T, SIZE > &  rhs  )

inline

Performs component-wise subtraction of two numeric arrays.

Parameters

[in]	lhs	numeric array instance on the left-hand side.
[in]	rhs	numeric array instance on the right-hand side.

Returns

C resulting numeric array from component-wise subtraction.

operator==() [3/4]

template<typename T , int NDIMS>

bool axom::primal::operator==	(	const OrientedBoundingBox< T, NDIMS > &	lhs,
		const OrientedBoundingBox< T, NDIMS > &	rhs
	)

Equality comparison operator for oriented bounding boxes.

Free functions implementing comparison and arithmetic operators.

Two oriented bounding boxes are equal when they have the same bounds and axes.

References axom::primal::OrientedBoundingBox< T, NDIMS >::contains().

operator!=() [3/4]

template<typename T , int NDIMS>

bool axom::primal::operator!=	(	const OrientedBoundingBox< T, NDIMS > &	lhs,
		const OrientedBoundingBox< T, NDIMS > &	rhs
	)

Inequality comparison operator for oriented bounding boxes.

Two bounding boxes are unequal when they have different bounds

operator<<() [3/12]

template<typename T , int NDIMS>

std::ostream & axom::primal::operator<<	(	std::ostream &	os,
		const OrientedBoundingBox< T, NDIMS > &	box
	)

Overloaded output operator for oriented bounding boxes.

References axom::primal::OrientedBoundingBox< T, NDIMS >::print().

compute_oriented_bounding_box()

template<typename T , int NDIMS>

OrientedBoundingBox< T, NDIMS > axom::primal::compute_oriented_bounding_box	(	const Point< T, NDIMS > *	pts,
		int	n
	)

Creates a bounding box which contains the collection of passed in points.

Call srand() to initialize the random number generator before using this function.

Parameters

[in]	pts	array of points
[in]	n	number of points

Note

if n <= 0, invokes default constructor

operator<<() [4/12]

template<typename T , int NDIMS>

std::ostream & axom::primal::operator<<	(	std::ostream &	os,
		const Plane< T, NDIMS > &	p
	)

operator==() [4/4]

template<typename T , int NDIMS>

bool axom::primal::operator==	(	const Point< T, NDIMS > &	lhs,
		const Point< T, NDIMS > &	rhs
	)

Equality comparison operator for points.

operator!=() [4/4]

template<typename T , int NDIMS>

bool axom::primal::operator!=	(	const Point< T, NDIMS > &	lhs,
		const Point< T, NDIMS > &	rhs
	)

Inequality comparison operator for points.

Inequality operator for points.

References A, AXOM_HOST_DEVICE, B, axom::primal::Point< T, NDIMS >::lerp(), axom::primal::Point< T, NDIMS >::make_point(), axom::primal::Point< T, NDIMS >::midpoint(), and axom::primal::Point< T, NDIMS >::print().

operator<<() [5/12]

template<typename T , int NDIMS>

std::ostream & axom::primal::operator<<	(	std::ostream &	os,
		const Point< T, NDIMS > &	pt
	)

Overloaded output operator for points.

Free functions implementing Point's operators.

operator<<() [6/12]

template<typename T , int NDIMS>

std::ostream & axom::primal::operator<<	(	std::ostream &	os,
		const Polygon< T, NDIMS > &	poly
	)

Overloaded output operator for polygons.

Free functions implementing Polygon's operators.

operator<<() [7/12]

template<typename T , int NDIMS>

std::ostream & axom::primal::operator<<	(	std::ostream &	os,
		const Ray< T, NDIMS > &	ray
	)

Overloaded output operator for rays.

Free functions implementing Ray's operators.

References axom::primal::Vector< T, NDIMS >::print().

operator<<() [8/12]

template<typename T , int NDIMS>

std::ostream & axom::primal::operator<<	(	std::ostream &	os,
		const Segment< T, NDIMS > &	seg
	)

Overloaded output operator for Segment.

Free functions implementing Segments's operators.

operator<<() [9/12]

template<typename T , int NDIMS>

std::ostream & axom::primal::operator<<	(	std::ostream &	os,
		const Sphere< T, NDIMS > &	s
	)

operator<<() [10/12]

template<typename T , int NDIMS>

std::ostream& axom::primal::operator<<	(	std::ostream &	os,
		const Tetrahedron< T, NDIMS > &	tet
	)

Free functions implementing Tetrahedron's operators.

References axom::primal::Point< T, NDIMS >::print().

operator<<() [11/12]

template<typename T , int NDIMS>

std::ostream & axom::primal::operator<<	(	std::ostream &	os,
		const Triangle< T, NDIMS > &	tri
	)

Overloaded output operator for triangles.

Free functions implementing Triangle's operators.

operator+() [2/2]

template<typename T , int NDIMS>

Vector< T, NDIMS > axom::primal::operator+ ( const Vector< T, NDIMS > &  A, const Vector< T, NDIMS > &  B  )

inline

Adds vectors A, B and stores the result into a new vector C.

Parameters

[in]	A	vector on the left-hand side.
[in]	B	vector on the right-hand side.

Returns

C resulting vector, \( C_i = A_i + B_i \forall i \)

operator-() [4/5]

template<typename T , int NDIMS>

Vector< T, NDIMS > axom::primal::operator- ( const Vector< T, NDIMS > &  A, const Vector< T, NDIMS > &  B  )

inline

Subtracts vectors A, B and stores the result into a new vector C.

Parameters

[in]	A	vector on the left-hand side.
[in]	B	vector on the right-hand side.

Returns

C resulting vector, \( C_i = A_i - B_i \forall i \)

operator-() [5/5]

template<typename T , int NDIMS>

Vector< T, NDIMS > axom::primal::operator- ( const Vector< T, NDIMS > &  vec1 )

inline

Unary negation of a vector instance.

Parameters

[in]

vec1

vector instance to negate.

Returns

C resulting vector from unary negation.

References axom::primal::Vector< T, NDIMS >::negate().

operator*() [4/5]

template<typename T , int NDIMS>

Vector< T, NDIMS > axom::primal::operator* ( const Vector< T, NDIMS > &  vec, const T  scalar  )

inline

Scalar multiplication of vector; Scalar on rhs.

Free functions involving vectors.

Parameters

[in]	vec	vector instance.
[in]	scalar	user-supplied scalar.

Returns

C resulting vector, \( \ni: C_i = scalar*vec_i, \forall i\)

operator*() [5/5]

template<typename T , int NDIMS>

Vector< T, NDIMS > axom::primal::operator* ( const T  scalar, const Vector< T, NDIMS > &  vec  )

inline

Scalar multiplication of vector; Scalar on lhs.

Parameters

[in]	scalar	user-supplied scalar.
[in]	vec	vector instance.

Returns

C resulting vector, \( \ni: C_i = scalar*vec_i, \forall i\)

operator/() [3/3]

template<typename T , int NDIMS>

Vector< T, NDIMS > axom::primal::operator/ ( const Vector< T, NDIMS > &  vec, const T  scalar  )

inline

Scalar division of vector; Scalar on rhs.

Parameters

[in]	vec	vector instance
[in]

operator<<() [12/12]

template<typename T , int NDIMS>

std::ostream & axom::primal::operator<<	(	std::ostream &	os,
		const Vector< T, NDIMS > &	vec
	)

Overloaded output operator for vectors.

Parameters

[in]	os	C++ output stream
[in]	vec	vector instance

References axom::primal::Vector< T, NDIMS >::print().

clip()

template<typename T >

Polygon<T, 3> axom::primal::clip	(	const Triangle< T, 3 > &	tri,
		const BoundingBox< T, 3 > &	bbox
	)

Clips a 3D triangle against an axis-aligned bounding box in 3D.

Parameters

[in]	tri	The triangle to clip
[in]	bbox	The bounding box to clip against

Returns

A planar polygon of the triangle clipped against the bounding box. If the triangle is completely outside the bounding box, the returned polygon is empty (i.e. it has no vertices).

Note

Using a specialization of the Sutherland-Hodgeman clipping algorithm for axis aligned planes

References axom::primal::Polygon< T, NDIMS >::addVertex(), axom::primal::BoundingBox< T, NDIMS >::contains(), axom::primal::BoundingBox< T, NDIMS >::getMax(), axom::primal::BoundingBox< T, NDIMS >::getMin(), axom::primal::BoundingBox< T, NDIMS >::intersectsWith(), and axom::utilities::swap().

Referenced by axom::quest::InOutOctree< DIM >::generateIndex().

closest_point() [1/2]

template<typename T , int NDIMS>

Point<T, NDIMS> axom::primal::closest_point ( const Point< T, NDIMS > &  P, const Triangle< T, NDIMS > &  tri, int *  loc = nullptr  )

inline

Computes the closest point from a point, P, to a given triangle.

Parameters

[in]	P	the query point
[in]	tri	user-supplied triangle.
[out]	loc	int pointer to store location of closest point (optional).

Returns

cp the closest point from a point P and a triangle.

Note

If the optional int pointer is supplied for loc, the method returns the location of the closest point, which is illustrated in the schematic diagram below and encoded as follows:

• loc \( \in [0,2] \), loc corresponds to the triangle node index
• loc \( \in [-3,-1] \), abs(loc) corresponds to an edge
• loc >= 3, loc is on a triangle face

*
*            2
*           /\
*    (-3)--/  \--(-2)
*         /    \
*        /_ _ _ \
*       0   |    1
*           |
*         (-1)
*
* 

Precondition

NDIMS==2 || NDIMS==3

Note

Implementation is based on "Real Time Collision Detection, Chapter 5.1.5 Closest Point on Triangle to Point".

References A, B, C, and axom::primal::Vector< T, NDIMS >::dot_product().

Referenced by axom::quest::SignedDistance< NDIMS >::computeDistance(), and squared_distance().

closest_point() [2/2]

template<typename T , int NDIMS>

Point<T, NDIMS> axom::primal::closest_point ( const Point< T, NDIMS > &  pt, const OrientedBoundingBox< T, NDIMS > &  obb  )

inline

Computes the closest point from a point to a given OBB.

Parameters

[in]	pt	the query pt.
[in]	obb	user-supplied oriented bounding box.

Returns

cp the closest point from a point pt and an OBB.

References axom::primal::OrientedBoundingBox< T, NDIMS >::getAxes(), axom::primal::OrientedBoundingBox< T, NDIMS >::getCentroid(), axom::primal::OrientedBoundingBox< T, NDIMS >::getExtents(), and axom::primal::OrientedBoundingBox< T, NDIMS >::toLocal().

merge_boxes() [1/2]

template<typename T , int NDIMS>

OrientedBoundingBox<T, NDIMS> axom::primal::merge_boxes	(	const OrientedBoundingBox< T, NDIMS > &	l,
		const OrientedBoundingBox< T, NDIMS > &	r
	)

Creates an oriented bounding box which contains the passed in OBBs.

Call srand() to initialize the random number generator before using this function.

Parameters

[in]	l	left obb
[in]	r	right obb

References axom::primal::OrientedBoundingBox< T, NDIMS >::contains(), and axom::primal::OrientedBoundingBox< T, NDIMS >::vertices().

merge_boxes() [2/2]

template<typename T , int NDIMS>

BoundingBox<T, NDIMS> axom::primal::merge_boxes	(	const BoundingBox< T, NDIMS > &	l,
		const BoundingBox< T, NDIMS > &	r
	)

Constructor. Creates a bounding box which contains the passed in bounding boxes.

Parameters

[in]	l	left bb
[in]	r	right bb

References axom::primal::BoundingBox< T, NDIMS >::addBox().

compute_bounding_box()

template<typename T , int NDIMS>

AXOM_HOST_DEVICE BoundingBox<T, NDIMS> axom::primal::compute_bounding_box

(

const Triangle< T, NDIMS > &

tri

)

Creates a bounding box around a Triangle.

Parameters

[in]

tri

The Triangle

References axom::primal::BoundingBox< T, NDIMS >::addPoint().

Referenced by axom::quest::getMeshTriangle().

in_sphere() [1/2]

template<typename T >

bool axom::primal::in_sphere ( const Point< T, 2 > &  q, const Point< T, 2 > &  p0, const Point< T, 2 > &  p1, const Point< T, 2 > &  p2  )

inline

Tests whether a query point lies inside a 2D triangle's circumcircle.

A triangle's circumscircle is the unique circle (i.e. a 2-sphere) that passes through each of its three vertices.

Parameters

[in]	q	the query point
[in]	p0	the first vertex of the triangle
[in]	p1	the second vertex of the triangle
[in]	p2	the third vertex of the triangle

Returns

true if the point is inside the circumcircle, false if it is on the circle's boundary or outside the circle

References axom::numerics::determinant().

in_sphere() [2/2]

template<typename T >

bool axom::primal::in_sphere ( const Point< T, 3 > &  q, const Point< T, 3 > &  p0, const Point< T, 3 > &  p1, const Point< T, 3 > &  p2, const Point< T, 3 > &  p3  )

inline

Tests whether a query point lies inside a 3D tetrahedron's circumsphere.

A tetrahedron's circumsphere is the unique sphere that passes through each of its four vertices.

Parameters

[in]	q	the query point
[in]	p0	the first vertex of the tetrahedron
[in]	p1	the second vertex of the tetrahedron
[in]	p2	the third vertex of the tetrahedron
[in]	p3	the fourth vertex of the tetrahedron

Returns

true if the point is inside the circumsphere, false if it is on the sphere's boundary or outside the sphere

References axom::numerics::determinant().

intersect() [1/18]

template<typename T >

AXOM_HOST_DEVICE bool axom::primal::intersect	(	const Triangle< T, 3 > &	t1,
		const Triangle< T, 3 > &	t2,
		bool	includeBoundary = false,
		double	EPS = 1E-08
	)

Tests if 3D Triangles t1 and t2 intersect.

Parameters

[in]	t1	The first triangle
[in]	t2	The second triangle
[in]	includeBoundary	Indicates if boundaries should be considered when detecting intersections (default: false)
[in]	EPS	Tolerance for determining intersections (default: 1E-8)

Returns

status true iff t1 intersects with t2, otherwise, false.

If parameter includeBoundary is false (default), this function will return true if the interior of t1 intersects the interior of t2. To include triangle boundaries in intersections, specify includeBoundary as true.

Referenced by axom::quest::detail::InOutOctreeValidator< DIM >::checkBlockIndexingConsistency(), axom::quest::InOutOctree< DIM >::generateIndex(), and axom::quest::getMeshTriangle().

intersect() [2/18]

template<typename T >

bool axom::primal::intersect	(	const Triangle< T, 2 > &	t1,
		const Triangle< T, 2 > &	t2,
		bool	includeBoundary = false,
		double	EPS = 1E-08
	)

Tests if 2D Triangles t1 and t2 intersect.

Parameters

[in]	t1	The first triangle
[in]	t2	The second triangle
[in]	includeBoundary	Indicates if boundaries should be considered when detecting intersections (default: false)
[in]	EPS	Tolerance for determining intersections (default: 1E-8)

Returns

status true iff t1 intersects with t2, otherwise, false.

If parameter includeBoundary is false (default), this function will return true if the interior of t1 intersects the interior of t2. To include triangle boundaries in intersections, specify includeBoundary as true.

intersect() [3/18]

template<typename T >

bool axom::primal::intersect	(	const Triangle< T, 3 > &	tri,
		const BoundingBox< T, 3 > &	bb
	)

Determines if a triangle and a bounding box intersect.

Parameters

[in]	tri	user-supplied triangle (with three vertices).
[in]	bb	user-supplied axis aligned bounding box.

Returns

true iff tri intersects with bb, otherwise, false.

intersect() [4/18]

template<typename T >

bool axom::primal::intersect	(	const Triangle< T, 3 > &	tri,
		const Ray< T, 3 > &	ray
	)

Determines if a 3D triangle intersects a 3D ray.

Parameters

[in]	tri	A 3D triangle
[in]	ray	A 3D ray

Returns

true iff tri intersects with ray, otherwise, false.

intersect() [5/18]

template<typename T >

bool axom::primal::intersect	(	const Triangle< T, 3 > &	tri,
		const Ray< T, 3 > &	ray,
		T &	t
	)

Determines if a 3D triangle intersects a 3D ray.

Parameters

[in]	tri	A 3D triangle
[in]	ray	A 3D ray
[out]	t	Intersection point of tri and R, w.r.t. parametrization of R

Note

If there is an intersection, the intersection point is: R.at(t)

Returns

true iff tri intersects with ray, otherwise, false.

intersect() [6/18]

template<typename T >

bool axom::primal::intersect	(	const Triangle< T, 3 > &	tri,
		const Ray< T, 3 > &	ray,
		T &	t,
		Point< double, 3 > &	p
	)

Determines if a 3D triangle intersects a 3D ray.

Parameters

[in]	tri	A 3D triangle
[in]	ray	A 3D ray
[out]	t	Intersection point of tri and R, w.r.t. parametrization of R
[out]	p	Intersection point of tri and R, in barycentric coordinates relative to tri.

Note

If there is an intersection, the intersection point is: R.at(t)

Returns

true iff tri intersects with ray, otherwise, false.

Note

t and p only valid when function returns true

References axom::primal::Point< T, NDIMS >::array().

intersect() [7/18]

template<typename T >

bool axom::primal::intersect	(	const Triangle< T, 3 > &	tri,
		const Segment< T, 3 > &	seg
	)

Determines if a 3D triangle intersects a 3D segment.

Parameters

[in]	tri	A 3D triangle
[in]	seg	A 3D line segment

Returns

true iff tri intersects with seg, otherwise, false.

intersect() [8/18]

template<typename T >

bool axom::primal::intersect	(	const Triangle< T, 3 > &	tri,
		const Segment< T, 3 > &	seg,
		T &	t
	)

Determines if a 3D triangle intersects a 3D segment.

Parameters

[in]	tri	A 3D triangle
[in]	seg	A 3D line segment
[out]	t	Intersection point of tri and seg, w.r.t. seg's parametrization

Returns

true iff tri intersects with seg, otherwise, false.

intersect() [9/18]

template<typename T >

bool axom::primal::intersect	(	const Triangle< T, 3 > &	tri,
		const Segment< T, 3 > &	seg,
		T &	t,
		Point< double, 3 > &	p
	)

Determines if a 3D triangle intersects a 3D segment.

Parameters

[in]	tri	A 3D triangle
[in]	seg	A 3D line segment
[out]	t	Intersection point of tri and seg, w.r.t. seg's parametrization
[out]	p	Intersection point of tri and R, in barycentric coordinates relative to tri.

Note

If there is an intersection, the intersection point pt is: pt = seg.source() + t * ( seg.dest() - seg.target() )

Returns

true iff tri intersects with seg, otherwise, false.

Note

t and p only valid when function returns true

References axom::primal::Point< T, NDIMS >::array().

intersect() [10/18]

template<typename T >

bool axom::primal::intersect	(	const Ray< T, 2 > &	R,
		const Segment< T, 2 > &	S,
		Point< T, 2 > &	ip
	)

Computes the intersection of the given ray, R, with the segment, S.

Parameters

[in]	R	the specified ray
[in]	S	the segment to check
[out]	ip	the intersection point on S, valid only if status=true.

Returns

status true iff R intersects with S, otherwise, false.

See also

primal::Ray

primal::Segment

primal::Point

intersect() [11/18]

template<typename T , int DIM>

bool axom::primal::intersect	(	const Ray< T, DIM > &	R,
		const BoundingBox< T, DIM > &	bb,
		Point< T, DIM > &	ip
	)

Computes the intersection of the given ray, R, with the Box, bb.

Parameters

[in]	R	the specified ray
[in]	bb	the user-supplied axis-aligned bounding box
[out]	ip	the intersection point where R intersects bb.

Returns

status true iff bb intersects with R, otherwise, false.

See also

primal::Ray

primal::Segment

primal::BoundingBox

Note

Computes Ray Box intersection using the slab method from pg 180 of Real Time Collision Detection by Christer Ericson.

intersect() [12/18]

template<typename T , int DIM>

bool axom::primal::intersect	(	const Segment< T, DIM > &	S,
		const BoundingBox< T, DIM > &	bb,
		Point< T, DIM > &	ip
	)

Computes the intersection of the given segment, S, with the Box, bb. If an intersection is found, output parameter ip contains an intersection point.

Returns

status true iff bb intersects with S, otherwise, false.

Note

The intersection between segment S and box bb intersect, will, in general, be a along a (1D) subset of segment S. This function returns a single point of the intersection of S and bb found while determining if there is a valid intersection.

Computes Segment-Box intersection using the slab method from pg 180 of Real Time Collision Detection by Christer Ericson.

intersect() [13/18]

template<typename T , int DIM>

bool axom::primal::intersect	(	const BoundingBox< T, DIM > &	bb1,
		const BoundingBox< T, DIM > &	bb2
	)

Determines if two axis aligned bounding boxes intersect.

Parameters

[in]	bb1	user-supplied axis aligned bounding box.
[in]	bb2	user-supplied axis aligned bounding box.

Returns

true iff bb1 intersects with bb2, otherwise, false.

References axom::primal::BoundingBox< T, NDIMS >::intersectsWith().

intersect() [14/18]

template<typename T , int DIM>

bool axom::primal::intersect	(	const Sphere< T, DIM > &	s1,
		const Sphere< T, DIM > &	s2,
		double	TOL = 1.e-9
	)

Determines if two spheres intersect.

Parameters

[in]	s1	user-supplied sphere object to check for intersection.
[in]	s2	user-supplied sphere object to check for intersection.
[in]	TOL	tolerance used for intersection check (optional)

Note

If TOL is not supplied, the default is 1.e-9.

Returns

status true iff s1 intersects with s2, otherwise, false.

References axom::primal::Sphere< T, NDIMS >::intersectsWith().

intersect() [15/18]

template<typename T >

bool axom::primal::intersect	(	const OrientedBoundingBox< T, 1 > &	b1,
		const OrientedBoundingBox< T, 1 > &	b2
	)

intersect() [16/18]

template<typename T >

bool axom::primal::intersect	(	const OrientedBoundingBox< T, 2 > &	b1,
		const OrientedBoundingBox< T, 2 > &	b2
	)

Determines if a 2D OBB intersects a 2D OBB.

Parameters

[in]	b1	A 2D OrientedBoundingBox
[in]	b2	A 2D OrientedBoundingBox

Returns

true iff b1 intersects with b2, otherwise, false.

intersect() [17/18]

template<typename T >

bool axom::primal::intersect	(	const OrientedBoundingBox< T, 3 > &	b1,
		const OrientedBoundingBox< T, 3 > &	b2,
		double	EPS = 1E-4
	)

Determines if a 3D OBB intersects a 3D OBB.

Parameters

[in]	b1	A 3D OrientedBoundingBox
[in]	b2	A 3D OrientedBoundingBox
[in]	EPS	error tolerance for intersection

Returns

true iff b1 intersects with b2, otherwise, false.

intersect() [18/18]

template<typename T , int NDIMS>

bool axom::primal::intersect	(	const BezierCurve< T, NDIMS > &	c1,
		const BezierCurve< T, NDIMS > &	c2,
		std::vector< T > &	sp,
		std::vector< T > &	tp,
		double	tol = 1E-8
	)

Tests if two Bezier Curves c1 and c2 intersect.

Returns

status true iff c1 intersects c2, otherwise false.

Parameters

[in]	c1	the first BezierCurve, parametrized in [0,1)
[in]	c2	the second BezierCurve, parametrized in [0,1)
[out]	sp	vector of parameter space intersection points for c1
[out]	tp	vector of parameter space intersection points for c2
[in]	tol	tolerance parameter for determining if a curve can be approximated by a line segment.

Returns

True if the curves intersect, false otherwise. Intersection parameters are stored in sp and tp

Finds all intersection points between the two curves.

Note

This function assumes two dimensional curves in a plane.

This function assumes that the curves are in general position. Specifically, we assume that all intersections are at points and that the curves don't overlap.

This function assumes the all intersections have multiplicity one, i.e. there are no points at which the curves and their derivatives both intersect. Thus, the function does not find tangencies.

This function assumes that the curves are half-open, i.e. they contain their first endpoint, but not their last endpoint. Thus, the curves do not intersect at \( s==1 \) or at \( t==1 \).

orientation() [1/2]

template<typename T >

int axom::primal::orientation ( const Point< T, 3 > &  p, const Triangle< T, 3 > &  tri  )

inline

Computes the orientation of the given point, p, with respect to a supplied oriented triangle.

Parameters

[in]	p	the query point.
[in]	tri	oriented triangle.

Returns

The orientation of the point with respect to the given triangle.

Note

The triangle lies in a plane that divides space into the positive half-space and the negative half-space. The triangle's normal vector points into the positive half-space. The return value of this routine can be one of the following:

• ON_BOUNDARY, when the point is coplanar with the triangle
• ON_POSITIVE_SIDE, when the point lies in the positive half-space
• ON_NEGATIVE_SIDE, when the point lies in the negative half-space

References axom::numerics::determinant(), axom::utilities::isNearlyEqual(), ON_BOUNDARY, ON_NEGATIVE_SIDE, and ON_POSITIVE_SIDE.

orientation() [2/2]

template<typename T >

int axom::primal::orientation ( const Point< T, 2 > &  p, const Segment< T, 2 > &  seg  )

inline

Computes the orientation of the given point, p, with respect to a supplied oriented segment.

Parameters

[in]	p	the query point.
[in]	seg	the user-supplied segment.

Returns

The orientation of the point with respect to the given segment.

Note

The return value can be one of the following:

• ON_BOUNDARY, when the point is collinear with the points that define the segment.
• ON_POSITIVE_SIDE, when the point is clockwise, i.e., to the right of the directed segment.
• ON_NEGATIVE_SIDE, when the point is counter-clockwise, i.e., to the left of the directed segment.

References axom::numerics::determinant(), axom::utilities::isNearlyEqual(), ON_BOUNDARY, ON_NEGATIVE_SIDE, ON_POSITIVE_SIDE, axom::primal::Segment< T, NDIMS >::source(), and axom::primal::Segment< T, NDIMS >::target().

squared_distance() [1/5]

double axom::primal::squared_distance ( const double *  A, const double *  B, int  N  )

inline

Computes the squared distance from point A to point B, represented by arrays of length N.

Parameters

[in]	A	source point
[in]	B	end point.
[in]	N	length of A and B.

Returns

d the distance from point A to point B. If N < 1, return 0.

Precondition

A and B are arrays of at least length N.

Referenced by axom::quest::SignedDistance< NDIMS >::computeDistance(), axom::spin::BVHTree< T, NDIMS >::contains(), axom::quest::InOutOctree< DIM >::generateIndex(), and squared_distance().

squared_distance() [2/5]

template<typename T , int NDIMS>

double axom::primal::squared_distance ( const Point< T, NDIMS > &  A, const Point< T, NDIMS > &  B  )

inline

Computes the squared distance from point A to point B.

Parameters

[in]	A	source point
[in]	B	end point.

Returns

d the distance from point A to point B.

References axom::primal::Vector< T, NDIMS >::squared_norm().

squared_distance() [3/5]

template<typename T , int NDIMS>

double axom::primal::squared_distance ( const Point< T, NDIMS > &  P, const BoundingBox< T, NDIMS > &  B  )

inline

Computes the minimum squared distance from a query point, P, to a given axis-aligned bounding box B.

Parameters

[in]	P	the query point.
[in]	B	the axis-aligned bounding box.

Returns

d the signed distance from P to the closest point on B.

References axom::primal::BoundingBox< T, NDIMS >::contains(), axom::primal::BoundingBox< T, NDIMS >::getMax(), axom::primal::BoundingBox< T, NDIMS >::getMin(), and squared_distance().

squared_distance() [4/5]

template<typename T , int NDIMS>

double axom::primal::squared_distance ( const Point< T, NDIMS > &  P, const Segment< T, NDIMS > &  S  )

inline

Computes the minimum squared distance from a query point, P, to the given segment, S.

Parameters

[in]	P	the query point.
[in]	S	the input segment.

Returns

d the minimum distance from P on the

References axom::primal::Vector< T, NDIMS >::dot_product(), axom::primal::Segment< T, NDIMS >::source(), axom::primal::Vector< T, NDIMS >::squared_norm(), and axom::primal::Segment< T, NDIMS >::target().

squared_distance() [5/5]

template<typename T , int NDIMS>

double axom::primal::squared_distance ( const Point< T, NDIMS > &  P, const Triangle< T, NDIMS > &  tri  )

inline

Computes the minimum squared distance from a query point, P, to the closest point on the given triangle.

Parameters

[in]	P	the query point.
[in]	tri	the supplied triangle.

Returns

d the distance from Q to the closest point on the triangle T.

References closest_point(), and squared_distance().