[petsc-dev] API changes in MatIS

[Stefano Zampini]

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.

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