public void solve_for_diagonal()
{
float[][] value =
{
new float[] { 2.1f, 3.1f },
new float[] { 1.6f, 4.2f },
};
float[] rhs = { 6.1f, 4.3f };
float[][] expected =
{
new float[] { 6.63730669f, -3.45336843f },
new float[] { -2.528498f, 2.3393786f }
};
var target = new JaggedLuDecompositionF(value);
float[][] actual = target.SolveForDiagonal(rhs);
Assert.IsTrue(Matrix.IsEqual(expected, actual, 1e-6));
Assert.IsTrue(Matrix.IsEqual(value, target.Reverse(), 1e-6));
}
