Add 2nd order flux function
Put the diffusion functions in the filter format