## Revision 1110 trunk/src/df1b2-separable/expm.cpp

expm.cpp (revision 1110)
1 1
/*
2 2
* $Id$
3 3
*
4
* Authors: Anders Nielsen <anders@nielsensweb.org> and Casper W. Berg <cbe@aqua.dtu.dk>
5
* Copyright (c) 2008-2012 Regents of the University of California
4
* Authors: Anders Nielsen <anders@nielsensweb.org> and Casper W. Berg <cbe@aqua.dtu.dk>
5
* Copyright (c) 2008-2012 Regents of the University of California
6 6
*/
7 7
/**
8 8
* \file
9
* Matrix exponential function for df1b2_matrix.
9
* Matrix exponential function for df1b2_matrix.
10 10
*/
11 11

12 12
#include <df1b2fun.h>
......
195 195

196 196
/**
197 197
\ingroup matop
198
Matrix exponential.
199

198
Matrix exponential.
199

200 200
The matrix exponential is calculated using the Pade approximation adapted from Moler, Cleve; Van Loan, Charles F. (2003), "Nineteen Dubious Ways to Compute the Exponential of a Matrix, Twenty-Five Years Later"
201 201

202

203
The main use of the matrix exponential is to solve linear ordinary differential equation (ODE) systems:
202
The main use of the matrix exponential is to solve linear ordinary differential equation (ODE) systems:
204 203
\f[
205 204
\frac{d}{dt}y(t) = Ay(t)\ , \ \mbox{with } y(0) = y_0
206
\f]
205
\f]
207 206
\item then the solution becomes
208 207
\f[
209 208
y(t) = e^{At}y_0

Also available in: Unified diff