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

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

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

Python / NumPy / SciPy: матрицы с d-типом float128 дают собственные векторы в float64.

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

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