We've been using Lanyon to discover (and prove correctness of) new state-of-the-art discontinuous Galerkin algorithms for solving PDEs.
Just select an order-of-accuracy and (optionally) a polynomial basis, and have Lanyon search the space of possible modal DG algorithms until it finds one that provably works for your problem. Then have it generate an end-to-end proof of correctness, and an optimized implementation, all from the same spec. All in a few minutes, all with negligible inference cost.
It turns out: automated algorithm discovery is a lot of fun!
Come and read more about it in our updated Research Notes, or by visiting our GitHub. 👇
We've been having a lot of fun using Lanyon to perform automated discovery of new discontinuous Galerkin (DG) algorithms for discretizing PDEs. We're starting with simple linear and nonlinear advection problems, with some more advanced (and exciting!) examples on the way.
The workflow is simple: have Lanyon autonomously explore the space of possible modal DG discretizations (of any order of accuracy), while autonomously checking its work by attempt to prove various correctness properties of the modal reconstruction polynomials, surface terms, and volume integrals, as well as entropy stability properties of the underlying fluxes. Once Lanyon finds a new discretization that works, it generates an optimized implementation in C and starts autonomously running numerical simulations with it. It turns out that you can very quickly discover all kinds of weird and wonderful (but provably correct!) numerical algorithms this way, many of which go beyond the current state of the art. And we're only just scratching the surface of what is possible with formally verified numerics.
Even for these simple cases, the optimized implementations (for p=1 and p=2, i.e. second-order and third-order accuracy) span around 30,000 lines of generated C code and around 40,000 lines of generated Lean proof. But Lanyon completes the whole algorithm discovery task in a few minutes, and with less than a dollar in total inference cost. We've updated our Research Notes with some of these new algorithmic results, but stay tuned for more discoveries still to come. GitHub links and Research Notes below. 👇