科学Python cod索引器

Question

我一直在做一些代码，集成，一些操作，然后更多的集成。 这是code（谢谢@JRichardSnape！）。 基本上，这段代码解决了一个矩阵方程，这就是mesolve所做的。它需要一个哈密顿量（一个关键的物理矩阵），一个初始密度矩阵，和一个时间列表来计算rho（t）at（我们要求解的）和折叠算子。然后我提取结果并乘以另外两个数组并绘制它。你知道吗

我使用qutip量子力学模块，因为他们有我需要的解算器：mesolve。qutip模块要求将所有矩阵转换成一个量子对象，这是由Qobj（x）完成的。你知道吗

我在collapse操作符定义之后添加了一行定义rho0=L1。发生这种情况时，它会给我一个索引错误：

IndexError: index 0 is out of bounds for axis 0 with size 0

我已经试着找出问题所在。它在定义f\u t时出错，问题是它试图索引数组（n.data是一个复数128的单元素数组）。怎么回事？你知道吗

Answer 1

这再现了您的错误：

In [34]: data = np.zeros((0,10))    
In [35]: data
Out[35]: array([], shape=(0, 10), dtype=float64)

In [36]: data[0]
                                     -
IndexError                                Traceback (most recent call last)
<ipython-input-36-88cae4a5bf31> in <module>()
  > 1 data[0]

IndexError: index 0 is out of bounds for axis 0 with size 0

在不深入研究链接代码的情况下，很明显，所讨论的数组具有长度为0的第一维。x.shape[0]是0。你知道吗

你可能会得到这样一个数组，用一个空列表索引另一个数组

In [44]: data=np.ones((3,4))

In [45]: data[[],...]
Out[45]: array([], shape=(0, 4), dtype=float64)

由于你提供的信息有限，很难说得更具体。检查所有可疑阵列的形状。你知道吗

---

所以Qobj记录在http://qutip.org/docs/2.2.0/apidoc/classes.html

和http://qutip.org/docs/2.2.0/apidoc/functions.html#qutip.mesolve.mesolve中的mesolve

而rho0预计是rho0 : qutip.qobj。你知道吗

ground的底层数组是(7,1)形状，对于L1（问题是rho0？）是(7,7)和除[0,0]之外的所有0。你知道吗

看起来这是Integration not successful in Python QuTiP的衍生产品