public void SolveMatrix10()
    {
      int i, j;
      float e, me;
      ComplexFloatLevinson cfl = new ComplexFloatLevinson(LC10, TR10);

      // check inverse
      ComplexFloatMatrix I = cfl.Solve(ComplexFloatMatrix.CreateIdentity(10));
      me = 0.0f;
      for (i = 0; i < cfl.Order; i++)
      {
        for (j = 0; j < cfl.Order; j++)
        {
          e = ComplexMath.Absolute((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()
 {
   ComplexFloatLevinson cfl = new ComplexFloatLevinson(LC10, TR10);
   ComplexFloatMatrix X = cfl.Solve(I5);
 }
 public void NullParameterTestforSolveMatrix()
 {
   ComplexFloatLevinson cfl = new ComplexFloatLevinson(LC10, TR10);
   ComplexFloatMatrix X = cfl.Solve(null as ComplexFloatMatrix);
 }
 public void SolveVector10()
 {
   int i;
   float e, me;
   ComplexFloatLevinson cfl = new ComplexFloatLevinson(LC10, TR10);
   ComplexFloatVector X = cfl.Solve(Y10);
   
   // determine the maximum error
   me = 0.0f;
   for (i = 0; i < cfl.Order; i++)
   {
     e = ComplexMath.Absolute((X10[i] - X[i]) / X10[i]);
     if (e > me)
     {
       me = e;
     }
   }
   Assert.IsTrue(me < Tolerance10, "Maximum Error = " + me.ToString());
 }
 public void MismatchRowsTestforSolveVector()
 {
   ComplexFloatLevinson cfl = new ComplexFloatLevinson(LC10, TR10);
   ComplexFloatVector X = cfl.Solve(X5);
 }
 public void NullParameterTestforSolveVector()
 {
   ComplexFloatLevinson cfl = new ComplexFloatLevinson(LC10, TR10);
   ComplexFloatVector X = cfl.Solve(null as ComplexFloatVector);
 }