本文整理汇总了C++中Matrix3x3::QLAlgorithm方法的典型用法代码示例。如果您正苦于以下问题:C++ Matrix3x3::QLAlgorithm方法的具体用法?C++ Matrix3x3::QLAlgorithm怎么用?C++ Matrix3x3::QLAlgorithm使用的例子?那么, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类Matrix3x3
的用法示例。
在下文中一共展示了Matrix3x3::QLAlgorithm方法的1个代码示例,这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞,您的评价将有助于系统推荐出更棒的C++代码示例。
示例1:
//-----------------------------------------------------------------------
void Matrix3x3::EigenSolveSymmetric (double afEigenvalue[3],
Vector3 akEigenvector[3]) const
{
Matrix3x3 kMatrix = *this;
double afSubDiag[3];
kMatrix.Tridiagonal(afEigenvalue,afSubDiag);
kMatrix.QLAlgorithm(afEigenvalue,afSubDiag);
for (size_t i = 0; i < 3; i++)
{
akEigenvector[i][0] = kMatrix[0][i];
akEigenvector[i][1] = kMatrix[1][i];
akEigenvector[i][2] = kMatrix[2][i];
}
// make eigenvectors form a right--handed system
Vector3 kCross = akEigenvector[1].CrossProduct(akEigenvector[2]);
double fDet = akEigenvector[0].DotProduct(kCross);
if ( fDet < 0.0 )
{
akEigenvector[2][0] = - akEigenvector[2][0];
akEigenvector[2][1] = - akEigenvector[2][1];
akEigenvector[2][2] = - akEigenvector[2][2];
}
}