Фаворит / Лучший AX = B Солвер, который поддерживает длинные пары?

Question

Я надеюсь повторно решить проблему X = B (A ^ -1). То есть решить линейную систему. Для C ++, какие числовые решатели поддерживают 128-битные двойные числа (квадраты)?

Использование массивов в стиле C является основным плюсом, так как все мои 2D данные хранятся в виде одного std :: vector.

Я надеялся скомпилировать код на linux с GCC или ICC.

Answer 1

Многие библиотеки линейной алгебры C ++ основаны на шаблонах, включая NT2, Boost.uBLAS, Eigen (см. Каковы наиболее широко используемые библиотеки C ++ для векторной / матричной математики / линейной алгебры, а также их соотношение цены и качества? 1002 * для ссылок). Таким образом, они должны иметь возможность поддерживать квадраторы, если ваш компилятор / библиотека может выполнять математику с квадраторами. Например, в Eigen тип Eigen::Matrix<long double, Dynamic, Dynamic> обозначает матрицу произвольного размера, содержащую длинные двойники, и вы можете использовать стандартные функции для решения с такими матрицами.