Numpy Toeplitz矩陣

Question

我嘗試在Python中使用scipy.linalg.toeplitz(c, r=None)構建toeplitz矩陣。雖然我成功了，但我無法維持Fortran訂購。我需要確保構造的toeplitz數組保持Fortran排序，因爲我正在調用BLAS函數。

有沒有一種方法可以做到這一點？

Answer 1

在Fortran排序中創建副本最簡單。 toeplitz創建一個新的數組，但不能控制排序。

E.g.

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

如果您想節約內存，也可以在原地進行此操作。例如。

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