Problem with rectangular matrices
Matthew Knepley
knepley at gmail.com
Tue May 26 14:26:50 CDT 2009
On Tue, May 26, 2009 at 2:25 PM, Lisandro Dalcin <dalcinl at gmail.com> wrote:
> On Tue, May 26, 2009 at 4:18 PM, Matthew Knepley <knepley at gmail.com>
> 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. I guess
> this
> > is better
> > than changing the code unless I can think of a case where is plainly does
> > not apply.
> >
>
> An then, going back to my former observations, we should ensure that
> colmaps are properly initialized after preallocation, right
Yes, this is done in MPIAIJ, but might be missing elsewhere.
Matt
>
> --
> Lisandro Dalcín
> ---------------
> Centro Internacional de Métodos Computacionales en Ingeniería (CIMEC)
> Instituto de Desarrollo Tecnológico para la Industria Química (INTEC)
> Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET)
> PTLC - Güemes 3450, (3000) Santa Fe, Argentina
> Tel/Fax: +54-(0)342-451.1594
>
--
What most experimenters take for granted before they begin their experiments
is infinitely more interesting than any results to which their experiments
lead.
-- Norbert Wiener
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20090526/19ad6267/attachment.html>
More information about the petsc-dev
mailing list