ILU(n) works (as far as I could test it). Unfortunately only

n=0 and n=1 are practically useful. n=2 is already a direct
solver. (May I have still something wrong ...)

[[Imported from SVN: r76]]