Bug in csd2by1 functions
Posted: Tue Oct 14, 2014 5:45 amI 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.
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.