- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content

Hi all,

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).

I run Microsoft Visual Studio 2010 Version 10.040219.1 SP1Rel on Windows 7 SP1, with an Intel Xeon X5650.

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!

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content

Link Copied

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content

*Does PARDISO have problems if a complex and unsymmetric matrix happens to have specific symmetries?*To answer that would lead one to indulge oneself in idle speculation.

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]

- Subscribe to RSS Feed
- Mark Topic as New
- Mark Topic as Read
- Float this Topic for Current User
- Bookmark
- Subscribe
- Printer Friendly Page