> On Feb 12, 2018, at 8:36 AM, Jed Brown <j...@jedbrown.org> wrote: > > Karl Rupp <r...@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.
This is feasible. I will give it a try when I have time. Another option that works now is to use two matrices with different Mat types, e.g. aij for the preconditioner and sell for Jacobian. Hong (Mr.)