Implement damping using the directional restriction.

Now the "coarse correction" solver is an inexact
CG with ILU0 preconditioner. We only use a small number
of CG steps such that CGs rate will dominate the nonlinear
iteration.