Hand over the FE Basis type as a template parameter

This is much cleaner than handing over the Grid only, and then trying
to guess the basis using the preprocessor.  After all, the
RiemannianTrustRegionSolver calls the assembler, so it may know the
function space basis.

[[Imported from SVN: r10035]]