Error propagation in adaptive Simpson algorithm
07-31-2018, 07:53 PM
Post: #6
 Albert Chan Senior Member Posts: 2,566 Joined: Jul 2018
RE: Error propagation in adaptive Simpson algorithm
Hi, Claudio

The problem with the sin integral is not cancellation, but its periodic nature.
Simpson's Rule is also based on periodic sampling (all equally spaced).

If the sampling and periodic function were in sync, the samples are going to be biased.
For sin(x) 0 to 200, the first few iterations, none of sin(x) samples were above zero.

A non-linear transformed function can fix this, by scrambling the sample points.
Have you tried the non-linear transformed sin integral (same thread, post #13) ?

---

Regarding "excessive" accuracy, it will happen even if tolerance is good estimater for accuracy.
Not all integral have the same convergence rate.

Say, tolerance of 1e-8 somehow guaranteed 7 digits accuracy (it does not)
What is the chance of really getting 7 digits accuracy ? Almost zero.

With guaranteed minimum accuracy, average accuracy is going to be higher, say, 10 digits.
Iterations that not quite make it to tolerance will doubled the points, "wasting" accuracy.
 « Next Oldest | Next Newest »

 Messages In This Thread Error propagation in adaptive Simpson algorithm - Claudio L. - 07-30-2018, 06:42 PM RE: Error propagation in adaptive Simpson algorithm - Dieter - 07-31-2018, 09:35 AM RE: Error propagation in adaptive Simpson algorithm - Albert Chan - 07-31-2018, 11:42 AM RE: Error propagation in adaptive Simpson algorithm - Claudio L. - 07-31-2018, 05:02 PM RE: Error propagation in adaptive Simpson algorithm - Claudio L. - 07-31-2018, 06:37 PM RE: Error propagation in adaptive Simpson algorithm - Albert Chan - 07-31-2018 07:53 PM RE: Error propagation in adaptive Simpson algorithm - Claudio L. - 07-31-2018, 08:40 PM RE: Error propagation in adaptive Simpson algorithm - Vtile - 07-31-2018, 09:06 PM RE: Error propagation in adaptive Simpson algorithm - Dieter - 08-01-2018, 08:02 AM RE: Error propagation in adaptive Simpson algorithm - Albert Chan - 08-01-2018, 02:23 PM RE: Error propagation in adaptive Simpson algorithm - Claudio L. - 08-01-2018, 04:28 PM RE: Error propagation in adaptive Simpson algorithm - Albert Chan - 08-01-2018, 05:53 PM

User(s) browsing this thread: 1 Guest(s)