我注意到这篇文章，推荐在 Matlab 中高效组装的位置spalloc和位置。sparse我个人对简单的情况使用稀疏组装。

然而，当涉及到耦合 PDE 的情况，比如 3-PDE 耦合时，标量未知数变成$3\times 3$张量。在这种情况下，我想不出一种利用的方法sparse()，不得不使用for-loops 来组装。

当最终组装的稀疏矩阵大到$30\text{k}\times 30\text{k}$（在 PC 上），组装过程变得非常慢（~10 分钟），而最终的矩阵求解步骤仍然很快（~不到 3 秒）。

有什么建议吗？任何通用解决方案都值得赞赏（不一定是 Matlab）