Cleanup

Author: lkdvos

Project.toml

@@ -4,7 +4,6 @@ authors = ["Maarten Van Damme", "Lukas Devos", "Gertian Roose", "Klaas Gunst"]
 version = "0.4.0"
 
 [deps]
-BlockTensorKit = "5f87ffc2-9cf1-4a46-8172-465d160bd8cd"
 LinearAlgebra = "37e2e46d-f89d-539d-b4ee-838fcccc9c8e"
 MPSKit = "bb1c41ca-d63c-52ed-829e-0820dda26502"
 MacroTools = "1914dd2f-81c6-5fcd-8719-6d5c9610ff09"

src/MPSKitModels.jl

@@ -1,8 +1,8 @@
 module MPSKitModels
 
-using TensorKit, MPSKit, BlockTensorKit
+using TensorKit, MPSKit
 using MacroTools: @capture, postwalk
-using MPSKit: @plansor, _lastspace, _firstspace
+using MPSKit: @plansor
 using TensorOperations
 using TupleTools
 

src/lattices/latticepoints.jl

@@ -61,5 +61,7 @@ end
 Base.:-(i::Int, j::LatticePoint{1}) = LatticePoint(i .- j, j.lattice)
 
 Base.isless(i::L, j::L) where {L<:LatticePoint} = linearize_index(i) < linearize_index(j)
-Base.isfinite(::Union{LatticePoint{N,G}, Type{<:LatticePoint{N,G}}}) where {N,G} = isfinite(G)
+function Base.isfinite(::Union{LatticePoint{N,G},Type{<:LatticePoint{N,G}}}) where {N,G}
+    return isfinite(G)
+end
 latticetype(::Union{LatticePoint{N,G},Type{<:LatticePoint{N,G}}}) where {N,G} = G

src/models/transfermatrices.jl

@@ -30,7 +30,7 @@ function classical_ising(elt::Type{<:Number}=ComplexF64, ::Type{Trivial}=Trivial
 
     @tensor o[-1 -2; -3 -4] := O[1 2; 3 4] * nt[-1; 1] * nt[-2; 2] * nt[-3; 3] * nt[-4; 4]
 
-    return InfiniteMPO(TensorMap(o, ℂ^2 * ℂ^2, ℂ^2 * ℂ^2))
+    return InfiniteMPO([TensorMap(o, ℂ^2 * ℂ^2, ℂ^2 * ℂ^2)])
 end
 
 function classical_ising(elt::Type{<:Number}, ::Type{Z2Irrep}; beta=log(1 + sqrt(2)) / 2)
@@ -42,7 +42,7 @@ function classical_ising(elt::Type{<:Number}, ::Type{Z2Irrep}; beta=log(1 + sqrt
     block(mpo, Irrep[ℤ₂](0)) .= [2x^2 2x*y; 2x*y 2y^2]
     block(mpo, Irrep[ℤ₂](1)) .= [2x*y 2x*y; 2x*y 2x*y]
 
-    return InfiniteMPO(mpo)
+    return InfiniteMPO([mpo])
 end
 
 #===========================================================================================
@@ -63,23 +63,23 @@ function sixvertex(elt::Type{<:Number}=ComplexF64, ::Type{Trivial}=Trivial; a=1.
             0 c b 0
             0 b c 0
             0 0 0 a]
-    return InfiniteMPO(permute(TensorMap(d, ℂ^2 ⊗ ℂ^2, ℂ^2 ⊗ ℂ^2), ((1, 2), (4, 3))))
+    return InfiniteMPO([permute(TensorMap(d, ℂ^2 ⊗ ℂ^2, ℂ^2 ⊗ ℂ^2), ((1, 2), (4, 3)))])
 end
 function sixvertex(elt::Type{<:Number}, ::Type{U1Irrep}; a=1.0, b=1.0, c=1.0)
     pspace = U1Space(-1 // 2 => 1, 1 // 2 => 1)
     mpo = zeros(elt, pspace ⊗ pspace, pspace ⊗ pspace)
-    blocks(mpo, Irrep[U₁](0)) .= [b c; c b]
+    block(mpo, Irrep[U₁](0)) .= [b c; c b]
     block(mpo, Irrep[U₁](1)) .= reshape([a], (1, 1))
     block(mpo, Irrep[U₁](-1)) .= reshape([a], (1, 1))
-    return InfiniteMPO(permute(mpo, ((1, 2), (4, 3))))
+    return InfiniteMPO([permute(mpo, ((1, 2), (4, 3)))])
 end
 function sixvertex(elt::Type{<:Number}, ::Type{CU1Irrep}; a=1.0, b=1.0, c=1.0)
     pspace = CU1Space(1 // 2 => 1)
     mpo = zeros(elt, pspace ⊗ pspace, pspace ⊗ pspace)
     block(mpo, Irrep[CU₁](0, 0)) .= reshape([b + c], (1, 1))
     block(mpo, Irrep[CU₁](0, 1)) .= reshape([-b + c], (1, 1))
     block(mpo, Irrep[CU₁](1, 2)) .= reshape([a], (1, 1))
-    return InfiniteMPO(permute(mpo, ((1, 2), (4, 3))))
+    return InfiniteMPO([permute(mpo, ((1, 2), (4, 3)))])
 end
 
 #===========================================================================================
@@ -95,7 +95,7 @@ function hard_hexagon(elt::Type{<:Number}=ComplexF64)
     P = Vect[FibonacciAnyon](:τ => 1)
     O = ones(elt, P ⊗ P ← P ⊗ P)
     block(O, FibonacciAnyon(:I)) .*= 0
-    return InfiniteMPO(O)
+    return InfiniteMPO([O])
 end
 
 #===========================================================================================
@@ -116,5 +116,5 @@ function qstate_clock(elt::Type{<:Number}=ComplexF64, ::Type{Trivial}=Trivial;
                             (comega(i - j) + comega(j - k) + comega(k - l) + comega(l - i)))
     end
 
-    return InfiniteMPO(TensorMap(O, ℂ^q * ℂ^q, ℂ^q * ℂ^q))
+    return InfiniteMPO([TensorMap(O, ℂ^q * ℂ^q, ℂ^q * ℂ^q)])
 end

test/sixvertex.jl

@@ -16,21 +16,21 @@ end
 
 @testset "U₁" begin
     mpo = @inferred sixvertex(ComplexF64, U1Irrep)
-    mpo2 = MPOMultiline([mpo.opp mpo.opp; mpo.opp mpo.opp])
+    mpo2 = repeat(mpo, 2, 2)
     vspaces = [U1Space(0 => 20, 1 => 10, -1 => 10, 2 => 5, -2 => 5),
                U1Space(1 // 2 => 15, -1 // 2 => 15, 3 // 2 => 5, -3 // 2 => 5)]
-    ψ = MPSMultiline(repeat(space.(mpo.opp, 2), 2, 2), [vspaces circshift(vspaces, 1)])
+    ψ = MultilineMPS(repeat(physicalspace(mpo), 2, 2), [vspaces circshift(vspaces, 1)])
     ψ, _ = leading_boundary(ψ, mpo2, alg)
     F = prod(expectation_value(ψ, mpo2))
     @test abs(F)^(1 / 4) ≈ F₀ atol = 1e-2
 end
 
 @testset "CU₁" begin
     mpo = @inferred sixvertex(ComplexF64, CU1Irrep)
-    mpo2 = MPOMultiline([mpo.opp mpo.opp; mpo.opp mpo.opp])
+    mpo2 = repeat(mpo, 2, 2)
     vspaces = [CU1Space((0, 0) => 10, (0, 1) => 10, (1, 2) => 5, (2, 2) => 5),
                CU1Space((1 // 2, 2) => 15, (3 // 2, 2) => 5)]
-    ψ = MPSMultiline(repeat(space.(mpo.opp, 2), 2, 2), [vspaces circshift(vspaces, 1)])
+    ψ = MultilineMPS(repeat(physicalspace(mpo), 2, 2), [vspaces circshift(vspaces, 1)])
     ψ, _ = leading_boundary(ψ, mpo2, alg)
     F = prod(expectation_value(ψ, mpo2))
     @test abs(F)^(1 / 4) ≈ F₀ atol = 1e-2