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