09-Nov-2018 07:32:15

ASA007_TEST
  MATLAB version
  Test the ASA007 library.

TEST01:
  SYMINV computes the inverse of a positive
  definite symmetric matrix.
  A compressed storage format is used.

  Here we look at the matrix A which is
  N+1 on the diagonal and
  N   on the off diagonals.

  Matrix order N = 1
  Maxtrix nullity NULLTY = 0
  RMS ( C * A - I ) = 1.110223e-16

  Matrix order N = 2
  Maxtrix nullity NULLTY = 0
  RMS ( C * A - I ) = 4.577567e-16

  Matrix order N = 3
  Maxtrix nullity NULLTY = 0
  RMS ( C * A - I ) = 4.742875e-16

  Matrix order N = 4
  Maxtrix nullity NULLTY = 0
  RMS ( C * A - I ) = 1.472877e-15

  Matrix order N = 5
  Maxtrix nullity NULLTY = 0
  RMS ( C * A - I ) = 2.045268e-15

  Matrix order N = 6
  Maxtrix nullity NULLTY = 0
  RMS ( C * A - I ) = 2.430159e-15

  Matrix order N = 7
  Maxtrix nullity NULLTY = 0
  RMS ( C * A - I ) = 5.028281e-15

  Matrix order N = 8
  Maxtrix nullity NULLTY = 0
  RMS ( C * A - I ) = 5.773160e-15

  Matrix order N = 9
  Maxtrix nullity NULLTY = 0
  RMS ( C * A - I ) = 3.804970e-15

  Matrix order N = 10
  Maxtrix nullity NULLTY = 0
  RMS ( C * A - I ) = 1.081835e-14

  Matrix order N = 11
  Maxtrix nullity NULLTY = 0
  RMS ( C * A - I ) = 2.068084e-14

  Matrix order N = 12
  Maxtrix nullity NULLTY = 0
  RMS ( C * A - I ) = 1.407139e-14

  Matrix order N = 13
  Maxtrix nullity NULLTY = 0
  RMS ( C * A - I ) = 2.061927e-14

  Matrix order N = 14
  Maxtrix nullity NULLTY = 0
  RMS ( C * A - I ) = 2.508551e-14

  Matrix order N = 15
  Maxtrix nullity NULLTY = 0
  RMS ( C * A - I ) = 2.622392e-14

TEST02:
  SYMINV computes the inverse of a positive
  definite symmetric matrix.
  A compressed storage format is used.

  Here we look at the Hilbert matrix
  A(I,J) = 1 / ( I + J - 1 )

  We expect errors to grow quickly with N.

  Matrix order N = 1
  Maxtrix nullity NULLTY = 0
  RMS ( C * A - I ) = 0.000000e+00

  Matrix order N = 2
  Maxtrix nullity NULLTY = 0
  RMS ( C * A - I ) = 4.937408e-16

  Matrix order N = 3
  Maxtrix nullity NULLTY = 0
  RMS ( C * A - I ) = 7.351357e-15

  Matrix order N = 4
  Maxtrix nullity NULLTY = 0
  RMS ( C * A - I ) = 3.339551e-13

  Matrix order N = 5
  Maxtrix nullity NULLTY = 0
  RMS ( C * A - I ) = 6.812212e-12

  Matrix order N = 6
  Maxtrix nullity NULLTY = 0
  RMS ( C * A - I ) = 9.882934e-11

  Matrix order N = 7
  Maxtrix nullity NULLTY = 0
  RMS ( C * A - I ) = 4.323682e-09

  Matrix order N = 8
  Maxtrix nullity NULLTY = 0
  RMS ( C * A - I ) = 2.151565e-07

  Matrix order N = 9
  Maxtrix nullity NULLTY = 0
  RMS ( C * A - I ) = 6.084266e-06

  Matrix order N = 10
  Maxtrix nullity NULLTY = 1
  RMS ( C * A - I ) = 1.408078e+01

  Matrix order N = 11
  Maxtrix nullity NULLTY = 1
  RMS ( C * A - I ) = 3.782687e+00

  Matrix order N = 12
  Maxtrix nullity NULLTY = 1
  RMS ( C * A - I ) = 3.925892e+00

  Matrix order N = 13
  Maxtrix nullity NULLTY = 1
  RMS ( C * A - I ) = 4.062417e+00

  Matrix order N = 14
  Maxtrix nullity NULLTY = 1
  RMS ( C * A - I ) = 4.193206e+00

  Matrix order N = 15
  Maxtrix nullity NULLTY = 1
  RMS ( C * A - I ) = 4.318995e+00

ASA007_TEST
  Normal end of execution.

09-Nov-2018 07:32:16