I’m not sure classical methods can do any better here (long time spans), but we do keep trying to optimize this more and more. saveat at longer intervals probably helps, since I assume you don’t need every step in there?
BTW, I think for estimation, directly using the PDE formulation might be more efficient if you’re fitting via distributional information. That might be something to look into.