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 <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
>
>
>
> --
> Stefano
>



-- 
Stefano

Reply via email to