The rule for SVD is here. It would probably be very simple to add one for svdvals
, perhaps calling the same function to do the work svd_rev(nt, NoTangent(), S̄, NoTangent())
or perhaps just keeping the few lines which are nonzero there.
1 Like