@@ -883,32 +883,32 @@ end
883883@testset " Finite temperature methods" begin
884884 L = 6
885885 H = transverse_field_ising (; L)
886-
886+
887887 beta = 0.1
888-
888+
889889 # exact diagonalization
890890 H_dense = convert (TensorMap, H)
891- Z_dense_1 = tr (exp (- beta * H_dense))^ (1 / N)
892- Z_dense_2 = tr (exp (- 2 beta * H_dense))^ (1 / N)
891+ Z_dense_1 = tr (exp (- beta * H_dense))^ (1 / N)
892+ Z_dense_2 = tr (exp (- 2 beta * H_dense))^ (1 / N)
893893
894894 # taylor cluster
895895 rho_taylor_1 = make_time_mpo (H, - im * beta, TaylorCluster (; N= 2 ))
896- Z_taylor_1 = tr (rho_taylor_1)^ (1 / N)
896+ Z_taylor_1 = tr (rho_taylor_1)^ (1 / N)
897897 @test Z_taylor_1 ≈ Z_dense_1 atol = 1e-2
898- Z_taylor_2 = real (dot (rho_taylor_1, rho_taylor_1))^ (1 / N)
898+ Z_taylor_2 = real (dot (rho_taylor_1, rho_taylor_1))^ (1 / N)
899899 @test Z_taylor_2 ≈ Z_dense_2 atol = 1e-2
900-
900+
901901 # MPO multiplication
902902 rho_mps = convert (FiniteMPS, rho_taylor_1)
903903 rho_mps, = approximate (rho_taylor_1, (rho_taylor_1, rho_mps), DMRG (; verbosity= 0 ))
904- Z_mpomul = tr (rho_mps)^ (1 / N)
904+ Z_mpomul = tr (rho_mps)^ (1 / N)
905905 @test Z_mpomul ≈ Z_dense_2 atol = 1e-2
906-
906+
907907 # TDVP
908908 rho_0 = infinite_temperature_density_matrix (H)
909909 rho_0_mps = convert (FiniteMPS, rho_0)
910910 rho_mps = timestep (rho_0_mps, H, 0.0 , - im * beta, TDVP ())
911- Z_tdvp = tr (rho_mps)^ (1 / N)
911+ Z_tdvp = tr (rho_mps)^ (1 / N)
912912 @test Z_tdvp ≈ Z_dense_1 atol = 1e-2
913913end
914914
0 commit comments