This implements a discretized diffusion operator
The discretization is a consistent finite element discretization. It uses interior quadrature and evaluates D and C at the quadrature points within an element. This operator must be used for spatial adaptation.
Parameters: D C
D is either a field, function, or variable
C is a field or function