Is there any documentation or examples on how to simulate with this code? I also have a Firedrake acoustic wave propagation package for seismic inversions and perhaps we could merge some of this functionality together.