Mercurial > octave
view test/bug-49904.tst @ 31248:8b75954a4670
delaunayn: adjust node ordering for positive outward normal vectors (bug #53397)
* delaunayn.m: Check sign of simplex volume, flip node order for negative
volumes to ensure positive (outward-pointing) normal vectors. Add BISTs to
check for positive volumes.
* etc/News.8.md: Append function improvement note to delaunayn change
paragraph under General Improvements.
| author | Nicholas R. Jankowski <jankowski.nicholas@gmail.com> |
|---|---|
| date | Thu, 29 Sep 2022 23:09:05 -0400 |
| parents | 9080316864bf |
| children |
line wrap: on
line source
%!shared A %! A = [ ... %! 1.00003+0.00000i, -0.68292+0.26406i, -0.57670+0.25291i, 0.34986+0.03335i; %! -0.68292-0.26406i, 0.94539+0.00000i, 0.69006-0.03967i, -0.43847+0.38619i; %! -0.57670-0.25291i, 0.69006+0.03967i, 0.95260+0.00000i, -0.34734+0.40728i; %! 0.34986-0.03335i, -0.43847-0.38619i, -0.34734-0.40728i, 0.98356+0.00000i]; %! %!test <*49904> %! [Ainv, rcond] = inv (A); %! assert (ishermitian (Ainv)); %! assert (A * Ainv, eye (4, 'double'), eps ('double') / rcond) %!test <*49904> %! [Ainv, rcond] = inv (single (A)); %! assert (ishermitian (Ainv)); %! assert (A * Ainv, eye (4, 'single'), eps ('single') / rcond)