Recently due to my research project, I want to calculate the inverse of an matrix in higher precision. I think it is best to use quadruple precision.
I learned that the ifort can define DOUBLE PRECISION declaration as REAL(KIND=16) ,i.e. in quadruple precision.
It seems that the source code of Lapack and my code which calls the subroutine in Lapack, can be compiled using -double-size 128 option.
However, I do not think Intel MKL offers quadruple precision version of Lapack.
Can this be done with MKL? If not, what is the possible way?
Link Copied
yuriisig wrote:Thanks very much! If you like, could you share your code with me? I have PM you my email account.I have the fast program of diagonalization of quad-precision accuracy which is partially based on BLAS and LAPACK (x64). For example, my qgemm only in 13 times more slowly dgemm Intel MKL.
For more complete information about compiler optimizations, see our Optimization Notice.