We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
There was an error while loading. Please reload this page.
1 parent 335ff0b commit 8d378c6Copy full SHA for 8d378c6
examples/windowmps.jl
@@ -2,8 +2,8 @@ using MPSKit, MPSKitModels, TensorKit, Plots
2
3
let
4
#defining the hamiltonian
5
- th = nonsym_ising_ham(; lambda = 0.3)
6
- sx, sy, sz = nonsym_spintensors(1 // 2)
+ th = transverse_field_ising(; g = 0.3)
+ sx, sy, sz = σˣ(ComplexF64), σʸ(ComplexF64), σᶻ(ComplexF64)
7
8
#initilizing a random mps
9
ts = InfiniteMPS([ℂ^2], [ℂ^12])
@@ -23,11 +23,12 @@ let
23
24
envs = environments(mpco, th)
25
26
- szdat = [expectation_value(mpco, sz)]
+ szdat = [expectation_value(mpco, i => sz) for i in 1:length(mpco)]
27
+ szdat = [szdat]
28
29
for i in 1:(totaltime / deltat)
- mpco, envs = timestep(mpco, th, deltat, TDVP2(; trscheme = truncrank(20)), envs)
30
- push!(szdat, expectation_value(mpco, sz))
+ mpco, envs = timestep(mpco, th, 0, deltat, TDVP2(; trscheme = truncbelow(10^(-8)) & truncdim(25)), envs)
31
+ push!(szdat, [expectation_value(mpco, i => sz) for i in 1:length(mpco)])
32
end
33
34
display(heatmap(real.(reduce((a, b) -> [a b], szdat))))
0 commit comments