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