public void MultiplyMatrix() { float[] values = new float[] { 12, 23, 45, 56, 67, 89, 90, 12, 43, 65, 87, 43, 34, -12, 84, 44 }; MatrixF m = new MatrixF(4, 4, values, MatrixOrder.RowMajor); Assert.AreEqual(new MatrixF(4, 4), MatrixF.Multiply(m, new MatrixF(4, 4))); Assert.AreEqual(new MatrixF(4, 4), MatrixF.Multiply(new MatrixF(4, 4), m)); Assert.AreEqual(m, MatrixF.Multiply(m, MatrixF.CreateIdentity(4, 4))); Assert.AreEqual(m, MatrixF.Multiply(MatrixF.CreateIdentity(4, 4), m)); Assert.IsTrue(MatrixF.AreNumericallyEqual(MatrixF.CreateIdentity(4, 4), MatrixF.Multiply(m, m.Inverse))); Assert.IsTrue(MatrixF.AreNumericallyEqual(MatrixF.CreateIdentity(4, 4), MatrixF.Multiply(m.Inverse, m))); MatrixF m1 = new MatrixF(3, 4, rowMajor, MatrixOrder.RowMajor); MatrixF m2 = new MatrixF(4, 4, values, MatrixOrder.RowMajor); MatrixF result = MatrixF.Multiply(m1, m2); Assert.AreEqual(3, result.NumberOfRows); Assert.AreEqual(4, result.NumberOfColumns); for (int column = 0; column < 4; column++) for (int row = 0; row < 3; row++) Assert.AreEqual(VectorF.Dot(m1.GetRow(row), m2.GetColumn(column)), result[row, column]); }
public void GetRowException2() { MatrixF m = new MatrixF(3, 4, rowMajor, MatrixOrder.RowMajor); m.GetRow(3); }
public void SetRow() { MatrixF m = new MatrixF(3, 4, rowMajor, MatrixOrder.RowMajor); m.SetRow(0, new VectorF(new float[]{ 0.1f, 0.2f, 0.3f, 0.4f })); Assert.AreEqual(new VectorF(new float[] { 0.1f, 0.2f, 0.3f, 0.4f }), m.GetRow(0)); Assert.AreEqual(new VectorF(new float[] { 5.0f, 6.0f, 7.0f, 8.0f }), m.GetRow(1)); Assert.AreEqual(new VectorF(new float[] { 9.0f, 10.0f, 11.0f, 12.0f }), m.GetRow(2)); m.SetRow(1, new VectorF(new float[] { 0.4f, 0.5f, 0.6f, 0.7f })); Assert.AreEqual(new VectorF(new float[] { 0.1f, 0.2f, 0.3f, 0.4f }), m.GetRow(0)); Assert.AreEqual(new VectorF(new float[] { 0.4f, 0.5f, 0.6f, 0.7f }), m.GetRow(1)); Assert.AreEqual(new VectorF(new float[] { 9.0f, 10.0f, 11.0f, 12.0f }), m.GetRow(2)); m.SetRow(2, new VectorF(new float[] { 0.7f, 0.8f, 0.9f, 1.0f })); Assert.AreEqual(new VectorF(new float[] { 0.1f, 0.2f, 0.3f, 0.4f }), m.GetRow(0)); Assert.AreEqual(new VectorF(new float[] { 0.4f, 0.5f, 0.6f, 0.7f }), m.GetRow(1)); Assert.AreEqual(new VectorF(new float[] { 0.7f, 0.8f, 0.9f, 1.0f }), m.GetRow(2)); }
public void GetRow() { MatrixF m = new MatrixF(3, 4, rowMajor, MatrixOrder.RowMajor); Assert.AreEqual(new VectorF(new float[] { 1.0f, 2.0f, 3.0f, 4.0f }), m.GetRow(0)); Assert.AreEqual(new VectorF(new float[] { 5.0f, 6.0f, 7.0f, 8.0f }), m.GetRow(1)); Assert.AreEqual(new VectorF(new float[] { 9.0f, 10.0f, 11.0f, 12.0f }), m.GetRow(2)); }
public void MultiplyVectorOperator() { VectorF v = new VectorF(new float[] { 2.34f, 3.45f, 4.56f, 23.4f }); Assert.AreEqual(v, MatrixF.CreateIdentity(4, 4) * v); Assert.AreEqual(new VectorF(4), new MatrixF(4, 4) * v); float[] values = new float[] { 12, 23, 45, 56, 67, 89, 90, 12, 43, 65, 87, 43, 34, -12, 84, 44 }; MatrixF m = new MatrixF(4, 4, values, MatrixOrder.RowMajor); Assert.IsTrue(VectorF.AreNumericallyEqual(v, m * m.Inverse * v)); for (int i = 0; i < 4; i++) Assert.IsTrue(Numeric.AreEqual(VectorF.Dot(m.GetRow(i), v), (m * v)[i])); }