Some improvements that I want to do when I have time: - [ ] use higher-order finite elements - [ ] use higher-order time integrator