public void MatrixMultiplyTest2()
{
var a = MatrixOperations.FromArray(4, 4, new double[]{
1, 2, 3, 4,
2, 4, 4, 2,
8, 6, 4, 1,
0, 0, 0, 1
});
var b = MatrixOperations.FromArray(4, 1, new double[]{
1, 2, 3, 1
});
var c = MatrixOperations.Multiply(a, b);
var expected = MatrixOperations.FromArray(4, 1, new double[]{
18,
24,
33,
1
});
Assert.Equal(expected, c);
output.WriteLine(c.ToString());
}
public void MatrixMultiplyTest1()
{
var a = MatrixOperations.FromArray(4, 4, new double[]{
1, 2, 3, 4,
5, 6, 7, 8,
9, 8, 7, 6,
5, 4, 3, 2
});
var b = MatrixOperations.FromArray(4, 4, new double[]{
-2, 1, 2, 3,
3, 2, 1, -1,
4, 3, 6, 5,
1, 2, 7, 8
});
var c = MatrixOperations.Multiply(a, b);
var expected = MatrixOperations.FromArray(4, 4, new double[]{
20, 22, 50, 48,
44, 54, 114, 108,
40, 58, 110, 102,
16, 26, 46, 42
});
Assert.Equal(expected, c);
output.WriteLine(c.ToString());
}
public void MatrixDeterminantTest1()
{
var a = MatrixOperations.FromArray(2, 2, new double[]{
1, 5,
-3, 2,
});
Assert.True(Constants.EpsilonCompare(17, MatrixOperations.Determinant(a, MatrixOperation.Cofactor)));
output.WriteLine(MatrixOperations.Determinant(a, MatrixOperation.Cofactor).ToString());
}
public void MatrixDeterminantTest2()
{
var a = MatrixOperations.FromArray(3, 3, new double[]{
1, 2, 6,
-5, 8, -4,
2, 6, 4
});
Assert.True(Constants.EpsilonCompare(-196, MatrixOperations.Determinant(a, MatrixOperation.Cofactor)));
Assert.True(Constants.EpsilonCompare(-196, MatrixOperations.Determinant(a, MatrixOperation.Gauss)));
output.WriteLine(MatrixOperations.Determinant(a, MatrixOperation.Cofactor).ToString());
}
public void MatrixDeterminantTest3()
{
var a = MatrixOperations.FromArray(4, 4, new double[]{
-2, -8, 3, 5,
-3, 1, 7, 3,
1, 2, -9, 6,
-6, 7, 7, -9
});
Assert.True(Constants.EpsilonCompare(-4071, MatrixOperations.Determinant(a, MatrixOperation.Gauss)));
Assert.True(Constants.EpsilonCompare(-4071, MatrixOperations.Determinant(a, MatrixOperation.Cofactor)));
output.WriteLine(MatrixOperations.Determinant(a, MatrixOperation.Cofactor).ToString());
}
public void MatrixInvertTest1()
{
var a = MatrixOperations.FromArray(4, 4, new double[]{
-2, -8, 3, 5,
-3, 1, 7, 3,
1, 2, -9, 6,
-6, 7, 7, -9
});
var a11 = MatrixOperations.Invert(a, MatrixOperation.Gauss);
var a12 = MatrixOperations.Invert(a, MatrixOperation.Cofactor);
var r1 = MatrixOperations.Multiply(a, a11);
var r2 = MatrixOperations.Multiply(a, a12);
var i = MatrixOperations.Identity(4);
Assert.True(a11.Equals(a12));
Assert.True(i.Equals(r1));
Assert.True(i.Equals(r2));
output.WriteLine(a11.ToString());
}
public void MatrixPowerTest1()
{
var a = MatrixOperations.FromArray(4, 4, new double[]{
-2, -8, 3, 5,
-3, 1, 7, 3,
1, 2, -9, 6,
-6, 7, 7, -9
});
// NOTE: For this input starts to fail when p = 15
for (int p=0; p<15; ++p)
{
var b = MatrixOperations.Identity(4);
for (int i=0; i<p; ++i)
{
b = MatrixOperations.Multiply(b, a);
}
var c = MatrixOperations.Power(a, (uint)p);
Assert.Equal(c, b);
}
}