geodynamics/aspect — pull request #132 by tjhei (2014-06-11 15:21:20+00:00)
Title: fix direct solver with nonlinear schemes

We need to compute and return the correct initial_residual in solve_stokes() when using a direct solver for the nonlinear methods like iterated IMPES to work.