Implement LDG and interior penalty methods for 1d and 2d heat equation. Some of this is explained in the DG slides, see the references for some test cases.