We have many 2 and 3 dimensional arrays which are allocated based on the size of the problem. We are running into memory overflow errors with large problems. These arrays are very sparse so I only need to store the non-zero elements. Is there a sparse matrix handler in Intel Fortran I could use, preferably with the same syntax, eg:
a(i,j,k) = 123.456
or xyz = a(i,j,k)
to save replacing the occurrences in the code?
Adrian
