到目前为止,我已经使用GSL例程来计算特征值/向量并用C ++求解常微分方程(ODEs)。
现在我必须使用GNU多精度库(GMP),但GSL不支持其他数据类型而不是double。
有没有人知道计算特征值和解决与GMP数据类型(mpf_class)兼容的ODE的好方法?
目前我正在考虑使用Numerical Recipes中提供的例程。有谁知道这对我的任务是否是一个好主意?
祝你好运 多米尼克
答案 0 :(得分:1)
虽然显然不是您所希望的答案,但我认为您可以通过修改GSL库来使其与GMP库兼容。两者都是开源的,所以至少对于许可证来说没有问题,顺便说一句,我想很多人会非常感激。由于问题仅在于数据类型,因此扫描GSL库的所有代码并逐个替换例程就足够了。这不会是一项快速的工作,但这将是一项简单的工作。