_jitcdde.py 文件源码

def integrate_blindly(self, target_time, step=None):
        """
        Like `jitcdde`’s `integrate_blindly`, except for orthonormalising the separation functions after each step and the output being analogous to `jitcdde_lyap`’s `integrate`.
        """

        dt,number,total_integration_time = self._prepare_blind_int(target_time, step)

        instantaneous_lyaps = []

        for _ in range(number):
            self.DDE.get_next_step(dt)
            self.DDE.accept_step()
            self.DDE.forget(self.max_delay)
            norms = self.DDE.orthonormalise(self._n_lyap, self.max_delay)
            instantaneous_lyaps.append(np.log(norms)/dt)

        lyaps = np.average(instantaneous_lyaps, axis=0)

        return self.DDE.get_current_state()[:self.n_basic], lyaps, total_integration_time