We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
There was an error while loading. Please reload this page.
2 parents d21a949 + 05c1993 commit 1bd8e1aCopy full SHA for 1bd8e1a
.gitignore
@@ -19,3 +19,7 @@ dist/**
19
20
# Virtual environment
21
venv*/**
22
+.venv
23
+
24
+# Pyenv
25
+.python-version
pydmd/bopdmd.py
@@ -624,8 +624,8 @@ def compute_B(alpha):
624
q_out, djac_out, j_pvt = qr(
625
djac_matrix, mode="economic", pivoting=True
626
)
627
- ij_pvt = np.arange(IA)
628
- ij_pvt = ij_pvt[j_pvt]
+ ij_pvt = np.zeros(IA).astype(int)
+ ij_pvt[j_pvt] = np.arange(IA).astype(int)
629
rjac[:IA] = np.triu(djac_out[:IA])
630
rhs_top = q_out.conj().T.dot(rhs_temp)
631
scales_pvt = scales[j_pvt[:IA]]
0 commit comments