[Rd] svd For Large Matrix

[Dario Strbenac]

Good day,

I have a real scenario involving 45 million biological cells (samples) and 60 proteins (variables) which leads to a segmentation fault for svd. I thought this might be a good example of why it might benefit from a long vector upgrade.

test <- matrix(rnorm(45000000*60), ncol = 60)
testSVD <- svd(test)

 *** caught segfault ***
address 0x7fe93514d618, cause 'memory not mapped'

Traceback:
 1: La.svd(x, nu, nv)
 2: svd(test)

--------------------------------------
Dario Strbenac
University of Sydney
Camperdown NSW 2050
Australia