B. W. Lewis

Hi! Sorry you ran into trouble, it does look like a bug :(. I'm looking in to this and will let you know as soon as i have a workaround...

Hi Aaron, This is great! I'll try to help review and test your code; I'm not great with C++ but I can at least help with testing, etc. On 6/29/21,...

p.s. would you consider being also a co-maintainer of the original irlba package since I am so negligent in updates sometimes? You're into this code quite a bit and have...

Sure no problem about the package; there are some subtle ideas there. I like the idea of a write up of a basic intro to projected/partial SVDs (IRLBA for dummies)....

Hi, sorry about this behavior! The good news is it should be easy to diagnose. I'm investigating... On Mon, Nov 2, 2020 at 6:02 PM Elio Campitelli wrote: > I'm...

Sorry, can you give a code example and I will investigate? Thanks.

Sorry I should be more specific, I don't follow which example you're referring to. There are several examples there and none of them have exactly 200 singular values (some many...

Hmmmm. I wrote this test: ```{r} ibrary(irlba) f = function(seed) { d = c(1.2010e+05,3.3975e+04,1.5586e+04,8.3248e+03,7.5574e+03,7.0686e+03,6.2953e+03,6.1647e+03,5.5930e+03,5.3994e+03,4.9431e+03,4.6991e+03,4.3388e+03,4.1773e+03,3.9298e+03,3.8293e+03,3.7712e+03,3.7303e+03,3.6607e+03,3.5395e+03,3.4908e+03,3.3796e+03,3.2971e+03,3.2209e+03,3.1451e+03,3.0029e+03,2.8943e+03,2.8009e+03,2.7219e+03,2.6431e+03,9.7071e-10,6.6343e-10,2.8837e-10,1.8375e-10,1.0578e-10,1.0274e-10,7.8556e-11,7.1549e-11,5.6877e-11,4.8820e-11,4.0724e-11,3.5984e-11,3.4707e-11,2.8522e-11,2.5319e-11,2.5076e-11,2.3933e-11,2.3539e-11,2.1689e-11,2.1335e-11,2.0271e-11,1.9699e-11,1.9426e-11,1.8862e-11,1.8385e-11,1.8283e-11,1.7684e-11,1.7657e-11,1.6216e-11,1.5041e-11,1.3888e-11,1.3781e-11,1.2863e-11,1.2862e-11,1.2635e-11,1.2570e-11,1.2540e-11,1.1974e-11,1.1850e-11,1.1368e-11,1.1355e-11,1.1192e-11,1.0941e-11,1.0780e-11,1.0678e-11,1.0511e-11,9.8443e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12,9.8073e-12) set.seed(seed) N = 3 tol = 1e-15 m = length(d) u = qr.Q(qr(matrix(rnorm(m ^ 2), m)))...

Thanks for this suggestion, will put this on the todo list for the next update. On 4/29/21, leahrosen ***@***.***> wrote: > Hi, I've rewritten my code to replace prcomp with...

That looks like a bug. Can you tell me which version of the package (CRAN or GitHub) and the output of `sessionInfo()`?