c - 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()