MatGetSubMatrix question

[Matthew Knepley]

On Tue, May 26, 2009 at 12:11 PM, Lisandro Dalcin <dalcinl at gmail.com> wrote:

> On Mon, May 25, 2009 at 10:20 AM, Jed Brown <jed at 59a2.org> wrote:
> >
> > Is it really necessary that iscol be a sequential index set that covers
> > all columns, with csize selecting the local columns?  I realize this
> > might be natural/easy for AIJ formats, but it feels clumsy (not to
> > mention not scalable due to the ISAllGather).
> >
>
> Your are right, this is not scalable, but...
>
> > I'm guessing this is a relic from when IS was not a real parallel
> > object.  Am I missing some deep reason for this interface, or could
> > iscol be parallel and csize be dropped?
> >
>
> unless things have changed very recently (after 3.0.0 release), IS is
> far from being a real parallel object.


I agree. I am trying to think of a scenario is which you would want to
create
a distributed IS, which would then gather those columns to each process
without gathering the integers. Theoretically, we might want this, but is it
worth any time to implement?

  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/24575ec8/attachment.html>