B. W. Lewis

Do you have any NA's in your matrix (see https://github.com/bwlewis/irlba/issues/14)? (Also, can you try the current CRAN version, 2.3.5?).

thankd for this... will investigate. On Fri, Apr 27, 2018, 18:07 ilyakorsunsky wrote: > Hi, for my larger datasets (250,000 x 2000) run with the R code > (fastpath=FALSE), I...

Yes indeed, I can replicate these behaviors with badly scaled data due to floating point overflow. For example: ```{r} x = rep(sqrt(.Machine$double.xmax) * 10, 2) # now its 2-norm: sqrt(drop(crossprod(x)))...

Anyone want to help with this?

thanks for this report, working on it...

Aaron's bug fixes may have fixed this for you...Sorry it took so long!

Thanks for this! I agree a speedup is desirable in the sparse case, and the order alignment with prcomp is a good idea. I'm carefully considering the @LTLA 's comments....

Hi Jan, Aaron: Sorry I have been so negligent in maintenance here -- it has been a bit of a crazy year (indeed for the whole world). Jan, I have...

Thanks, investigating now.

I get: ```{r} library(irlba) packageVersion('irlba') ## [1] ‘2.3.3’ Xh