https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/What-is-the-type-of-variable-pt-required-by-pardiso/m-p/1223983#M30258 <P>As you can see from the examples in the folder<STRONG> ...\mkl\examples\solverf</STRONG>, or in the include file <STRONG>mkl_pardiso.fi,</STRONG> the type of pt is&nbsp;<STRONG>MKL_PARDISO_HANDLE</STRONG> . Since this is not an intrinsic type, you will need a definition for it in the subprogram that calls Pardiso (for example, by inserting&nbsp;<STRONG>USE mkl_pardiso</STRONG>).</P> Thu, 29 Oct 2020 23:54:03 GMT mecej4 2020-10-29T23:54:03Z https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/What-is-the-type-of-variable-pt-required-by-pardiso/m-p/1223969#M30257 <P>I am attempting to invert a matrix using pardiso in intel MKL, Fortran. (ifort 2017)</P> <P>&nbsp;</P> <LI-CODE lang="fortran">subroutine test use vars use vars_parallel use lapack95, only: DGBSV_F95 implicit none include 'mkl.fi' integer ii,im,in,ip,ix,iy,ik,ir,ic,ixr,ixc,ipr,ipc,ikr,ikc,iter,niter,i,j,n,kl,ku,ldab, ipiv, ldb integer is,l_success,l_ss,info,i1,i2,ndiag,IROWB,irow,icol,icount,nband integer maxfct,mnum,mtype,phase,nrhs,iparm(1:64),msglvl,error,irow_previous integer*8 pt(1:64) !Band storage: One diagonal + kl lowers + ku uppers kl=5+2 ku=5+2 n=nx*np ldab=2*kl+ku+1 nband=kl+ku+1 pt=0 ia(n+1)=icount+1 maxfct=1 mnum=1 mtype=1 nrhs=1 msglvl=0 call pardisoinit(pt, mtype, iparm) call pardiso_d(pt,maxfct,mnum,mtype,phase,n,Apimp,ia,ja,perm,nrhs,iparm,msglvl,Bimp,xsol,error) end</LI-CODE> <P>The compiler complains about the variable "pt". It says "The type of the actual argument differs from the type of the dummy argument".</P> <P>According to the documentation, pt is an integer array on 32-bit systems and an integer*8 on 64-bit systems. My system is 64-bit. I have also tried declaring pt as an integer (rather than integer*8).</P> <P>How can I find out what type ot shoudl really be?</P> <P>Thanks.</P> Thu, 29 Oct 2020 23:07:43 GMT https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/What-is-the-type-of-variable-pt-required-by-pardiso/m-p/1223969#M30257 mark10 2020-10-29T23:07:43Z https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/What-is-the-type-of-variable-pt-required-by-pardiso/m-p/1223983#M30258 <P>As you can see from the examples in the folder<STRONG> ...\mkl\examples\solverf</STRONG>, or in the include file <STRONG>mkl_pardiso.fi,</STRONG> the type of pt is&nbsp;<STRONG>MKL_PARDISO_HANDLE</STRONG> . Since this is not an intrinsic type, you will need a definition for it in the subprogram that calls Pardiso (for example, by inserting&nbsp;<STRONG>USE mkl_pardiso</STRONG>).</P> Thu, 29 Oct 2020 23:54:03 GMT https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/What-is-the-type-of-variable-pt-required-by-pardiso/m-p/1223983#M30258 mecej4 2020-10-29T23:54:03Z https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/What-is-the-type-of-variable-pt-required-by-pardiso/m-p/1224281#M30263 <P>Thanks for the speedy reply, Mecej, problem solved. Now I'll sift through the runtime errors...</P> Sat, 31 Oct 2020 01:44:22 GMT https://community.intel.com/t5/Intel-oneAPI-Math-Kernel-Library/What-is-the-type-of-variable-pt-required-by-pardiso/m-p/1224281#M30263 mark10 2020-10-31T01:44:22Z