[petsc-users] custom sparse matrix

[Jed Brown]

2011/3/21 Barry Smith <bsmith at mcs.anl.gov>

>
> On Mar 21, 2011, at 10:14 AM, Alejandro Marcos Aragón wrote:
>
> > Hi Barry,
> >
> > Thanks for your answer. I'll look into those functions, but instead of
> copying, do you know if there is a way to give the address of the arrays?
>
> MatCreateSeqAIJWithArrays() will not copy the arrays but just use them. In
> parallel MatCreateMPIAIJWithSplitArrays() will not copy the arrays, just use
> them but you need to have the correct format for the arrays (in parallel we
> use a "non-standard" parallel compressed sparse row format that likely you
> don't use but check the manual page
> http://www.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-3.0.0/docs/manualpages/Mat/MatCreateMPIAIJWithSplitArrays.html


There is slightly more info in the more recent man page

http://www.mcs.anl.gov/petsc/petsc-as/snapshots/petsc-dev/docs/manualpages/Mat/MatCreateMPIAIJWithSplitArrays.html
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110321/76308273/attachment.htm>