Dear LAPACK team,

I am sorry that I can't send in a report that uses the LAPACK routines directly, but I am not sure how to do this.

I can report two equivalent reproducible examples, one in numpy and one in R. Both are using LAPACK, to do the SVD within the pseudoinverse calculation:

See the attached code for the reproducible examples.

I am running all of this on a machine with 128Gb or RAM. These matrices require around 40Gb to process. I have tried this with OPENMP num threads of 1 or higher, and

the result is the same.

Do you have any large-memory machines you can use to try and reproduce this? Is there any reason why it does not like large matrices?

Thanks so much,

Cyrus