Hi, 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). Thanks, ACAR
For more complete information about compiler optimizations, see our Optimization Notice.