enabled DUNE-ISTL::AMG and DUNE-ISTL::CGSolver in CPRPreconditioner.

The default behavior is the same as before, i.e. DUNE-ISTL::ILU0 and DUNE-ISTL::BiCGSolver are used 
for solving the elliptic problem. 

For some setting the AMG/CG combination can really speed up the solver, for others it doesn't. I still need to investigate more on the solver part. Maybe Markus can have  a look on the AMG args.
I was hoping that keeping the right hand side from the previous iteration would speed up the solver, but in all tested cases it didn't matter. So ve = 0 is the better choice, because the 0 is always in the Krylov subspace.