<div dir="ltr">This is the branch I'm currently working on <a href="https://bitbucket.org/petsc/petsc/branch/stefano_zampini/feature-factor-sell">https://bitbucket.org/petsc/petsc/branch/stefano_zampini/feature-factor-sell</a><div>So far, it provides support for MUMPS and SUPERLU</div></div><div class="gmail_extra"><br><div class="gmail_quote">2018-02-12 19:25 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">Barry,<div><br></div><div>for sure Amat,Pmat is the right approach; however, with complicated user codes, we are not always in control of having a different Jacobian matrix.</div><div>Since Mat*SELL does not currently support any preconditioning except PCSOR and PCJACOBI, we ask the user to put codes like </div><div><br></div><div>if (type is SELL)</div><div> create two matrices (and maybe modify the code in many other parts)</div><div>else</div><div>  ok with the previous code</div><div><br></div><div>Just my two cents.</div><div><br></div></div><div class="gmail_extra"><div><div class="h5"><br><div class="gmail_quote">2018-02-12 19:10 GMT+03:00 Smith, Barry F. <span dir="ltr"><<a href="mailto:bsmith@mcs.anl.gov" target="_blank">bsmith@mcs.anl.gov</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span><br>
<br>
> On Feb 12, 2018, at 9:59 AM, Stefano Zampini <<a href="mailto:stefano.zampini@gmail.com" target="_blank">stefano.zampini@gmail.com</a>> wrote:<br>
><br>
> FYI, I just checked and MatSOR_*SELL does not use any vectorized instruction.<br>
> Why just not converting to SeqAIJ, factor and then use the AIJ implementation for MatSolve for the moment?<br>
<br>
</span>  Why not use the mat, pmat feature of the solvers to pass in both matrices and have the solvers handle using two formats simultaneously instead of burdening the MatSELL code with tons of special code for automatically converting to AIJ for solvers etc?<br>
<div class="m_1622390754512268798HOEnZb"><div class="m_1622390754512268798h5"><br>
<br>
><br>
> 2018-02-12 18:06 GMT+03:00 Stefano Zampini <<a href="mailto:stefano.zampini@gmail.com" target="_blank">stefano.zampini@gmail.com</a>>:<br>
><br>
><br>
> 2018-02-12 17:36 GMT+03:00 Jed Brown <<a href="mailto:jed@jedbrown.org" target="_blank">jed@jedbrown.org</a>>:<br>
> 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>
> 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>
><br>
> Yes, this was my idea. Today I have started coding something. I'll push the branch whenever I have anything working<br>
><br>
><br>
><br>
> --<br>
> Stefano<br>
><br>
><br>
><br>
> --<br>
> Stefano<br>
<br>
</div></div></blockquote></div><br><br clear="all"><div><br></div></div></div><span class="HOEnZb"><font color="#888888">-- <br><div class="m_1622390754512268798gmail_signature" data-smartmail="gmail_signature">Stefano</div>
</font></span></div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature">Stefano</div>
</div>