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?