public void Simple3x3()
{
NDArray <double> np1 = new NDArray <double>().Zeros(3, 3);
np1[0, 0] = 5;
np1[0, 1] = 1;
np1[0, 2] = 2;
np1[1, 0] = 1;
np1[1, 1] = 0;
np1[1, 2] = 1;
np1[2, 0] = 1;
np1[2, 1] = 1;
np1[2, 2] = 0;
NDArray <double> np1Inv = np1.Inv();
var OncesMatrix = np1.Dot(np1Inv);
Assert.IsTrue(Math.Abs(OncesMatrix[0, 0]) < 1.000001);
Assert.IsTrue(Math.Abs(OncesMatrix[1, 1]) < 1.000001);
Assert.IsTrue(Math.Abs(OncesMatrix[2, 2]) < 1.000001);
Assert.IsTrue(Math.Abs(OncesMatrix[0, 1]) < 0.000001);
Assert.IsTrue(Math.Abs(OncesMatrix[0, 2]) < 0.000001);
Assert.IsTrue(Math.Abs(OncesMatrix[1, 0]) < 0.000001);
Assert.IsTrue(Math.Abs(OncesMatrix[1, 2]) < 0.000001);
Assert.IsTrue(Math.Abs(OncesMatrix[2, 0]) < 0.000001);
Assert.IsTrue(Math.Abs(OncesMatrix[2, 1]) < 0.000001);
}
public void MatrixMultiplyComplex()
{
NDArray <Complex> matrix1 = new NDArray <Complex>();
matrix1.Data = new Complex[] { new Complex(1, -1), new Complex(2, -2), new Complex(3, 0), new Complex(4, 0), 5, 6 };
matrix1.Shape = new int[] { 3, 2 };
NDArray <Complex> matrix2 = new NDArray <Complex>();
matrix2.Data = new Complex[] { 7, 8, 9, new Complex(10, -10), 11, new Complex(12, -12) };
matrix2.Shape = new int[] { 2, 3 };
var matrix3 = matrix1.Dot(matrix2);
var matrix4 = new Complex[9];
matrix4[0] = new Complex(7, -47);
matrix4[1] = new Complex(30, -30);
matrix4[2] = new Complex(9, -57);
matrix4[3] = new Complex(61, -40);
matrix4[4] = new Complex(68, 0);
matrix4[5] = new Complex(75, -48);
matrix4[6] = new Complex(95, -60);
matrix4[7] = new Complex(106, 0);
matrix4[8] = new Complex(117, -72);
Assert.IsTrue(Enumerable.SequenceEqual(matrix4, matrix3.Data));
}
public void DotTwo1DDouble()
{
var series1 = new NDArray <double>().ARange(4, 1, 1);
var series2 = new NDArray <double>();
series2.Data = new double[] { 0, 1, 0.5 };
var innerProduct = series1.Dot(series2);
Assert.IsTrue(innerProduct.Data[0] == 3.5);
}
public void Simple2x2()
{
NDArray <double> np1 = new NDArray <double>().ARange(4).ReShape(2, 2);
NDArray <double> np1Inv = np1.Inv();
var OncesMatrix = np1.Dot(np1Inv);
Assert.IsTrue(OncesMatrix[0, 0] == 1);
Assert.IsTrue(OncesMatrix[1, 1] == 1);
Assert.IsTrue(OncesMatrix[1, 0] == 0);
Assert.IsTrue(OncesMatrix[0, 1] == 0);
}
public void MatrixMultiplyDouble()
{
NDArray <double> matrix1 = new NDArray <double>().ARange(7, 1, 1).ReShape(3, 2);
Assert.IsTrue(matrix1[0, 0] == 1);
Assert.IsTrue(matrix1[0, 1] == 2);
Assert.IsTrue(matrix1[1, 0] == 3);
Assert.IsTrue(matrix1[2, 0] == 5);
NDArray <double> matrix2 = new NDArray <double>().ARange(13, 7, 1).ReShape(2, 3);
Assert.IsTrue(matrix2[0, 0] == 7);
Assert.IsTrue(matrix2[0, 1] == 8);
Assert.IsTrue(matrix2[1, 0] == 10);
Assert.IsTrue(matrix2[0, 2] == 9);
var matrix3 = matrix1.Dot(matrix2);
var matrix4 = new NDArray <double>().Zeros(3, 3);
for (int idx = 0; idx < 3; idx++)
{
for (int jdx = 0; jdx < 3; jdx++)
{
matrix4[idx, jdx] = 0;
for (int kdx = 0; kdx < 2; kdx++)
{
matrix4[idx, jdx] += matrix1[idx, kdx] * matrix2[kdx, jdx];
}
}
}
Assert.IsTrue(matrix4[0, 0] == matrix3[0, 0]);
Assert.IsTrue(matrix4[0, 1] == matrix3[0, 1]);
Assert.IsTrue(matrix4[0, 2] == matrix3[0, 2]);
Assert.IsTrue(matrix4[1, 0] == matrix3[1, 0]);
Assert.IsTrue(matrix4[1, 1] == matrix3[1, 1]);
Assert.IsTrue(matrix4[1, 2] == matrix3[1, 2]);
Assert.IsTrue(matrix4[2, 0] == matrix3[2, 0]);
Assert.IsTrue(matrix4[2, 1] == matrix3[2, 1]);
Assert.IsTrue(matrix4[2, 2] == matrix3[2, 2]);
}