@@ -103,3 +103,35 @@ dmrg(H, ψr; nsweeps = 10, maxdim = [10, 20, 40, 60], cutoff = 1.0e-12)
103103
104104println (" \n Run dmrg with free fermion starting state" )
105105dmrg (H, ψ0; nsweeps = 4 , maxdim = 60 , cutoff = 1.0e-12 )
106+
107+ # This will output something like:
108+ #=
109+ ```
110+ (N, Nf) = (20, 10)
111+ U = 1.0
112+
113+ Free fermion starting energy
114+ inner(ψ0', H, ψ0) = -2.3812770621299357
115+
116+ Random state starting energy
117+ inner(ψr', H, ψr) = 10.0
118+
119+ Run dmrg with random starting state
120+ After sweep 1 energy=6.261701784151 maxlinkdim=2 time=0.041
121+ After sweep 2 energy=2.844954346204 maxlinkdim=5 time=0.056
122+ After sweep 3 energy=0.245282430911 maxlinkdim=14 time=0.071
123+ After sweep 4 energy=-1.439072132586 maxlinkdim=32 time=0.098
124+ After sweep 5 energy=-2.220202191945 maxlinkdim=59 time=0.148
125+ After sweep 6 energy=-2.376787647893 maxlinkdim=60 time=0.186
126+ After sweep 7 energy=-2.381484153892 maxlinkdim=60 time=0.167
127+ After sweep 8 energy=-2.381489999291 maxlinkdim=57 time=0.233
128+ After sweep 9 energy=-2.381489999595 maxlinkdim=49 time=0.175
129+ After sweep 10 energy=-2.381489999595 maxlinkdim=49 time=0.172
130+
131+ Run dmrg with free fermion starting state
132+ After sweep 1 energy=-2.381489929965 maxlinkdim=49 time=0.139
133+ After sweep 2 energy=-2.381489999588 maxlinkdim=49 time=0.165
134+ After sweep 3 energy=-2.381489999594 maxlinkdim=48 time=0.161
135+ After sweep 4 energy=-2.381489999594 maxlinkdim=48 time=0.169
136+ ```
137+ =#
0 commit comments