accelerInt
v0.1
|
File containing functions for best rational approximation to matrix exponential. More...
#include <stdlib.h>
#include <math.h>
#include <string.h>
#include "solver_options.h"
#include "linear-algebra.h"
#include <complex.h>
#include <fftw3.h>
Go to the source code of this file.
Macros | |
#define | M_PI 4 * atan(1) |
Functions | |
void | cf (int n, double *poles_r, double *poles_i, double *res_r, double *res_i) |
Function that calculates the poles and residuals of best rational (partial fraction) approximant to the matrix exponential. More... | |
File containing functions for best rational approximation to matrix exponential.
Contains main and linear algebra functions for Carathéodory-Fejér method of Rational Approximants
Definition in file cf.c.
void cf | ( | int | n, |
double * | poles_r, | ||
double * | poles_i, | ||
double * | res_r, | ||
double * | res_i | ||
) |
Function that calculates the poles and residuals of best rational (partial fraction) approximant to the matrix exponential.
Complex math Fast Fourier tranform functionsUses the Carathéodory-Fejér method; based on the MATLAB code in L.N. Trefethen, J.A.C. Weideman, T. Schmelzer, "Talbot quadratures and rational approximations," BIT Numer. Math. 46 (2006) 653–670.
[in] | n | size of approximation (n, n) |
[out] | poles_r | array with real parts of poles, size n |
[out] | poles_i | array with imaginary parts of poles, size n |
[out] | res_r | array with real parts of residuals, size n |
[out] | res_i | array with imaginary parts of residuals, size n |