I am trying to optimize a function that calculates a model that I am fitting with emcee. I profiled it and found the part that was taking the most time, and did a fairly radical re-engineering of that section which makes that function run over twice as fast in my standalone tests.