我有一个通过 ODE 制定的物理问题。现在我喜欢使用 Pythons scipy.integrate 和其中的 complex_ode 对其进行数值求解。我弄清楚了它是如何工作的，但现在我想优化我的代码，我遇到了一个问题：

我需要大量的网格点才能使解决方案收敛。所以我尝试使用步长控制（dop853）设置另一种积分器方法。我不太熟悉该功能，但它也能正常工作，而且比我的方法快得多。这是一个简单的 MWE，它不反映我复杂的 ODE，但设置的结构（复杂值等）是相同的：

from scipy import *
from scipy.integrate import *
from pylab import *

grid = linspace(0, 8, 10e3)
dgrid = (grid[1]-grid[0])*ones(100e3)

def RHS(t, x):
    return -x

y = zeros(len(grid), dtype=complex)
y[0] = 1.0

solver = complex_ode(RHS)
solver.set_initial_value(y[0], grid[0]).set_integrator('dop853')

for t in range(len(grid)):
    solver.integrate(solver.t+dgrid[t])
    y[t] = solver.y[0]

plot(grid, y.real)
show()

但现在的问题是：我现在当然想将时间网格与解决方案联系起来。什么电网现在连接到解决方案？它仍然是我引入的网格（网格）还是步长控件在不同的时间步长计算并且我需要另一个时间数组？