Mercurial > octave
changeset 31945:2a75c771c451
maint: Merge stable to default.
author | Markus Mützel <markus.muetzel@gmx.de> |
---|---|
date | Fri, 31 Mar 2023 17:43:55 +0200 |
parents | c90de146a9ed (current diff) c687e0ba6606 (diff) |
children | 89850bb5eb31 |
files | |
diffstat | 1 files changed, 10 insertions(+), 0 deletions(-) [+] |
line wrap: on
line diff
--- a/libinterp/corefcn/lu.cc Tue Mar 28 11:56:03 2023 -0400 +++ b/libinterp/corefcn/lu.cc Fri Mar 31 17:43:55 2023 +0200 @@ -539,6 +539,16 @@ %! assert (u, single ([5, 6; 0, 4/5]), sqrt (eps ("single"))); %! assert (p(:,:), single ([0, 0, 1; 1, 0, 0; 0 1 0]), sqrt (eps ("single"))); +# complex matrix input +%!test +%! C = [1, 0, 1, 0; +%! 1i, 1/3, -1i, 1/3; +%! 1, 2i/3, 1, -2i/3; +%! 1i, -1/3, -1i, -1/3]; +%! [L, U, P] = lu (C); +%! assert (rcond (C), 1/8, eps); +%! assert (norm (P'*L*U - C, Inf) < eps); + %!testif HAVE_UMFPACK %! Bi = [1 2 3 4 5 2 3 6 7 8 4 5 7 8 9]; %! Bj = [1 3 4 5 6 7 8 9 11 12 13 14 15 16 17];