Constructing a PETSc MatMPIAIJ matrix from partitioned arrays

[Mads Hoel]

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