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); }