[sage-devel] Re: Matrix exponentials of diagonal matrices.

2019-04-16 Thread Nils Bruin
On Tuesday, April 16, 2019 at 8:35:27 PM UTC-7, John H Palmieri wrote: > > > I think it's any sparse matrix, not just sparse matrices over SR. > > sage: m = matrix(3, [1, 0, 0, 0, 2, 0, 0, 0, 3], sparse=True) > sage: exp(m) # hits the bug > > Indeed. That's because m.exp() first changes the ring

[sage-devel] Re: Matrix exponentials of diagonal matrices.

2019-04-16 Thread John H Palmieri
On Tuesday, April 16, 2019 at 8:30:13 PM UTC-7, Nils Bruin wrote: > > On Tuesday, April 16, 2019 at 8:06:06 PM UTC-7, Brian Fitzpatrick wrote: >> >> >> Is this a bug? >> >> Looks like it. The problem seems to be that "sparse matrices over SR" > aren't special cased the way dense matrices are

[sage-devel] Re: Matrix exponentials of diagonal matrices.

2019-04-16 Thread Nils Bruin
On Tuesday, April 16, 2019 at 8:06:06 PM UTC-7, Brian Fitzpatrick wrote: > > > Is this a bug? > > Looks like it. The problem seems to be that "sparse matrices over SR" aren't special cased the way dense matrices are (and a diagonal matrix is constructed as sparse: sage: D=diagonal_matrix([1, 2,