I am using the MethodOfLines package for solving a system of first order pde’s. Now I want to introduce callback functions to handle discrete events. It is unclear to me how to call the integrator, because the pde system has been changed by discretization into an ODE system. How are the fields of the ODE system obtained by discretization called?
ODESystem had this capability, but I’m not sure PDRSystem does yet. Open an issue.