|
7 | 7 | from varipeps.peps import PEPS_Unit_Cell |
8 | 8 |
|
9 | 9 |
|
10 | | -def calculate_triangular_correlation_length(unitcell: PEPS_Unit_Cell, num_eigvals: int = 8): |
| 10 | +def calculate_triangular_correlation_length( |
| 11 | + unitcell: PEPS_Unit_Cell, num_eigvals: int = 8 |
| 12 | +): |
11 | 13 | if num_eigvals < 2: |
12 | | - raise ValueError("Number of eigenvalues must be at least two to compute the correlation length.") |
| 14 | + raise ValueError( |
| 15 | + "Number of eigenvalues must be at least two to compute the correlation length." |
| 16 | + ) |
13 | 17 |
|
14 | 18 | initial_vector_60 = apply_contraction_jitted( |
15 | 19 | "triangular_ctmrg_corrlen_vec_60", |
@@ -111,14 +115,18 @@ def matvec_180(vec): |
111 | 115 |
|
112 | 116 | corr_len_60 = -1 / np.log(np.abs(eig_60[1])) |
113 | 117 |
|
114 | | - eig_120, eigvec_120 = eigs(lin_op_120, k=num_eigvals, v0=initial_vector_120, which="LM") |
| 118 | + eig_120, eigvec_120 = eigs( |
| 119 | + lin_op_120, k=num_eigvals, v0=initial_vector_120, which="LM" |
| 120 | + ) |
115 | 121 |
|
116 | 122 | eig_120 = eig_120[np.argsort(np.abs(eig_120))[::-1]] |
117 | 123 | eig_120 /= np.abs(eig_120[0]) |
118 | 124 |
|
119 | 125 | corr_len_120 = -1 / np.log(np.abs(eig_120[1])) |
120 | 126 |
|
121 | | - eig_180, eigvec_180 = eigs(lin_op_180, k=num_eigvals, v0=initial_vector_180, which="LM") |
| 127 | + eig_180, eigvec_180 = eigs( |
| 128 | + lin_op_180, k=num_eigvals, v0=initial_vector_180, which="LM" |
| 129 | + ) |
122 | 130 |
|
123 | 131 | eig_180 = eig_180[np.argsort(np.abs(eig_180))[::-1]] |
124 | 132 | eig_180 /= np.abs(eig_180[0]) |
|
0 commit comments