def test_warp_q_gamma(self):
M = 101
q1 = np.sin(np.linspace(0, 2 * np.pi, M))
timet = np.linspace(0, 1, M)
gam = fs.optimum_reparam(q1, timet, q1)
q1a = fs.warp_q_gamma(timet, q1, gam)
self.assertAlmostEqual(sum(q1 - q1a), 0)
def test_reparm(self):
M = 101
q1 = np.sin(np.linspace(0, 2 * np.pi, M))
timet = np.linspace(0, 1, M)
gam = fs.optimum_reparam(q1, timet, q1)
self.assertAlmostEqual(sum(gam - timet), 0)
