(PC-12xx~14xx) qthsh Tanh-Sinh quadrature

[emece67]

Hi,

These are the results for my Python script:

Code:

   TNFE err_r err_t
 1 0047 8e-12 0e+00
 2 0049 5e-12 2e-15
 3 0049 1e-08 8e-15
 4 0049 2e-13 2e-15
 5 0049 7e-12 3e-15
 6 0049 5e-08 6e-08
 7 0049 2e-12 2e-12
 8 0025 5e-07 2e-08
 9 0025 5e-07 2e-08
10 0197 6e-04 6e-03
11 0197 6e-04 5e-03
12 0049 8e-08 1e-07
13 0197 5e-04 4e-03
14 0025 4e-07 2e-08
15 0025 4e-07 2e-08
16 0137 2e-15 0e+00
17 0049 7e-10 1e-11
18 0049 2e-12 3e-15
19 0197 3e-03 2e-03
20 0197 1e-02 1e-02
21 0197 1e-02 1e-02

where err_r is the reported error and err_t is the true error. TNFE is the total number of function evaluations. The only adjustable parameter is mp.dps (set to 15, seems to work, at least, with values 8~128). It worked as (I) expected, except for functions 10~11... they hurt!

The script is suffering now from some refactoring, so I expect to upload it to GitHub in a few days.

Regards.


EDIT: results modified as a result of code refactoring.