Sparse support for AbstractQuMatrix.

Author: amitjamadagni

src/arrays/arraymath.jl

@@ -153,15 +153,23 @@ normalize(qarr::AbstractQuArray) = normalize!(copy(qarr))
 # matrix operations returning an array
 # sparse to dense
 function Base.full(qarr::AbstractQuMatrix)
-    fc = full(rawcoeffs(qarr))
+    fc = full(coeffs(qarr))
     QAT = similar_type(qarr)
     return QAT(fc, bases(qarr))
 end
+
+# convert to a sparse matrix
+function Base.sparse(qarr::AbstractQuMatrix)
+    sc = sparse(coeffs(qarr))
+    QAT = similar_type(qarr)
+    return QAT(sc, bases(qarr))
+end
+
 #Base.full(ct::CTranspose) = full(ct.qarr)'
 
 # exponential of dense matrix
 function Base.expm(qarr::AbstractQuMatrix)
-    fc = expm(full(rawcoeffs(qarr)))
+    fc = expm(full(coeffs(qarr)))
     QAT = similar_type(qarr)
     return QAT(fc, bases(qarr))
 end

test/multest.jl

@@ -55,3 +55,6 @@ qv1 = normalize!(QuArray(v1))
 
 # Vectorize
 @assert vec(qv) == vec(qv')
+
+# Sparse matrix
+@assert sparse(sigmax)*sparse(sigmax) == sparse(QuArray(eye(2)))