2018-02-12 17:36 GMT+03:00 Jed Brown <[email protected]>: > Karl Rupp <[email protected]> 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
