问题描述
我要解决大量一维非线性不动点问题,最有效的数值求解器是什么?我目前正在使用scipy.optimize.fixed_point,大约需要17秒才能运行1000个任务。感谢您的任何建议。
解决方法
如果这些都是一维的,则可以使用fixed_point源, https://github.com/scipy/scipy/blob/v1.5.2/scipy/optimize/minpack.py#L876 简化它(可以一次决定加速策略,不需要_lazywhere等),然后用cython或numba进行编译。