This implements a discretized diffusion operator
Parameters: D C
D is either a variable, field, or function
C is a field or function