ValueError ("Contour levels must be increasing") about corner.py HOT 6 CLOSED

Assuming you're running the example block under "Usage" in the README.rst, I can't reproduce this error. I'm on Ubuntu, using conda, python 2.7.11, numpy=1.10.2, matplotlib=1.5.0.

from corner.py.

Le 06/01/2016 20:04, Kyle Barbary a écrit :

Assuming you're running the example block under "Usage" in the
README.rst, I can't reproduce this error. I'm on Ubuntu, using conda,
python 2.7.11, numpy=1.10.2, matplotlib=1.5.0.

Well, thank you for trying. I am indeed doing:

import numpy as np
import corner

ndim, nsamples = 5, 10000
samples = np.random.randn(ndim * nsamples).reshape([nsamples, ndim])
figure = corner.corner(samples)
figure.savefig("corner.png")

to which I add that

corner.version
u'1.0.0'

np.version
'1.10.1'
Any other suggestions are welcome, David.

from corner.py.

Progress! Two things...

1. I execute to completion without errors, but with contour-less plots, by changing the line to
   corner.corner(samples, plot_contours=False)

2. Towards the end of corner.py, the contours array V is created:
   V = np.empty(len(levels))
   for i, v0 in enumerate(levels):
   try:
   V[i] = Hflat[sm <= v0][-1]
   except:
   V[i] = Hflat[0]

A print of "levels" shows that it is ascending. But V is descending!
I added a few lines:
print (" das corner VVVVV: ", V)
V = V[::-1]
print (" das corner VVVVV: ", V)
and here is a sample of what goes to my screen:
das corner VVVVV: [ 196. 140. 71. 29.]
das corner VVVVV: [ 29. 71. 140. 196.]

With this reverse of V,
corner.corner(samples, plot_contours=True)
now works nicely.

I will leave it to you experts to understand why V is okay on some systems but not on mine!
Thank you for listening, David.

from corner.py.

This error was added in matplotlib 1.5.1 because of a bug in all versions of contourf prior to 1.5.0. See #65 for the relevant issue. It's not quite as simple as just sorting the levels.

from corner.py.

This should be fixed in version 1.0.1 (see #70). Give it a shot.

from corner.py.

Yes! my problem is indeed resolved. Thanks again to all, David.

from corner.py.