Extracts overlap and Fock matrices from LCAO/CRYSTAL output files. Handles spin channels (alpha-alpha, beta-beta, alpha-beta, beta-alpha) and constructs spin-block matrices. Includes Hermiticity ...