@@ -74,7 +74,7 @@ implemented_symmetries = [(Trivial, Trivial),
7474 end
7575end
7676
77- function hamiltonian (particle_symmetry, spin_symmetry; t, J, mu, L, slave_fermion)
77+ function tjhamiltonian (particle_symmetry, spin_symmetry; t, J, mu, L, slave_fermion)
7878 num = tJ. e_number (particle_symmetry, spin_symmetry; slave_fermion)
7979 hop_heis = (- t) * (tJ. e_plusmin (particle_symmetry, spin_symmetry; slave_fermion) +
8080 tJ. e_minplus (particle_symmetry, spin_symmetry; slave_fermion)) +
9898 mu = randn ()
9999
100100 for slave_fermion in (false , true )
101- H_triv = hamiltonian (Trivial, Trivial; t, J, mu, L, slave_fermion)
101+ H_triv = tjhamiltonian (Trivial, Trivial; t, J, mu, L, slave_fermion)
102102 vals_triv = mapreduce (vcat, eigvals (H_triv)) do (c, v)
103103 return repeat (real .(v), dim (c))
104104 end
108108 if (particle_symmetry, spin_symmetry) == (Trivial, Trivial)
109109 continue
110110 end
111- H_symm = hamiltonian (particle_symmetry, spin_symmetry; t, J, mu, L,
112- slave_fermion)
111+ H_symm = tjhamiltonian (particle_symmetry, spin_symmetry; t, J, mu, L,
112+ slave_fermion)
113113 vals_symm = mapreduce (vcat, eigvals (H_symm)) do (c, v)
114114 return repeat (real .(v), dim (c))
115115 end
0 commit comments