Returns the default number of components to compute in a truncated SVD
for the given input matrix. Used by the CLAMP solvers when svd_k is
not provided explicitly. Other SVD contexts use their own heuristics.
select_svd_k(Y)An integer: max(2, floor((min(nrow(Y), ncol(Y)) - 1) / 4)).
select_svd_k(matrix(0, nrow = 100, ncol = 20))
#> [1] 4