Kernel experiments with solution + gradient
Rewrite cdefastdg.hh
using sum factorized kernel for solution + gradient and do some polishing of the code.
Matrix free and matrix based methods are directed to one base method since they differ only by constants such as the right-hand side and boundary conditions.