Rb stuff needs a non-const reference to dof_indices.
Use accessor API for FEMContext.
Update solvers to use updated FEMContext API.
Need non-const versions of these accessors. I'm not convinced this
is the right way to go here, but rolling with it for now.
Updating reduced_basis to use new FEMContext API's.
Assert on the private data structures.
Update to use new FEMContext API. Also needed a
non-const Elem&.
Return a reference to the QBase instead of a pointer.
Now make everything protected/private. The library builds for me,
but the examples still need to be updated so likely make check
will fail at this point.
Updated adjoints_ex1 to use the new FEMContext API.
Update adjoint_ex2 for new FEMContext API.
Update adjoints_ex3 to use new FEMContext API.
Updated adjoint_ex4 to use updated FEMContext API.
Updated adjoints_ex5 to use updated FEMContext API.
Updated fem_system examples to use updated
FEMContext API.
Updated reduced_basis examples to conform to new
FEMContext API.
UPdated vector_fe examples to conform to new FEMContext API.
Patch from @dknez to fixed reduced_basis_ex2 for
FEMContext API updated.
Patch from @dknez to fix reduced_basis_ex7 for
FEMContext API changes.
Revert "Now make everything protected/private. The library builds for me,"

This reverts commit 248f1facdb574d5142081d8c69c74c4dc656370d.
Will reapply after the 0.9.2 release.