I'm having problems with the "Pardiso" subroutine in Intel Visual Fortran Composer XE 2011 Update 8.
I get an "access violation reading location" during the computation of the "non-zeros for LL^T factorization," at the very beginning of the analysis phase.
I guess the problem is related either to my MKL configuration or to the way I'm declaring some variables.
I'm attaching a sample project (a simplified version of the actual code).
The options I chose in the "Properties" window of my project are in the txt file "Properties".
The other text files are the data used by the project (the size of the matrices involved is also very small for practical reasons).
Thank you for your help!
Your data initializations do not correctly describe the packed column form of the diagonal matrix. The MKL Reference should be consulted regarding compact columns storage conventions.
Try [fortran] ia(1:5) = (/ 1, 2, 3, 4, 5 /); ia(6:n+1)=nnz+1 ja(1:5) = (/ 1, 2, 3, 4, 5 /) a(1:5) = (/ 1d0, 1d0, 1d0, 1d0, 1d0 /); a(6:nnz)=0d0 [/fortran]