test_jitcode.py 文件源码

def tearDown(self):
        self.initialise_integrator()
        times = range(10,100000,10)
        data = np.vstack( self.ODE.integrate(time)[1] for time in times )
        result = np.average(data[1000:], axis=0)
        margin = standard_error(data[1000:], axis=0)
        print(data,result,margin)
        self.assertLess( np.max(margin), 0.003 )
        for i in range(self.n):
            self.assertLess( result[i]-lyaps[i], 3*margin[i] )