On Fri, Feb 10, 2012 at 00:05, Dmitry Karpeev <karpeev at mcs.anl.gov> wrote:
> I don't think so. Preallocating doesn't set any values and seems to leave > the matrix marked !assembled. > MatDuplicate for such a matrix will fail. Assemblying it before setting > values (just to force an assembled flag) > will squeeze out the extra values, won't it? I think it would just be > reasonable to allow to duplicate unassembled > matrices, or, better yet, have a matrix be "assembled" by default until > MatSetValues has been called. > But I'm not sure whether either solution will break something else. > Are you calling DMSetMatrixPreallocateOnly()? -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120210/ebb80ce8/attachment.html>
