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)