Problem with rectangular matrices

Jed Brown jed at
Tue May 26 14:29:47 CDT 2009

Matthew Knepley wrote:

> I am fine with this too, which is option 1 in my original
> mail. However, this means that when users preallocate, they must
> understand the column layout.

They are going to multiply with a vector, the layout of which must be
compatible with the column layout of the matrix.  Might be good to have
a nonsquare example in MatMPIAIJSetPreallocation or the manual.


-------------- next part --------------
A non-text attachment was scrubbed...
Name: signature.asc
Type: application/pgp-signature
Size: 260 bytes
Desc: OpenPGP digital signature
URL: <>

More information about the petsc-dev mailing list