update spin tests

qojulia · Mar 14, 2024 · ea4015b · ea4015b
1 parent 77e176b
commit ea4015b
Showing 1 changed file with 17 additions and 13 deletions.
diff --git a/test/test_spin.jl b/test/test_spin.jl
@@ -186,19 +186,23 @@ sol_cs1 = solve(prob_cs1,Tsit5(),abstol=1e-8,reltol=1e-8)
 # sol_cs1[Sz(1)] # TODO: error
 
 
-# ### collective spin definition (factor 1/2 in sum of atoms?)
-# Ncs1_ = 200
-# using QuantumOptics
-# b = SpinBasis(Ncs1_/2)
-# ψ0 = spinup(b)
-# expect(sigmax(b)/2, ψ0)
-# expect(sigmay(b)/2, ψ0)
-# expect(sigmaz(b)/2, ψ0)
-# expect(sigmax(b)/2*sigmax(b)/2, ψ0)
-# expect(sigmay(b)/2*sigmay(b)/2, ψ0)
-# expect(sigmax(b)/2*sigmay(b)/2, ψ0)
-# expect(sigmax(b)/2*sigmaz(b)/2, ψ0)
-# expect(sigmaz(b)/2*sigmaz(b)/2, ψ0)
+### collective spin definition (factor 1/2 in sum of atoms?)
+Ncs1_ = 200
+using QuantumOptics
+b = SpinBasis(Ncs1_/2)
+ψ0 = spinup(b)
+expect(sigmax(b)/2, ψ0)
+expect(sigmay(b)/2, ψ0)
+expect(sigmaz(b)/2, ψ0)
+expect(sigmax(b)/2*sigmax(b)/2, ψ0)
+expect(sigmay(b)/2*sigmay(b)/2, ψ0)
+expect(sigmax(b)/2*sigmay(b)/2, ψ0)
+expect(sigmax(b)/2*sigmaz(b)/2, ψ0)
+expect(sigmaz(b)/2*sigmaz(b)/2, ψ0)
+expect(sigmaz(b)/2*sigmaz(b)/2*sigmaz(b)/2, ψ0)
+expect(sigmax(b)/2*sigmax(b)/2*sigmaz(b)/2, ψ0)
+expect(sigmax(b)/2*sigmax(b)/2*sigmax(b)/2, ψ0)
+
 
 ### TODO: 
 # time evolution: simple example, superradiant decay, farokh paper,