- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content
I'm using MKL versiion 10.2. I'm trying a program that calls Pardiso multiple times. Each time has multiple right hand sides. PARDISO works good for the first call and it gaves access violation error for second call. As the error happens, the program will freeze the computer and I have to push the power button to turn off the computer. I set PT()=0 after the first call. I follow the procedure in pardiso_sym_f.f to call Pardiso. Any help?
WaterDog
WaterDog
Link Copied
2 Replies
- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content
Update: the problem is solved. PT()=0 is actually needed after releasing memory (phase =-1).
- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Report Inappropriate Content
Hi,
Sorry for delay. In MKL manual pt describe in next way: "ptis an integer array with 64 entries. It is very important that the pointer is initialized with zero at the first call of...". So you need to initialize pt by zero before call to pardiso, so if you call pardiso several times you need to set it by zero before each initialization step.
With best regards,
Alexander Kalinkin

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