exponential matrix exp(A)
What is the maximum size of a square matrix A to calculate exp(A)
asked 2012-03-08 09:51:07 +0100
This post is a wiki. Anyone with karma >750 is welcome to improve it.
What is the maximum size of a square matrix A to calculate exp(A)
answered 2012-03-08 11:58:25 +0100
This post is a wiki. Anyone with karma >750 is welcome to improve it.
I don't think there are any hard-coded size limits.
answered 2012-03-08 13:30:35 +0100
This post is a wiki. Anyone with karma >750 is welcome to improve it.
At the same time, experimentation with random_matrix(ZZ,i,i)
suggests that Maxima hits a wall with this already when i=4
. RR
doesn't look much better, and
sage: m = random_matrix(RDF,4,4)
sage: m.exp?
basically says as much. But
sage: m = random_matrix(RDF,4,4)
sage: exp(m)
[ 0.360479628838 -0.145794589097 1.4781378783 0.316593839994]
[ 0.317770329736 0.0984841855079 -0.733966866354 -0.820812967827]
[ -1.58402218787 0.783288361066 0.965206766202 0.919221986577]
[-0.254553742946 0.837705408602 -2.47262343564 1.63085510321]
sage: m = random_matrix(RDF,10,10)
sage: exp(m)
looks much better. So the answer to your question depends on what you are exactly looking for.
Indeed, using RDF or CDF will probably scale *much* better.
Please start posting anonymously - your entry will be published after you log in or create a new account.
Asked: 2012-03-08 09:51:07 +0100
Seen: 2,382 times
Last updated: Mar 08 '12
Is it possible to find the exponential of a symbolic matrix using sage?
how to obtain the real part of a complex function
evaluating the $U^\dagger U$ of an unitary matrix does not work
how to get the diagonal of a matrix?
complex exponential/trigonometric