I have to invert large, but densely populated matrices with higher precision arithmetic. Therefore I am looking for something like the PLASMA library, which can do Cholesky or LU factorization in parallel on shared memory machines. Unfortunately it does support only classical float or double arithmetic. Therefore I am looking for a templated library with the same features, but the ability to use arbitrary precision data types.
Thanks a lot!