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>
