From 8bc26f1150b47f472bc7846fef83ceb7334c6ac1 Mon Sep 17 00:00:00 2001
From: tl <thibaut.lacroix@gmail.com>
Date: Thu, 11 Apr 2024 11:15:00 +0200
Subject: [PATCH] Add a couple of comment lines

---
 examples/sbm_Htimedependent.jl | 17 +++++++++--------
 1 file changed, 9 insertions(+), 8 deletions(-)

diff --git a/examples/sbm_Htimedependent.jl b/examples/sbm_Htimedependent.jl
index 29a3267..8d7896c 100644
--- a/examples/sbm_Htimedependent.jl
+++ b/examples/sbm_Htimedependent.jl
@@ -35,6 +35,7 @@ Trabi = 30.0 # Rabi period of the drive
 ωdrive = ω0 # Frequency of the drive
 
 Ndrive = 1 #Number of the site on which the drive is applied
+
 #-----------------------
 # Simulation parameters
 #-----------------------
@@ -51,14 +52,14 @@ D = [6] # MPS bond dimension
 # MPO and initial state MPS
 #---------------------------
 
-numsteps = length(collect(0:dt:tfinal))-1
-timelist = [(i-1)*dt for i=1:numsteps+1]
+timelist = collect(0:dt:tfinal)
+numsteps = length(timelist)-1
 
-Ht = [ϵ*sx*sin(ωdrive*tstep) for tstep in timelist] # Time dependent Hamiltonian term MPO
+Ht = [ϵ*sx*sin(ωdrive*tstep) for tstep in timelist] # Time-dependent Hamiltonian term
 
 H = spinbosonmpo(ω0, Δ, d, N, cpars) # MPO representation of the Hamiltonian
 
-ψ = unitcol(2,2) # Initial low-z system state 
+ψ = unitcol(2,2) # Initial down-z system state 
 
 A = productstatemps(physdims(H), state=[ψ, fill(unitcol(1,d), N)...]) # MPS representation of |ψ>|Vacuum>
 
@@ -83,9 +84,9 @@ A, dat = runsim(dt, tfinal, A, H;
                 convobs = [ob1],
                 params = @LogParams(N, d, α, Δ, ω0, s),
                 convparams = D,
-                timedep = true,
-                Ndrive = Ndrive,
-                Htime = Ht,
+                timedep = true, # the Hamiltonian is time dependent
+                Ndrive = Ndrive, # the first site of the MPS/MPO (i.e. the system) is concerned
+                Htime = Ht, # list of time-dependent terms
                 verbose = false,
                 save = true,
                 plot = true,
@@ -95,4 +96,4 @@ A, dat = runsim(dt, tfinal, A, H;
 # Plots
 #----------
 
-plot(dat["data/times"], dat["data/sz"], label=["Dmax = $D"], xlabel=L"t",ylabel=L"\sigma_z", title="")
\ No newline at end of file
+plot(dat["data/times"], dat["data/sz"], label="Dmax = $(D...)", xlabel=L"t",ylabel=L"\sigma_z", title="")