Hi all, I have a partitioned matrix on CRS format across several processors. I have 4 arrays describing the parallel matrix; the typical CRS arrays irow, jcol and values, and an array describing the mapping of local nodes/rows to global nodes/rows. jcol is in local indices.
1. Does Petsc support overlap? 2. What is the optimal way of creating a PETSc MatMPIAIJ from these arrays? I looked at MatInsertValues, but would not that create overhead in re-partitioning/distributing the arrays? Thanks in advance. -- Mads Hoel
