@@ -41,7 +41,7 @@ peps₀ = InfinitePEPS(rand, ComplexF64, ComplexSpace(2), ComplexSpace(2))
4141````
4242
4343````
44- InfinitePEPS{TensorMap{ComplexF64, ComplexSpace, 1, 4, Vector{ComplexF64}}}(TensorMap{ComplexF64, ComplexSpace, 1, 4, Vector{ComplexF64}}[TensorMap(ℂ^2 ← (ℂ^2 ⊗ ℂ^2 ⊗ (ℂ^2)' ⊗ (ℂ^2)')):
44+ InfinitePEPS{TensorKit. TensorMap{ComplexF64, TensorKit. ComplexSpace, 1, 4, Vector{ComplexF64}}}(TensorKit. TensorMap{ComplexF64, TensorKit. ComplexSpace, 1, 4, Vector{ComplexF64}}[TensorMap(ℂ^2 ← (ℂ^2 ⊗ ℂ^2 ⊗ (ℂ^2)' ⊗ (ℂ^2)')):
4545[:, :, 1, 1, 1] =
4646 0.8343040072662887 + 0.15425705836788395im 0.4612746978522435 + 0.7411151918989216im
4747 0.6640771294125087 + 0.4428356798799721im 0.9163597170532635 + 0.24145695415210522im
@@ -99,7 +99,7 @@ T = InfiniteTransferPEPS(peps₀, dir, row)
9999````
100100
101101````
102- single site MPSKit.InfiniteMPO{Tuple{TensorMap{ComplexF64, ComplexSpace, 1, 4, Vector{ComplexF64}}, TensorMap{ComplexF64, ComplexSpace, 1, 4, Vector{ComplexF64}}}}:
102+ single site MPSKit.InfiniteMPO{Tuple{TensorKit. TensorMap{ComplexF64, TensorKit. ComplexSpace, 1, 4, Vector{ComplexF64}}, TensorKit. TensorMap{ComplexF64, TensorKit. ComplexSpace, 1, 4, Vector{ComplexF64}}}}:
103103╷ ⋮
104104┼ O[1]: (TensorMap(ℂ^2 ← (ℂ^2 ⊗ ℂ^2 ⊗ (ℂ^2)' ⊗ (ℂ^2)')), TensorMap(ℂ^2 ← (ℂ^2 ⊗ ℂ^2 ⊗ (ℂ^2)' ⊗ (ℂ^2)')))
105105╵ ⋮
@@ -137,7 +137,7 @@ mps, env, ϵ = leading_boundary(mps₀, T, VUMPS(; tol=1e-6, verbosity=2));
137137
138138````
139139[ Info: VUMPS init: obj = +5.052950412844e+00 +1.493192627823e-02im err = 8.4684e-01
140- [ Info: VUMPS conv 4: obj = +1.744071149764e+01 +2.162482215518e -08im err = 2.0274771864e -07 time = 2.78 sec
140+ [ Info: VUMPS conv 4: obj = +1.744071149764e+01 +2.162482270345e -08im err = 2.0274771848e -07 time = 4.39 sec
141141
142142````
143143
@@ -149,7 +149,7 @@ norm_vumps = abs(prod(expectation_value(mps, T)))
149149````
150150
151151````
152- 17.440711497641615
152+ 17.440711497641633
153153````
154154
155155This can be compared to the result obtained using CTMRG, where we see that the results match:
@@ -164,8 +164,8 @@ norm_ctmrg = abs(norm(peps₀, env_ctmrg))
164164
165165````
166166[ Info: CTMRG init: obj = -5.556349490423e-01 +1.605938670370e+00im err = 1.0000e+00
167- [ Info: CTMRG conv 37: obj = +1.744071151099e+01 err = 3.2013507863e -07 time = 0.15 sec
168- abs(norm_vumps - norm_ctmrg) / norm_vumps = 7.653157885484634e -10
167+ [ Info: CTMRG conv 37: obj = +1.744071151099e+01 err = 3.2046524013e -07 time = 0.19 sec
168+ abs(norm_vumps - norm_ctmrg) / norm_vumps = 7.653151774414538e -10
169169
170170````
171171
@@ -206,10 +206,11 @@ norm_2x2_ctmrg = abs(norm(peps₀_2x2, env_ctmrg_2x2))
206206
207207````
208208[ Info: VUMPS init: obj = +6.668046237341e+02 -1.267878277078e+01im err = 8.7901e-01
209- [ Info: VUMPS conv 101: obj = +9.723958853252e+04 +1.933478098677e-03im err = 4.1930714637e-07 time = 2.92 sec
209+ ┌ Warning: VUMPS cancel 200: obj = +9.713547617328e+04 -3.270767578613e+02im err = 4.7663743999e-02 time = 8.25 sec
210+ └ @ MPSKit ~/.julia/packages/MPSKit/XpTWn/src/algorithms/groundstate/vumps.jl:73
210211[ Info: CTMRG init: obj = +1.074898090007e+03 -2.096255594496e+02im err = 1.0000e+00
211- [ Info: CTMRG conv 41: obj = +9.723959008610e+04 err = 6.0876980113e -07 time = 0.37 sec
212- abs(norm_2x2_vumps - norm_2x2_ctmrg) / norm_2x2_vumps = 1.5976813021479506e-8
212+ [ Info: CTMRG conv 41: obj = +9.723959008610e+04 err = 6.0700367621e -07 time = 0.66 sec
213+ abs(norm_2x2_vumps - norm_2x2_ctmrg) / norm_2x2_vumps = 0.0010661671847507165
213214
214215````
215216
@@ -252,7 +253,7 @@ transfer_pepo = InfiniteTransferPEPO(peps₀, pepo, 1, 1)
252253````
253254
254255````
255- single site MPSKit.InfiniteMPO{Tuple{TensorMap{ComplexF64, ComplexSpace, 1, 4, Vector{ComplexF64}}, TensorMap{ComplexF64, ComplexSpace, 1, 4, Vector{ComplexF64}}, TensorMap{ComplexF64, ComplexSpace, 2, 4, Vector{ComplexF64}}}}:
256+ single site MPSKit.InfiniteMPO{Tuple{TensorKit. TensorMap{ComplexF64, TensorKit. ComplexSpace, 1, 4, Vector{ComplexF64}}, TensorKit. TensorMap{ComplexF64, TensorKit. ComplexSpace, 1, 4, Vector{ComplexF64}}, TensorKit. TensorMap{ComplexF64, TensorKit. ComplexSpace, 2, 4, Vector{ComplexF64}}}}:
256257╷ ⋮
257258┼ O[1]: (TensorMap(ℂ^2 ← (ℂ^2 ⊗ ℂ^2 ⊗ (ℂ^2)' ⊗ (ℂ^2)')), TensorMap(ℂ^2 ← (ℂ^2 ⊗ ℂ^2 ⊗ (ℂ^2)' ⊗ (ℂ^2)')), TensorMap((ℂ^2 ⊗ (ℂ^2)') ← (ℂ^2 ⊗ ℂ^2 ⊗ (ℂ^2)' ⊗ (ℂ^2)')))
258259╵ ⋮
@@ -269,9 +270,9 @@ norm_pepo = abs(prod(expectation_value(mps_pepo, transfer_pepo)));
269270````
270271
271272````
272- [ Info: VUMPS init: obj = +3.309203535702e+01 -4.227375981213e -01im err = 9.3280e-01
273- [ Info: VUMPS conv 5: obj = +2.483696258643e+02 +2.405967620689e -07im err = 5.0410779189e -08 time = 2.32 sec
274- norm_pepo = 248.36962586426355
273+ [ Info: VUMPS init: obj = +3.309203535702e+01 -4.227375981212e -01im err = 9.3280e-01
274+ [ Info: VUMPS conv 5: obj = +2.483696258643e+02 +2.405968608492e -07im err = 5.0410779107e -08 time = 3.38 sec
275+ norm_pepo = 248.36962586426358
275276
276277````
277278
0 commit comments