Plotting super slow with ODE model

I am using PyMC3 to fit an ODE model. Sampling is very slow but I understand that this is a limitation of the DifferentialEquation model in PyMC whichis not yet mature.

What I don’t understand is why plot_posterior() and other plotting functions are so slow when using the traces sampled from the model. Shouldn’t this be as fast as with any other trace? Is there any way I can improve this?