|
6 | 6 |
|
7 | 7 | import numpy as np |
8 | 8 | import symengine |
9 | | -from numpy import hstack, log |
10 | 9 |
|
11 | 10 | from jitcxde_common import checker, jitcxde |
12 | 11 | from jitcxde_common.helpers import find_dependent_helpers, sort_helpers, sympify_helpers |
@@ -729,7 +728,7 @@ def set_initial_value(self, y, time=0.0): |
729 | 728 | for _ in range(self._n_lyap): |
730 | 729 | new_y.append(random_direction(self.n_basic)) |
731 | 730 |
|
732 | | - super().set_initial_value(hstack(new_y), time) |
| 731 | + super().set_initial_value(np.hstack(new_y), time) |
733 | 732 |
|
734 | 733 | def set_integrator(self,name,interpolate=None,**kwargs): |
735 | 734 | """ |
@@ -769,7 +768,7 @@ def integrate(self, *args, **kwargs): |
769 | 768 | super().integrate(*args, **kwargs) |
770 | 769 | delta_t = self.t-old_t |
771 | 770 | norms, tangent_vectors = self.norms() |
772 | | - lyaps = log(norms) / delta_t |
| 771 | + lyaps = np.log(norms) / delta_t |
773 | 772 | self._y[self.n_basic:] = tangent_vectors.flatten() |
774 | 773 | super().set_initial_value(self._y, self.t) |
775 | 774 |
|
@@ -910,7 +909,7 @@ def integrate(self, *args, **kwargs): |
910 | 909 | super().integrate(*args, **kwargs) |
911 | 910 | delta_t = self.t-old_t |
912 | 911 | norm = self.norm() |
913 | | - lyap = log(norm) / delta_t |
| 912 | + lyap = np.log(norm) / delta_t |
914 | 913 | super().set_initial_value(self._y, self.t) |
915 | 914 |
|
916 | 915 | return self._y[self.main_indices], lyap |
|
0 commit comments