def test_getNewTimes_with_half_phase_two_day_bin(): """ Test to make sure times are shifted properly if the phase optional argument is used with a two day bin. """ times = np.random.uniform(0, 10, 75) newtimes = wm.getNewTimes(times, 2.) newtimes2 = wm.getNewTimes(times, 2., phase=0.5) assert np.round((np.min(newtimes2) - np.min(newtimes)), 7) == 1.000
def test_getNewTimes_one_day_bin_ten_days(): """ Test to make sure `getNewTimes` returns ten elements when 75 elements are entered over a ten day span with one day binning in place. """ times = np.random.uniform(0, 10, 75) newtimes = wm.getNewTimes(times, 1.) print(len(newtimes)) assert len(newtimes) == 10
def test_big_gaps_getNewVals(): """ Ensure getNewVals routine can handle big gaps in times """ timebin = 1. times = np.concatenate((np.random.uniform(0, 10, 50), np.random.uniform(30, 40, 50))) newtimes = wm.getNewTimes(times, timebin) rvs = np.random.normal(loc=0, scale=5, size=100) uncs = np.random.normal(loc=1., scale=0.5, size=100) newRVs, newUncs = wm.getNewVals(newtimes, times, rvs, uncs, timebin=timebin) fins = np.where(np.isfinite(newUncs)) newRVs = newRVs[fins] newUncs = newUncs[fins] newtimes = newtimes[fins] assert np.median(newUncs) < np.median(uncs)