diff --git a/python/hccd/path_tracker.py b/python/hccd/path_tracker.py index 9635ff5..2fe61ad 100644 --- a/python/hccd/path_tracker.py +++ b/python/hccd/path_tracker.py @@ -3,10 +3,6 @@ import typing import scipy -def _sign(val): - return -1 * (val < 0) + 1 * (val >= 0) - - class PathTracker: """ Path trakcer for the homotopy continuation method. Uses a @@ -43,7 +39,6 @@ class PathTracker: raise RuntimeError("Newton corrector did not converge") - # TODO: Make sure the implementation with null_space instead of QR is correct def _perform_euler_predictor_step(self, y, step_size) -> typing.Tuple[np.ndarray, np.ndarray]: # Obtain y_prime @@ -52,12 +47,6 @@ class PathTracker: y_prime = ns[:, 0] * self._sigma - # Q, R = np.linalg.qr(np.transpose(DH), mode="complete") - # y_prime = Q[:, 2] - # - # if _sign(np.linalg.det(Q)*np.linalg.det(R[:2, :])) != _sign(self._sigma): - # y_prime = -y_prime - # Perform prediction y_hat = y + step_size*y_prime