Numerical precision issue in rotate_vector

Personne n'a encore pris cette issue.

Hi,

I ran into a numerical precision issue when using the rotate_vector function with 32bit floats as input data. In this case, the lat/lon increments used to determine the relative rotation between map projection and the lat/lon grid is chosen too small. The result was a cristalline structure in the plotted wind barbs which patches of barbs pointing in the exact same direction.

https://github.com/matplotlib/basemap/blob/decfa95124dab76499734145d03b002b9db27477/lib/mpl_toolkits/basemap/__init__.py#L3118

Converting the input to 64bit floats solves the issue for me, as did an increase of the increment to 1.0e-4.

Best regards,
Clemens

Commencez dans lib/mpl_toolkits/basemap/__init__.py, autour de l’implémentation de rotate_vector référencée dans l’issue. Reproduisez l’artefact des barbules de vent avec des entrées à virgule flottante sur 32 bits, comparez-le avec des entrées sur 64 bits et avec l’incrément plus grand mentionné, et vérifiez que le motif cristallin est éliminé sans modifier le comportement de rotation attendu.

Rédigé par le modèle d'indexation à partir du texte de l'issue.