Thanks for the paper, I will read it. There are some SMC-ABC methods specially tailored to ODEs. It seems these methods are capable of generating better proposal distributions. I have not read those papers, but it seems the implementation in PyMC3 could accommodate those methods.