Implements Singular Value Decomposition for generic number types, such as
Quaternions. It internally overloads several Base functions such that existing methods (
svdvals) should work directly.
It uses a Golub-Kahan 2-stage algorithm of bidiagonalization with Householder reflections, followed by an implicit QR with shift.
Based on initial code by Andreas Noack.
12 days ago