746 lines
14 KiB
C++
746 lines
14 KiB
C++
///////////////////////////////////////////////////////////////////////////////////////////////////
|
|
// OpenGL Mathematics Copyright (c) 2005 - 2010 G-Truc Creation (www.g-truc.net)
|
|
///////////////////////////////////////////////////////////////////////////////////////////////////
|
|
// Created : 2008-08-03
|
|
// Updated : 2008-09-14
|
|
// Licence : This source is under MIT License
|
|
// File : glm/core/func_trigonometric.inl
|
|
///////////////////////////////////////////////////////////////////////////////////////////////////
|
|
|
|
namespace glm
|
|
{
|
|
namespace core{
|
|
namespace function{
|
|
namespace trigonometric{
|
|
|
|
// radians
|
|
template <typename genType>
|
|
inline genType radians
|
|
(
|
|
genType const & degrees
|
|
)
|
|
{
|
|
GLM_STATIC_ASSERT(detail::type<genType>::is_float);
|
|
|
|
const genType pi = genType(3.1415926535897932384626433832795);
|
|
return degrees * (pi / genType(180));
|
|
}
|
|
|
|
template <typename T>
|
|
inline detail::tvec2<T> radians
|
|
(
|
|
detail::tvec2<T> const & degrees
|
|
)
|
|
{
|
|
return detail::tvec2<T>(
|
|
radians(degrees.x),
|
|
radians(degrees.y));
|
|
}
|
|
|
|
template <typename T>
|
|
inline detail::tvec3<T> radians
|
|
(
|
|
detail::tvec3<T> const & degrees
|
|
)
|
|
{
|
|
return detail::tvec3<T>(
|
|
radians(degrees.x),
|
|
radians(degrees.y),
|
|
radians(degrees.z));
|
|
}
|
|
|
|
template <typename T>
|
|
inline detail::tvec4<T> radians
|
|
(
|
|
detail::tvec4<T> const & degrees
|
|
)
|
|
{
|
|
return detail::tvec4<T>(
|
|
radians(degrees.x),
|
|
radians(degrees.y),
|
|
radians(degrees.z),
|
|
radians(degrees.w));
|
|
}
|
|
|
|
// degrees
|
|
template <typename genType>
|
|
inline genType degrees
|
|
(
|
|
genType const & radians
|
|
)
|
|
{
|
|
GLM_STATIC_ASSERT(detail::type<genType>::is_float);
|
|
|
|
const genType pi = genType(3.1415926535897932384626433832795);
|
|
return radians * (genType(180) / pi);
|
|
}
|
|
|
|
template <typename T>
|
|
inline detail::tvec2<T> degrees
|
|
(
|
|
detail::tvec2<T> const & radians
|
|
)
|
|
{
|
|
return detail::tvec2<T>(
|
|
degrees(radians.x),
|
|
degrees(radians.y));
|
|
}
|
|
|
|
template <typename T>
|
|
inline detail::tvec3<T> degrees
|
|
(
|
|
detail::tvec3<T> const & radians
|
|
)
|
|
{
|
|
return detail::tvec3<T>(
|
|
degrees(radians.x),
|
|
degrees(radians.y),
|
|
degrees(radians.z));
|
|
}
|
|
|
|
template <typename T>
|
|
inline detail::tvec4<T> degrees
|
|
(
|
|
detail::tvec4<T> const & radians
|
|
)
|
|
{
|
|
return detail::tvec4<T>(
|
|
degrees(radians.x),
|
|
degrees(radians.y),
|
|
degrees(radians.z),
|
|
degrees(radians.w));
|
|
}
|
|
|
|
// sin
|
|
template <typename genType>
|
|
inline genType sin
|
|
(
|
|
genType const & angle
|
|
)
|
|
{
|
|
GLM_STATIC_ASSERT(detail::type<genType>::is_float);
|
|
|
|
return ::std::sin(angle);
|
|
}
|
|
|
|
template <typename T>
|
|
inline detail::tvec2<T> sin
|
|
(
|
|
detail::tvec2<T> const & angle
|
|
)
|
|
{
|
|
return detail::tvec2<T>(
|
|
sin(angle.x),
|
|
sin(angle.y));
|
|
}
|
|
|
|
template <typename T>
|
|
inline detail::tvec3<T> sin
|
|
(
|
|
detail::tvec3<T> const & angle
|
|
)
|
|
{
|
|
return detail::tvec3<T>(
|
|
sin(angle.x),
|
|
sin(angle.y),
|
|
sin(angle.z));
|
|
}
|
|
|
|
template <typename T>
|
|
inline detail::tvec4<T> sin
|
|
(
|
|
detail::tvec4<T> const & angle
|
|
)
|
|
{
|
|
return detail::tvec4<T>(
|
|
sin(angle.x),
|
|
sin(angle.y),
|
|
sin(angle.z),
|
|
sin(angle.w));
|
|
}
|
|
|
|
// cos
|
|
template <typename genType>
|
|
inline genType cos(genType const & angle)
|
|
{
|
|
GLM_STATIC_ASSERT(detail::type<genType>::is_float);
|
|
|
|
return ::std::cos(angle);
|
|
}
|
|
|
|
template <typename T>
|
|
inline detail::tvec2<T> cos
|
|
(
|
|
detail::tvec2<T> const & angle
|
|
)
|
|
{
|
|
return detail::tvec2<T>(
|
|
cos(angle.x),
|
|
cos(angle.y));
|
|
}
|
|
|
|
template <typename T>
|
|
inline detail::tvec3<T> cos
|
|
(
|
|
detail::tvec3<T> const & angle
|
|
)
|
|
{
|
|
return detail::tvec3<T>(
|
|
cos(angle.x),
|
|
cos(angle.y),
|
|
cos(angle.z));
|
|
}
|
|
|
|
template <typename T>
|
|
inline detail::tvec4<T> cos
|
|
(
|
|
detail::tvec4<T> const & angle
|
|
)
|
|
{
|
|
return detail::tvec4<T>(
|
|
cos(angle.x),
|
|
cos(angle.y),
|
|
cos(angle.z),
|
|
cos(angle.w));
|
|
}
|
|
|
|
// tan
|
|
template <typename genType>
|
|
inline genType tan
|
|
(
|
|
genType const & angle
|
|
)
|
|
{
|
|
GLM_STATIC_ASSERT(detail::type<genType>::is_float);
|
|
|
|
return ::std::tan(angle);
|
|
}
|
|
|
|
template <typename T>
|
|
inline detail::tvec2<T> tan
|
|
(
|
|
detail::tvec2<T> const & angle
|
|
)
|
|
{
|
|
return detail::tvec2<T>(
|
|
tan(angle.x),
|
|
tan(angle.y));
|
|
}
|
|
|
|
template <typename T>
|
|
inline detail::tvec3<T> tan
|
|
(
|
|
detail::tvec3<T> const & angle
|
|
)
|
|
{
|
|
return detail::tvec3<T>(
|
|
tan(angle.x),
|
|
tan(angle.y),
|
|
tan(angle.z));
|
|
}
|
|
|
|
template <typename T>
|
|
inline detail::tvec4<T> tan
|
|
(
|
|
detail::tvec4<T> const & angle
|
|
)
|
|
{
|
|
return detail::tvec4<T>(
|
|
tan(angle.x),
|
|
tan(angle.y),
|
|
tan(angle.z),
|
|
tan(angle.w));
|
|
}
|
|
|
|
// asin
|
|
template <typename genType>
|
|
inline genType asin
|
|
(
|
|
genType const & x
|
|
)
|
|
{
|
|
GLM_STATIC_ASSERT(detail::type<genType>::is_float);
|
|
|
|
return ::std::asin(x);
|
|
}
|
|
|
|
template <typename T>
|
|
inline detail::tvec2<T> asin
|
|
(
|
|
detail::tvec2<T> const & x
|
|
)
|
|
{
|
|
return detail::tvec2<T>(
|
|
asin(x.x),
|
|
asin(x.y));
|
|
}
|
|
|
|
template <typename T>
|
|
inline detail::tvec3<T> asin
|
|
(
|
|
detail::tvec3<T> const & x
|
|
)
|
|
{
|
|
return detail::tvec3<T>(
|
|
asin(x.x),
|
|
asin(x.y),
|
|
asin(x.z));
|
|
}
|
|
|
|
template <typename T>
|
|
inline detail::tvec4<T> asin
|
|
(
|
|
detail::tvec4<T> const & x
|
|
)
|
|
{
|
|
return detail::tvec4<T>(
|
|
asin(x.x),
|
|
asin(x.y),
|
|
asin(x.z),
|
|
asin(x.w));
|
|
}
|
|
|
|
// acos
|
|
template <typename genType>
|
|
inline genType acos
|
|
(
|
|
genType const & x
|
|
)
|
|
{
|
|
GLM_STATIC_ASSERT(detail::type<genType>::is_float);
|
|
|
|
return ::std::acos(x);
|
|
}
|
|
|
|
template <typename T>
|
|
inline detail::tvec2<T> acos
|
|
(
|
|
detail::tvec2<T> const & x
|
|
)
|
|
{
|
|
return detail::tvec2<T>(
|
|
acos(x.x),
|
|
acos(x.y));
|
|
}
|
|
|
|
template <typename T>
|
|
inline detail::tvec3<T> acos
|
|
(
|
|
detail::tvec3<T> const & x
|
|
)
|
|
{
|
|
return detail::tvec3<T>(
|
|
acos(x.x),
|
|
acos(x.y),
|
|
acos(x.z));
|
|
}
|
|
|
|
template <typename T>
|
|
inline detail::tvec4<T> acos
|
|
(
|
|
detail::tvec4<T> const & x
|
|
)
|
|
{
|
|
return detail::tvec4<T>(
|
|
acos(x.x),
|
|
acos(x.y),
|
|
acos(x.z),
|
|
acos(x.w));
|
|
}
|
|
|
|
// atan
|
|
template <typename genType>
|
|
inline genType atan
|
|
(
|
|
genType const & y,
|
|
genType const & x
|
|
)
|
|
{
|
|
GLM_STATIC_ASSERT(detail::type<genType>::is_float);
|
|
|
|
return ::std::atan2(y, x);
|
|
}
|
|
|
|
template <typename T>
|
|
inline detail::tvec2<T> atan
|
|
(
|
|
detail::tvec2<T> const & y,
|
|
detail::tvec2<T> const & x
|
|
)
|
|
{
|
|
return detail::tvec2<T>(
|
|
atan(y.x, x.x),
|
|
atan(y.y, x.y));
|
|
}
|
|
|
|
template <typename T>
|
|
inline detail::tvec3<T> atan
|
|
(
|
|
detail::tvec3<T> const & y,
|
|
detail::tvec3<T> const & x
|
|
)
|
|
{
|
|
return detail::tvec3<T>(
|
|
atan(y.x, x.x),
|
|
atan(y.y, x.y),
|
|
atan(y.z, x.z));
|
|
}
|
|
|
|
template <typename T>
|
|
inline detail::tvec4<T> atan
|
|
(
|
|
detail::tvec4<T> const & y,
|
|
detail::tvec4<T> const & x
|
|
)
|
|
{
|
|
return detail::tvec4<T>(
|
|
atan(y.x, x.x),
|
|
atan(y.y, x.y),
|
|
atan(y.z, x.z),
|
|
atan(y.w, x.w));
|
|
}
|
|
|
|
template <typename genType>
|
|
inline genType atan
|
|
(
|
|
genType const & x
|
|
)
|
|
{
|
|
GLM_STATIC_ASSERT(detail::type<genType>::is_float);
|
|
|
|
return ::std::atan(x);
|
|
}
|
|
|
|
template <typename T>
|
|
inline detail::tvec2<T> atan
|
|
(
|
|
detail::tvec2<T> const & x
|
|
)
|
|
{
|
|
return detail::tvec2<T>(
|
|
atan(x.x),
|
|
atan(x.y));
|
|
}
|
|
|
|
template <typename T>
|
|
inline detail::tvec3<T> atan
|
|
(
|
|
detail::tvec3<T> const & x
|
|
)
|
|
{
|
|
return detail::tvec3<T>(
|
|
atan(x.x),
|
|
atan(x.y),
|
|
atan(x.z));
|
|
}
|
|
|
|
template <typename T>
|
|
inline detail::tvec4<T> atan
|
|
(
|
|
detail::tvec4<T> const & x
|
|
)
|
|
{
|
|
return detail::tvec4<T>(
|
|
atan(x.x),
|
|
atan(x.y),
|
|
atan(x.z),
|
|
atan(x.w));
|
|
}
|
|
|
|
// sinh
|
|
template <typename genType>
|
|
inline genType sinh
|
|
(
|
|
genType const & angle
|
|
)
|
|
{
|
|
GLM_STATIC_ASSERT(detail::type<genType>::is_float);
|
|
|
|
return std::sinh(angle);
|
|
}
|
|
|
|
template <typename T>
|
|
inline detail::tvec2<T> sinh
|
|
(
|
|
detail::tvec2<T> const & angle
|
|
)
|
|
{
|
|
return detail::tvec2<T>(
|
|
sinh(angle.x),
|
|
sinh(angle.y));
|
|
}
|
|
|
|
template <typename T>
|
|
inline detail::tvec3<T> sinh
|
|
(
|
|
detail::tvec3<T> const & angle
|
|
)
|
|
{
|
|
return detail::tvec3<T>(
|
|
sinh(angle.x),
|
|
sinh(angle.y),
|
|
sinh(angle.z));
|
|
}
|
|
|
|
template <typename T>
|
|
inline detail::tvec4<T> sinh
|
|
(
|
|
detail::tvec4<T> const & angle
|
|
)
|
|
{
|
|
return detail::tvec4<T>(
|
|
sinh(angle.x),
|
|
sinh(angle.y),
|
|
sinh(angle.z),
|
|
sinh(angle.w));
|
|
}
|
|
|
|
// cosh
|
|
template <typename genType>
|
|
inline genType cosh
|
|
(
|
|
genType const & angle
|
|
)
|
|
{
|
|
GLM_STATIC_ASSERT(detail::type<genType>::is_float);
|
|
|
|
return std::cosh(angle);
|
|
}
|
|
|
|
template <typename T>
|
|
inline detail::tvec2<T> cosh
|
|
(
|
|
detail::tvec2<T> const & angle
|
|
)
|
|
{
|
|
return detail::tvec2<T>(
|
|
cosh(angle.x),
|
|
cosh(angle.y));
|
|
}
|
|
|
|
template <typename T>
|
|
inline detail::tvec3<T> cosh
|
|
(
|
|
detail::tvec3<T> const & angle
|
|
)
|
|
{
|
|
return detail::tvec3<T>(
|
|
cosh(angle.x),
|
|
cosh(angle.y),
|
|
cosh(angle.z));
|
|
}
|
|
|
|
template <typename T>
|
|
inline detail::tvec4<T> cosh
|
|
(
|
|
detail::tvec4<T> const & angle
|
|
)
|
|
{
|
|
return detail::tvec4<T>(
|
|
cosh(angle.x),
|
|
cosh(angle.y),
|
|
cosh(angle.z),
|
|
cosh(angle.w));
|
|
}
|
|
|
|
// tanh
|
|
template <typename genType>
|
|
inline genType tanh
|
|
(
|
|
genType const & angle
|
|
)
|
|
{
|
|
GLM_STATIC_ASSERT(detail::type<genType>::is_float);
|
|
|
|
return std::tanh(angle);
|
|
}
|
|
|
|
template <typename T>
|
|
inline detail::tvec2<T> tanh
|
|
(
|
|
detail::tvec2<T> const & angle
|
|
)
|
|
{
|
|
return detail::tvec2<T>(
|
|
tanh(angle.x),
|
|
tanh(angle.y));
|
|
}
|
|
|
|
template <typename T>
|
|
inline detail::tvec3<T> tanh
|
|
(
|
|
detail::tvec3<T> const & angle
|
|
)
|
|
{
|
|
return detail::tvec3<T>(
|
|
tanh(angle.x),
|
|
tanh(angle.y),
|
|
tanh(angle.z));
|
|
}
|
|
|
|
template <typename T>
|
|
inline detail::tvec4<T> tanh
|
|
(
|
|
detail::tvec4<T> const & angle
|
|
)
|
|
{
|
|
return detail::tvec4<T>(
|
|
tanh(angle.x),
|
|
tanh(angle.y),
|
|
tanh(angle.z),
|
|
tanh(angle.w));
|
|
}
|
|
|
|
// asinh
|
|
template <typename genType>
|
|
inline genType asinh
|
|
(
|
|
genType const & x
|
|
)
|
|
{
|
|
GLM_STATIC_ASSERT(detail::type<genType>::is_float);
|
|
|
|
return (x < genType(0) ? genType(-1) : (x > genType(0) ? genType(1) : genType(0))) * log(abs(x) + sqrt(genType(1) + x * x));
|
|
}
|
|
|
|
template <typename T>
|
|
inline detail::tvec2<T> asinh
|
|
(
|
|
detail::tvec2<T> const & x
|
|
)
|
|
{
|
|
return detail::tvec2<T>(
|
|
asinh(x.x),
|
|
asinh(x.y));
|
|
}
|
|
|
|
template <typename T>
|
|
inline detail::tvec3<T> asinh
|
|
(
|
|
detail::tvec3<T> const & x
|
|
)
|
|
{
|
|
return detail::tvec3<T>(
|
|
asinh(x.x),
|
|
asinh(x.y),
|
|
asinh(x.z));
|
|
}
|
|
|
|
template <typename T>
|
|
inline detail::tvec4<T> asinh
|
|
(
|
|
detail::tvec4<T> const & x
|
|
)
|
|
{
|
|
return detail::tvec4<T>(
|
|
asinh(x.x),
|
|
asinh(x.y),
|
|
asinh(x.z),
|
|
asinh(x.w));
|
|
}
|
|
|
|
// acosh
|
|
template <typename genType>
|
|
inline genType acosh
|
|
(
|
|
genType const & x
|
|
)
|
|
{
|
|
GLM_STATIC_ASSERT(detail::type<genType>::is_float);
|
|
|
|
if(x < genType(1))
|
|
return genType(0);
|
|
return log(x + sqrt(x * x - genType(1)));
|
|
}
|
|
|
|
template <typename T>
|
|
inline detail::tvec2<T> acosh
|
|
(
|
|
detail::tvec2<T> const & x
|
|
)
|
|
{
|
|
return detail::tvec2<T>(
|
|
acosh(x.x),
|
|
acosh(x.y));
|
|
}
|
|
|
|
template <typename T>
|
|
inline detail::tvec3<T> acosh
|
|
(
|
|
detail::tvec3<T> const & x
|
|
)
|
|
{
|
|
return detail::tvec3<T>(
|
|
acosh(x.x),
|
|
acosh(x.y),
|
|
acosh(x.z));
|
|
}
|
|
|
|
template <typename T>
|
|
inline detail::tvec4<T> acosh
|
|
(
|
|
detail::tvec4<T> const & x
|
|
)
|
|
{
|
|
return detail::tvec4<T>(
|
|
acosh(x.x),
|
|
acosh(x.y),
|
|
acosh(x.z),
|
|
acosh(x.w));
|
|
}
|
|
|
|
// atanh
|
|
template <typename genType>
|
|
inline genType atanh
|
|
(
|
|
genType const & x
|
|
)
|
|
{
|
|
GLM_STATIC_ASSERT(detail::type<genType>::is_float);
|
|
|
|
if(abs(x) >= genType(1))
|
|
return 0;
|
|
return genType(0.5) * log((genType(1) + x) / (genType(1) - x));
|
|
}
|
|
|
|
template <typename T>
|
|
inline detail::tvec2<T> atanh
|
|
(
|
|
detail::tvec2<T> const & x
|
|
)
|
|
{
|
|
return detail::tvec2<T>(
|
|
atanh(x.x),
|
|
atanh(x.y));
|
|
}
|
|
|
|
template <typename T>
|
|
inline detail::tvec3<T> atanh
|
|
(
|
|
detail::tvec3<T> const & x
|
|
)
|
|
{
|
|
return detail::tvec3<T>(
|
|
atanh(x.x),
|
|
atanh(x.y),
|
|
atanh(x.z));
|
|
}
|
|
|
|
template <typename T>
|
|
inline detail::tvec4<T> atanh
|
|
(
|
|
detail::tvec4<T> const & x
|
|
)
|
|
{
|
|
return detail::tvec4<T>(
|
|
atanh(x.x),
|
|
atanh(x.y),
|
|
atanh(x.z),
|
|
atanh(x.w));
|
|
}
|
|
|
|
}//namespace trigonometric
|
|
}//namespace function
|
|
}//namespace core
|
|
}//namespace glm
|