phase-space modeling

I spent all weekend madly coding up my likelihood or Bayesian moment modeling and full phase-space modeling for one-dimensional dynamical systems. I can't see any easy, simple ways to make the code fast so it is super-slow. I have a method for inferring the gravitational potential using the position distribution and second moments of the velocity distribution and then I have one using the full position-and-velocity two-dimensional distribution function. Although there is a lot of computation involved, the problem is technically easy because all one-dimensional potentials are integrable. That makes it a nice toy system for learning about inference with different noise models and different selection functions.