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