Work on examples. The following still does not work 100%:

DG -> minor bug, probably in values and sub_idx checking
Neutronics -> ??? (Milan)
Runge-Kutta for richards
N-S bearing, driven-cavity (circular-obstacle works somehow) -> probably Dirichlet
Schroedinger with no adaptivity (it does work for adaptivity).