So I decided to implement some test routines for my math library and stumbled into something that I do not understand.
I'am comparing the results from my library against the counterparts in the Direct3D library and the code below actually signals that my matrix is not equal to the d3d one.
BOOST_AUTO_TEST_CASE(Matrix44_RotationX)
{
{
float angle = 0.0f;
for(unsigned int i = 0; i < 100; ++i)
{
angle = (float)math::Random(-100.0f, 100.0f);
D3DXMATRIX dxA;
D3DXMatrixRotationX(&dxA, angle);
math::Matrix44f A;
A.LoadRotateX(angle);
BOOST_CHECK(A==dxA);
}
}
}
The thing is that they are almost identical so it is a precision issue somewhere. My rotation Matrix looks like this:
template <typename T>
void Matrix44T<T>::LoadRotateX(const T& amountRad)
{
E[ 0] = (T)1.0; E[ 4] = (T)0.0; E[ 8] = (T)0.0; E[12] = (T)0.0;
E[ 1] = (T)0.0; E[ 5] = Cos(amountRad); E[ 9] = -Sin(amountRad); E[13] = (T)0.0;
E[ 2] = (T)0.0; E[ 6] = Sin(amountRad); E[10] = Cos(amountRad); E[14] = (T)0.0;
E[ 3] = (T)0.0; E[ 7] = (T)0.0; E[11] = (T)0.0; E[15] = (T)1.0;
}
[font="Arial"]So it must be something about cos and sin because the rest of the values are just 1 and 0.[/font]
What is going on here?
I can use a tolerance when comparing the matrices but I want to understand why this is happening first.