Собственные значения, собственные векторы и логарифмы матриц с четверной точностью

Я пытаюсь диагонализировать матрицы с четырехкратной точностью и взять их логарифмы. Есть ли язык, на котором я могу сделать это с помощью встроенных функций?

Примечание, языки/пакеты в тегах недостаточны, страдают от следующих недостатков:

Matlab: Не поддерживает квадратичную точность.

Python/NumPy/SciPy: Матрицы с dtype float128 дают собственные векторы в float64.

Sage: Интерфейс через GP/PARI выдает загадочные сообщения об ошибках.

Кто-нибудь выполнял диагонализацию и логарифмы матриц с квадратичной точностью, и если да, то как?