I am hoping to repeatedly hammer an X=B(A^-1) problem. That is, solve a linear system. For C++, which numerical solvers have support for 128bit long doubles (quads)?
Using C style arrays is a major plus as all my 2D data is stored as a single std::vector.
I was hoping to compile the code on linux with either GCC or ICC.