Testing Piessens algorithm implemented in ply_legser_module with the data provided in his description in [1] Communications of the ACM, January 1974, Volume 17, Number 1, page 25. Used is the expansion of the function f(x) = 1 / (2-x).
Type | Attributes | Name | Initial | |||
---|---|---|---|---|---|---|
real(kind=rk) | :: | chebymodes(0:20) | ||||
real(kind=rk) | :: | legmodes(0:20) | ||||
real(kind=rk) | :: | refleg(0:20) | ||||
real(kind=rk), | parameter | :: | rt3fourth | = | sqrt(0.75_rk) | |
real(kind=rk) | :: | cterm | ||||
real(kind=rk) | :: | ci | ||||
real(kind=rk) | :: | maxdiff | ||||
integer | :: | i |