Skip to content

Conversation

@basilib
Copy link
Contributor

@basilib basilib commented Feb 14, 2024

No description provided.

@basilib basilib requested a review from obackhouse February 14, 2024 19:11
@@ -225,8 +225,13 @@ def gradient(self, x, fock=None, out=None):
h1 = -np.dot(c[gf.nphys :, gf_occ.naux :].conj().T, c[gf.nphys :, : gf_occ.naux])
Copy link
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

we probably also need some control over hermiticity here

@obackhouse
Copy link
Collaborator

before merging this we should probably make sure that the behaviour is sound regarding complex eigenvalues -- should AufbauPrinciple throw an exception for complex eigenvalues? Perhaps it can throw away small imaginary parts?

I also don't think that states are necessarily sorted at any point, so the filling might be in the wrong order even just by the real parts -- I guess we should write some unit tests with complex states.

@basilib basilib marked this pull request as draft February 28, 2024 18:42
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

3 participants