Replace InfiniteITensorSum with a more generic parametrized type InfiniteSum{T} (#49)

Author: mtfishman

examples/vumps/transfer_matrix_spectrum.jl

@@ -39,7 +39,7 @@ initstate(n) = "↑"
 ψ = InfMPS(s, initstate)
 
 # Form the Hamiltonian
-H = InfiniteITensorSum(model, s; model_params...)
+H = InfiniteSum{MPO}(model, s; model_params...)
 
 # Check translational invariance
 @show norm(contract(ψ.AL[1:N]..., ψ.C[N]) - contract(ψ.C[0], ψ.AR[1:N]...))

examples/vumps/vumps_from_dmrg.jl

@@ -60,7 +60,7 @@ end
 
 s¹ = siteinds(ψ∞, Cell(1))
 
-Σ∞h = InfiniteITensorSum(model, s¹; J=J, h=h)
+Σ∞h = InfiniteSum{MPO}(model, s¹; J=J, h=h)
 
 ψ∞_opt = vumps(Σ∞h, ψ∞)
 

examples/vumps/vumps_heisenberg.jl

@@ -27,7 +27,7 @@ initstate(n) = isodd(n) ? "↑" : "↓"
 model = Model("heisenberg")
 
 # Form the Hamiltonian
-H = InfiniteITensorSum(model, s)
+H = InfiniteSum{MPO}(model, s)
 
 # Check translational invariance
 @show norm(contract(ψ.AL[1:N]..., ψ.C[N]) - contract(ψ.C[0], ψ.AR[1:N]...))

examples/vumps/vumps_hubbard_extended.jl

@@ -36,7 +36,7 @@ model = Model"hubbard"()
 @show model, model_params
 
 # Form the Hamiltonian
-H = InfiniteITensorSum(model, s; model_params...)
+H = InfiniteSum{MPO}(model, s; model_params...)
 
 # Check translational invariance
 println("\nCheck translational invariance of initial infinite MPS")

examples/vumps/vumps_ising.jl

@@ -42,7 +42,7 @@ initstate(n) = "↑"
 ψ = InfMPS(s, initstate)
 
 # Form the Hamiltonian
-H = InfiniteITensorSum(model, s; model_params...)
+H = InfiniteSum{MPO}(model, s; model_params...)
 
 # Check translational invariance
 @show norm(contract(ψ.AL[1:nsite]..., ψ.C[nsite]) - contract(ψ.C[0], ψ.AR[1:nsite]...))

examples/vumps/vumps_ising_extended.jl

@@ -29,7 +29,7 @@ initstate(n) = "↑"
 ψ = InfMPS(s, initstate);
 
 model = Model("ising_extended");
-H = InfiniteITensorSum(model, s; J=J, J₂=J₂, h=h);
+H = InfiniteSum{MPO}(model, s; J=J, J₂=J₂, h=h);
 
 @show norm(contract(ψ.AL[1:N]..., ψ.C[N]) - contract(ψ.C[0], ψ.AR[1:N]...));
 

examples/vumps/vumps_ising_noncontiguous.jl

@@ -28,7 +28,7 @@ initstate(n) = "↑"
 ψ = InfMPS(s, initstate);
 
 model = Model("ising");
-H = InfiniteITensorSum(model, s; J=J, h=h);
+H = InfiniteSum{MPO}(model, s; J=J, h=h);
 #to test the case where the range is larger than the unit cell size
 for x in 1:N
   temp = MPO(3)

examples/vumps/vumps_localham.jl

@@ -22,7 +22,7 @@ h = 1.5
 model = Model(:ising)
 
 # Form the Hamiltonian
-Σ∞h = InfiniteITensorSum(model, s; J=J, h=h)
+Σ∞h = InfiniteSum{MPO}(model, s; J=J, h=h)
 
 χ = 6
 ψ = InfiniteMPS(s; space=χ)

examples/vumps/vumps_localham_qns.jl

@@ -31,7 +31,7 @@ energy_finite, ψfinite = dmrg(Hfinite, ψfinite, sweeps)
 @show energy_finite / Nfinite
 
 # Form the Hamiltonian
-Σ∞h = InfiniteITensorSum(model, s; model_kwargs...)
+Σ∞h = InfiniteSum{MPO}(model, s; model_kwargs...)
 
 d = 1
 χ = [QN("SzParity", 1, 2) => d, QN("SzParity", 0, 2) => d]

src/ITensorInfiniteMPS.jl

@@ -50,7 +50,7 @@ export Cell,
   InfiniteMPS,
   InfiniteCanonicalMPS,
   InfMPS,
-  InfiniteITensorSum,
+  InfiniteSum,
   InfiniteMPO,
   InfiniteSumLocalOps,
   ITensorMap,