This implements a discretized diffusion operator
The discretization decouples obtuse elements at the assembled diffusion matrix level and uses nodal values for D and C. This technique only conserves dose if the diffusivity is constant over a region.
Parameters: D C
D is a variable
C is a field or function