Tests added for displaceop.

JuliaAttic · Apr 19, 2015 · be7e3ce · be7e3ce
1 parent de8d52b
commit be7e3ce
Show file tree

Hide file tree

Showing 2 changed files with 13 additions and 7 deletions.
diff --git a/src/arrays/ladderops.jl b/src/arrays/ladderops.jl
@@ -41,7 +41,7 @@ squeeze_construct(a::AbstractQuMatrix,b::AbstractQuMatrix, z::Number) = scale!(0
 squeezingop(a::AbstractQuMatrix, b::AbstractQuMatrix, z::Number) = expm(squeeze_construct(a,b,z))
 squeezingop(n::Int,z::Number) = squeezingop(lowerop(n),lowerop(n),z)
 
-displaceop(n::Int,alpha::Number) = expm(scale!(alpha,lowerop(n)')-scale!(alpha',lowerop(n)))
+displaceop(n::Int,alpha::Number) = expm(scale(alpha,lowerop(n)')-scale(alpha',lowerop(n)))
 
 export raiseop,
     lowerop,

diff --git a/test/operatortest.jl b/test/operatortest.jl
@@ -10,19 +10,25 @@
 
 @assert coeffs(commutator(sigma_x, sigma_x)) == spzeros(2,2)
 
-######################################
-# Position & Momentum Operators Test #
-######################################
+####################################################
+# Position, Displacement & Momentum Operators Test #
+####################################################
 
 p = positionop(2)
 m = momentumop(2)
 @assert coeffs(commutator(sigma_x, p)) == spzeros(2,2)
 @assert coeffs(commutator(sigma_y, m)) == spzeros(2,2)
 
-#####################################
-# Squeeze, Squeezing Operators Test #
-#####################################
+sv1 = [0.5403023058681398,0.8414709848078965]
+sv2 = [0.5403023058681398 + 0.0im,0.0 + 0.8414709848078965im]
+@assert displaceop(2,1)*statevec(1,FiniteBasis(2)) == QuArray(sv1)
+@assert displaceop(2,1*im)*statevec(1,FiniteBasis(2)) == QuArray(sv2)
+
+############################
+# Squeezing Operators Test #
+############################
 
 @assert squeezingop(2,1.0)== QuArray(eye(2))
 @assert squeezingop(lowerop(2), lowerop(2), 1.0) == squeezingop(2,1)
 
+