This implements a discretized diffusion operator
The discretization decouples obtuse elements at the element level and uses nodal values for D and C.
Parameters: D C
D is either a field, function, or variable
C is a field or function
