[petsc-dev] plans for preconditioners for SeqSELL

[Stefano Zampini]

2018-02-12 17:36 GMT+03:00 Jed Brown <jed at jedbrown.org>:

> Karl Rupp <rupp at iue.tuwien.ac.at> writes:
>
> > Hi Stefano,
> >
> >> Is there any plan to write code for native ILU/ICC etc for SeqSELL, at
> least to have BJACOBI in parallel?
> >
> > (imho) ILU/ICC is a pain to do with SeqSELL. Point-Jacobi should be
> > possible, yes. SELL is really just tailored to MatMults and a pain for
> > anything that is not very similar to a MatMult...
>
> There is already MatSOR_*SELL.  MatSolve_SeqSELL wouldn't be any harder.
> I think it would be acceptable to convert to SeqAIJ, factor, and convert
> the factors back to SELL.
>

Yes, this was my idea. Today I have started coding something. I'll push the
branch whenever I have anything working



-- 
Stefano
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20180212/4c33a299/attachment.html>