Hello, this is just to say that pardiso works very nicely and is very fast ... but it is a DIRECT SOLVER!!
ILU is for solving really large systems, which can not be solved directly. I do understand that it is extremely difficult to write an efficient ILU preconditioner, that's exactly the reason for which i orginally started this thread
I want to use iterative methods (Bi-CG Stab, TFQMR and other) for solving dense linear system. I'm doing at the beginning prefiltration dense matrix to obtain a sparse preconditioner, and then use ILU0 or ILUT algorithm for decomposition.
Since the last post on this topic has been a couple of years ago, I'd like to ask, what is the current state of parallelism for ilu0 preconditioners? Is such a feature on the horizon (or implemented already?).