Skip to content

Commit 671e1a7

Browse files
lkdvoslkdvos
committed
Update examples
1 parent 97fb433 commit 671e1a7

File tree

22 files changed

+1214
-3388
lines changed

22 files changed

+1214
-3388
lines changed

‎docs/src/examples/classic2d/1.hard-hexagon/index.md‎

Lines changed: 16 additions & 68 deletions
Large diffs are not rendered by default.

‎docs/src/examples/classic2d/1.hard-hexagon/main.ipynb‎

Lines changed: 17 additions & 15 deletions
Original file line numberDiff line numberDiff line change
@@ -45,7 +45,7 @@
4545
" return Vect[FibonacciAnyon](sector => _D for sector in (:I, :τ))\n",
4646
"end\n",
4747
"\n",
48-
"@assert isapprox(dim(virtual_space(100)), 100; atol=3)"
48+
"@assert isapprox(dim(virtual_space(100)), 100; atol = 3)"
4949
],
5050
"metadata": {},
5151
"execution_count": null
@@ -69,10 +69,10 @@
6969
"D = 10\n",
7070
"V = virtual_space(D)\n",
7171
"ψ₀ = InfiniteMPS([P], [V])\n",
72-
"ψ, envs, = leading_boundary(ψ₀, mpo,\n",
73-
" VUMPS(; verbosity=0,\n",
74-
" alg_eigsolve=MPSKit.Defaults.alg_eigsolve(;\n",
75-
" ishermitian=false))) # use non-hermitian eigensolver\n",
72+
"ψ, envs, = leading_boundary(\n",
73+
" ψ₀, mpo,\n",
74+
" VUMPS(; verbosity = 0, alg_eigsolve = MPSKit.Defaults.alg_eigsolve(; ishermitian = false))\n",
75+
") # use non-hermitian eigensolver\n",
7676
"F = real(expectation_value(ψ, mpo))\n",
7777
"S = real(first(entropy(ψ)))\n",
7878
"ξ = correlation_length(ψ)\n",
@@ -97,8 +97,10 @@
9797
"outputs": [],
9898
"cell_type": "code",
9999
"source": [
100-
"function scaling_simulations(ψ₀, mpo, Ds; verbosity=0, tol=1e-6,\n",
101-
" alg_eigsolve=MPSKit.Defaults.alg_eigsolve(; ishermitian=false))\n",
100+
"function scaling_simulations(\n",
101+
" ψ₀, mpo, Ds; verbosity = 0, tol = 1.0e-6,\n",
102+
" alg_eigsolve = MPSKit.Defaults.alg_eigsolve(; ishermitian = false)\n",
103+
" )\n",
102104
" entropies = similar(Ds, Float64)\n",
103105
" correlations = similar(Ds, Float64)\n",
104106
" alg = VUMPS(; verbosity, tol, alg_eigsolve)\n",
@@ -108,7 +110,7 @@
108110
" correlations[1] = correlation_length(ψ)\n",
109111
"\n",
110112
" for (i, d) in enumerate(diff(Ds))\n",
111-
" ψ, envs = changebonds(ψ, mpo, OptimalExpand(; trscheme=truncdim(d)), envs)\n",
113+
" ψ, envs = changebonds(ψ, mpo, OptimalExpand(; trscheme = truncrank(d)), envs)\n",
112114
" ψ, envs, = leading_boundary(ψ, mpo, alg, envs)\n",
113115
" entropies[i + 1] = real(entropy(ψ)[1])\n",
114116
" correlations[i + 1] = correlation_length(ψ)\n",
@@ -122,9 +124,9 @@
122124
"\n",
123125
"f = fit(log.(ξs), 6 * Ss, 1)\n",
124126
"c = f.coeffs[2]\n",
125-
"p = plot(; xlabel=\"logarithmic correlation length\", ylabel=\"entanglement entropy\")\n",
126-
"p = plot(log.(ξs), Ss; seriestype=:scatter, label=nothing)\n",
127-
"plot!(p, ξ -> f(ξ) / 6; label=\"fit\")"
127+
"p = plot(; xlabel = \"logarithmic correlation length\", ylabel = \"entanglement entropy\")\n",
128+
"p = plot(log.(ξs), Ss; seriestype = :scatter, label = nothing)\n",
129+
"plot!(p, ξ -> f(ξ) / 6; label = \"fit\")"
128130
],
129131
"metadata": {},
130132
"execution_count": null
@@ -145,13 +147,13 @@
145147
"file_extension": ".jl",
146148
"mimetype": "application/julia",
147149
"name": "julia",
148-
"version": "1.10.4"
150+
"version": "1.12.0"
149151
},
150152
"kernelspec": {
151-
"name": "julia-1.10",
152-
"display_name": "Julia 1.10.4",
153+
"name": "julia-1.12",
154+
"display_name": "Julia 1.12.0",
153155
"language": "julia"
154156
}
155157
},
156158
"nbformat": 4
157-
}
159+
}

‎docs/src/examples/quantum1d/1.ising-cft/index.md‎

Lines changed: 87 additions & 288 deletions
Large diffs are not rendered by default.

‎docs/src/examples/quantum1d/1.ising-cft/main.ipynb‎

Lines changed: 24 additions & 21 deletions
Original file line numberDiff line numberDiff line change
@@ -56,12 +56,11 @@
5656
"outputs": [],
5757
"cell_type": "code",
5858
"source": [
59-
"energies, states = exact_diagonalization(H; num=18, alg=Lanczos(; krylovdim=200));\n",
60-
"plot(real.(energies);\n",
61-
" seriestype=:scatter,\n",
62-
" legend=false,\n",
63-
" ylabel=\"energy\",\n",
64-
" xlabel=\"#eigenvalue\")"
59+
"energies, states = exact_diagonalization(H; num = 18, alg = Lanczos(; krylovdim = 200));\n",
60+
"plot(\n",
61+
" real.(energies);\n",
62+
" seriestype = :scatter, legend = false, ylabel = \"energy\", xlabel = \"#eigenvalue\"\n",
63+
")"
6564
],
6665
"metadata": {},
6766
"execution_count": null
@@ -167,11 +166,13 @@
167166
"Δ = real.(energies[1:18] .- energies[1]) ./ (2π * v / L)\n",
168167
"S = momenta ./ (2π / L)\n",
169168
"\n",
170-
"p = plot(S, real.(Δ);\n",
171-
" seriestype=:scatter, xlabel=\"conformal spin (S)\", ylabel=\"scaling dimension (Δ)\",\n",
172-
" legend=false)\n",
173-
"vline!(p, -3:3; color=\"gray\", linestyle=:dash)\n",
174-
"hline!(p, [0, 1 / 8, 1, 9 / 8, 2, 17 / 8]; color=\"gray\", linestyle=:dash)\n",
169+
"p = plot(\n",
170+
" S, real.(Δ);\n",
171+
" seriestype = :scatter, xlabel = \"conformal spin (S)\", ylabel = \"scaling dimension (Δ)\",\n",
172+
" legend = false\n",
173+
")\n",
174+
"vline!(p, -3:3; color = \"gray\", linestyle = :dash)\n",
175+
"hline!(p, [0, 1 / 8, 1, 9 / 8, 2, 17 / 8]; color = \"gray\", linestyle = :dash)\n",
175176
"p"
176177
],
177178
"metadata": {},
@@ -212,7 +213,7 @@
212213
"outputs": [],
213214
"cell_type": "code",
214215
"source": [
215-
"E_ex, qps = excitations(H_mps, QuasiparticleAnsatz(), ψ, envs; num=18)\n",
216+
"E_ex, qps = excitations(H_mps, QuasiparticleAnsatz(), ψ, envs; num = 18)\n",
216217
"states_mps = vcat(ψ, map(qp -> convert(FiniteMPS, qp), qps))\n",
217218
"energies_mps = map(x -> expectation_value(x, H_mps), states_mps)\n",
218219
"\n",
@@ -231,11 +232,13 @@
231232
"Δ_mps = real.(energies_mps[1:18] .- energies_mps[1]) ./ (2π * v / L_mps)\n",
232233
"S_mps = momenta_mps ./ (2π / L_mps)\n",
233234
"\n",
234-
"p_mps = plot(S_mps, real.(Δ_mps);\n",
235-
" seriestype=:scatter, xlabel=\"conformal spin (S)\",\n",
236-
" ylabel=\"scaling dimension (Δ)\", legend=false)\n",
237-
"vline!(p_mps, -3:3; color=\"gray\", linestyle=:dash)\n",
238-
"hline!(p_mps, [0, 1 / 8, 1, 9 / 8, 2, 17 / 8]; color=\"gray\", linestyle=:dash)\n",
235+
"p_mps = plot(\n",
236+
" S_mps, real.(Δ_mps);\n",
237+
" seriestype = :scatter, xlabel = \"conformal spin (S)\",\n",
238+
" ylabel = \"scaling dimension (Δ)\", legend = false\n",
239+
")\n",
240+
"vline!(p_mps, -3:3; color = \"gray\", linestyle = :dash)\n",
241+
"hline!(p_mps, [0, 1 / 8, 1, 9 / 8, 2, 17 / 8]; color = \"gray\", linestyle = :dash)\n",
239242
"p_mps"
240243
],
241244
"metadata": {},
@@ -257,13 +260,13 @@
257260
"file_extension": ".jl",
258261
"mimetype": "application/julia",
259262
"name": "julia",
260-
"version": "1.10.4"
263+
"version": "1.12.0"
261264
},
262265
"kernelspec": {
263-
"name": "julia-1.10",
264-
"display_name": "Julia 1.10.4",
266+
"name": "julia-1.12",
267+
"display_name": "Julia 1.12.0",
265268
"language": "julia"
266269
}
267270
},
268271
"nbformat": 4
269-
}
272+
}

0 commit comments

Comments
 (0)