Ticket #12101: trac_12101-sparse-matrix-exp.patch

sage/matrix/matrix2.pyx

@@ -11527 +11527 @@
             sage: a.change_ring(RDF).exp()
             [42748127.3153 7368259.24416]
             [234538976.138 40426191.4516]
+
+        This works for sparse matrices, too; see :trac:`12101`::
+
+            sage: D = matrix(SR,[1],sparse=True)
+            sage: D.exp()
+            [e]
+            sage: D
+            [1]
+            sage: exp(diagonal_matrix([1, 2, 3]))
+            [  e   0   0]
+            [  0 e^2   0]
+            [  0   0 e^3]
         """
         from sage.symbolic.ring import SR
-        return self.change_ring(SR).exp()
+        if not self.is_sparse():
+            return self.change_ring(SR).exp()
+        return self.change_ring(SR).dense_matrix().exp()
 
     def elementary_divisors(self):
         r"""