<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Tue, Aug 4, 2015 at 12:20 PM, Jed Brown <span dir="ltr"><<a href="mailto:jed@jedbrown.org" target="_blank">jed@jedbrown.org</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Cong Li <<a href="mailto:solvercorleone@gmail.com">solvercorleone@gmail.com</a>> writes:<br>
<br>
> Thanks very much for your suggestions.<br>
><br>
> Actually I am also considering using communication-avoiding matrix power<br>
> kernel (CA-MPK) to do SPMM. However, the communication pattern of CA-MPK<br>
> depends on the sparsity pattern. So the implementation could be very<br>
> complex for some of problems.The efficient implementation of CA-MPK is<br>
> actually one of problems I want to solve during my PhD course.<br>
<br>
You can get the pattern with MatIncreaseOverlap.<br>
<br>
Lack of useful preconditioning all but destroys the practical utility of<br>
these matrix powers kernels for solvers, even if the surface area/volume<br>
ratio is favorable (an odd corner of the relevant problem space).  I<br>
consider it a fashion that has already gotten more attention than it<br>
deserves and should be on its way out.  I would recommend a different<br>
thesis topic if you want to do work with a tangible impact on practical<br>
computational science and engineering.  If you're really excited about<br>
this niche, by all means, have fun.<br>
</blockquote></div><br>Totally true if you are doing this for solvers, particularly elliptic solvers. If you are</div><div class="gmail_extra">planning to use it for something like Tall-Skinny QR, or maybe some graph problems</div><div class="gmail_extra">it could make more sense.</div><div class="gmail_extra"><br></div><div class="gmail_extra">  Matt<br clear="all"><div><br></div>-- <br><div class="gmail_signature">What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div>
</div></div>