test_lyaps.py 文件源码

def tearDown(self):
        lyaps = []
        weights = []
        for T in np.arange(self.DDE.t, self.DDE.t+1000, 10):
            _, lyap, weight = self.DDE.integrate(T)
            lyaps.append(lyap)
            weights.append(weight)
        lyaps = np.vstack(lyaps)

        lyap_start = 40
        for i,lyap_control in enumerate(lyap_controls):
            lyap = np.average(lyaps[lyap_start:,i], weights=weights[lyap_start:])
            stderr = sem(lyaps[lyap_start:,i])
            print(lyap,stderr)
            self.assertAlmostEqual(lyap_control, lyap, delta=3*stderr)