Numpy Toeplitz矩阵

Question

我尝试使用scipy.linalg.toeplitz(c, r=None)在Python语言中构建toeplitz矩阵。虽然我成功了，但我不能保持Fortran的顺序。我需要确保正在构造的toeplitz数组保持Fortran顺序，因为我调用的是BLAS函数。

有什么办法可以做到这一点吗？

Answer 1

在Fortran排序中制作副本是最简单的方法。toeplitz会创建一个新的数组，但不会让您控制排序。

例如。

x = scipy.linalg.toeplitz([1, 2, 3, 4])
x = np.asfortranarray(x)

如果您想节省内存，也可以就地执行此操作。例如。

x[:] = x.T
x = x.T