[petsc-users] Recursive field split with MatNest

Jed Brown jed at jedbrown.org
Mon Jan 12 11:15:00 CST 2015

Colin McAuliffe <cjm2176 at columbia.edu> writes:

> I'd be interested in implementing this, I think I should have some time to
> spend on it in a few weeks. I suppose the ISSubset function should
> determine which ISs were used to make up a third IS? Can this be done by
> having ISSum attach a list of pointers to the constituent ISs to the summed
> IS for retrieval by ISSubset?

It could and that would probably be fairly harmless in terms of peak
memory usage.  An alternative would be to compute a cryptographic hash
of the index sets and only store that.  This would allow comparison
without needing to flatten the index sets to arrays of integers (see
current implementation of ISEqual).
-------------- next part --------------
A non-text attachment was scrubbed...
Name: signature.asc
Type: application/pgp-signature
Size: 818 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20150112/ab55b595/attachment.pgp>

More information about the petsc-users mailing list