Non-uniform 2D mesh questions

[Berend van Wachem]

Hi,

I guess my question might be superfluous - but I've read a number of 
contradictory things.

You would like multiple DA vectors to share one matrix? If so, this is 
exactly what I am doing and if you want, I can send you parts of code 
how to do so.

Berend.

> Barry-
> So far I only thought of having a single large sparse matrix.
>
> Yaron
>
>
> 	-------Original Message-------
>   From: Barry Smith
>   Subject: Re: Non-uniform 2D mesh questions
>   Sent: 30 Jan '07 10:58
>
>
>    Yaron,
>
>    Do you want to end up generating a single large sparse matrix?
> Like a MPIAIJ
>   matrix? Or do you want to somehow not store the entire huge matrix
> but still
>   be able to solve with the composed matrix? Or both?
>
>    Barry
>
>
>   On Mon, 29 Jan 2007, [LINK:
> http://webmail.oak-research.com/compose.php?to=yaron@oak-research.co
>m]
>
> yaron at oak-research.com wrote:
>   > Barry-
>   > Yes, each block is a rectangular portion of the domain. Not so
>   > small though (more like 100 x 100 nodes)
>   >
>   > Yaron
>   >
>   >
>   > 	 -------Original Message-------
>   > From: Barry Smith
>   > Subject: Re: Non-uniform 2D mesh questions
>   > Sent: 29 Jan '07 19:40
>   >
>   >
>   > Yaron,
>   >
>   > Is each one of these "blocks" a small rectangular part of the
>   > domain (like a 4 by 5 set of nodes)? I don't understand what you
>   > want to do.
>   >
>   > Barry
>   >
>   >
>   > On Mon, 29 Jan 2007, [LINK:
>   > [LINK:
>
> http://webmail.oak-research.com/compose.php?to=yaron@oak-research.co
>m]
> http://webmail.oak-research.com/compose.php?to=yaron@oak-research.co
>m]
>
>   > [LINK:
>
> http://webmail.oak-research.com/compose.php?to=yaron@oak-research.co
>m]
>
> yaron at oak-research.com wrote:
>   > > Hi all
>   > > I have a laplace-type problem that's physically built from
>   > > repeating instances of the same block.
>   > > I'm creaing matrices for the individual blocks, and I'd like
>   > > to
>
> reuse
>
>   > > the individual block matrices in order to compose the complete
>   >
>   > problem.
>   >
>   > > (i.e if there 10K instances of 20 blocks, I'd like to build 20
>   >
>   > matrices,
>   >
>   > > then use them to compose the large complete matrix)
>   > > Is a 2D DA the right object to do that? And if so, where can I
>   > > find
>
> a
>
>   > > small example of building the DA object in parallel, then
>   > > using the different (for every instance) mappings of local
>   > > nodes to global
>
> nodes
>
>   > in
>   >
>   > > order to build the complete matrix?
>   > >
>   > >
>   > > Thanks
>   > > Yaron