本文整理汇总了C++中mat3::At方法的典型用法代码示例。如果您正苦于以下问题：C++ mat3::At方法的具体用法？C++ mat3::At怎么用？C++ mat3::At使用的例子？那么, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类mat3的用法示例。
在下文中一共展示了mat3::At方法的1个代码示例，这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞，您的评价将有助于系统推荐出更棒的C++代码示例。
示例1: SymmetricEigenSystem
bool SymmetricEigenSystem(const mat3& m, vec3& value, mat3& r)
{
const int max_sweep = 32;
const float epsilon = 1.0e-10f;
float m11 = m.At(0,0);
float m12 = m.At(0,1);
float m13 = m.At(0,2);
float m22 = m.At(1,1);
float m23 = m.At(1,2);
float m33 = m.At(2,2);
r.Identity();
for (int a = 0; a < max_sweep; ++a)
{
if ((Abs(m12) < epsilon) && (Abs(m13) < epsilon) && (Abs(m23) < epsilon))
break;
// annihilate 1 2
if (m12 != 0.0f)
{
float u = (m22 - m11) * 0.5f / m12;
float u2 = u*u;
float u2p1 = u2 + 1.0f;
float t = (u2p1 != u2) ? ((u < 0.0f) ? -1.0f : 1.0f) * (sqrt(u2p1) - fabs(u)) : 0.5f / u;
float c = 1.0f / sqrt(t*t + 1.0f);
float s = c * t;
m11 -= t * m12;
m22 += t * m12;
m12 = 0.0f;
float temp = c*m13 - s * m23;
m23 = s * m13 + c * m23;
m13 = temp;
for (int i = 0; i < 3; ++i)
{
float temp = c * r.At(i, 0) - s * r.At(i,1);
r.At(i,1) = s * r.At(i,0) + c * r.At(i,1);
r.At(i,0) = temp;
}
}
// annihilate 1 3
if (m13 != 0.0f)
{
float u = (m33 - m11) * 0.5f / m13;
float u2 = u*u;
float u2p1 = u2 + 1.0f;
float t = (u2p1 != u2) ? ((u < 0.0f) ? -1.0f : 1.0f) * (sqrt(u2p1) - fabs(u)) : 0.5f / u;
float c = 1.0f / sqrt(t*t + 1.0f);
float s = c * t;
m11 -= t * m13;
m33 += t * m13;
m13 = 0.0f;
float temp = c * m12 - s * m23;
m23 = s * m12 + c * m23;
m12 = temp;
for (int i = 0; i < 3; ++i)
{
float temp = c * r.At(i, 0) - s * r.At(i, 2);
r.At(i, 2) = s * r.At(i, 0) + c * r.At(i, 2);
r.At(i, 0) = temp;
}
}
// annihilate 2 3
if (m23 != 0.0f)
{
float u = (m33 - m22) * 0.5f / m23;
float u2 = u*u;
float u2p1 = u2 + 1.0f;
float t = (u2p1 != u2) ? ((u < 0.0f) ? -1.0f : 1.0f) * (sqrt(u2p1) - fabs(u)) : 0.5f / u;
float c = 1.0f / sqrt(t*t + 1.0f);
float s = c * t;
m22 -= t * m23;
m33 += t * m23;
m23 = 0.0f;
float temp = c * m12 - s * m13;
m13 = s * m12 + c * m13;
m12 = temp;
for (int i = 0; i < 3; ++i)
{
float temp = c * r.At(i, 1) - s * r.At(i, 2);
r.At(i, 2) = s * r.At(i, 1) + c * r.At(i, 2);
r.At(i, 1) = temp;
}
}
}
value.Set(m11, m22, m33);
return true;
}