public void DivideByDiagonalTest()
{
double[,] a =
{
{ 1, 0, 5 },
{ 1, 2, 1 },
{ 0, 6, 1 },
{ 2, 6, 5 },
{ 2, 1, 1 },
{ 5, 1, 1 }
};
double[] b = { 2, 6, 1 };
double[,] result = new double[6, 3];
Matrix.DivideByDiagonal(a, b, result);
double[,] expected = Matrix.Divide(a, Matrix.Diagonal(b));
Assert.IsTrue(expected.IsEqual(result, 1e-6));
result = Matrix.DivideByDiagonal(a.ToJagged(), b).ToMatrix();
expected = Matrix.Divide(a.ToJagged(), Jagged.Diagonal(b)).ToMatrix();
Assert.IsTrue(expected.IsEqual(result, 1e-6));
}
#endregion