ekmett / ekmett/integration

The two last results are always identical

Open
#2 1 comment 0 reactions 0 assignees View on GitHub
Dominant language
Haskell
Stars
12
Forks
2
PR merge metrics
No merged PRs in 30d

Description

Hello,
Using `everywhere` or `nonNegative`, I've observed that the two last results are always identical. That sounds strange.
```
*Gauss> everywhere trap (\x -> exp(-x*x))
[Result {result = 1.7576521878389109, errorEstimate = 0.25846444760875, evaluations = 13},Result {result = 1.772497391257178, errorEstimate = 1.4845203418266999e-2, evaluations = 25},Result {result = 1.7724538509479826, errorEstimate = 4.354030919522532e-5, evaluations = 49},Result {result = 1.7724538509055163, errorEstimate = 4.246636850407919e-11, evaluations = 97},Result {result = 1.7724538509055163, errorEstimate = 1.743934249004316e-16, evaluations = 193},Result {result = 1.7724538509055159, errorEstimate = 3.487868498008632e-16, evaluations = 385},Result {result = 1.7724538509055159, errorEstimate = 3.487868498008632e-16, evaluations = 769}]
*Gauss> nonNegative trap (\x -> exp(-x*x))
[Result {result = 0.8857864863947141, errorEstimate = 8.891160064614978e-2, evaluations = 13},Result {result = 0.8862152677982271, errorEstimate = 4.287814035129345e-4, evaluations = 25},Result {result = 0.8862269252508314, errorEstimate = 1.1657452604230123e-5, evaluations = 49},Result {result = 0.8862269254527348, errorEstimate = 2.0190354943494526e-10, evaluations = 97},Result {result = 0.8862269254527361, errorEstimate = 1.2212453270876722e-15, evaluations = 193},Result {result = 0.8862269254527364, errorEstimate = 3.885780586188048e-16, evaluations = 385},Result {result = 0.8862269254527364, errorEstimate = 3.885780586188048e-16, evaluations = 769}]
*Gauss> nonNegative trap (\x -> exp(-x))
[Result {result = 0.9997341755432854, errorEstimate = 2.1288525221500487e-2, evaluations = 13},Result {result = 0.9999999702600391, errorEstimate = 2.657947167536978e-4, evaluations = 25},Result {result = 0.9999999999999781, errorEstimate = 2.9739939066963927e-8, evaluations = 49},Result {result = 0.9999999999999769, errorEstimate = 1.2212453270876722e-15, evaluations = 97},Result {result = 0.9999999999999779, errorEstimate = 9.992007221626409e-16, evaluations = 193},Result {result = 0.999999999999978, errorEstimate = 1.6653345369377348e-16, evaluations = 385},Result {result = 0.999999999999978, errorEstimate = 1.6653345369377348e-16, evaluations = 769}]
*Gauss> nonNegative trap (\x -> 1/(x+1))
[Result {result = 31.62342541843333, errorEstimate = 0.47682097989158123, evaluations = 13},Result {result = 31.506682390464704, errorEstimate = 0.11674302796862968, evaluations = 25},Result {result = 31.480003758402475, errorEstimate = 2.6678632062230534e-2, evaluations = 49},Result {result = 31.473026453375567, errorEstimate = 6.977305026909519e-3, evaluations = 97},Result {result = 31.47126525542729, errorEstimate = 1.7611979482783369e-3, evaluations = 193},Result {result = 31.47091823638368, errorEstimate = 3.4701904360900926e-4, evaluations = 385},Result {result = 31.47091823638368, errorEstimate = 3.4701904360900926e-4, evaluations = 769}]
```

Contributor guide

No contributing guide indexed for this repository

Assessment

This issue has not been assessed yet.

Get new issues in your inbox

A short digest of beginner-friendly GitHub issues.