![]() |
Eigen-unsupported
5.0.1-dev+284dcc12
|
#include <unsupported/Eigen/src/MatrixFunctions/MatrixExponential.h>
Compute the (17,17)-Padé approximant to the exponential.
After exit, \( (V+U)(V-U)^{-1} \) is the Padé approximant of \( \exp(A) \) around \( A = 0 \).
This function activates only if your long double is double-double or quadruple.
Static Public Member Functions | |
| static void | run (const MatrixType &arg, MatrixType &U, MatrixType &V, int &squarings) |
| Compute Padé approximant to the exponential. | |
|
static |
Compute Padé approximant to the exponential.
Computes U, V and squarings such that \( (V+U)(V-U)^{-1} \) is a Padé approximant of \( \exp(2^{-\mbox{squarings}}M) \) around \( M = 0 \), where \( M \) denotes the matrix arg. The degree of the Padé approximant and the value of squarings are chosen such that the approximation error is no more than the round-off error.