core dump segmentation fault

[neuronphysics]

I have been trying to use this library to plot my mcmc results. I have 26 free parameters that have been measured by emcee but when I tried to plot them just before the script return fig I got the core dump error message (I tested where the code would crash). How could I solve the problem. Wouldn't it be more convenient and lighter if you use gs = gridspec.GridSpec(K, K) instead of fig, axes = pl.subplots(K, K, figsize=(dim, dim)) to make subplots?

Thanks in advance.