public void SolveMatrix10()
{
int i, j;
double e, me;
DoubleSymmetricLevinson dsl = new DoubleSymmetricLevinson(T10);
// check inverse
DoubleMatrix I = dsl.Solve(DoubleMatrix.CreateIdentity(10));
me = 0.0;
for (i = 0; i < dsl.Order; i++)
{
for (j = 0; j < dsl.Order; j++)
{
e = System.Math.Abs((I10[i, j] - I[i, j]) / I10[i, j]);
if (e > me)
{
me = e;
}
}
}
Assert.IsTrue(me < Tolerance10, "Maximum Error = " + me.ToString());
}
public void MismatchRowsTestforSolveMatrix()
{
DoubleSymmetricLevinson dsl = new DoubleSymmetricLevinson(T10);
DoubleMatrix X = dsl.Solve(I5);
}
public void NullParameterTestforSolveMatrix()
{
DoubleSymmetricLevinson dsl = new DoubleSymmetricLevinson(T10);
DoubleMatrix X = dsl.Solve(null as DoubleMatrix);
}
public void SolveVector10()
{
int i;
double e, me;
DoubleSymmetricLevinson dsl = new DoubleSymmetricLevinson(T10);
DoubleVector X = dsl.Solve(Y10);
// determine the maximum error
me = 0.0;
for (i = 0; i < dsl.Order; i++)
{
e = System.Math.Abs((X10[i] - X[i]) / X10[i]);
if (e > me)
{
me = e;
}
}
Assert.IsTrue(me < Tolerance10, "Maximum Error = " + me.ToString());
}
public void MismatchRowsTestforSolveVector()
{
DoubleSymmetricLevinson dsl = new DoubleSymmetricLevinson(T10);
DoubleVector X = dsl.Solve(X5);
}
public void NullParameterTestforSolveVector()
{
DoubleSymmetricLevinson dsl = new DoubleSymmetricLevinson(T10);
DoubleVector X = dsl.Solve(null as DoubleVector);
}