public void MultiplyMatrix()
{
Matrix33D m = new Matrix33D(12, 23, 45,
67, 89, 90,
43, 65, 87);
Assert.AreEqual(Matrix33D.Zero, Matrix33D.Multiply(m, Matrix33D.Zero));
Assert.AreEqual(Matrix33D.Zero, Matrix33D.Multiply(Matrix33D.Zero, m));
Assert.AreEqual(m, Matrix33D.Multiply(m, Matrix33D.Identity));
Assert.AreEqual(m, Matrix33D.Multiply(Matrix33D.Identity, m));
Assert.IsTrue(Matrix33D.AreNumericallyEqual(Matrix33D.Identity, Matrix33D.Multiply(m, m.Inverse)));
Assert.IsTrue(Matrix33D.AreNumericallyEqual(Matrix33D.Identity, Matrix33D.Multiply(m.Inverse, m)));
Matrix33D m1 = new Matrix33D(columnMajor, MatrixOrder.ColumnMajor);
Matrix33D m2 = new Matrix33D(12, 23, 45,
67, 89, 90,
43, 65, 87);
Matrix33D result = Matrix33D.Multiply(m1, m2);
for (int column = 0; column < 3; column++)
{
for (int row = 0; row < 3; row++)
{
Assert.AreEqual(Vector3D.Dot(m1.GetRow(row), m2.GetColumn(column)), result[row, column]);
}
}
}
public void Multiplication()
{
double s = 0.1234;
Matrix33D m = new Matrix33D(columnMajor, MatrixOrder.ColumnMajor);
m = Matrix33D.Multiply(s, m);
for (int i = 0; i < 9; i++)
Assert.AreEqual(rowMajor[i] * s, m[i]);
}
public void MultiplyVector()
{
Vector3D v = new Vector3D(2.34, 3.45, 4.56);
Assert.AreEqual(v, Matrix33D.Multiply(Matrix33D.Identity, v));
Assert.AreEqual(Vector3D.Zero, Matrix33D.Multiply(Matrix33D.Zero, v));
Matrix33D m = new Matrix33D(12, 23, 45,
67, 89, 90,
43, 65, 87);
Assert.IsTrue(Vector3D.AreNumericallyEqual(v, Matrix33D.Multiply(m * m.Inverse, v)));
for (int i = 0; i < 3; i++)
Assert.AreEqual(Vector3D.Dot(m.GetRow(i), v), Matrix33D.Multiply(m, v)[i]);
}
