Exact diagonalization of the S=1/2 Heisenberg chain
(using semi-momentum states and spin-inversion symmetry)
This program
 Source code for diagonalization hchain_mkpz.f90 dsyev.f (LAPAC diagonalization routine; real symmetric matrix)
 Source code for expectation values hchain_temp.f90
 Running instructions and example results Diagonalization instructions and results Thermodynamics instructions and results