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.

