Problem

Follow up on#21309.

The problem iscontour on discontinuous functions. Mpl makes the (sane) assumption that, if a quad's corners have the values 0.0 and 0.78, the level contour of 0.51 (for example) has to pass somewhere through the quad, and draws it. This is incorrect where the function isactually discontinuous, but of course mpl has no way of knowing that.

Seethis stackoverflow post orhere.

Proposed solution

A simple solution for this would be to allow for a user-supplied float parameter,max_jump ordiscontinuity_threshold or something like that, defaulting tonp.inf. If within a quad the function values differ bymore than this threshold, a discontinuity is assumed and the level contour is assumednot to pass through the quad.

Additional context and prior art

No response