Merge branch 'master' of github.com:quantumghent/PEPSKit.jl

Author: pbrehmer

.github/workflows/Tests.yml

@@ -41,7 +41,7 @@ jobs:
       group: "${{ matrix.group }}"
       julia-version: "${{ matrix.version }}"
       os: "${{ matrix.os }}"
-      timeout-minutes: 75
+      timeout-minutes: 120
     secrets:
       CODECOV_TOKEN: ${{ secrets.CODECOV_TOKEN }}
 

src/algorithms/time_evolution/simpleupdate.jl

@@ -70,6 +70,11 @@ function su_iter(
     if bipartite
         @assert Nr == Nc == 2
     end
+    (Nr >= 2 && Nc >= 2) || throw(
+        ArgumentError(
+            "iPEPS unit cell size for simple update should be no smaller than (2, 2)."
+        ),
+    )
     # TODO: make algorithm independent on the choice of dual in the network
     for (r, c) in Iterators.product(1:Nr, 1:Nc)
         @assert [isdual(space(peps.vertices[r, c], ax)) for ax in 1:5] == [0, 1, 1, 0, 0]
@@ -132,10 +137,6 @@ function simpleupdate(
     check_interval::Int=500,
 )
     time_start = time()
-    Nr, Nc = size(peps)
-    if bipartite
-        @assert Nr == Nc == 2
-    end
     # exponentiating the 2-site Hamiltonian gate
     gate = get_gate(alg.dt, ham)
     wtdiff = 1.0

src/states/infiniteweightpeps.jl

@@ -49,6 +49,17 @@ function compare_weights(wts1::SUWeight, wts2::SUWeight)
     return sum(_singular_value_distance, zip(wts1.data, wts2.data)) / length(wts1)
 end
 
+function Base.show(io::IO, ::MIME"text/plain", wts::SUWeight)
+    println(io, typeof(wts))
+    for idx in CartesianIndices(wts.data)
+        println(io, Tuple(idx), ":")
+        for (k, b) in blocks(wts.data[idx])
+            println(io, k, " = ", diag(b))
+        end
+    end
+    return nothing
+end
+
 """
     struct InfiniteWeightPEPS{T<:PEPSTensor,E<:PEPSWeight}