this is wrong, it needs to be e - ((shape[0] - (npy_uintp)d) % sizeof(npy_uint64)) % bits128)
+ tests which cover these cases, see testing.utils.gen_alignment_data used for the SIMD code