[examples.linearelliptic] add preliminary thermalblock demo

Warning: Does not work with inhomogeneous Dirichlet boundary
conditions. Also, the error estimator does not work as there
is currently no way to obtain the vector belonging to a linear
functional.