Thanks a lot, this function does a great job! The sptrisolve function is really much faster than the \. I observed a speedup factor of 7 for my ILU-preconditioned BiCG method and a factor of 4.7 in case of ILU-preconditioned BiCGSTAB. Now preconditioning not only reduces the number of iterates, but also computing time.