Estoy intentando diagonalizar las matrices con precisión cuádruple y tomar sus logaritmos. ¿Hay algún lenguaje en el que pueda lograr esto usando las funciones integradas?

Nota, los idiomas / paquetes en las etiquetas son insuficientes y sufren las siguientes deficiencias:

Matlab: no admite precisión cuádruple.

Python / NumPy / SciPy: las matrices con dtype float128 producen vectores propios en float64.

Sage: la interfaz a través de GP / PARI genera mensajes de error crípticos.

¿Alguien ha realizado diagonalización y logaritmos de matriz con precisión cuádruple, y si es así, cómo?