FYI, I just checked and MatSOR_*SELL does not use any vectorized
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 <stefano.zamp...@gmail.com>:
> 2018-02-12 17:36 GMT+03:00 Jed Brown <j...@jedbrown.org>:
>> 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.
> Yes, this was my idea. Today I have started coding something. I'll push
> the branch whenever I have anything working