I'm trying to modify the Singular Poisson demo found at
http://fenicsproject.org/documentation/dolfin/dev/python/demo/documented/singular-poisson/python/documentation.html
so that I can use the AdaptiveLinearVariationalSolver, the sole reason being so that I can call solver.summary() and get an estimation of the error. I want everything else about the solving procedure to be the same. For example, I still want to assemble the matrix A and vector b and specify the null space of A. The main trouble I'm having is I don't know how to use the AdaptiveLinearVariationalSolver when I've already assembled A and b. Any help with this would be appreciated.