Matrix and vector functionality. More...

Collaboration diagram for Matrix/Vector Math:

Macros
#define	bn_mat_zero(_m)

#define	bn_mat_idn(_m)

#define	bn_mat_copy(_d, _s)

Functions
void	bn_mat_print (const char *title, const mat_t m)

void	bn_mat_print_guts (const char title, const mat_t m, char buf, int buflen)

void	bn_mat_print_vls (const char title, const mat_t m, struct bu_vls vls)

double	bn_atan2 (double x, double y)

void	bn_mat_mul (mat_t o, const mat_t a, const mat_t b)

void	bn_mat_mul2 (const mat_t i, mat_t o)

void	bn_mat_mul3 (mat_t o, const mat_t a, const mat_t b, const mat_t c)

void	bn_mat_mul4 (mat_t o, const mat_t a, const mat_t b, const mat_t c, const mat_t d)

void	bn_matXvec (hvect_t ov, const mat_t im, const hvect_t iv)

void	bn_mat_inv (mat_t output, const mat_t input)

int	bn_mat_inverse (mat_t output, const mat_t input)

void	bn_vtoh_move (hvect_t h, const vect_t v)

void	bn_htov_move (vect_t v, const hvect_t h)

void	bn_mat_trn (mat_t om, const mat_t im)

void	bn_mat_ae (mat_t m, double azimuth, double elev)

void	bn_ae_vec (fastf_t azp, fastf_t elp, const vect_t v)

void	bn_aet_vec (fastf_t az, fastf_t el, fastf_t *twist, vect_t vec_ae, vect_t vec_twist, fastf_t accuracy)

void	bn_vec_ae (vect_t vec, fastf_t az, fastf_t el)

void	bn_vec_aed (vect_t vec, fastf_t az, fastf_t el, fastf_t dist)

void	bn_mat_angles (mat_t mat, double alpha, double beta, double ggamma)

void	bn_mat_angles_rad (mat_t mat, double alpha, double beta, double ggamma)

void	bn_eigen2x2 (fastf_t val1, fastf_t val2, vect_t vec1, vect_t vec2, fastf_t a, fastf_t b, fastf_t c)

void	bn_vec_perp (vect_t new_vec, const vect_t old_vec)

void	bn_mat_fromto (mat_t m, const fastf_t from, const fastf_t to, const struct bn_tol *tol)

void	bn_mat_xrot (mat_t m, double sinx, double cosx)

void	bn_mat_yrot (mat_t m, double siny, double cosy)

void	bn_mat_zrot (mat_t m, double sinz, double cosz)

void	bn_mat_lookat (mat_t rot, const vect_t dir, int yflip)

void	bn_vec_ortho (vect_t out, const vect_t in)

int	bn_mat_scale_about_pnt (mat_t mat, const point_t pnt, const double scale)

void	bn_mat_xform_about_pnt (mat_t mat, const mat_t xform, const point_t pnt)

int	bn_mat_is_equal (const mat_t a, const mat_t b, const struct bn_tol *tol)

int	bn_mat_is_identity (const mat_t m)

void	bn_mat_arb_rot (mat_t m, const point_t pt, const vect_t dir, const fastf_t ang)

matp_t	bn_mat_dup (const mat_t in)

int	bn_mat_ck (const char *title, const mat_t m)

fastf_t	bn_mat_det3 (const mat_t m)

fastf_t	bn_mat_determinant (const mat_t m)

int	bn_mat_is_non_unif (const mat_t m)

void	bn_wrt_point_direc (mat_t out, const mat_t change, const mat_t in, const point_t point, const vect_t direc, const struct bn_tol *tol)

void	persp_mat (mat_t m, fastf_t fovy, fastf_t aspect, fastf_t near1, fastf_t far1, fastf_t backoff)

void	mike_persp_mat (fastf_t pmat, const fastf_t eye)

void	deering_persp_mat (fastf_t m, const fastf_t l, const fastf_t h, const fastf_t eye)

int	bn_opt_mat (struct bu_vls msg, size_t argc, const char argv, void set_var)

Variables
const mat_t	bn_mat_identity

Detailed Description

Matrix and vector functionality.

Macro Definition Documentation

◆ bn_mat_zero

#define bn_mat_zero ( _m )

Value:

        { \
	bu_log("%s:%d bn_mat_zero() is deprecated, use MAT_ZERO()\n", \
               __FILE__, __LINE__); \
        (_m)[0] = (_m)[1] = (_m)[2] = (_m)[3] = \
        (_m)[4] = (_m)[5] = (_m)[6] = (_m)[7] = \
        (_m)[8] = (_m)[9] = (_m)[10] = (_m)[11] = \
        (_m)[12] = (_m)[13] = (_m)[14] = (_m)[15] = 0.0; }

Definition at line 67 of file mat.h.

◆ bn_mat_idn

#define bn_mat_idn ( _m )

Value:

        { \
	bu_log("%s:%d bn_mat_idn() is deprecated, use MAT_IDN()\n", \
               __FILE__, __LINE__); \
        (_m)[1] = (_m)[2] = (_m)[3] = (_m)[4] = \
        (_m)[6] = (_m)[7] = (_m)[8] = (_m)[9] = \
        (_m)[11] = (_m)[12] = (_m)[13] = (_m)[14] = 0.0; \
        (_m)[0] = (_m)[5] = (_m)[10] = (_m)[15] = 1.0; }

Definition at line 77 of file mat.h.

◆ bn_mat_copy

#define bn_mat_copy	(	_d,
		_s
	)

Value:

        { \
	bu_log("%s:%d bn_mat_copy() is deprecated, use MAT_COPY()\n", \
               __FILE__, __LINE__); \
        (_d)[0] = (_s)[0];\
        (_d)[1] = (_s)[1];\
        (_d)[2] = (_s)[2];\
        (_d)[3] = (_s)[3];\
        (_d)[4] = (_s)[4];\
        (_d)[5] = (_s)[5];\
        (_d)[6] = (_s)[6];\
        (_d)[7] = (_s)[7];\
        (_d)[8] = (_s)[8];\
        (_d)[9] = (_s)[9];\
        (_d)[10] = (_s)[10];\
        (_d)[11] = (_s)[11];\
        (_d)[12] = (_s)[12];\
        (_d)[13] = (_s)[13];\
        (_d)[14] = (_s)[14];\
        (_d)[15] = (_s)[15]; }

Definition at line 88 of file mat.h.

Function Documentation

◆ bn_mat_print()

void bn_mat_print	(	const char *	title,
		const mat_t	m
	)

extern

◆ bn_mat_print_guts()

void bn_mat_print_guts	(	const char *	title,
		const mat_t	m,
		char *	buf,
		int	buflen
	)

extern

◆ bn_mat_print_vls()

void bn_mat_print_vls	(	const char *	title,
		const mat_t	m,
		struct bu_vls *	vls
	)

extern

◆ bn_atan2()

double bn_atan2	(	double	x,
		double	y
	)

extern

A wrapper for the system atan2(). On the Silicon Graphics, and perhaps on others, x==0 incorrectly returns infinity.

◆ bn_mat_mul()

void bn_mat_mul	(	mat_t	o,
		const mat_t	a,
		const mat_t	b
	)

extern

Multiply matrix "a" by "b" and store the result in "o".

This is different from multiplying "b" by "a" (most of the time!) Also, "o" must not be the same as either of the inputs.

◆ bn_mat_mul2()

void bn_mat_mul2	(	const mat_t	i,
		mat_t	o
	)

extern

o = i * o

A convenience wrapper for bn_mat_mul() to update a matrix in place. The argument ordering is confusing either way.

◆ bn_mat_mul3()

void bn_mat_mul3	(	mat_t	o,
		const mat_t	a,
		const mat_t	b,
		const mat_t	c
	)

extern

o = a * b * c

The output matrix may be the same as 'b' or 'c', but may not be 'a'.

◆ bn_mat_mul4()

void bn_mat_mul4	(	mat_t	o,
		const mat_t	a,
		const mat_t	b,
		const mat_t	c,
		const mat_t	d
	)

extern

o = a * b * c * d

The output matrix may be the same as any input matrix.

◆ bn_matXvec()

void bn_matXvec	(	hvect_t	ov,
		const mat_t	im,
		const hvect_t	iv
	)

extern

Multiply the matrix "im" by the vector "iv" and store the result in the vector "ov". Note this is post-multiply, and operates on 4-tuples. Use MAT4X3VEC() to operate on 3-tuples.

◆ bn_mat_inv()

void bn_mat_inv	(	mat_t	output,
		const mat_t	input
	)

extern

The matrix pointed at by "input" is inverted and stored in the area pointed at by "output".

Calls bu_bomb if matrix is singular.

◆ bn_mat_inverse()

int bn_mat_inverse	(	mat_t	output,
		const mat_t	input
	)

extern

The matrix pointed at by "input" is inverted and stored in the area pointed at by "output".

Invert a 4-by-4 matrix using Algorithm 120 from ACM. This is a modified Gauss-Jordan algorithm. Note: Inversion is done in place, with 3 work vectors.

Returns: 1 if OK.; 0 if matrix is singular.

◆ bn_vtoh_move()

void bn_vtoh_move	(	hvect_t	h,
		const vect_t	v
	)

extern

Takes a pointer to a [x, y, z] vector, and a pointer to space for a homogeneous vector [x, y, z, w], and builds [x, y, z, 1].

◆ bn_htov_move()

void bn_htov_move	(	vect_t	v,
		const hvect_t	h
	)

extern

Takes a pointer to [x, y, z, w], and converts it to an ordinary vector [x/w, y/w, z/w]. Optimization for the case of w==1 is performed.

Todo:: make tolerance configurable

◆ bn_mat_trn()

void bn_mat_trn	(	mat_t	om,
		const mat_t	im
	)

extern

◆ bn_mat_ae()

void bn_mat_ae	(	mat_t	m,
		double	azimuth,
		double	elev
	)

extern

Compute a 4x4 rotation matrix given Azimuth and Elevation.

Azimuth is +X, Elevation is +Z, both in degrees.

Formula due to Doug Gwyn, BRL.

◆ bn_ae_vec()

void bn_ae_vec	(	fastf_t *	azp,
		fastf_t *	elp,
		const vect_t	v
	)

extern

Find the azimuth and elevation angles that correspond to the direction (not including twist) given by a direction vector.

◆ bn_aet_vec()

void bn_aet_vec	(	fastf_t *	az,
		fastf_t *	el,
		fastf_t *	twist,
		vect_t	vec_ae,
		vect_t	vec_twist,
		fastf_t	accuracy
	)

extern

Find the azimuth, elevation, and twist from two vectors. Vec_ae is in the direction of view (+z in mged view) and vec_twist points to the viewers right (+x in mged view). Accuracy (degrees) is used to stabilize flutter between equivalent extremes of atan2(), and to snap twist to zero when elevation is near +/- 90

◆ bn_vec_ae()

void bn_vec_ae	(	vect_t	vec,
		fastf_t	az,
		fastf_t	el
	)

extern

Find a unit vector from the origin given azimuth and elevation.

◆ bn_vec_aed()

void bn_vec_aed	(	vect_t	vec,
		fastf_t	az,
		fastf_t	el,
		fastf_t	dist
	)

extern

Find a vector from the origin given azimuth, elevation, and distance.

◆ bn_mat_angles()

void bn_mat_angles	(	mat_t	mat,
		double	alpha,
		double	beta,
		double	ggamma
	)

extern

This routine builds a Homogeneous rotation matrix, given alpha, beta, and gamma as angles of rotation, in degrees.

Alpha is angle of rotation about the X axis, and is done third. Beta is angle of rotation about the Y axis, and is done second. Gamma is angle of rotation about Z axis, and is done first.

Todo:: make tolerance configurable

◆ bn_mat_angles_rad()

void bn_mat_angles_rad	(	mat_t	mat,
		double	alpha,
		double	beta,
		double	ggamma
	)

extern

This routine builds a Homogeneous rotation matrix, given alpha, beta, and gamma as angles of rotation, in radians.

Alpha is angle of rotation about the X axis, and is done third. Beta is angle of rotation about the Y axis, and is done second. Gamma is angle of rotation about Z axis, and is done first.

Todo:: make tolerance configurable

◆ bn_eigen2x2()

void bn_eigen2x2	(	fastf_t *	val1,
		fastf_t *	val2,
		vect_t	vec1,
		vect_t	vec2,
		fastf_t	a,
		fastf_t	b,
		fastf_t	c
	)

extern

Find the eigenvalues and eigenvectors of a symmetric 2x2 matrix. (a b) (b c)

The eigenvalue with the smallest absolute value is returned in val1, with its eigenvector in vec1.

◆ bn_vec_perp()

void bn_vec_perp	(	vect_t	new_vec,
		const vect_t	old_vec
	)

extern

Given a vector, create another vector which is perpendicular to it. The output vector will have unit length only if the input vector did.

Todo:: make tolerance configurable

◆ bn_mat_fromto()

void bn_mat_fromto	(	mat_t	m,
		const fastf_t *	from,
		const fastf_t *	to,
		const struct bn_tol *	tol
	)

extern

Given two vectors, compute a rotation matrix that will transform space by the angle between the two. There are many candidate matrices.

The input 'from' and 'to' vectors need not be unit length. MAT4X3VEC(to, m, from) is the identity that is created.

◆ bn_mat_xrot()

void bn_mat_xrot	(	mat_t	m,
		double	sinx,
		double	cosx
	)

extern

Given the sin and cos of an X rotation angle, produce the rotation matrix.

◆ bn_mat_yrot()

void bn_mat_yrot	(	mat_t	m,
		double	siny,
		double	cosy
	)

extern

Given the sin and cos of a Y rotation angle, produce the rotation matrix.

◆ bn_mat_zrot()

void bn_mat_zrot	(	mat_t	m,
		double	sinz,
		double	cosz
	)

extern

Given the sin and cos of a Z rotation angle, produce the rotation matrix.

◆ bn_mat_lookat()

void bn_mat_lookat	(	mat_t	rot,
		const vect_t	dir,
		int	yflip
	)

extern

Given a direction vector D of unit length, product a matrix which rotates that vector D onto the -Z axis. This matrix will be suitable for use as a "model2view" matrix.

XXX This routine will fail if the vector is already more or less aligned with the Z axis.

This is done in several steps.

1) Rotate D about Z to match +X axis.  Azimuth adjustment.
2) Rotate D about Y to match -Y axis.  Elevation adjustment.
3) Rotate D about Z to make projection of X axis again point
in the +X direction.  Twist adjustment.
4) Optionally, flip sign on Y axis if original Z becomes inverted.
This can be nice for static frames, but is astonishing when
used in animation.

◆ bn_vec_ortho()

void bn_vec_ortho	(	vect_t	out,
		const vect_t	in
	)

extern

Given a vector, create another vector which is perpendicular to it, and with unit length. This algorithm taken from Gift's arvec.f; a faster algorithm may be possible.

Todo:: make tolerance configurable

◆ bn_mat_scale_about_pnt()

int bn_mat_scale_about_pnt	(	mat_t	mat,
		const point_t	pnt,
		const double	scale
	)

extern

Build a matrix to scale uniformly around a given point.

Returns: -1 if scale is too small.; 0 if OK.

Todo:: make tolerance configurable

◆ bn_mat_xform_about_pnt()

void bn_mat_xform_about_pnt	(	mat_t	mat,
		const mat_t	xform,
		const point_t	pnt
	)

extern

Build a matrix to apply arbitrary 4x4 transformation around a given point.

◆ bn_mat_is_equal()

int bn_mat_is_equal	(	const mat_t	a,
		const mat_t	b,
		const struct bn_tol *	tol
	)

extern

Returns: 0 When matrices are not equal; 1 When matrices are equal

◆ bn_mat_is_identity()

int bn_mat_is_identity ( const mat_t m )

extern

This routine is intended for detecting identity matrices read in from ascii or binary files, where the numbers are pure ones or zeros. This routine is not intended for tolerance-based "near-zero" comparisons; as such, it shouldn't be used on matrices which are the result of calculation.

Returns: 0 non-identity; 1 a perfect identity matrix

◆ bn_mat_arb_rot()

void bn_mat_arb_rot	(	mat_t	m,
		const point_t	pt,
		const vect_t	dir,
		const fastf_t	ang
	)

extern

Construct a transformation matrix for rotation about an arbitrary axis. The axis is defined by a point (pt) and a unit direction vector (dir). The angle of rotation is "ang"

Todo:: make tolerance configurable

◆ bn_mat_dup()

matp_t bn_mat_dup ( const mat_t in )

extern

Return a pointer to a copy of the matrix in dynamically allocated memory.

◆ bn_mat_ck()

int bn_mat_ck	(	const char *	title,
		const mat_t	m
	)

extern

Check to ensure that a rotation matrix preserves axis perpendicularity. Note that not all matrices are rotation matrices.

Returns: -1 FAIL; 0 OK

◆ bn_mat_det3()

fastf_t bn_mat_det3 ( const mat_t m )

extern

Calculates the determinant of the 3X3 "rotation" part of the passed matrix.

◆ bn_mat_determinant()

fastf_t bn_mat_determinant ( const mat_t m )

extern

Calculates the determinant of the 4X4 matrix

◆ bn_mat_is_non_unif()

int bn_mat_is_non_unif ( const mat_t m )

extern

Todo:: make tolerance configurable

◆ bn_wrt_point_direc()

void bn_wrt_point_direc	(	mat_t	out,
		const mat_t	change,
		const mat_t	in,
		const point_t	point,
		const vect_t	direc,
		const struct bn_tol *	tol
	)

extern

Given a model-space transformation matrix "change", return a matrix which applies the change with-respect-to given "point" and "direc".

◆ persp_mat()

void persp_mat	(	mat_t	m,
		fastf_t	fovy,
		fastf_t	aspect,
		fastf_t	near1,
		fastf_t	far1,
		fastf_t	backoff
	)

extern

◆ mike_persp_mat()

void mike_persp_mat	(	fastf_t *	pmat,
		const fastf_t *	eye
	)

extern

◆ deering_persp_mat()

void deering_persp_mat	(	fastf_t *	m,
		const fastf_t *	l,
		const fastf_t *	h,
		const fastf_t *	eye
	)

extern

◆ bn_opt_mat()

int bn_opt_mat	(	struct bu_vls *	msg,
		size_t	argc,
		const char **	argv,
		void *	set_var
	)

extern

Variable Documentation

◆ bn_mat_identity

const mat_t bn_mat_identity

extern

Macros

Functions

Variables

Detailed Description

Macro Definition Documentation

◆ bn_mat_zero

◆ bn_mat_idn

◆ bn_mat_copy

Function Documentation

◆ bn_mat_print()

◆ bn_mat_print_guts()

◆ bn_mat_print_vls()

◆ bn_atan2()

◆ bn_mat_mul()

◆ bn_mat_mul2()

◆ bn_mat_mul3()

◆ bn_mat_mul4()

◆ bn_matXvec()

◆ bn_mat_inv()

◆ bn_mat_inverse()

◆ bn_vtoh_move()

◆ bn_htov_move()

◆ bn_mat_trn()

◆ bn_mat_ae()

◆ bn_ae_vec()

◆ bn_aet_vec()

◆ bn_vec_ae()

◆ bn_vec_aed()

◆ bn_mat_angles()

◆ bn_mat_angles_rad()

◆ bn_eigen2x2()

◆ bn_vec_perp()

◆ bn_mat_fromto()

◆ bn_mat_xrot()

◆ bn_mat_yrot()

◆ bn_mat_zrot()

◆ bn_mat_lookat()

◆ bn_vec_ortho()

◆ bn_mat_scale_about_pnt()

◆ bn_mat_xform_about_pnt()

◆ bn_mat_is_equal()

◆ bn_mat_is_identity()

◆ bn_mat_arb_rot()

◆ bn_mat_dup()

◆ bn_mat_ck()

◆ bn_mat_det3()

◆ bn_mat_determinant()

◆ bn_mat_is_non_unif()

◆ bn_wrt_point_direc()

◆ persp_mat()

◆ mike_persp_mat()

◆ deering_persp_mat()

◆ bn_opt_mat()

Variable Documentation

◆ bn_mat_identity