MatGetArray() not supporting Mat type mpiaij.

Matthew Knepley knepley at gmail.com
Tue Jun 3 09:28:37 CDT 2008


On Tue, Jun 3, 2008 at 8:11 AM, Stephen R Ball <Stephen.R.Ball at awe.co.uk> wrote:
> Hi
>
> I have been trying to extract an array containing the local matrix
> values using MatGetArray() via the Fortran interface but get the error
> message that Mat type mpiaij is not supported with this routine. All I
> want to do is to extract the local matrix values so that I can output
> them to file in the format I want rather than via use of the MatView()
> routine. Can you suggest a way of how I can go about extracting the
> local matrix values?

This is no "local matrix". The Mat interface is supposed to be data structure
neutral so we can optimize for different architectures. If you want the values
directly, I would use MatGetRow() for each row.

  Matt

> Thanks
>
> Stephen
>
>



-- 
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which
their experiments lead.
-- Norbert Wiener




More information about the petsc-users mailing list