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