[petsc-dev] API changes in MatIS

Stefano Zampini stefano.zampini at gmail.com
Mon Jun 11 08:35:10 CDT 2012

I missed your answer...my e-mail client hid your last email.

>> Do you mean assembling the matrix on preselected vertices during
>> MatAssemblyBegin/End? Note that this will imply that standard
>> Neumann-Neumann methods will not work (they need the unassembled matrix to
>> solve for the local Schur complements).
> I'm not too concerned about that since I consider the classic N-N and
> original FETI methods to be rather special-purpose compared to the newer
> generation. I would like to limit the number of copies of a matrix to
> control peak memory usage.

Duplicate vertices don't concur so much to memory usage as faces and edge
nodes do, so why limit your matrix class? In principle the new class should
handle _every_ nonoverlapping method, either belonging to the old or the
new generation.

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20120611/857556e5/attachment.html>

More information about the petsc-dev mailing list