I have code linked to CLAPACK 3.2.1. I'm running in to a problem with dlarfp. I've seen
and I'm hoping someone can confirm the problem I'm seeing is related to the bug (or another bug documented somewhere else) so I can determine if I should port to something newer (lapacke or magma) or try to understand how to avoid the situation if it's something I'll run in to anyway. (Obviously I *should* upgrade at some point, but I also hope to graduate this May...)
The issue is that I'm doing a least squares regression using dgels. On occasion the code does not return but computes for hours. When I attach to the process with GDB I'm able to 'next' down a few stack frames until I get to dlarfp_(), but the code never returns from dlarfp. as shown below (sorry, I lost the rest of the info above the scrollback).
If there's info I could get to help I can try to dig it out, but the algorithm is a genetic algorithm and so the B array is randomly generated.
In short, is the bug I mentioned, or another one I'm not aware of, one which would cause computation to carry on without stoping?
Thanks for the help.