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>
&gt; Barry-<br>
&gt;   Yes, each block is a rectangular portion  of the domain. Not so small<br>
&gt; though (more like 100 x 100 nodes)<br>
&gt; <br>
&gt;   Yaron<br>
&gt; <br>
&gt; <br>
&gt;            -------Original Message-------<br>
&gt;   From: Barry Smith<br>
&gt;   Subject: Re: Non-uniform 2D mesh questions<br>
&gt;   Sent: 29 Jan '07 19:40<br>
&gt; <br>
&gt; <br>
&gt;    Yaron,<br>
&gt; <br>
&gt;    Is each one of these "blocks" a small rectangular part of the<br>
&gt;   domain (like a 4 by 5 set of nodes)? I don't understand what you<br>
&gt;   want to do.<br>
&gt; <br>
&gt;    Barry<br>
&gt; <br>
&gt; <br>
&gt;   On Mon, 29 Jan 2007, [LINK:<br>
&gt; <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>
&gt; <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>
&gt; <br>
&gt;   &gt; Hi all<br>
&gt;   &gt; I have a laplace-type problem that's physically built from repeating<br>
&gt;   &gt; instances of the same block.<br>
&gt;   &gt; I'm creaing matrices for the individual blocks, and I'd like to reuse<br>
&gt;   &gt; the individual block matrices in order to compose the complete<br>
&gt; problem.<br>
&gt;   &gt; (i.e if there 10K instances of 20 blocks, I'd like to build 20<br>
&gt; matrices,<br>
&gt;   &gt; then use them to compose the large complete matrix)<br>
&gt;   &gt; Is a 2D DA the right object to do that? And if so, where can I find a<br>
&gt;   &gt; small example of building the DA object in parallel, then using the<br>
&gt;   &gt; different (for every instance) mappings of local nodes to global nodes<br>
&gt; in<br>
&gt;   &gt; order to build the complete matrix?<br>
&gt;   &gt;<br>
&gt;   &gt;<br>
&gt;   &gt; Thanks<br>
&gt;   &gt; Yaron<br>
&gt;   &gt;<br>
&gt; <br>
<br>
</bsmith@mcs.anl.gov></td></tr></tbody></table>