Sparse support for AbstractQuMatrix.

JuliaAttic · Oct 2, 2015 · 9f77b9c · 9f77b9c
1 parent c8fdb89
commit 9f77b9c
Showing 1 changed file with 10 additions and 2 deletions.
diff --git a/src/arrays/arraymath.jl b/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