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