Python 不指定时间步长数组的galpy中的动态观察积分

Python 不指定时间步长数组的galpy中的动态观察积分,python,astropy,Python,Astropy,我刚刚发现了galpy,花了一整天的时间浏览了文档和示例,但我仍然对在不指定时间集的情况下集成轨道的可能性有点困惑 文档中的一个简单示例如下: from galpy import potential import numpy ts= numpy.linspace(0,100,10000) o= Orbit([1.,0.1,1.1,0.,0.1,0.]) o.integrate(ts,potential.MWPotential2014) 但是我想问,是否可以使用integrate方法,而不指定整

我刚刚发现了galpy,花了一整天的时间浏览了文档和示例,但我仍然对在不指定时间集的情况下集成轨道的可能性有点困惑

文档中的一个简单示例如下:

from galpy import potential
import numpy
ts= numpy.linspace(0,100,10000)
o= Orbit([1.,0.1,1.1,0.,0.1,0.])
o.integrate(ts,potential.MWPotential2014)
但是我想问,是否可以使用integrate方法,而不指定整个时间数组,而只指定以适当单位表示的时间量。例如:

t = 100
o.integrate(t, potential.MWPotential2014)

谢谢。

如果文档中没有说明存在这样的接口,您可以假设它不存在,尽管快速查看代码可以确认这一点。在任何情况下,即使存在这样一个快捷方式,积分器函数本身仍然需要一系列时间步来进行积分,因此这只是一个快捷方式

即使您只想指定“时间量”,您仍然需要提供一个时间步长,尽管您也可以有一个默认的时间步长。您可以编写自己的快捷方式函数,如:

def积分_轨道(轨道,pot,t,dt=0.01):
t=np.arange(0,t,dt)
返回轨道积分(t,pot)