#ifndef INCLUDED_IMATHMATRIXALGO_H
#define INCLUDED_IMATHMATRIXALGO_H
#include "ImathMatrix.h"
#include "ImathQuat.h"
#include "ImathEuler.h"
#include "ImathExc.h"
#include "ImathVec.h"
#include "ImathLimits.h"
#include <math.h>
#ifdef OPENEXR_DLL
#ifdef IMATH_EXPORTS
#define IMATH_EXPORT_CONST extern __declspec(dllexport)
#else
#define IMATH_EXPORT_CONST extern __declspec(dllimport)
#endif
#else
#define IMATH_EXPORT_CONST extern const
#endif
namespace Imath {
IMATH_EXPORT_CONST M33f identity33f;
IMATH_EXPORT_CONST M44f identity44f;
IMATH_EXPORT_CONST M33d identity33d;
IMATH_EXPORT_CONST M44d identity44d;
template <class T> bool extractScaling
(const Matrix44<T> &mat,
Vec3<T> &scl,
bool exc = true);
template <class T> Matrix44<T> sansScaling (const Matrix44<T> &mat,
bool exc = true);
template <class T> bool removeScaling
(Matrix44<T> &mat,
bool exc = true);
template <class T> bool extractScalingAndShear
(const Matrix44<T> &mat,
Vec3<T> &scl,
Vec3<T> &shr,
bool exc = true);
template <class T> Matrix44<T> sansScalingAndShear
(const Matrix44<T> &mat,
bool exc = true);
template <class T> void sansScalingAndShear
(Matrix44<T> &result,
const Matrix44<T> &mat,
bool exc = true);
template <class T> bool removeScalingAndShear
(Matrix44<T> &mat,
bool exc = true);
template <class T> bool extractAndRemoveScalingAndShear
(Matrix44<T> &mat,
Vec3<T> &scl,
Vec3<T> &shr,
bool exc = true);
template <class T> void extractEulerXYZ
(const Matrix44<T> &mat,
Vec3<T> &rot);
template <class T> void extractEulerZYX
(const Matrix44<T> &mat,
Vec3<T> &rot);
template <class T> Quat<T> extractQuat (const Matrix44<T> &mat);
template <class T> bool extractSHRT
(const Matrix44<T> &mat,
Vec3<T> &s,
Vec3<T> &h,
Vec3<T> &r,
Vec3<T> &t,
bool exc ,
typename Euler<T>::Order rOrder);
template <class T> bool extractSHRT
(const Matrix44<T> &mat,
Vec3<T> &s,
Vec3<T> &h,
Vec3<T> &r,
Vec3<T> &t,
bool exc = true);
template <class T> bool extractSHRT
(const Matrix44<T> &mat,
Vec3<T> &s,
Vec3<T> &h,
Euler<T> &r,
Vec3<T> &t,
bool exc = true);
template <class T> bool checkForZeroScaleInRow
(const T &scl,
const Vec3<T> &row,
bool exc = true);
template <class T> Matrix44<T> outerProduct
( const Vec4<T> &a,
const Vec4<T> &b);
template <class T> Matrix44<T> rotationMatrix (const Vec3<T> &fromDirection,
const Vec3<T> &toDirection);
template <class T> Matrix44<T> rotationMatrixWithUpDir
(const Vec3<T> &fromDir,
const Vec3<T> &toDir,
const Vec3<T> &upDir);
template <class T> void alignZAxisWithTargetDir
(Matrix44<T> &result,
Vec3<T> targetDir,
Vec3<T> upDir);
template <class T> Matrix44<T> computeLocalFrame( const Vec3<T>& p,
const Vec3<T>& xDir,
const Vec3<T>& normal);
template <class T> Matrix44<T> addOffset( const Matrix44<T>& inMat,
const Vec3<T>& tOffset,
const Vec3<T>& rOffset,
const Vec3<T>& sOffset,
const Vec3<T>& ref);
template <class T> Matrix44<T> computeRSMatrix( bool keepRotateA,
bool keepScaleA,
const Matrix44<T>& A,
const Matrix44<T>& B);
template <class T> bool extractScaling
(const Matrix33<T> &mat,
Vec2<T> &scl,
bool exc = true);
template <class T> Matrix33<T> sansScaling (const Matrix33<T> &mat,
bool exc = true);
template <class T> bool removeScaling
(Matrix33<T> &mat,
bool exc = true);
template <class T> bool extractScalingAndShear
(const Matrix33<T> &mat,
Vec2<T> &scl,
T &h,
bool exc = true);
template <class T> Matrix33<T> sansScalingAndShear
(const Matrix33<T> &mat,
bool exc = true);
template <class T> bool removeScalingAndShear
(Matrix33<T> &mat,
bool exc = true);
template <class T> bool extractAndRemoveScalingAndShear
(Matrix33<T> &mat,
Vec2<T> &scl,
T &shr,
bool exc = true);
template <class T> void extractEuler
(const Matrix33<T> &mat,
T &rot);
template <class T> bool extractSHRT (const Matrix33<T> &mat,
Vec2<T> &s,
T &h,
T &r,
Vec2<T> &t,
bool exc = true);
template <class T> bool checkForZeroScaleInRow
(const T &scl,
const Vec2<T> &row,
bool exc = true);
template <class T> Matrix33<T> outerProduct
( const Vec3<T> &a,
const Vec3<T> &b);
template <class T>
bool
extractScaling (const Matrix44<T> &mat, Vec3<T> &scl, bool exc)
{
Vec3<T> shr;
Matrix44<T> M (mat);
if (! extractAndRemoveScalingAndShear (M, scl, shr, exc))
return false;
return true;
}
template <class T>
Matrix44<T>
sansScaling (const Matrix44<T> &mat, bool exc)
{
Vec3<T> scl;
Vec3<T> shr;
Vec3<T> rot;
Vec3<T> tran;
if (! extractSHRT (mat, scl, shr, rot, tran, exc))
return mat;
Matrix44<T> M;
M.translate (tran);
M.rotate (rot);
M.shear (shr);
return M;
}
template <class T>
bool
removeScaling (Matrix44<T> &mat, bool exc)
{
Vec3<T> scl;
Vec3<T> shr;
Vec3<T> rot;
Vec3<T> tran;
if (! extractSHRT (mat, scl, shr, rot, tran, exc))
return false;
mat.makeIdentity ();
mat.translate (tran);
mat.rotate (rot);
mat.shear (shr);
return true;
}
template <class T>
bool
extractScalingAndShear (const Matrix44<T> &mat,
Vec3<T> &scl, Vec3<T> &shr, bool exc)
{
Matrix44<T> M (mat);
if (! extractAndRemoveScalingAndShear (M, scl, shr, exc))
return false;
return true;
}
template <class T>
Matrix44<T>
sansScalingAndShear (const Matrix44<T> &mat, bool exc)
{
Vec3<T> scl;
Vec3<T> shr;
Matrix44<T> M (mat);
if (! extractAndRemoveScalingAndShear (M, scl, shr, exc))
return mat;
return M;
}
template <class T>
void
sansScalingAndShear (Matrix44<T> &result, const Matrix44<T> &mat, bool exc)
{
Vec3<T> scl;
Vec3<T> shr;
if (! extractAndRemoveScalingAndShear (result, scl, shr, exc))
result = mat;
}
template <class T>
bool
removeScalingAndShear (Matrix44<T> &mat, bool exc)
{
Vec3<T> scl;
Vec3<T> shr;
if (! extractAndRemoveScalingAndShear (mat, scl, shr, exc))
return false;
return true;
}
template <class T>
bool
extractAndRemoveScalingAndShear (Matrix44<T> &mat,
Vec3<T> &scl, Vec3<T> &shr, bool exc)
{
Vec3<T> row[3];
row[0] = Vec3<T> (mat[0][0], mat[0][1], mat[0][2]);
row[1] = Vec3<T> (mat[1][0], mat[1][1], mat[1][2]);
row[2] = Vec3<T> (mat[2][0], mat[2][1], mat[2][2]);
T maxVal = 0;
for (int i=0; i < 3; i++)
for (int j=0; j < 3; j++)
if (Imath::abs (row[i][j]) > maxVal)
maxVal = Imath::abs (row[i][j]);
if (maxVal != 0)
{
for (int i=0; i < 3; i++)
if (! checkForZeroScaleInRow (maxVal, row[i], exc))
return false;
else
row[i] /= maxVal;
}
scl.x = row[0].length ();
if (! checkForZeroScaleInRow (scl.x, row[0], exc))
return false;
row[0] /= scl.x;
shr[0] = row[0].dot (row[1]);
row[1] -= shr[0] * row[0];
scl.y = row[1].length ();
if (! checkForZeroScaleInRow (scl.y, row[1], exc))
return false;
row[1] /= scl.y;
shr[0] /= scl.y;
shr[1] = row[0].dot (row[2]);
row[2] -= shr[1] * row[0];
shr[2] = row[1].dot (row[2]);
row[2] -= shr[2] * row[1];
scl.z = row[2].length ();
if (! checkForZeroScaleInRow (scl.z, row[2], exc))
return false;
row[2] /= scl.z;
shr[1] /= scl.z;
shr[2] /= scl.z;
if (row[0].dot (row[1].cross (row[2])) < 0)
for (int i=0; i < 3; i++)
{
scl[i] *= -1;
row[i] *= -1;
}
for (int i=0; i < 3; i++)
{
mat[i][0] = row[i][0];
mat[i][1] = row[i][1];
mat[i][2] = row[i][2];
}
scl *= maxVal;
return true;
}
template <class T>
void
extractEulerXYZ (const Matrix44<T> &mat, Vec3<T> &rot)
{
Vec3<T> i (mat[0][0], mat[0][1], mat[0][2]);
Vec3<T> j (mat[1][0], mat[1][1], mat[1][2]);
Vec3<T> k (mat[2][0], mat[2][1], mat[2][2]);
i.normalize();
j.normalize();
k.normalize();
Matrix44<T> M (i[0], i[1], i[2], 0,
j[0], j[1], j[2], 0,
k[0], k[1], k[2], 0,
0, 0, 0, 1);
rot.x = Math<T>::atan2 (M[1][2], M[2][2]);
Matrix44<T> N;
N.rotate (Vec3<T> (-rot.x, 0, 0));
N = N * M;
T cy = Math<T>::sqrt (N[0][0]*N[0][0] + N[0][1]*N[0][1]);
rot.y = Math<T>::atan2 (-N[0][2], cy);
rot.z = Math<T>::atan2 (-N[1][0], N[1][1]);
}
template <class T>
void
extractEulerZYX (const Matrix44<T> &mat, Vec3<T> &rot)
{
Vec3<T> i (mat[0][0], mat[0][1], mat[0][2]);
Vec3<T> j (mat[1][0], mat[1][1], mat[1][2]);
Vec3<T> k (mat[2][0], mat[2][1], mat[2][2]);
i.normalize();
j.normalize();
k.normalize();
Matrix44<T> M (i[0], i[1], i[2], 0,
j[0], j[1], j[2], 0,
k[0], k[1], k[2], 0,
0, 0, 0, 1);
rot.x = -Math<T>::atan2 (M[1][0], M[0][0]);
Matrix44<T> N;
N.rotate (Vec3<T> (0, 0, -rot.x));
N = N * M;
T cy = Math<T>::sqrt (N[2][2]*N[2][2] + N[2][1]*N[2][1]);
rot.y = -Math<T>::atan2 (-N[2][0], cy);
rot.z = -Math<T>::atan2 (-N[1][2], N[1][1]);
}
template <class T>
Quat<T>
extractQuat (const Matrix44<T> &mat)
{
Matrix44<T> rot;
T tr, s;
T q[4];
int i, j, k;
Quat<T> quat;
int nxt[3] = {1, 2, 0};
tr = mat[0][0] + mat[1][1] + mat[2][2];
if (tr > 0.0) {
s = Math<T>::sqrt (tr + T(1.0));
quat.r = s / T(2.0);
s = T(0.5) / s;
quat.v.x = (mat[1][2] - mat[2][1]) * s;
quat.v.y = (mat[2][0] - mat[0][2]) * s;
quat.v.z = (mat[0][1] - mat[1][0]) * s;
}
else {
i = 0;
if (mat[1][1] > mat[0][0])
i=1;
if (mat[2][2] > mat[i][i])
i=2;
j = nxt[i];
k = nxt[j];
s = Math<T>::sqrt ((mat[i][i] - (mat[j][j] + mat[k][k])) + T(1.0));
q[i] = s * T(0.5);
if (s != T(0.0))
s = T(0.5) / s;
q[3] = (mat[j][k] - mat[k][j]) * s;
q[j] = (mat[i][j] + mat[j][i]) * s;
q[k] = (mat[i][k] + mat[k][i]) * s;
quat.v.x = q[0];
quat.v.y = q[1];
quat.v.z = q[2];
quat.r = q[3];
}
return quat;
}
template <class T>
bool
extractSHRT (const Matrix44<T> &mat,
Vec3<T> &s,
Vec3<T> &h,
Vec3<T> &r,
Vec3<T> &t,
bool exc ,
typename Euler<T>::Order rOrder )
{
Matrix44<T> rot;
rot = mat;
if (! extractAndRemoveScalingAndShear (rot, s, h, exc))
return false;
extractEulerXYZ (rot, r);
t.x = mat[3][0];
t.y = mat[3][1];
t.z = mat[3][2];
if (rOrder != Euler<T>::XYZ)
{
Imath::Euler<T> eXYZ (r, Imath::Euler<T>::XYZ);
Imath::Euler<T> e (eXYZ, rOrder);
r = e.toXYZVector ();
}
return true;
}
template <class T>
bool
extractSHRT (const Matrix44<T> &mat,
Vec3<T> &s,
Vec3<T> &h,
Vec3<T> &r,
Vec3<T> &t,
bool exc)
{
return extractSHRT(mat, s, h, r, t, exc, Imath::Euler<T>::XYZ);
}
template <class T>
bool
extractSHRT (const Matrix44<T> &mat,
Vec3<T> &s,
Vec3<T> &h,
Euler<T> &r,
Vec3<T> &t,
bool exc )
{
return extractSHRT (mat, s, h, r, t, exc, r.order ());
}
template <class T>
bool
checkForZeroScaleInRow (const T& scl,
const Vec3<T> &row,
bool exc )
{
for (int i = 0; i < 3; i++)
{
if ((abs (scl) < 1 && abs (row[i]) >= limits<T>::max() * abs (scl)))
{
if (exc)
throw Imath::ZeroScaleExc ("Cannot remove zero scaling "
"from matrix.");
else
return false;
}
}
return true;
}
template <class T>
Matrix44<T>
outerProduct (const Vec4<T> &a, const Vec4<T> &b )
{
return Matrix44<T> (a.x*b.x, a.x*b.y, a.x*b.z, a.x*b.w,
a.y*b.x, a.y*b.y, a.y*b.z, a.x*b.w,
a.z*b.x, a.z*b.y, a.z*b.z, a.x*b.w,
a.w*b.x, a.w*b.y, a.w*b.z, a.w*b.w);
}
template <class T>
Matrix44<T>
rotationMatrix (const Vec3<T> &from, const Vec3<T> &to)
{
Quat<T> q;
q.setRotation(from, to);
return q.toMatrix44();
}
template <class T>
Matrix44<T>
rotationMatrixWithUpDir (const Vec3<T> &fromDir,
const Vec3<T> &toDir,
const Vec3<T> &upDir)
{
if (fromDir.length () == 0)
return Matrix44<T> ();
else
{
Matrix44<T> zAxis2FromDir( Imath::UNINITIALIZED );
alignZAxisWithTargetDir (zAxis2FromDir, fromDir, Vec3<T> (0, 1, 0));
Matrix44<T> fromDir2zAxis = zAxis2FromDir.transposed ();
Matrix44<T> zAxis2ToDir( Imath::UNINITIALIZED );
alignZAxisWithTargetDir (zAxis2ToDir, toDir, upDir);
return fromDir2zAxis * zAxis2ToDir;
}
}
template <class T>
void
alignZAxisWithTargetDir (Matrix44<T> &result, Vec3<T> targetDir, Vec3<T> upDir)
{
if ( targetDir.length () == 0 )
targetDir = Vec3<T> (0, 0, 1);
if ( upDir.length () == 0 )
upDir = Vec3<T> (0, 1, 0);
if (upDir.cross (targetDir).length () == 0)
{
upDir = targetDir.cross (Vec3<T> (1, 0, 0));
if (upDir.length() == 0)
upDir = targetDir.cross(Vec3<T> (0, 0, 1));
}
Vec3<T> targetPerpDir = upDir.cross (targetDir);
Vec3<T> targetUpDir = targetDir.cross (targetPerpDir);
Vec3<T> row[3];
row[0] = targetPerpDir.normalized ();
row[1] = targetUpDir .normalized ();
row[2] = targetDir .normalized ();
result.x[0][0] = row[0][0];
result.x[0][1] = row[0][1];
result.x[0][2] = row[0][2];
result.x[0][3] = (T)0;
result.x[1][0] = row[1][0];
result.x[1][1] = row[1][1];
result.x[1][2] = row[1][2];
result.x[1][3] = (T)0;
result.x[2][0] = row[2][0];
result.x[2][1] = row[2][1];
result.x[2][2] = row[2][2];
result.x[2][3] = (T)0;
result.x[3][0] = (T)0;
result.x[3][1] = (T)0;
result.x[3][2] = (T)0;
result.x[3][3] = (T)1;
}
template <class T>
Matrix44<T>
computeLocalFrame( const Vec3<T>& p,
const Vec3<T>& xDir,
const Vec3<T>& normal)
{
Vec3<T> _xDir(xDir);
Vec3<T> x = _xDir.normalize();
Vec3<T> y = (normal % x).normalize();
Vec3<T> z = (x % y).normalize();
Matrix44<T> L;
L[0][0] = x[0];
L[0][1] = x[1];
L[0][2] = x[2];
L[0][3] = 0.0;
L[1][0] = y[0];
L[1][1] = y[1];
L[1][2] = y[2];
L[1][3] = 0.0;
L[2][0] = z[0];
L[2][1] = z[1];
L[2][2] = z[2];
L[2][3] = 0.0;
L[3][0] = p[0];
L[3][1] = p[1];
L[3][2] = p[2];
L[3][3] = 1.0;
return L;
}
template <class T>
Matrix44<T>
addOffset( const Matrix44<T>& inMat,
const Vec3<T>& tOffset,
const Vec3<T>& rOffset,
const Vec3<T>& sOffset,
const Matrix44<T>& ref)
{
Matrix44<T> O;
Vec3<T> _rOffset(rOffset);
_rOffset *= M_PI / 180.0;
O.rotate (_rOffset);
O[3][0] = tOffset[0];
O[3][1] = tOffset[1];
O[3][2] = tOffset[2];
Matrix44<T> S;
S.scale (sOffset);
Matrix44<T> X = S * O * inMat * ref;
return X;
}
template <class T>
Matrix44<T>
computeRSMatrix( bool keepRotateA,
bool keepScaleA,
const Matrix44<T>& A,
const Matrix44<T>& B)
{
Vec3<T> as, ah, ar, at;
extractSHRT (A, as, ah, ar, at);
Vec3<T> bs, bh, br, bt;
extractSHRT (B, bs, bh, br, bt);
if (!keepRotateA)
ar = br;
if (!keepScaleA)
as = bs;
Matrix44<T> mat;
mat.makeIdentity();
mat.translate (at);
mat.rotate (ar);
mat.scale (as);
return mat;
}
template <class T>
bool
extractScaling (const Matrix33<T> &mat, Vec2<T> &scl, bool exc)
{
T shr;
Matrix33<T> M (mat);
if (! extractAndRemoveScalingAndShear (M, scl, shr, exc))
return false;
return true;
}
template <class T>
Matrix33<T>
sansScaling (const Matrix33<T> &mat, bool exc)
{
Vec2<T> scl;
T shr;
T rot;
Vec2<T> tran;
if (! extractSHRT (mat, scl, shr, rot, tran, exc))
return mat;
Matrix33<T> M;
M.translate (tran);
M.rotate (rot);
M.shear (shr);
return M;
}
template <class T>
bool
removeScaling (Matrix33<T> &mat, bool exc)
{
Vec2<T> scl;
T shr;
T rot;
Vec2<T> tran;
if (! extractSHRT (mat, scl, shr, rot, tran, exc))
return false;
mat.makeIdentity ();
mat.translate (tran);
mat.rotate (rot);
mat.shear (shr);
return true;
}
template <class T>
bool
extractScalingAndShear (const Matrix33<T> &mat, Vec2<T> &scl, T &shr, bool exc)
{
Matrix33<T> M (mat);
if (! extractAndRemoveScalingAndShear (M, scl, shr, exc))
return false;
return true;
}
template <class T>
Matrix33<T>
sansScalingAndShear (const Matrix33<T> &mat, bool exc)
{
Vec2<T> scl;
T shr;
Matrix33<T> M (mat);
if (! extractAndRemoveScalingAndShear (M, scl, shr, exc))
return mat;
return M;
}
template <class T>
bool
removeScalingAndShear (Matrix33<T> &mat, bool exc)
{
Vec2<T> scl;
T shr;
if (! extractAndRemoveScalingAndShear (mat, scl, shr, exc))
return false;
return true;
}
template <class T>
bool
extractAndRemoveScalingAndShear (Matrix33<T> &mat,
Vec2<T> &scl, T &shr, bool exc)
{
Vec2<T> row[2];
row[0] = Vec2<T> (mat[0][0], mat[0][1]);
row[1] = Vec2<T> (mat[1][0], mat[1][1]);
T maxVal = 0;
for (int i=0; i < 2; i++)
for (int j=0; j < 2; j++)
if (Imath::abs (row[i][j]) > maxVal)
maxVal = Imath::abs (row[i][j]);
if (maxVal != 0)
{
for (int i=0; i < 2; i++)
if (! checkForZeroScaleInRow (maxVal, row[i], exc))
return false;
else
row[i] /= maxVal;
}
scl.x = row[0].length ();
if (! checkForZeroScaleInRow (scl.x, row[0], exc))
return false;
row[0] /= scl.x;
shr = row[0].dot (row[1]);
row[1] -= shr * row[0];
scl.y = row[1].length ();
if (! checkForZeroScaleInRow (scl.y, row[1], exc))
return false;
row[1] /= scl.y;
shr /= scl.y;
if (row[0][0] * row[1][1] - row[0][1] * row[1][0] < 0)
{
row[1][0] *= -1;
row[1][1] *= -1;
scl[1] *= -1;
shr *= -1;
}
for (int i=0; i < 2; i++)
{
mat[i][0] = row[i][0];
mat[i][1] = row[i][1];
}
scl *= maxVal;
return true;
}
template <class T>
void
extractEuler (const Matrix33<T> &mat, T &rot)
{
Vec2<T> i (mat[0][0], mat[0][1]);
Vec2<T> j (mat[1][0], mat[1][1]);
i.normalize();
j.normalize();
rot = - Math<T>::atan2 (j[0], i[0]);
}
template <class T>
bool
extractSHRT (const Matrix33<T> &mat,
Vec2<T> &s,
T &h,
T &r,
Vec2<T> &t,
bool exc)
{
Matrix33<T> rot;
rot = mat;
if (! extractAndRemoveScalingAndShear (rot, s, h, exc))
return false;
extractEuler (rot, r);
t.x = mat[2][0];
t.y = mat[2][1];
return true;
}
template <class T>
bool
checkForZeroScaleInRow (const T& scl,
const Vec2<T> &row,
bool exc )
{
for (int i = 0; i < 2; i++)
{
if ((abs (scl) < 1 && abs (row[i]) >= limits<T>::max() * abs (scl)))
{
if (exc)
throw Imath::ZeroScaleExc ("Cannot remove zero scaling "
"from matrix.");
else
return false;
}
}
return true;
}
template <class T>
Matrix33<T>
outerProduct (const Vec3<T> &a, const Vec3<T> &b )
{
return Matrix33<T> (a.x*b.x, a.x*b.y, a.x*b.z,
a.y*b.x, a.y*b.y, a.y*b.z,
a.z*b.x, a.z*b.y, a.z*b.z );
}
template <typename T>
Imath::M44d
procrustesRotationAndTranslation (const Imath::Vec3<T>* A,
const Imath::Vec3<T>* B,
const T* weights,
const size_t numPoints,
const bool doScaling = false);
template <typename T>
Imath::M44d
procrustesRotationAndTranslation (const Imath::Vec3<T>* A,
const Imath::Vec3<T>* B,
const size_t numPoints,
const bool doScaling = false);
template <typename T>
void
jacobiSVD (const Imath::Matrix33<T>& A,
Imath::Matrix33<T>& U,
Imath::Vec3<T>& S,
Imath::Matrix33<T>& V,
const T tol = Imath::limits<T>::epsilon(),
const bool forcePositiveDeterminant = false);
template <typename T>
void
jacobiSVD (const Imath::Matrix44<T>& A,
Imath::Matrix44<T>& U,
Imath::Vec4<T>& S,
Imath::Matrix44<T>& V,
const T tol = Imath::limits<T>::epsilon(),
const bool forcePositiveDeterminant = false);
template <typename T>
void
jacobiEigenSolver (Matrix33<T>& A,
Vec3<T>& S,
Matrix33<T>& V,
const T tol);
template <typename T>
inline
void
jacobiEigenSolver (Matrix33<T>& A,
Vec3<T>& S,
Matrix33<T>& V)
{
jacobiEigenSolver(A,S,V,limits<T>::epsilon());
}
template <typename T>
void
jacobiEigenSolver (Matrix44<T>& A,
Vec4<T>& S,
Matrix44<T>& V,
const T tol);
template <typename T>
inline
void
jacobiEigenSolver (Matrix44<T>& A,
Vec4<T>& S,
Matrix44<T>& V)
{
jacobiEigenSolver(A,S,V,limits<T>::epsilon());
}
template <typename TM, typename TV>
void
maxEigenVector (TM& A, TV& S);
template <typename TM, typename TV>
void
minEigenVector (TM& A, TV& S);
}
#endif