vmml Namespace Reference

heavily inspired by boost::enable_if http://www.boost.org, file: boost/utility/enable_if.hpp, Copyright 2003 Jaakko Järvi, Jeremiah Willcock, Andrew Lumsdaine More...

Classes
class	AABB
	An axis-aligned bounding box. More...
struct	enable_if
struct	enable_if< false, T >
struct	except_here
class	exception
class	frustum
class	frustum_culler
	Helper class for OpenGL view frustum culling. More...
class	intersection
class	lowpass_filter
class	matrix
struct	set_to_zero_functor
struct	set_to_identity_functor
struct	matrix_equals_allow_inverted_rows
struct	matrix_equals_allow_inverted_columns
class	quaternion
class	vector

Typedefs
typedef AABB< float >	AABBf
typedef AABB< double >	AABBd
typedef frustum< float >	Frustumf
typedef frustum< double >	Frustumd
typedef frustum_culler< float >	frustum_cullerf
typedef frustum_culler< double >	frustum_cullerd
typedef vmml::matrix < 3, 3, double >	Matrix3d
	A 3x3 double matrix.
typedef vmml::matrix < 4, 4, double >	Matrix4d
	A 4x4 double matrix.
typedef vmml::matrix< 3, 3, float >	Matrix3f
	A 3x3 float matrix.
typedef vmml::matrix< 4, 4, float >	Matrix4f
	A 4x4 float matrix.
typedef quaternion< float >	quaternionf
typedef quaternion< double >	quaterniond
typedef vmml::vector< 2, int >	Vector2i
typedef vmml::vector< 3, int >	Vector3i
typedef vmml::vector< 4, int >	Vector4i
typedef vmml::vector< 2, unsigned >	Vector2ui
typedef vmml::vector< 3, unsigned >	Vector3ui
typedef vmml::vector< 4, unsigned >	Vector4ui
typedef vmml::vector< 3, double >	Vector3d
typedef vmml::vector< 4, double >	Vector4d
typedef vmml::vector< 2, float >	Vector2f
typedef vmml::vector< 3, float >	Vector3f
typedef vmml::vector< 4, float >	Vector4f
typedef vmml::vector< 3, uint8_t >	Vector3ub
typedef vmml::vector< 4, uint8_t >	Vector4ub

Enumerations
enum	Visibility { VISIBILITY_NONE = 0, VISIBILITY_PARTIAL = 1, VISIBILITY_FULL = 2 }

Functions
template<typename T >
std::ostream &	operator<< (std::ostream &os, const AABB< T > &aabb)
void	error_noexcept (const std::string &desc, const except_here &here)
template<size_t M, size_t N, typename T >
bool	equals (const matrix< M, N, T > &m0, const matrix< M, N, T > &m1, const T tolerance=std::numeric_limits< T >::epsilon())
template<size_t M, size_t N, typename T >
void	multiply (const matrix< M, N, T > &left, const matrix< M, N, T > &right, matrix< M, N, T > &result)
template<size_t M, size_t N, size_t P, typename T >
void	multiply (const matrix< M, P, T > &left, const matrix< P, N, T > &right, matrix< M, N, T > &result)
template<size_t M, size_t N, typename T >
enable_if< M==N >::type *	identity (matrix< M, N, T > &m)
template<typename T >
T	compute_determinant (const matrix< 1, 1, T > &matrix_)
template<typename T >
T	compute_determinant (const matrix< 2, 2, T > &matrix_)
template<typename T >
T	compute_determinant (const matrix< 3, 3, T > &m_)
template<typename T >
T	compute_determinant (const matrix< 4, 4, T > &m)
template<typename T >
bool	compute_inverse (const matrix< 2, 2, T > &m_, matrix< 2, 2, T > &inverse_, T tolerance_=std::numeric_limits< T >::epsilon())
template<typename T >
bool	compute_inverse (const matrix< 3, 3, T > &m_, matrix< 3, 3, T > &inverse_, T tolerance_=std::numeric_limits< T >::epsilon())
template<typename T >
bool	compute_inverse (const matrix< 4, 4, T > &m_, matrix< 4, 4, T > &inv_, T tolerance_=std::numeric_limits< T >::epsilon())
template<size_t M, size_t N, typename T >
matrix< N, M, T >	transpose (const matrix< M, N, T > &matrix_)
template<size_t M, size_t N, typename T >
bool	is_positive_definite (const matrix< M, N, T > &matrix_, const T limit=-1e12, typename enable_if< M==N &&M<=3 >::type *=0)
template<size_t M, size_t N, typename T >
void	qr_decompose_gram_schmidt (const matrix< M, N, T > &A_, matrix< M, M, T > &Q, matrix< N, N, T > &R)
template<typename T , typename U >
bool	equals (T arg0, U arg1)
template<size_t M, size_t N, typename T >
matrix< M, N, T >	create_translation (const vector< M-1, T > &arg)
template<typename T >
matrix< 4, 4, T >	create_translation (const vector< 3, T > &arg)
template<size_t M, size_t N, typename T >
void	apply_translation (matrix< M, N, T > &m, T t0, T t1, T t2)
template<size_t M, size_t N, typename T >
matrix< M, N, T >	create_rotation (T angle, const vector< M-1, T > &axis)
template<typename T >
matrix< 4, 4, T >	create_rotation (T angle, const vector< 3, T > &axis)
template<size_t M, size_t N, typename T >
void	apply_rotation (matrix< M, N, T > &m, T angle, T t0, T t1, T t2)
template<size_t M, size_t N, typename T >
matrix< M, N, T >	create_scaling (const vector< N-1, T > &arg)
template<typename T >
matrix< 4, 4, T > const	create_scaling (const vector< 3, T > &arg)
template<typename T >
matrix< 4, 4, T >	create_scaling (T arg)
template<size_t M, typename T >
void	zero (vector< M, T > &arg)
template<size_t M, typename T >
vector< M, T >	min (const vector< M, T > &arg0, const vector< M, T > &arg1)
template<size_t M, typename T >
vector< M, T >	max (const vector< M, T > &arg0, const vector< M, T > &arg1)
template<size_t M, typename T >
T	manhattan (const vector< M, T > &arg)
template<size_t M, typename T >
bool	equals (const vector< M, T > &a, const vector< M, T > &b)
template<size_t M, typename T >
static vector< M, T >	operator* (T factor, const vector< M, T > &vector_)
template<size_t M, typename T >
T	dot (const vector< M, T > &first, const vector< M, T > &second)
template<size_t M, typename T >
vector< M, T >	cross (const vector< 3, T > &a, const vector< 3, T > &b)
template<size_t M, typename T >
vector< M, T >	normalize (const vector< M, T > &vector_)
template<typename T >
vector< 4, T >	compute_plane (const vector< 3, T > &a, const vector< 3, T > &b, const vector< 3, T > &c)
std::ostream &	operator<< (std::ostream &os, const Visibility &v)

Variables
template<size_t O, size_t P>
enable_if< O<=M &&P<=N >::type *matrix< M, N, T > ::get_sub_matrix(matrix< O, P, T > &result, size_t row_offset, size_t col_offset) const {for(size_t row=0;row< O;++row){for(size_t col=0;col< P;++col){result.at(row, col)=at(row_offset+row, col_offset+col);}}return 0;}template< size_t M, size_t N, typename T >template < size_t O, size_t P >typename enable_if< O<=M &&P<=N >::type *matrix< M, N, T > ::set_sub_matrix(const matrix < O, P, T > &sub_matrix, size_t row_offset, size_t col_offset){for(size_t row=0;row < O;++row){for(size_t col=0;col < P;++col){at(row_offset+row, col_offset+col)=sub_matrix.at(row, col);}}return 0;}template < size_t M, size_t N, typename T >inline Tmatrix< M, N, T > ::det() const {return compute_determinant(*this);}template < size_t M, size_t N, typename T >template< size_t O, size_t P, typename TT >inline bool matrix< M, N, T >::inverse(matrix < O, P, TT > &inverse_, T tolerance, typename enable_if < M==N &&O==P &&O==M &&M > =2 &&M<=4, TT >::type *) const {return compute_inverse(*this, inverse_, tolerance);}template < size_t M, size_t N, typename T >template< size_t O, size_t P >typename enable_if< O==P &&M==N &&O==M &&M >=2 >::type *matrix< M, N, T > ::get_adjugate(matrix< O, P, T > &adjugate) const {get_cofactors(adjugate);adjugate=transpose(adjugate);return 0;}template< size_t M, size_t N, typename T >template < size_t O, size_t P >typename enable_if< O==P &&M==N &&O==M &&M >=2 >::type *matrix< M, N, T >::get_cofactors(matrix< O, P, T > &cofactors) const {matrix< M-1, N-1, T > minor_;const size_t _negate=1u;for(size_t row_index=0;row_index< M;++row_index){for(size_t col_index=0;col_index< N;++col_index){if((row_index+col_index)&_negate) cofactors(row_index, col_index)=-get_minor(minor_, row_index, col_index);else cofactors(row_index, col_index)=get_minor(minor_, row_index, col_index);}}return 0;}template< size_t M, size_t N, typename T >template< size_t O, size_t P >Tmatrix< M, N, T >::get_minor(matrix< O, P, T > &minor_, size_t row_to_cut, size_t col_to_cut, typename enable_if< O==M-1 &&P==N-1 &&M==N &&M >=2 >::type *) const {size_t row_offset=0;size_t col_offset=0;for(size_t row_index=0;row_index< M;++row_index){if(row_index==row_to_cut) row_offset=-1;else{for(size_t col_index=0;col_index< M;++col_index){if(col_index==col_to_cut) col_offset=-1;else minor_.at(row_index+row_offset, col_index+col_offset)=at(row_index, col_index);}col_offset=0;}}return compute_determinant(minor_);}template< size_t M, size_t N, typename T >template< typename TT >matrix< M, N, T > &matrix< M, N, T >::rotate(const TT angle_, const vector< M-1, T > &axis, typename enable_if< M==N &&M==4, TT >::type *){const T angle=static_cast< T > angle_);const T sine=std::sin(angle);const T cosine=std::cos(angle);const T _zero=0.0;const T one=1.0;const T two=2.0;array[0]=cosine+(one-cosine)*std::pow(axis.array[0], two);array[1]=(one-cosine)*axis.array[0]*axis.array[1]+sine *axis.array[2];array[2]=(one-cosine)*axis.array[0]*axis.array[2]-sine *axis.array[1];array[3]=_zero;array[4]=(one-cosine)*axis.array[0]*axis.array[1]-sine *axis.array[2];array[5]=cosine+(one-cosine)*std::pow(axis.array[1], two);array[6]=(one-cosine)*axis.array[1]*axis.array[2]+sine *axis.array[0];array[7]=_zero;array[8]=(one-cosine)*axis.array[0]*axis.array[2]+sine *axis.array[1];array[9]=(one-cosine)*axis.array[1]*axis.array[2]-sine *axis.array[0];array[10]=cosine+(one-cosine)*std::pow(axis.array[2], two);array[11]=_zero;array[12]=_zero;array[13]=_zero;array[14]=_zero;array[15]=one;return *this;}template< size_t M, size_t N, typename T >template< typename TT >matrix< M, N, T > &matrix< M, N, T >::rotate_x(const TT angle_, typename enable_if< M==N &&M==4, TT >::type *){const T angle=static_cast< T > angle_);const T sine=std::sin(angle);const T cosine=std::cos(angle);T tmp;tmp=array[4]*cosine+array[8]*sine;array[8]=-array[4]*sine+array[8]*cosine;array[4]=tmp;tmp=array[5]*cosine+array[9]*sine;array[9]=-array[5]*sine+array[9]*cosine;array[5]=tmp;tmp=array[6]*cosine+array[10]*sine;array[10]=-array[6]*sine+array[10]*cosine;array[6]=tmp;tmp=array[7]*cosine+array[11]*sine;array[11]=-array[7]*sine+array[11]*cosine;array[7]=tmp;return *this;}template< size_t M, size_t N, typename T >template< typename TT >matrix< M, N, T > &matrix< M, N, T >::rotate_y(const TT angle_, typename enable_if< M==N &&M==4, TT >::type *){const T angle=static_cast< T > angle_);const T sine=std::sin(angle);const T cosine=std::cos(angle);T tmp;tmp=array[0]*cosine-array[8]*sine;array[8]=array[0]*sine+array[8]*cosine;array[0]=tmp;tmp=array[1]*cosine-array[9]*sine;array[9]=array[1]*sine+array[9]*cosine;array[1]=tmp;tmp=array[2]*cosine-array[10]*sine;array[10]=array[2]*sine+array[10]*cosine;array[2]=tmp;tmp=array[3]*cosine-array[11]*sine;array[11]=array[3]*sine+array[11]*cosine;array[3]=tmp;return *this;}template< size_t M, size_t N, typename T >template< typename TT >matrix< M, N, T > &matrix< M, N, T >::rotate_z(const TT angle_, typename enable_if< M==N &&M==4, TT >::type *){const T angle=static_cast< T > angle_);const T sine=std::sin(angle);const T cosine=std::cos(angle);T tmp;tmp=array[0]*cosine+array[4]*sine;array[4]=-array[0]*sine+array[4]*cosine;array[0]=tmp;tmp=array[1]*cosine+array[5]*sine;array[5]=-array[1]*sine+array[5]*cosine;array[1]=tmp;tmp=array[2]*cosine+array[6]*sine;array[6]=-array[2]*sine+array[6]*cosine;array[2]=tmp;tmp=array[3]*cosine+array[7]*sine;array[7]=-array[3]*sine+array[7]*cosine;array[3]=tmp;return *this;}template< size_t M, size_t N, typename T >template< typename TT >matrix< M, N, T > &matrix< M, N, T >::pre_rotate_x(const TT angle_, typename enable_if< M==N &&M==4, TT >::type *){const T angle=static_cast< T > angle_);const T sine=std::sin(angle);const T cosine=std::cos(angle);T tmp;tmp=array[1];array[1]=array[1]*cosine+array[2]*sine;array[2]=tmp *-sine+array[2]*cosine;tmp=array[5];array[5]=array[5]*cosine+array[6]*sine;array[6]=tmp *-sine+array[6]*cosine;tmp=array[9];array[9]=array[9]*cosine+array[10]*sine;array[10]=tmp *-sine+array[10]*cosine;tmp=array[13];array[13]=array[13]*cosine+array[14]*sine;array[14]=tmp *-sine+array[14]*cosine;return *this;}template< size_t M, size_t N, typename T >template< typename TT >matrix< M, N, T > &matrix< M, N, T >::pre_rotate_y(const TT angle_, typename enable_if< M==N &&M==4, TT >::type *){const T angle=static_cast< T > angle_);const T sine=std::sin(angle);const T cosine=std::cos(angle);T tmp;tmp=array[0];array[0]=array[0]*cosine-array[2]*sine;array[2]=tmp *sine+array[2]*cosine;tmp=array[4];array[4]=array[4]*cosine-array[6]*sine;array[6]=tmp *sine+array[6]*cosine;tmp=array[8];array[8]=array[8]*cosine-array[10]*sine;array[10]=tmp *sine+array[10]*cosine;tmp=array[12];array[12]=array[12]*cosine-array[14]*sine;array[14]=tmp *sine+array[14]*cosine;return *this;}template< size_t M, size_t N, typename T >template< typename TT >matrix< M, N, T > &matrix< M, N, T >::pre_rotate_z(const TT angle_, typename enable_if< M==N &&M==4, TT >::type *){const T angle=static_cast< T > angle_);const T sine=std::sin(angle);const T cosine=std::cos(angle);T tmp;tmp=array[0];array[0]=array[0]*cosine+array[1]*sine;array[1]=tmp *-sine+array[1]*cosine;tmp=array[4];array[4]=array[4]*cosine+array[5]*sine;array[5]=tmp *-sine+array[5]*cosine;tmp=array[8];array[8]=array[8]*cosine+array[9]*sine;array[9]=tmp *-sine+array[9]*cosine;tmp=array[12];array[12]=array[12]*cosine+array[13]*sine;array[13]=tmp *-sine+array[13]*cosine;return *this;}template< size_t M, size_t N, typename T >template< typename TT >matrix< M, N, T > &matrix< M, N, T >::scale(const TT _scale[3], typename enable_if< M==N &&M==4, TT >::type *){const T scale0=static_cast< T > _scale[0]);const T scale1=static_cast< T > _scale[1]);const T scale2=static_cast< T > _scale[2]);array[0]*=scale0;array[1]*=scale0;array[2]*=scale0;array[3]*=scale0;array[4]*=scale1;array[5]*=scale1;array[6]*=scale1;array[7]*=scale1;array[8]*=scale2;array[9]*=scale2;array[10]*=scale2;array[11]*=scale2;return *this;}template< size_t M, size_t N, typename T >template< typename TT >matrix< M, N, T > &matrix< M, N, T >::scale(const TT x_, const T y_, const T z_, typename enable_if< M==N &&M==4, TT >::type *){const T _x=static_cast< T > x_);array[0]*=_x;array[1]*=_x;array[2]*=_x;array[3]*=_x;array[4]*=y_;array[5]*=y_;array[6]*=y_;array[7]*=y_;array[8]*=z_;array[9]*=z_;array[10]*=z_;array[11]*=z_;return *this;}template< size_t M, size_t N, typename T >template< typename TT >inline matrix< M, N, T > &matrix< M, N, T >::scale(const vector< 3, TT > &scale_, typename enable_if< M==N &&M==4, TT >::type *){return scale(scale_.array);}template< size_t M, size_t N, typename T >template< typename TT >matrix< M, N, T > &matrix< M, N, T >::scale_translation(const TT scale_[3], typename enable_if< M==N &&M==4, TT >::type *){array[12]*=static_cast< T > scale_[0]);array[13]*=static_cast< T > scale_[1]);array[14]*=static_cast< T > scale_[2]);return *this;}template< size_t M, size_t N, typename T >template< typename TT >inline matrix< M, N, T > &matrix< M, N, T >::scale_translation(const vector< 3, TT > &scale_, typename enable_if< M==N &&M==4, TT >::type *){return scale_translation(scale_.array);}template< size_t M, size_t N, typename T >template< typename TT >inline matrix< M, N, T > &matrix< M, N, T >::set_translation(const TT x_, const TT y_, const TT z_, typename enable_if< M==N &&M==4, TT >::type *){array[12]=static_cast< T > x_);array[13]=static_cast< T > y_);array[14]=static_cast< T > z_);return *this;}template< size_t M, size_t N, typename T >template< typename TT >inline matrix< M, N, T > &matrix< M, N, T >::set_translation(const TT trans[3], typename enable_if< M==N &&M==4, TT >::type *){array[12]=static_cast< T > trans[0]);array[13]=static_cast< T > trans[1]);array[14]=static_cast< T > trans[2]);return *this;}template< size_t M, size_t N, typename T >template< typename TT >inline matrix< M, N, T > &matrix< M, N, T >::set_translation(const vector< 3, TT > &translation_, typename enable_if< M==N &&M==4, TT >::type *){return set_translation(translation_.array);}template< size_t M, size_t N, typename T > template< typename TT > inlinevoid matrix< M, N, T >::get_translation(vector< N-1, TT > &translation_) const {for(size_t i=0;i< N-1;++i) translation_.array[i]=array[i+M *(N-1)];}template< size_t M, size_t N, typename T >inline vector< N-1, T > matrix< M, N, T >::get_translation() const {vector< N-1, T > result;get_translation(result);return result;}template< size_t M, size_t N, typename T >size_tmatrix< M, N, T >::size() const {return M *N;}template< size_t M, size_t N, typename T >typename matrix< M, N, T >::iteratormatrix< M, N, T >::begin(){return array;}template< size_t M, size_t N, typename T >typename matrix< M, N, T >::iteratormatrix< M, N, T >::end(){return array+size();}template< size_t M, size_t N, typename T >typename matrix< M, N, T >::const_iteratormatrix< M, N, T >::begin() const {return array;}template< size_t M, size_t N, typename T >typename matrix< M, N, T >::const_iteratormatrix< M, N, T >::end() const {return array+size();}template< size_t M, size_t N, typename T >typename matrix< M, N, T >::reverse_iteratormatrix< M, N, T >::rbegin(){return array+size()-1;}template< size_t M, size_t N, typename T >typename matrix< M, N, T >::reverse_iteratormatrix< M, N, T >::rend(){return array-1;}template< size_t M, size_t N, typename T >typename matrix< M, N, T >::const_reverse_iteratormatrix< M, N, T >::rbegin() const {return array+size()-1;}template< size_t M, size_t N, typename T >typename matrix< M, N, T >::const_reverse_iteratormatrix< M, N, T >::rend() const {return array-1;}template< size_t M, size_t N, typename T > template< typename init_functor_t >const matrix< M, N, T > matrix< M, N, T >::get_initialized_matrix(){matrix< M, N, T > matrix_;init_functor_t()(matrix_);return matrix_;}template< size_t M, size_t N, typename T >const matrix< M, N, T >matrix< M, N, T >::IDENTITY(matrix< M, N, T >::get_initialized_matrix< set_to_identity_functor< matrix< M, N, T > > >));template< size_t M, size_t N, typename T >const matrix< M, N, T >matrix< M, N, T >::ZERO(matrix< M, N, T >::get_initialized_matrix< set_to_zero_functor< matrix< M, N, T > > >));template< size_t M, size_t N, typename T >doublematrix< M, N, T >::frobenius_norm() const {double norm=0.0;const_iterator it=begin(), it_end=end();for(;it!=it_end;++it){norm+=*it **it;}return std::sqrt(norm);}template< size_t M, size_t N, typename T >doublematrix< M, N, T >::p_norm(double p) const {double norm=0.0;const_iterator it=begin(), it_end=end();for(;it!=it_end;++it){norm+=std::pow(*it, p);}return std::pow(norm, 1./p);}template< size_t M, size_t N, typename T >template< size_t O >voidmatrix< M, N, T >::khatri_rao_product(const matrix< O, N, T > &right_, matrix< M *O, N, T > &prod_) const {for(size_t col=0;col< N;++col){for(size_t m=0;m< M;++m){for(size_t o=0;o< O;++o){prod_.at(O *m+o, col)=at(m, col)*right_.at(o, col);}}}}template< size_t M, size_t N, typename T >template< size_t O, size_t P >voidmatrix< M, N, T >::kronecker_product(const matrix< O, P, T > &right_, matrix< M *O, N *P, T > &result_) const {for(size_t m=0;m< M;++m){for(size_t n=0;n< N;++n){for(size_t o=0;o< O;++o){for(size_t p=0;p< P;++p){result_.at(O *m+o, P *n+p)=at(m, n)*right_.at(o, p);}}}}}template< size_t M, size_t N, typename T >template< typename TT >voidmatrix< M, N, T >::cast_from(const matrix< M, N, TT > &other){typedef vmml::matrix< M, N, TT > matrix_tt_type;typedef typename matrix_tt_type::const_iterator tt_const_iterator;iterator it=begin(), it_end=end();tt_const_iterator other_it=other.begin();for(;it!=it_end;++it,++other_it){*it=static_cast< T > *	other_it

Detailed Description

heavily inspired by boost::enable_if http://www.boost.org, file: boost/utility/enable_if.hpp, Copyright 2003 Jaakko Järvi, Jeremiah Willcock, Andrew Lumsdaine