PETSc中具有不规则形状的结构化网格（DM上下文）

Question

我在PETSc的结构化网格中存在一个有限的差异问题，DM上下文帮助我创建矩阵非常简单，DM给我一个矩形域的整个矩阵，对我好，因为我使用这个参考系统，但我的问题有一个不规则的形状，这意味着我没有在这个矩阵中使用很多条目。如何在矩阵汇编中分配此空间并告诉KSP上下文不使用此条目来解决问题？

我的代码是这样的：

DMDACreate2d()
DMCreateMatrix()
DMDAGetCorners()

! loop along the local matrix
DO i,1,width_y 
    DO i,1,width_y
        IF (IsInsideProblemDomian) THEN
            SetMatStencil()
            MatSetValuesStencil()
        END IF
    END DO
END DO

MatAssemblyBegin()
MatAssemblyBegin()