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 variable, field, or function
C is a field or function
