are the clapack functions implemented in the ATLAS libraries supposed to work with large matrices (i.e. larger than 65536x65536)? Calling ATLAS function clapack_dposv() and clapack_dgesv() with larger sizes generates a segfault after running ~ 1 minute on both 3.8.4 stable and 3.9.51 devel here. The test was done on a 1TB machine. ATLAS was compiled with gcc 4.6.2 and -DATL_USE64BITS . I checked first that there was no problem on my side to allocate such large matrices and access matrix elements with size_t array indices on the same platform.
I read somewhere that the netlib Lapack itself is supposed to support matrix sizes in excess of 100,000x100,000. So is the limitation in the ATLAS Lapack subset or in the ATLAS implementation of BLAS?
Thanks in advance for your help and keep up the good work!