Add function to compute residual and to test linear convergence
fixed merge auto-conflict linear implicit system