我正在尝试在 Julia 中制作对角块矩阵。我有一个nxn数组，我想将其P复制为对角线上的块矩阵，而矩阵的其余部分是稀疏的。
如果arr是我的阵列，我知道我可以使用：
blockdiag(sparse(arr),sparse(arr))
创建一个块，其中 P=2 对角线上的数组副本。
但是，对于大 P，我如何使用变量以一般方式执行此操作P？
我尝试制作一个 ( nxnxP)数组，但BlockDiag()不接受 3D 数组。