R example:
require(MASS)
l=30000
w=30000
x=rpois(l*w,0.00126)
bigmatrix=matrix(x,nrow=l,ncol=w)
inverted=ginv(bigmatrix)
This causes R to crash inside the La.svd routine (which is not based on LAPACK 3.4.2). So I tried it in numpy:
Numpy Example in python:
import numpy as np
N = 30000
b = np.random.random_integers(0,3,size=(N,N))
c = numpy.linalg.pinv(b)
Segmentation fault (core dumped)
This also causes a crash.