When I am running this code for psi0 being the ground state of 64-site heisenberg model,
auto tau = 0.0001;
auto expH = toExpH(ampo,tau*Cplx_i);
printfln("%.16c", innerC(psi0, expH, psi0));
I'm expecting the imaginary part of the answer to be around the value of tau*E0, which is estimated by the taylor expansion. However it's half the value.
vN Entropy at center bond b=32 = 0.845274370251
Eigs at center bond b=32: 0.7615 0.0773 0.0773 0.0773 0.0018 0.0018 0.0018
Largest link dim during sweep 100/100 was 82
Largest truncation error: 8.2362e-11
Energy after sweep 100/100 is -28.175424840952
Sweep 100/100 CPU time = 0.433s (Wall time = 0.478s)
Ground State Energy = -28.1754248410
I don't know what's the wrong in this difference between my expectation and result.
hope for your best help.