不同处理器行索引不连续时，PETSc应选用哪种并行稀疏矩阵类型？

MATMPIAIJ并非强制要求行索引连续，它本身支持任意的行所有权分布，只是默认采用连续行拆分的方式。针对你需要的非连续行分配场景（处理器0处理行0、2，处理器1处理行1、3），无需更换矩阵类型，只需手动指定每个进程负责的行索引即可。

具体可参考两种常见实现方式：

• 使用MatCreateMPIAIJWithArrays函数，直接传入当前进程负责的行索引数组、对应列索引数组及数值数组完成矩阵初始化。
• 先通过MatCreateMPI创建空的MATMPIAIJ矩阵，再调用MatSetValues逐个添加属于当前进程的行数据，确保仅操作自身负责的非连续行索引。

如果你的矩阵是块稀疏结构（比如每行对应固定大小的块），也可以考虑MATMPISBAIJ（块AIJ并行矩阵类型），它同样支持非连续行分布，适配块结构化的稀疏矩阵场景。

内容的提问来源于stack exchange，提问作者Pentaery