Re: [Lapackpp-devel] Interesting crash in QR solve

[Michael P. <mr...@cm...>]

Christian,

Thanks for responding so quickly and thoughtfully.  As for negative 
submatrix increments, I'm not using them, but I also may have reported 
the wrong failed assertion; it may have been "assert(cI.end() < 
mat.size(0))" on line 273. 

For now I've only been doing overdetermined systems (for example: 50 to 
200 equations, 3 unknowns).  I think the change of linslv.cc:220 makes 
sense to correctly initialize Xtmp for either m >= n or m < n.  I wonder 
why that was changed to Xtmp = B?

Anyway, changing the work vector size for the QR solve routine to agree 
with the manual (as you suggested) eliminated the crash I was 
experiencing.  I will take a closer look at it later, since we will also 
be using the sparse LU solver and eigenvalue functions (which I 
understand we may need to debug somewhat).  I'll keep you posted.

Thanks
Michael