Documentation for module minimax_rpa

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.F
Loading...

public Subroutines/Functions:

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.

SUBROUTINEget_rpa_minimax_coeff(k, e_range, aw, ierr, print_warning)

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.

Arguments:
INTEGER
:: k ...
REAL(dp)
:: e_range ...
REAL(dp),
ALLOCATABLE
:: aw(:) ...
INTEGER
:: ierr ...
LOGICAL,
OPTIONAL
:: print_warning ...