Merge branch 'main' into ITensorMPS_random_mps

Author: mtfishman

src/lib/ITensorMPS/examples/finite_temperature/metts.jl

@@ -108,7 +108,7 @@ function main(; N=10, cutoff=1E-8, δτ=0.1, beta=2.0, NMETTS=3000, Nwarm=10)
       samp = sample!(psi)
       new_state = [samp[j] == 1 ? "Z+" : "Z-" for j in 1:N]
     end
-    psi = productMPS(s, new_state)
+    psi = MPS(s, new_state)
   end
 
   return nothing

src/lib/ITensorMPS/examples/gate_evolution/mpo_gate_evolution.jl

@@ -31,7 +31,7 @@ function main(; N=10, cutoff=1E-8, δt=0.1, ttotal=5.0)
   end
 
   # Initialize psi to be a product state (alternating up and down)
-  psi0 = productMPS(s, n -> isodd(n) ? "Up" : "Dn")
+  psi0 = MPS(s, n -> isodd(n) ? "Up" : "Dn")
 
   c = div(N, 2)
 

src/lib/ITensorMPS/examples/gate_evolution/quantum_simulator.jl

@@ -9,17 +9,17 @@ H = ops(s, [("H", n) for n in 1:N])
 CX = ops(s, [("CX", n, m) for n in 1:N, m in 1:N])
 
 # Start with the state |0000...⟩
-ψ0 = productMPS(s, "0")
+ψ0 = MPS(s, "0")
 
 # Change to the state |1010...⟩
 gates = [X[n] for n in 1:2:N]
 ψ = apply(gates, ψ0; cutoff=1e-15)
-@assert inner(ψ, productMPS(s, n -> isodd(n) ? "1" : "0")) ≈ 1
+@assert inner(ψ, MPS(s, n -> isodd(n) ? "1" : "0")) ≈ 1
 
 # Change to the state |10111011...⟩
 append!(gates, [CX[n, n + 3] for n in 1:4:(N - 3)])
 ψ = apply(gates, ψ0; cutoff=1e-15)
-@assert inner(ψ, productMPS(s, ["1", "0", "1", "1", "1", "0", "1", "1", "1", "0"])) ≈ 1
+@assert inner(ψ, MPS(s, ["1", "0", "1", "1", "1", "0", "1", "1", "1", "0"])) ≈ 1
 
 # Change the state |10111011...⟩ to the (|+⟩, |-⟩) basis
 append!(gates, [H[n] for n in 1:N])

src/lib/ITensorMPS/test/base/test_autompo.jl

@@ -912,15 +912,15 @@ end
     a3 += "Cdag", 1, "N", 2, "C", 3
     M3 = MPO(a3, s)
 
-    p011 = productMPS(s, [1, 2, 2, 1, 1])
-    p110 = productMPS(s, [2, 2, 1, 1, 1])
+    p011 = MPS(s, [1, 2, 2, 1, 1])
+    p110 = MPS(s, [2, 2, 1, 1, 1])
 
     @test inner(p110', M1, p011) ≈ -1.0
     @test inner(p110', M2, p011) ≈ -1.0
     @test inner(p110', M3, p011) ≈ -1.0
 
-    p001 = productMPS(s, [1, 1, 2, 1, 1])
-    p100 = productMPS(s, [2, 1, 1, 1, 1])
+    p001 = MPS(s, [1, 1, 2, 1, 1])
+    p100 = MPS(s, [2, 1, 1, 1, 1])
 
     @test inner(p100', M1, p001) ≈ +1.0
     @test inner(p100', M2, p001) ≈ +1.0
@@ -941,14 +941,14 @@ end
     a2 += -1, "Cdn", 3, "Cdagdn", 1
     M2 = MPO(a2, s)
 
-    p0uu = productMPS(s, [1, 2, 2, 1, 1])
-    puu0 = productMPS(s, [2, 2, 1, 1, 1])
-    p0ud = productMPS(s, [1, 2, 3, 1, 1])
-    pdu0 = productMPS(s, [3, 2, 1, 1, 1])
-    p00u = productMPS(s, [1, 1, 2, 1, 1])
-    pu00 = productMPS(s, [2, 1, 1, 1, 1])
-    p00d = productMPS(s, [1, 1, 3, 1, 1])
-    pd00 = productMPS(s, [3, 1, 1, 1, 1])
+    p0uu = MPS(s, [1, 2, 2, 1, 1])
+    puu0 = MPS(s, [2, 2, 1, 1, 1])
+    p0ud = MPS(s, [1, 2, 3, 1, 1])
+    pdu0 = MPS(s, [3, 2, 1, 1, 1])
+    p00u = MPS(s, [1, 1, 2, 1, 1])
+    pu00 = MPS(s, [2, 1, 1, 1, 1])
+    p00d = MPS(s, [1, 1, 3, 1, 1])
+    pd00 = MPS(s, [3, 1, 1, 1, 1])
 
     @test inner(puu0', M1, p0uu) ≈ -1.0
     @test inner(pdu0', M2, p0ud) ≈ -1.0
@@ -1039,8 +1039,8 @@ end
         os += -1im, "S-", i, "S+", i + 1
       end
       H = MPO(os, sites)
-      psiud = productMPS(sites, [1, 2, 1, 2])
-      psidu = productMPS(sites, [2, 1, 1, 2])
+      psiud = MPS(sites, [1, 2, 1, 2])
+      psidu = MPS(sites, [2, 1, 1, 2])
       @test inner(psiud', H, psidu) ≈ +1im
       @test inner(psidu', H, psiud) ≈ -1im
     end
@@ -1146,7 +1146,7 @@ end
       os += V2, "N", j, "N", j + 2
     end
     H = MPO(os, sites)
-    psi0 = productMPS(sites, n -> isodd(n) ? "0" : "1")
+    psi0 = MPS(sites, n -> isodd(n) ? "0" : "1")
     @test abs(inner(psi0', H, psi0) - 0.00018) < 1E-10
   end
 

src/lib/ITensorMPS/test/base/test_dmrg.jl

@@ -386,7 +386,7 @@ using ITensors.ITensorMPS: nsite, set_nsite!, site_range
     state[3] = 2
     state[5] = 2
     state[7] = 2
-    psi0 = productMPS(s, state)
+    psi0 = MPS(s, state)
 
     os = OpSum()
     for j in 1:(N - 1)

src/lib/ITensorMPS/test/base/test_mps.jl

@@ -61,7 +61,7 @@ include(joinpath(@__DIR__, "utils", "util.jl"))
     @test psi ⋅ psi ≈ *(dag(psi)..., psi...)[]
   end
 
-  @testset "productMPS" begin
+  @testset "MPS" begin
     @testset "vector of string input" begin
       sites = siteinds("S=1/2", 10)
       state = fill("", length(sites))
@@ -73,13 +73,13 @@ include(joinpath(@__DIR__, "utils", "util.jl"))
         sign = isodd(j) ? +1.0 : -1.0
         @test (psi[j] * op(sites, "Sz", j) * dag(prime(psi[j], "Site")))[] ≈ sign / 2
       end
-      psi = productMPS(sites, state)
+      psi = MPS(sites, state)
       for j in 1:length(psi)
         sign = isodd(j) ? +1.0 : -1.0
         @test (psi[j] * op(sites, "Sz", j) * dag(prime(psi[j], "Site")))[] ≈ sign / 2
       end
       @test_throws DimensionMismatch MPS(sites, fill("", length(psi) - 1))
-      @test_throws DimensionMismatch productMPS(sites, fill("", length(psi) - 1))
+      @test_throws DimensionMismatch MPS(sites, fill("", length(psi) - 1))
     end
 
     @testset "String input" begin
@@ -89,7 +89,7 @@ include(joinpath(@__DIR__, "utils", "util.jl"))
         sign = -1.0
         @test (psi[j] * op(sites, "Sz", j) * dag(prime(psi[j], "Site")))[] ≈ sign / 2
       end
-      psi = productMPS(sites, "Dn")
+      psi = MPS(sites, "Dn")
       for j in 1:length(psi)
         sign = -1.0
         @test (psi[j] * op(sites, "Sz", j) * dag(prime(psi[j], "Site")))[] ≈ sign / 2
@@ -108,7 +108,7 @@ include(joinpath(@__DIR__, "utils", "util.jl"))
         sign = -1.0
         @test (psi[j] * op(sites, "Sz", j) * dag(prime(psi[j], "Site")))[] ≈ sign / 2
       end
-      psi = productMPS(sites, 2)
+      psi = MPS(sites, 2)
       for j in 1:length(psi)
         sign = -1.0
         @test (psi[j] * op(sites, "Sz", j) * dag(prime(psi[j], "Site")))[] ≈ sign / 2
@@ -126,7 +126,7 @@ include(joinpath(@__DIR__, "utils", "util.jl"))
         sign = isodd(j) ? +1.0 : -1.0
         @test (psi[j] * op(sites, "Sz", j) * dag(prime(psi[j], "Site")))[] ≈ sign / 2
       end
-      psi = productMPS(sites, state)
+      psi = MPS(sites, state)
       for j in 1:length(psi)
         sign = isodd(j) ? +1.0 : -1.0
         @test (psi[j] * op(sites, "Sz", j) * dag(prime(psi[j], "Site")))[] ≈ sign / 2
@@ -145,7 +145,7 @@ include(joinpath(@__DIR__, "utils", "util.jl"))
         sign = isodd(j) ? +1.0 : -1.0
         @test (psi[j] * op(sites, "Sz", j) * dag(prime(psi[j], "Site")))[] ≈ sign / 2
       end
-      psi = productMPS(ivals)
+      psi = MPS(ivals)
       for j in 1:length(psi)
         sign = isodd(j) ? +1.0 : -1.0
         @test (psi[j] * op(sites, "Sz", j) * dag(prime(psi[j], "Site")))[] ≈ sign / 2
@@ -157,7 +157,7 @@ include(joinpath(@__DIR__, "utils", "util.jl"))
         for j in 1:length(psi)
           @test eltype(psi[j]) == ComplexF64
         end
-        psi = productMPS(ComplexF64, sites, fill(1, length(psi)))
+        psi = MPS(ComplexF64, sites, fill(1, length(psi)))
         for j in 1:length(psi)
           @test eltype(psi[j]) == ComplexF64
         end
@@ -171,13 +171,13 @@ include(joinpath(@__DIR__, "utils", "util.jl"))
       psi = MPS([site], [1])
       @test psi[1][1] ≈ 1.0
       @test psi[1][2] ≈ 0.0
-      psi = productMPS([site], [1])
+      psi = MPS([site], [1])
       @test psi[1][1] ≈ 1.0
       @test psi[1][2] ≈ 0.0
       psi = MPS([site], [2])
       @test psi[1][1] ≈ 0.0
       @test psi[1][2] ≈ 1.0
-      psi = productMPS([site], [2])
+      psi = MPS([site], [2])
       @test psi[1][1] ≈ 0.0
       @test psi[1][2] ≈ 1.0
     end

src/lib/ITensorMPS/test/base/test_threading.jl

@@ -44,7 +44,7 @@ end
     H = MPO(opsum, sites)
     Hsplit = splitblocks(linkinds, H)
     state = [isodd(n) ? "↑" : "↓" for n in 1:N]
-    ψ0 = productMPS(sites, state)
+    ψ0 = MPS(sites, state)
     enabled = ITensors.enable_threaded_blocksparse(true)
     energy, _ = dmrg(H, ψ0, sweeps; outputlevel=outputlevel)
     energy_split, _ = dmrg(Hsplit, ψ0, sweeps; outputlevel=outputlevel)