This is an interesting take on how to solve the problem and am very keen to see how it does once you get it to work. Like you said, it’s a different angle because you are going to get the derivatives with the filter instead of using it to draw samples outright. I’d like to work on your code and play around with it - which version of PyMC are you using currently?