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 at 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 at oak-research.co >m] > http://webmail.oak-research.com/compose.php?to=yaron at oak-research.co >m] > > > [LINK: > > http://webmail.oak-research.com/compose.php?to=yaron at 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
