`find_vertices()` undercounts for some float64 values of `lattice_spacing`
If lattice_spacing
is specified in float64 precision, some values (e.g. 244.91610425766507
) will produce a faulty count of vertices.
Note: Does not reproduce if spin_angle = 0.0
(or for e.g. SquareSpinIceClosed
).
from flatspin.model import PinwheelSpinIceDiamond
ls64 = np.float64(244.91610425766507) # 244.91610425766507
ls32 = np.float32(244.91610425766507) # 244.9161
asi64 = PinwheelSpinIceDiamond(lattice_spacing=ls64)
print(len(asi64.vertices())) # 1
asi32 = PinwheelSpinIceDiamond(lattice_spacing=ls32)
print(len(asi32.vertices())) # 25