<div dir="ltr">FYI, I just checked and MatSOR_*SELL does not use any vectorized instruction.<div>Why just not converting to SeqAIJ, factor and then use the AIJ implementation for MatSolve for the moment? </div></div><div class="gmail_extra"><br><div class="gmail_quote">2018-02-12 18:06 GMT+03:00 Stefano Zampini <span dir="ltr"><<a href="mailto:stefano.zampini@gmail.com" target="_blank">stefano.zampini@gmail.com</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><br><div class="gmail_extra"><br><div class="gmail_quote"><div><div class="h5">2018-02-12 17:36 GMT+03:00 Jed Brown <span dir="ltr"><<a href="mailto:jed@jedbrown.org" target="_blank">jed@jedbrown.org</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span>Karl Rupp <<a href="mailto:rupp@iue.tuwien.ac.at" target="_blank">rupp@iue.tuwien.ac.at</a>> writes:<br>
<br>
> Hi Stefano,<br>
><br>
>> Is there any plan to write code for native ILU/ICC etc for SeqSELL, at least to have BJACOBI in parallel?<br>
><br>
> (imho) ILU/ICC is a pain to do with SeqSELL. Point-Jacobi should be<br>
> possible, yes. SELL is really just tailored to MatMults and a pain for<br>
> anything that is not very similar to a MatMult...<br>
<br>
</span>There is already MatSOR_*SELL.  MatSolve_SeqSELL wouldn't be any harder.<br>
I think it would be acceptable to convert to SeqAIJ, factor, and convert<br>
the factors back to SELL.<br></blockquote><div><br></div></div></div><div>Yes, this was my idea. Today I have started coding something. I'll push the branch whenever I have anything working </div></div><span class="HOEnZb"><font color="#888888"><br><br clear="all"><div><br></div>-- <br><div class="m_2616199336360193854gmail_signature" data-smartmail="gmail_signature">Stefano</div>
</font></span></div></div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature">Stefano</div>
</div>