I am trying to use the new feature of building more-than-two-operator terms using AutoMPO with the sample code:
ampo += "Sz",1,"Sz",2,"Sz",3;
auto H = MPO(ampo);
auto E0 = dmrg(psi,H,sweeps,"Quiet");
And when running I get the error message
Only at most 2-operator terms allowed for exact AutoMPO conversion to MPO/IQMPO
What is the correct way of handling this?