I’m trying to use cusolver on fortran and I succeeded using a 3 by 3 matrix. However, when I try to use a larger, sparse matrix (348k by 348k), the program crashes with what I suspect was insufficent memory.
I tried with -stack=50000000
Would anyone please advise?
real, allocatable :: A(:,:),ATest(:,:) real, allocatable :: B(:,:) real, allocatable :: X(:,:) integer, allocatable::Ipiv(:) integer, target::Lwork real, allocatable :: IRN(:), JCN(:), val(:) ...... open(121, file='Input.dat') read(121, *) lda, ldb, nnz ! <--- 348780, 348780, 4492548 allocate(A(lda,ldb)) allocate(ATest(lda,ldb)) allocate(B(n,1)) allocate(X(n,1)) allocate(Ipiv(lda)) allocate(IRN(nnz)) allocate(JCN(nnz)) allocate(val(nnz)) A_size = SIZEOF(A) B_size = SIZEOF(B) X_size = SIZEOF(X) Ipiv_size = SIZEOF(Ipiv) devInfo_size = SIZEOF(devInfo) Lwork_size = SIZEOF(Lwork) do i=1, nnz read(121, *) IRN(i), JCN(i), val(i) col = IRN(i) row = JCN(i) A(row,col) = val(i) print *, row," ",col, " ", val(i) !<----- error happens when col > 8000 enddo close(121)