FYI, I just checked and MatSOR_*SELL does not use any vectorized instruction. Why just not converting to SeqAIJ, factor and then use the AIJ implementation for MatSolve for the moment?
2018-02-12 18:06 GMT+03:00 Stefano Zampini <[email protected]>: > > > 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 > -- Stefano
