在 MATLAB 中,使用如下命令创建五对角稀疏矩阵非常方便:
I = eye(m); % create identity matrix
e = ones(m,1); % create an array of all 1's
T = spdiags([e -4*e e],[-1 0 1],m,m);
S = spdiags([e e],[-1 1],m,m);
A = (kron(I,T) + kron(S,I))/hˆ2;
我想知道在 c/c++ 中是否有任何巧妙的技巧可以做同样的事情。