Lisandro Dalcin wrote: > If you do MatSetType(A, MATAIJ), you end-up having a MATSEQAIJ or a > MPIMPIAIJ... Vec does not have such "aliasing" on type names ...
How is this a problem? You aren't matching type names, you're just setting a function to be called after the sizes are known. Since this happens in MatCreate_{Seq,MPI}AIJ, it doesn't matter whether MatSetType was called with MATAIJ. Also note that we currently have mat->ops->setsizes (only implemented by SeqDense) so we should either use it or remove it, as long as we value one-way-to-do-it. Jed -------------- next part -------------- A non-text attachment was scrubbed... Name: signature.asc Type: application/pgp-signature Size: 260 bytes Desc: OpenPGP digital signature URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20090722/9135cfd4/attachment.pgp>