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

bfast/bfast/monitor/utils.py

Line 477 in 4fd0c17

def map_indices(dates):

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.