我正在使用 petsc 解决由有限元离散的线性弹性问题。初始网格由网格文件读取，每个处理器中的分布使用 METIS 完成。我仅使用 petsc 的基本例程来创建子矩阵和向量。解决后我希望每个处理器都在文件中编写解决方案。我的问题是，由于 petsc 在内部重新排序自由度（例如处理器 0 从 0-N 等）我如何从初始网格跟踪哪个自由度对应于从 petsc 解向量获得的解值？例如，处理器 1 具有来自初始网格的 {4,6,8,5,13,​​9} 节点，并且 petsc 解向量的全局范围从 5 到 11（使用 VecGetOwnershipRange） .如何知道为节点 5 分配了哪个值？