Open
Description
Description
The documentation of the blocked LQ decomposition methods ?gelqt contains a typo. The size of the work array should be mb * m instead of mb * n. If the latter value is used for matrices where m > n, all kinds of segmentation faults or memory corruption errors occur. These disappear completely when a work vector of length mb * m is used instead, independent of the value of mb.
Checklist
- I've included a minimal example to reproduce the issue
- I'd be willing to make a PR to solve this issue