The LAPACK forum has moved to https://github.com/Reference-LAPACK/lapack/discussions.

Bug in csd2by1 functions

Post here if you want to report a bug to the LAPACK team

Bug in csd2by1 functions

Postby chereshnev » Tue Oct 14, 2014 5:45 am

I found a bug in new functions for CS decomposition from LAPACK 3.5.0 (functions cuncsd2by1, dorcsd1by1, sorcsd2by1, zuncsd2by1)
Reproducer for bug is attached. It tests orthogonality of a matrix U2 after call of function dorcsd2by1.
The bug is caused by overwriting of work array. Proposed fix is attached as diff file.
Attachments
csd2by1_fix_diff.txt
(4.21 KiB) Downloaded 424 times
dorcsd2by1_bug.cpp
(2.67 KiB) Downloaded 445 times
chereshnev
 
Posts: 14
Joined: Thu Oct 09, 2014 6:17 am
Location: Intel Corp., Russia, Novosibirsk

Return to Bug report

Who is online

Users browsing this forum: No registered users and 1 guest