From d@tr7320 @end|ng |rom un|@@ydney@edu@@u  Sat Aug 14 02:00:04 2021
From: d@tr7320 @end|ng |rom un|@@ydney@edu@@u (Dario Strbenac)
Date: Sat, 14 Aug 2021 00:00:04 +0000
Subject: [Rd] [External]  svd For Large Matrix
In-Reply-To: <alpine.DEB.2.22.394.2108130958100.3514@luke-Latitude-7480>
References: <SYBPR01MB4761CA7FD01EB93E4CF8511ACDFA9@SYBPR01MB4761.ausprd01.prod.outlook.com>
 <alpine.DEB.2.22.394.2108130958100.3514@luke-Latitude-7480>
Message-ID: <SYBPR01MB4761D9277370E4FDF84CA821CDFA9@SYBPR01MB4761.ausprd01.prod.outlook.com>

Good day,

Ah, I was confident it wouldn't be environment-specific but it is. My environment is

R version 4.1.0 (2021-05-18)
Platform: x86_64-pc-linux-gnu (64-bit)
Running under: Debian GNU/Linux 10 (buster)

Matrix products: default
BLAS:   /usr/lib/x86_64-linux-gnu/atlas/libblas.so.3.10.3
LAPACK: /usr/lib/x86_64-linux-gnu/atlas/liblapack.so.3.10.3

It crashes at about 180 GB RAM usage. The server has 1024 GB physical RAM in it. Modestly downsampling to 30 million cells avoids the segmentation fault. The segmentation fault originates from BLAS

Program received signal SIGSEGV, Segmentation fault.
0x00007ffff7649c10 in ATL_dgecopy () from /usr/lib/x86_64-linux-gnu/libblas.so.3

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


[Rd] [External] svd For Large Matrix
Dario Strbenac