Re: [eigen] Matrix exponential |

[ Thread Index | Date Index | More lists.tuxfamily.org/eigen Archives ]

*To*: eigen@xxxxxxxxxxxxxxxxxxx*Subject*: Re: [eigen] Matrix exponential*From*: Jitse Niesen <jitse@xxxxxxxxxxxxxxxxx>*Date*: Wed, 6 May 2009 15:21:48 +0100 (BST)

On Wed, 6 May 2009, Benoit Jacob wrote:

The complexity of the present approach is said to be roughly 20n^3 which is huge, though, and in practice, the most common case is exponentiating a skew-adjoint matrix (think exp(itH) where H is a hamiltonian, think Stone's theorem...), whose diagonalization boils down to diagonalization of a selfadjoint matrix. So it'd be very important to make sure the selfadjointeigensolver is robust enough that it can always safely be used without convergence issues!

I don't know whether complex numbers make a difference here. Jitse

**Follow-Ups**:**Re: [eigen] Matrix exponential***From:*Hauke Heibel

**References**:**Re: [eigen] Matrix exponential***From:*Jitse Niesen

**Re: [eigen] Matrix exponential***From:*Benoit Jacob

**Re: [eigen] Matrix exponential***From:*Thomas Capricelli

**Re: [eigen] Matrix exponential***From:*Benoit Jacob

**Messages sorted by:**[ date | thread ]- Prev by Date:
**Re: [eigen] [patch] Missing function in StdVector** - Next by Date:
**Re: [eigen] [patch] Missing function in StdVector** - Previous by thread:
**Re: [eigen] Matrix exponential** - Next by thread:
**Re: [eigen] Matrix exponential**

Mail converted by MHonArc 2.6.19+ | http://listengine.tuxfamily.org/ |