On Wed, Apr 20, 2011 at 03:50, Debao Shao <Debao.Shao at brion.com> wrote:
> ierr = MatZeroEntries( M ); assert( ierr == 0); > ierr = MatDiagonalSet( M, vec.pv, INSERT_VALUES ); > ierr = MatCopy( ms->M, mStorage->M, DIFFERENT_NONZERO_PATTERN ); > Is M somehow related to ms->M or mStorage->M? What do you actually want to do? -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110420/4eb2391b/attachment-0001.htm>
