11using MPSKit
2+ using MPSKit: _transpose_front, _transpose_tail
23using MPSKit: GeometryStyle, InfiniteChainStyle, TransferMatrix
34using TensorKit
45using TensorKit: ℙ
@@ -9,10 +10,10 @@ using Adapt, CUDA
910 tol = Float64 (eps (real (elt)) * 100 )
1011
1112 ψ = adapt (CuArray, InfiniteMPS ([rand (elt, D * d, D), rand (elt, D * d, D)]; tol))
12- @test TensorKit. storagetype (ψ) == CuVector{ComplexF64}
13+ @test TensorKit. storagetype (ψ) == CuVector{ComplexF64, CUDA . DeviceMemory }
1314 @test eltype (ψ) == eltype (typeof (ψ))
1415
15- for i in 1 : length (ψ)
16+ #= for i in 1:length(ψ)
1617 @plansor difference[-1 -2; -3] := ψ.AL[i][-1 -2; 1] * ψ.C[i][1; -3] -
1718 ψ.C[i - 1][-1; 1] * ψ.AR[i][1 -2; -3]
1819 @test norm(difference, Inf) < tol * 10
@@ -26,19 +27,19 @@ using Adapt, CUDA
2627 @test TransferMatrix(ψ.AL[i], ψ.AR[i]) * r_LR(ψ, i) ≈ r_LR(ψ, i + 1)
2728 @test TransferMatrix(ψ.AR[i], ψ.AL[i]) * r_RL(ψ, i) ≈ r_RL(ψ, i + 1)
2829 @test TransferMatrix(ψ.AR[i], ψ.AR[i]) * r_RR(ψ, i) ≈ r_RR(ψ, i + 1)
29- end
30+ end=# # TODO
3031
3132 L = rand (3 : 20 )
3233 ψ = adapt (CuArray, FiniteMPS (rand, elt, L, d, D))
33- @test TensorKit. storagetype (ψ) == CuVector{ComplexF64}
34+ @test TensorKit. storagetype (ψ) == CuVector{ComplexF64, CUDA . DeviceMemory }
3435 @test eltype (ψ) == eltype (typeof (ψ))
35- ovl = dot (ψ, ψ)
36+ #= ovl = dot(ψ, ψ)
3637
3738 @test ovl ≈ norm(ψ.AC[1])^2
3839
3940 for i in 1:length(ψ)
4041 @test ψ.AC[i] ≈ ψ.AL[i] * ψ.C[i]
41- @test ψ. AC[i] ≈ _transpose_front (ψ. C[i - 1 ] * _transpose_tail (ψ. AR[i]))
42+ # @test ψ.AC[i] ≈ _transpose_front(ψ.C[i - 1] * _transpose_tail(ψ.AR[i])) # TODO
4243 end
4344
4445 @test ComplexF64 == scalartype(ψ)
@@ -47,5 +48,5 @@ using Adapt, CUDA
4748 ψ = 3 * ψ
4849 @test ovl * 9 * 9 ≈ norm(ψ)^2
4950
50- @test norm (2 * ψ + ψ - 3 * ψ) ≈ 0.0 atol = sqrt (eps (real (ComplexF64)))
51+ @test norm(2 * ψ + ψ - 3 * ψ) ≈ 0.0 atol = sqrt(eps(real(ComplexF64)))=# # TODO
5152end
0 commit comments