[petsc-users] Which matrix object type should I use?

[Klaus Burkart]

I am confused with matrix object types and which one to use in my case:


I want to solve a linear system using multiple processes. The steps I am struggling with is the import of the global matrix and global rhs vector into applicable PETSc objects.

The global matrix is currently stored in CSR format using 3 arrays. Array one contains all the row indices, array two all the column indices, array three all the values of the global matrix. In addition to that, there's a vector object containing the rhs vector.

I think for the vector object I should use VectorCreateMPI with m set to PETSC_DECIDE?


What's the matrix object I should use in my case which also takes care of the allocation local matrices to the processes at a later stage i.e. creates the local matrices in line with the rhs vector?
Klaus
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20180412/fe5d079b/attachment.html>