I am writing a finite element application using IVF and the PARDISO linear solver that comes with the MKL libraries. My application was working fine until 3rd January 2011 when I upgraded tow_mkl_10.3.1.127_ia32. Following the upgrade the solver failed at the numerical factorisation phase (PHASE=22). Instructions in the userguide tell you that to set the default parameters one should pass IPARM(1)=0 at the analysis phase (PHASE=11) however on inspection of IPARM following the analysis phase I found that some elements of IPARM had not been set to their default values, for example IPARM(28) which only admits values of 0 or 1 with 0 being the default remained unchanged at -1 which was the value present in the uninitialised array. In seeking a solution to this issue I replaced the initialisation IPARM(1)=0 with IPARM=0, i.e. setting all elements of IPARM to zero, prior to performing the analysis phase. This appeared to solve the problem.
Inspection of the results produced by PARDISO seems to indicate that they are fine with this work-around but I am now left with some uncertainty as to whether or not IPARM is being correctly initialised to appropriate default values (it could be the case that by setting IPARM=0 allowed PARDISO to run for my particular problem configuration but some of the elements of IPARM are not now at sensible default values).
Yes that is what I observed - the code was working correctly until I installed the upgrade for mkl and for IVF. That is not to say, of course, that the vector IPARM did not enter phase 1 with different values as I did not initialise it explicitly since according to the documentation this is not required. Unfortunately now that I have installed the upgrades I cannot confirm the outcome of the initialisation for the older version of the libraries.
Almost a week since I asked a fairly important question and still no response (admittedly running over a holiday period)! Is anyone taking note of the question?
Just rebooted my machine and tried to find the mkl forum via the intel.com website - it took me numerous attempts via the intel search engine to find anything to do with forum. This is a general problem that I have with the intel website - it is not easy to return to your previous settings. Now I have found my way in I shall record it for future use but really should there be an issue with shoving "forum" into the intel searh engine?