python - 无需 QuTip 工具箱即可重新实现 QuTip 功能

Question

有没有办法实现这个功能，它是 QuTip 函数（Python 中的 Quantum Toolbox）的一部分，而不需要安装 QuTip？

def destroy(N, offset=0):
    if not isinstance(N, (int, np.integer)):  # raise error if N not integer
        raise ValueError("Hilbert space dimension must be integer value")
    data = np.sqrt(np.arange(offset+1, N+offset, dtype=complex))
    ind = np.arange(1,N, dtype=np.int32)
    ptr = np.arange(N+1, dtype=np.int32)
    ptr[-1] = N-1
    return Qobj(fast_csr_matrix((data,ind,ptr),shape=(N,N)), isherm=False)

DESTROY 返回 N 维希尔伯特空间的湮没算子 out = destroy(N), N 是整数值 & N>0