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?
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.