在求解ODE时，我希望工作精度为40或更多小数位。 谁能告诉我是否有可能使用gsl ODE求解器获得更高的精度结果，也许使用libquadmath？从http://www.gnu.org/software/gsl/manual/gsl-ref.html#Driver处的gsl_odeiv2文档中可以看出它们只采用了双打，但情况确实如此？
如果没有人能提出替代方案吗？像'odeint'这样的东西支持高于双精度吗？我是否必须从头开始编写自己的文章？
（供参考我的具体问题是：Writing Runge-Kutta ODE solver using gsl，现在由于Daniel Fischer而修复了错误）