Mercurial > octave-nkf
changeset 13042:ca7aaf2689c3
codesprint: 8 tests for pinv.cc
| author | Jordi Gutiérrez Hermoso <jordigh@gmail.com> |
|---|---|
| date | Sat, 03 Sep 2011 10:09:50 -0500 |
| parents | e5a49b2f8225 |
| children | 098c8e696236 |
| files | src/DLD-FUNCTIONS/pinv.cc |
| diffstat | 1 files changed, 20 insertions(+), 0 deletions(-) [+] |
line wrap: on
line diff
--- a/src/DLD-FUNCTIONS/pinv.cc Sat Sep 03 07:58:26 2011 -0700 +++ b/src/DLD-FUNCTIONS/pinv.cc Sat Sep 03 10:09:50 2011 -0500 @@ -169,3 +169,23 @@ return retval; } + +/* +%!shared a, b, hitol, d, u, x, y, tolmult +%! a = reshape (rand*[1:16], 4, 4); ## Rank 2 matrix +%! b = pinv (a); +%! hitol = sqrt(eps); +%! d = diag ([rand, rand, hitol, hitol]); +%! u = rand (4); ## Could be singular by freak accident +%! x = inv (u)*d*u; +%! y = pinv (x, sqrt(eps)); +%! tolmult = 15; +%!assert(a*b*a, a, tolmult*eps); +%!assert(b*a*b, b, tolmult*eps); +%!assert((b*a)', b*a, tolmult*eps); +%!assert((a*b)', a*b, tolmult*eps); +%!assert(x*y*x, x, tolmult*hitol); +%!assert(y*x*y, y, tolmult*hitol); +%!assert((x*y)', x*y, tolmult*hitol); +%!assert((y*x)', y*x, tolmult*hitol); +*/