Just a quick comment: > 3/ One would need to compute the Newton rhs vector. This is simply the > residual -- should not be overly complicated. This would also be useful by itself and could be used/tested for the current iterated IMPES scheme. -- Timo Heister http://www.math.clemson.edu/~heister/