Routines to calculate the minimax coefficients in order to approximate 1/x as a sum over exponential functions 1/x ~ SUM_{i}^{K} w_i EXP(-a_i * x) for x belonging to [1:Rc]. This module contains coefficients for minimax approximations with 1 <= k <= 15.

source: minimax_exp_k15.FLoading...

Check that the range for the minimax approximation is not too small for the chosen number of integration point, or that the number of integration point is not exceeding the maximum allowed by the actual implementation.

Get minimax coefficients: k15 implementation (coefficients up to k=15 terms). All a_i and w_i have been fitted with a 12th order polynomial as a function of Rc for each value of K.

Sample numerical error and return its maximum.

Check that the range for the minimax approximation is not too small for the chosen number of integration point, or that the number of integration point is not exceeding the maximum allowed by the actual implementation.

INTEGER, |
INTENT(in) |
:: | k | ... | |

REAL(dp), |
INTENT(in) |
:: | rc | ... | |

INTEGER |
:: | ierr | ... |

Get minimax coefficients: k15 implementation (coefficients up to k=15 terms). All a_i and w_i have been fitted with a 12th order polynomial as a function of Rc for each value of K.

INTEGER, |
INTENT(in) |
:: | k | ... | |

REAL(dp), |
INTENT(in) |
:: | rc | ... | |

REAL(dp), |
INTENT(out) |
:: | aw(2*k) | ... | |

REAL(dp), |
INTENT(out), |
OPTIONAL |
:: | mm_error | ... |

Sample numerical error and return its maximum.

Return Value :: REAL(dp)

REAL(dp), |
INTENT(in) |
:: | rc | ... | |

REAL(dp), |
INTENT(in) |
:: | aw(:) | ... |