在Python 3中创建任意维稀疏矩阵

所以我使用Python 3创建一个表单的矩阵

L=[B 0 0
   0 B 0
   0 0 B]

哪里

B=[4 -1  0
  -1  4 -1
   0 -1  4]

B 我想做三次 N 次数(取决于输入)。我目前的拍摄如下

import numpy as np
import scipy.sparse as sp

one=np.ones(N)
four=4*np.ones(N)
data = np.array([-one, four, -one])
diags = np.array([-1,0, 1])
B=sp.spdiags(data, diags, N, N).toarray() # Create matrix B

L=np.kron(np.eye(N), B)

然而,当 N 是大的(这是必要的,因为这是为了解一个微分方程)。有没有更有效的方法?