mirror of
https://github.com/triqs/dft_tools
synced 2024-07-17 16:33:46 +02:00
remove debug output, revert to pre-PR
This commit is contained in:
parent
39f6af050c
commit
d36467402e
@ -416,9 +416,8 @@ class ProjectorGroup:
|
|||||||
overlap = np.dot(p_matrix, p_matrix.conj().T)
|
overlap = np.dot(p_matrix, p_matrix.conj().T)
|
||||||
# Calculate [O^{-1/2}]_{m m'}
|
# Calculate [O^{-1/2}]_{m m'}
|
||||||
eig, eigv = np.linalg.eigh(overlap)
|
eig, eigv = np.linalg.eigh(overlap)
|
||||||
eig = np.around(eig,10)
|
assert np.all(eig > 0.0), ("Negative eigenvalues of the overlap matrix:"
|
||||||
#assert np.all(eig > 0.0), ("Negative eigenvalues of the overlap matrix:"
|
"projectors are ill-defined")
|
||||||
# "projectors are ill-defined")
|
|
||||||
sqrt_eig = 1.0 / np.sqrt(eig)
|
sqrt_eig = 1.0 / np.sqrt(eig)
|
||||||
shalf = np.dot(eigv * sqrt_eig, eigv.conj().T)
|
shalf = np.dot(eigv * sqrt_eig, eigv.conj().T)
|
||||||
# Apply \tilde{P}_{m v} = \sum_{m'} [O^{-1/2}]_{m m'} P_{m' v}
|
# Apply \tilde{P}_{m v} = \sum_{m'} [O^{-1/2}]_{m m'} P_{m' v}
|
||||||
|
Loading…
Reference in New Issue
Block a user