在 python 中计算大量固定点的最快方法？

Question

我有大量的一维非线性不动点问题要解决，什么是最有效的数值求解器？我目前正在使用 scipy.optimize.fixed_point，我的 运行 1000 个任务大约需要 17 秒。感谢您的任何建议。

Answer 1

如果这些都是一维的，可以取fixed_point源码， https://github.com/scipy/scipy/blob/v1.5.2/scipy/optimize/minpack.py#L876 简化它（可以决定一次加速策略，不需要 _lazywhere 等）并用 cython 或 numba 编译它。