Routines to calculate the minimax coefficients for approximating 1/x as 1/x ~ 1/pi SUM_{i}^{K} w_i x^2/(x^2+a_i^2)^2 for x belonging to [1:Rc].

source: minimax_rpa.FLoading...

The a_i and w_i coefficient are stored in aw such that the first 1:K elements correspond to a_i and the K+1:2k correspond to w_i. All a_i and w_i have been fitted with a 12th order polynomial as a function of Rc for each value of K.

The a_i and w_i coefficient are stored in aw such that the first 1:K elements correspond to a_i and the K+1:2k correspond to w_i. 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 |
:: | k | ... | ||

REAL(dp) |
:: | e_range | ... | ||

REAL(dp), |
ALLOCATABLE |
:: | aw(:) | ... | |

INTEGER |
:: | ierr | ... | ||

LOGICAL, |
OPTIONAL |
:: | print_warning | ... |