PR summary

The lookup table generated bySegmentedBivarColormap is slightly inaccurate. Most notably, the darkest color in the lookup table is not as dark as the darkest color in the input, and the same goes for the brightest color. This has been glossed over in the tests due to the liberal use ofatol.

>>>frommatplotlib.cmimportbivar_cmaps>>>cmap=bivar_cmaps['BiOrangeBlue']>>>print(cmap.patch[0,0, :])[0.0.0.]>>>print(cmap.patch[-1,0, :])[1.0.50. ]>>>print(cmap.patch[-1,-1, :])[1.1.1.]>>>print(cmap.lut[0,0, :])[0.002441410.002422240.002441411.        ]>>>print(cmap.lut[-1,0, :])[0.998535160.500484370.002441411.        ]>>>print(cmap.lut[-1,-1, :])[0.998535160.998546750.998535161.        ]

SegmentedBivarColormap uses the image-plotting resampler (Agg) to perform bilinear interpolation, and these inaccuracies result from a bug in the transform definition plus a second bug that would be fixed by#30184. However, just fixing these two bugs is not a good solution. The Agg resampler – which is intended for producing good-enough images for visual display – is problematic because it divides each pixel into 256 subpixels, which means it internally calculates 256 * M - 1 intermediate color candidates between the end colors, but the color table wants 254 intermediate colors between the end colors (to get 256 colors in total). That means fixing the aforementioned two bugs results in a non-smooth color table because 254 does not divide into 256 * M - 1 cleanly.

So, this PR instead takes the approach of eschewing the Agg resampler entirely and performing the bilinear interpolation using straight NumPy. Now the extreme colors are exactly what they should be.I'll note that the tests still need to use a bit ofatol because the definition ofBiOrangeBlue uses input that is rounded to 3 decimal places. Update: now fixed.

>>>frommatplotlib.cmimportbivar_cmaps>>>cmap=bivar_cmaps['BiOrangeBlue']>>>print(cmap.patch[0,0, :])[0.0.0.]>>>print(cmap.patch[-1,0, :])[1.0.50. ]>>>print(cmap.patch[-1,-1, :])[1.1.1.]>>>print(cmap.lut[0,0, :])[0.0.0.1.]>>>print(cmap.lut[-1,0, :])[1.0.50.1. ]>>>print(cmap.lut[-1,-1, :])[1.1.1.1.]

PR checklist

• [] "closes #0000" is in the body of the PR description tolink the related issue
• new and changed code istested
• []Plotting related features are demonstrated in anexample
• []New Features andAPI Changes are noted with adirective and release note
• [] Documentation complies withgeneral anddocstring guidelines