Sorting of times should be stable?
I think the sorting of observed times done by merge_series
should use a stable algorithm, using
ksort = np.argsort(full_t, kind='stable')
I didn't fully test the effects of not doing a stable sort, but can it lead to mixing up observations from different time series?