Estou tentando diagonalizar matrizes com precisão quádrupla e usar seus logaritmos. Existe um idioma no qual eu possa fazer isso usando funções internas?

Observe que os idiomas / pacotes nas tags são insuficientes, sofrendo das seguintes deficiências:

Matlab: Não suporta precisão quad.

Python / NumPy / SciPy: matrizes com o tipo dtype float128 produzem autovetores em float6

Sage: A interface através do GP / PARI produz mensagens de erro enigmática

lguém já realizou logaritmos de diagonalização e matriz com precisão quad e, em caso afirmativo, com