Mercurial > octave
comparison scripts/general/quad2d.m @ 31205:b0e90ca8e679 stable
quad2d: fix unintended complex conjugate return (bug #62972)
quad2d: use .' instead of ' to avoid complex conjugate in q and qerr.
| author | Arun Giridhar <arungiridhar@gmail.com> |
|---|---|
| date | Sun, 28 Aug 2022 12:21:17 -0400 |
| parents | 796f54d4ddbf |
| children | 1c40fc2344f4 |
comparison
equal
deleted
inserted
replaced
| 31200:075443476dfb | 31205:b0e90ca8e679 |
|---|---|
| 405 x = ones (15, 1) * x; | 405 x = ones (15, 1) * x; |
| 406 tx = ones (15,1) * tx; | 406 tx = ones (15,1) * tx; |
| 407 ty = ty * ones (1, 15); | 407 ty = ty * ones (1, 15); |
| 408 | 408 |
| 409 z = yhalfwidth .* f (x, y) .* ztrans(tx, ty) .* xhalfwidth; | 409 z = yhalfwidth .* f (x, y) .* ztrans(tx, ty) .* xhalfwidth; |
| 410 q = weights15 * (weights15 * z)'; | 410 q = weights15 * (weights15 * z).'; |
| 411 qerr = abs (weights7 * (weights7 * z)' - q); | 411 qerr = abs (weights7 * (weights7 * z).' - q); |
| 412 | 412 |
| 413 endfunction | 413 endfunction |
| 414 | 414 |
| 415 | 415 |
| 416 %!shared f | 416 %!shared f |