浮動小数点の扱いは面倒。jismeshでメッシュコードが意図しないものに

緯度経度から地域メッシュコードを求める

を使っていたけど，スカラーと配列で値が違っていた。

from jismesh.utils import to_meshcode

lat0 = 34.875
lon0 = 135.1875

print(to_meshcode(lat0, lon0, 3))
print(to_meshcode(np.array([lat0, lat0]), np.array([lon0, lon0]), 3))

を実行すると、

52352155
[52352154 52352154]

と、同じ緯度経度のはずなのにメッシュコードが違う。

$34.875$　が そのまま扱われず $34.87499999999999$ と解釈され、
切り捨てられたんだろうなと思いつつ、どうしてやろうかと。。。