Skip to content

*chklq: fix test failures when NB=1 with M=N#1298

Open
jschueller wants to merge 1 commit into
Reference-LAPACK:masterfrom
jschueller:issue973
Open

*chklq: fix test failures when NB=1 with M=N#1298
jschueller wants to merge 1 commit into
Reference-LAPACK:masterfrom
jschueller:issue973

Conversation

@jschueller

Copy link
Copy Markdown
Contributor

The LQ test drivers (dchklq, schklq, cchklq, zchklq) called DGELS/SGELS which re-factorizes AF internally using GEQRF when M>=N. This corrupted the LQ factorization expected by subsequent LQT02/LQT03 calls on later K iterations, causing test(1) and test(2) failures.

Fix by re-copying A to AF and re-factorizing with GELQF after the GELS call, restoring the LQ factorization for the remaining tests.

Closes #973

The LQ test drivers (dchklq, schklq, cchklq, zchklq) called DGELS/SGELS
which re-factorizes AF internally using GEQRF when M>=N. This corrupted
the LQ factorization expected by subsequent LQT02/LQT03 calls on later
K iterations, causing test(1) and test(2) failures.

Fix by re-copying A to AF and re-factorizing with GELQF after the GELS
call, restoring the LQ factorization for the remaining tests.

Closes Reference-LAPACK#973
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

fails of tests for QR and LQ decompositions

1 participant