From d@tr7320 @end|ng |rom un|@@ydney@edu@@u  Fri Aug 13 11:00:05 2021
From: d@tr7320 @end|ng |rom un|@@ydney@edu@@u (Dario Strbenac)
Date: Fri, 13 Aug 2021 09:00:05 +0000
Subject: [Rd] svd For Large Matrix
Message-ID: <SYBPR01MB4761CA7FD01EB93E4CF8511ACDFA9@SYBPR01MB4761.ausprd01.prod.outlook.com>

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



[Rd] svd For Large Matrix
Dario Strbenac