有没有办法在 Python 中使用带有 scipy.integrate 的十进制数据类型具有 28 个精度数字？

问题描述

我目前正在尝试使用命令 solve_ivp 在 Python 中使用 scipy.integrate 解决 ODE 系统。有没有办法使用十进制数据类型集成这个系统？我需要超过 16 位的精度。例如，我需要 28 位数字的输出。

解决方法

Solve_ivp 在内部将 Python 例程与 Numpy 数组一起使用（线性求解和调用 Fortran livrary 的 LSODA 方法除外）。 Numpy 仅限于双精度，因此您将无法获得比这更高的精度。

我相信 Julia 语言能够在其所有算法中原生使用四精度。