IndexError when datetime is used with hour != 0

[oceanites]

Whenever one of the dates fed to model.fit() has the hour or minute set to something else then 0, the error IndexError: boolean index did not match indexed array along dimension 1; dimension is 151 but corresponding boolean dimension is 152 is thrown.

This works also for the example code, when one value is modified like:

i = 34
di = dates_example[i]
dates_example[i] = datetime(di.year, di.month, di.day, 0, 34)

The error occurs in map_indices https://github.com/diku-dk/bfast/blob/4fd0c1727ec6ac309966d0ecfeef67e3ef02776b/bfast/monitor/utils.py#L477 where an empty array with shape (0, ) is returned, because ~np.isnan(ts).to_numpy() is False for every element.

A fix workaround for me is to use date instead of datetime objects. I can also try to come up with a better fix, but that would need some time.