- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content
Hi,
I have been using MKL non-symmetric solver for my finite element code (up to MKL Version10), but have a probelm with updated version of MKL when call pardiso_unsym ( reordering):
phase = 11 ! only reordering and symbolic factorization
CALL pardiso (pt, maxfct, mnum, mtype, phase, n, a, ia, ja,
1 idum, nrhs, iparm, msglvl, ddum, ddum, error)
Errors:
forrtl: severe (174): SIGSEGV, segmentation fault occurred
Image PC Routine Line Source
libpthread.so.0 00002AAAAC8B754B Unknown Unknown Unknown
Stack trace terminated abnormally.
------
The compiler link I used are:
Currently Loaded Modulefiles:
1) Barrine/2.0 2) mpt/2.00 3) intel-cmkl-11/12.12.1.2.8.273 4) intel-fc-11/12.1.2.8.273
COM = ifort -m64 -w -fpp -I/sw/sdev/intel/ComposerXE/composer_xe_2011_sp1.8.273/mkl/include
LIBS = -L"/sw/sdev/intel/ComposerXE/composer_xe_2011_sp1.8.273/mkl/lib/intel64" -lmkl_intel_lp64 \
-lmkl_intel_thread \
-lmkl_core \
-liomp5 -lpthread -lm
In addtion, my code works well with the SuperLU, but not with MKL Pardiso at the moment. Any advices to fix it? Thanks in advance
Link Copied
- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content
Hi Henry,
Could you provide testcase or matrix elements with iparm array to reproduce issue on our side? It could help us to investigate it
With best regards,
Alexander Kalinkin
- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content
Did you follow a checklist such as
to track this down? No doubt, you're aware that the most common situation is a requirement to increase stack limit or try heap-arrays option.
- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content
Thanks Alexander and Tim.
I had a probelm to post it here and thus I emailed you directly this Monday. I used to check the code up to item 5) as listed by Tim.
My probelm is not solved yet and thus I attached the matrix element and solver code (iparm) below. Look forward to your further advice.
Thanks
Henry
- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content
Hi,
I've turn on matrix checker on your matrix (iparm(27)=1) and got next problem:
*** Error in PARDISO (incorrect input matrix ) error_num= 24
*** Input check: j=4908159, ja(j)=35835, ja(j+1)=35830 are incompatible
The following ERROR was detected: -1
Looks like your ja array is incorrect.
With best regards,
Alexander Kalinkin
- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content
did you check the input consistency with matrix checker option? iparm(27) =1
- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content
Thanks. My probelm is solved after reordering JA.

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