expH for time evolution

Question

Hi

I am facing a Segmentation fault error with the function expH (my goal is to simulate the dynamics of a spin boson hamiltonian with long-range interactions)

Please find below a minimal example which causes the crash:

It seems to be due to the interplay between the Complex argument Cplxi + the flip-flip terms Sm S_m..

Thanks for your help

Benoit

int main(int argc, char* argv[])
{
auto sites = SpinHalf(3);
auto ampo = AutoMPO(sites);
ampo += "Sm",1,"Sm",2;
ampo += "Sp",1,"Sp",2;
ampo += "Sp",2,"Sp",3;
ampo += "Sm",2,"Sm",3;
auto H = MPO(ampo);
PrintData(H);
auto tau = 0.1;
auto expH = toExpH(ampo,tau*Cplx_i);
PrintData(expH);
return 0;
}
~

Comments

Hi Benoit,

Which version of ITensor are you using? Using version 2.0.11 retrieved on May 20, I can get the code running, provided that I change the toExpH line to

auto expH = toExpH<ITensor>(ampo,tau*Cplx_i);

Hope that helps!

Best,
Chengshu

---

Hi Chengshu,

I had itensor 2.0.11 but I just installed the new version 2.1.0 and now it works.
Thanks a lot!

Benoit

Answer 1

Marking this as answered because of the discussion above. Thanks Chengshu!