<div dir="ltr"><div>On Sun, Feb 8, 2015 at 9:09 PM, Jed Brown <span dir="ltr"><<a href="mailto:jed@jedbrown.org" target="_blank">jed@jedbrown.org</a>></span> wrote:<br></div><div class="gmail_extra"><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex"><span class="">Matthew Knepley <<a href="mailto:knepley@gmail.com">knepley@gmail.com</a>> writes:<br>
<br>
> On Sun, Feb 8, 2015 at 7:20 PM, Ronal Celaya <<a href="mailto:ronalcelayavzla@gmail.com">ronalcelayavzla@gmail.com</a>><br>
> wrote:<br>
><br>
>> I know that. I want to have all the vector x replicated in all processes<br>
>> and update it in each iteration, so I don't need to communicate the vector<br>
>> x each time MatMult() is called.<br>
>> I'm not sure I'm making myself clear, sorry<br>
>><br>
><br>
> 1) This is not a scalable strategy.<br>
<br>
</span>Also note that you have to communicate many layers of overlap of A to<br>
run CG without neighbor communication.  The overhead is especially large<br>
if you have small subdomains (the case where communication latency is<br>
more important than bandwidth).<br></blockquote><div> </div><div> I need to explain this to my partners. Thank you</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">
<div class=""><div class="h5"><br>
> 2) If you know A and all the updates to x locally, why don't you just<br>
> compute y directly?<br></div></div></blockquote><div>  This is exactly what I'm doing in C implementation </div></div><br><div>Many thanks for your replies. They are very useful.</div><div>I am doing some tests with CG, comparing C and PETSc implementations and I need to explain the communication used in MatMult().</div><div><br></div><div>Once again, thanks for your help.</div><div><br></div><div>Best regards,</div><div><br></div>-- <br><div class="gmail_signature"><div dir="ltr">Ronal Celaya<br></div></div>
</div></div>