Problem with rectangular matrices

Jed Brown jed at 59A2.org
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.

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/20090526/e40fc310/attachment.sig>


More information about the petsc-dev mailing list