Support complex-Hermitian Green's functions in the multi-orbital solver - #1
Merged
Conversation
Allow complex off-diagonal G_ab(tau) / spectra (e.g. from spin-orbit coupling), so the semi-positive-definiteness of the full complex-Hermitian spectral matrix is enforced during continuation rather than dropping the imaginary part. - predict_Gtau: keep the dtype of rho_l so the complex off-diagonal G_ab(tau) is no longer truncated to its real part. - write_Gtau / write_rhol / write_rho: also write the imaginary part (0 for real data, backward compatible for readers of the first columns). Requires admmsolver with complex-Hermitian L1 / SPD support. Co-Authored-By: Claude Opus 4.8 (1M context) <noreply@anthropic.com> Claude-Session: https://claude.ai/code/session_01TMeDqXABvd876kMmVY7K34
Allow optional column_imag when loading G(tau), and write Re/Im for single- and multi-orbital spectrum/Gtau outputs so complex data is preserved end-to-end. Co-authored-by: Cursor <cursoragent@cursor.com>
This file contains hidden or bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
2 participants
Add this suggestion to a batch that can be applied as a single commit.This suggestion is invalid because no changes were made to the code.Suggestions cannot be applied while the pull request is closed.Suggestions cannot be applied while viewing a subset of changes.Only one suggestion per line can be applied in a batch.Add this suggestion to a batch that can be applied as a single commit.Applying suggestions on deleted lines is not supported.You must change the existing code in this line in order to create a valid suggestion.Outdated suggestions cannot be applied.This suggestion has been applied or marked resolved.Suggestions cannot be applied from pending reviews.Suggestions cannot be applied on multi-line comments.Suggestions cannot be applied while the pull request is queued to merge.Suggestion cannot be applied right now. Please check back later.
Summary
Allows complex off-diagonal
G_ab(τ)/ spectra (e.g. from spin–orbit coupling), so the semi-positive-definiteness of the full complex-Hermitian spectral matrix is enforced during continuation instead of dropping the imaginary part.Changes
predict_Gtau(multi-orbital solver): keep the dtype ofrho_lso the complex off-diagonalG_ab(τ)is no longer truncated to its real part.write_Gtau/write_rhol/write_rho: also write the imaginary part (0 for real data; readers of the leading columns are unaffected).Dependency
Requires
admmsolverwith complex-HermitianL1/ SPD support (SpM-lab/admmsolver#1).Validation
An artificial two-orbital spectrum rotated by a complex unitary
U(θ, φ)(φ ≠ 0) is continued end-to-end: the complex off-diagonal and the Hermitian positive-semidefiniteness are recovered (previously the imaginary part was lost). Existing tests pass.🤖 Generated with Claude Code
https://claude.ai/code/session_01TMeDqXABvd876kMmVY7K34