Closed
Description
Description
LAPACKE_?larfb
and LAPACKE_?larfb_work
routines compute an uplo value based on the input parameters in order to transpose the trapezoidal matrix V. This computation wrongly uses the parameter left
(side) instead of the correct forward
(direct).
In other words, instead of
uplo = ( ( left && col ) || !( left || col ) ) ? 'l' : 'u';
the code should do
uplo = ( ( forward && col ) || !( forward || col ) ) ? 'l' : 'u';
With this fix our Gonum tests (for Dlarfb) pass again.
Checklist
- I've included a minimal example to reproduce the issue
- I'd be willing to make a PR to solve this issue