Templates in the diffusion code
Use templates in 4th order pressure solver