Accuracy of Integral with epsilon

10082021, 09:57 PM
Post: #1




Accuracy of Integral with epsilon
I was surprised that my TI30X Pro MathPrint came up with a more accurate solution to the integral (2*e^x)/(cosh (x) + sinh (x)) from 10 to 10 using an epsilon of 10^4 or greater. I was naively thinking that the stock epsilon of 10^5 or smaller would yield a more accurate result. The actual answer is 40. My TI cx ii was accurate to 5 digits as is the TI30X Pro MathPrint and CASIO fx991X. The CASIO fxCG50 is accurate to 6 digits, the HP Prime in home mode thought about it a bit, and was accurate to 7 digits. The WP 34S on my SwissMicros in double precision quickly produced an answer accurate to 27 digits!


10082021, 10:17 PM
(This post was last modified: 10092021 12:51 AM by toml_12953.)
Post: #2




RE: Accuracy of Integral with epsilon
(10082021 09:57 PM)lrdheat Wrote: I was surprised that my TI30X Pro MathPrint came up with a more accurate solution to the integral (2*e^x)/(cosh (x) + sinh (x)) from 10 to 10 using an epsilon of 10^4 or greater. I was naively thinking that the stock epsilon of 10^5 or smaller would yield a more accurate result. The actual answer is 40. My TI cx ii was accurate to 5 digits as is the TI30X Pro MathPrint and CASIO fx991X. The CASIO fxCG50 is accurate to 6 digits, the HP Prime in home mode thought about it a bit, and was accurate to 7 digits. The WP 34S on my SwissMicros in double precision quickly produced an answer accurate to 27 digits! I don't understand what you mean about Prime being accurate to 7 digits. I get 40 exactly on Prime in home mode. Tom L Cui bono? 

10082021, 10:18 PM
Post: #3




RE: Accuracy of Integral with epsilon
I don't even understand what's supposed to be challenging about that integral, since that function should be the constant 2...


10082021, 10:28 PM
Post: #4




RE: Accuracy of Integral with epsilon
(10082021 10:18 PM)Thomas Okken Wrote: I don't even understand what's supposed to be challenging about that integral, since that function should be the constant 2... Maybe he was hoping we wouldn't notice. Tom L Cui bono? 

10082021, 11:28 PM
Post: #5




RE: Accuracy of Integral with epsilon
Did not think the integral was remarkable, but was fascinated by how different calculators handled it as written. My Prime in home shows 40.0000000143 (!). Impressed with the WP 34S.


10092021, 12:47 AM
Post: #6




RE: Accuracy of Integral with epsilon
(10082021 11:28 PM)lrdheat Wrote: Did not think the integral was remarkable, but was fascinated by how different calculators handled it as written. My Prime in home shows 40.0000000143 (!). Impressed with the WP 34S. What firmware version are you using? I use the 2021 10 01 version on a G2. I wonder why the big difference between yours and mine. Tom L Cui bono? 

10092021, 01:41 AM
Post: #7




RE: Accuracy of Integral with epsilon
(10082021 11:28 PM)lrdheat Wrote: Did not think the integral was remarkable, but was fascinated by how different calculators handled it as written. My Prime in home shows 40.0000000143 (!). Impressed with the WP 34S. This integral is efficiently integrated with the TanhSinh method (HP forum post). This method is particularly suitable for your function. I did a quick check with my implementation: for a specified accuracy of eps=10^4 TanhSinh only performs 27 integrand evaluations to converge to 40 exactly, or at least very close (depends on the floating point format, single/double IEEE 754 or BCD etc.) The WP 34S uses TanhSinh. That implementation can be improved though. TanhSinh is not considered a general integration method as it may perform quite poorly in certain (known) cases.  Rob "I count on old friends to remain rational" 

10092021, 02:22 AM
Post: #8




RE: Accuracy of Integral with epsilon
(10092021 01:41 AM)robve Wrote:(10082021 11:28 PM)lrdheat Wrote: Did not think the integral was remarkable, but was fascinated by how different calculators handled it as written. My Prime in home shows 40.0000000143 (!). Impressed with the WP 34S. Someone here is not paying attention. Are you saying that a complex stateoftheart TanhSinh method requires 27 evaluations of the integrand (the constant 2) to integrate the constant 2 from 10 to 10 ? Really !? What a lousy implementation ! V. All My Articles & other Materials here: Valentin Albillo's HP Collection 

10092021, 02:27 AM
Post: #9




RE: Accuracy of Integral with epsilon
(10092021 02:22 AM)Valentin Albillo Wrote: Are you saying that a complex stateoftheart TanhSinh method requires 27 evaluations of the integrand (the constant 2) to integrate the constant 2 from 10 to 10 ? Really !? It's neither complex nor state of the art. How many function evaluations should be performed? Clearly one evaluation is sufficient to get the correct answer here. Two might be prudent to verify that it is a constant function. Three perhaps? Pauli 

10092021, 03:13 AM
Post: #10




RE: Accuracy of Integral with epsilon
I’m on 2.1.14588
The Prime does not seem to directly know how to recognize that the integrand is a constant, even when simplify is applied to the integrand. The WP 34S was quick/accurate, regardless of number of iterations. The MathPrint and 991X took ~15 seconds. I messed with a few home settings, but still came up with the 40.0000000143 answer on the Prime G2 

10092021, 05:16 AM
Post: #11




RE: Accuracy of Integral with epsilon
What surprises me most is my CASIO Classpad readily simplifies the integrand to a constant (2), yet, when I integrate the integral with the unsimplified integrand, it fails to recognize the constant integrand, and produces 40.00000008


10092021, 05:44 AM
Post: #12




RE: Accuracy of Integral with epsilon
(10092021 03:13 AM)lrdheat Wrote: I’m on 2.1.14588 You should get the latest firmware, then: 2.1.14596 (2021 10 01) It looks like that has been fixed. Tom L Cui bono? 

10092021, 08:47 AM
Post: #13




RE: Accuracy of Integral with epsilon
(10092021 02:27 AM)Paul Dale Wrote: Clearly one evaluation is sufficient to get the correct answer here. Two might be prudent to verify that it is a constant function. Three perhaps? It seems what's really being tested here is how the integrator handles the slight noisiness of the function at the low end of the range, when sinh(x) and cosh(x) become large but with opposite signs, so numerical cancellation becomes an issue. 

10092021, 11:26 AM
Post: #14




RE: Accuracy of Integral with epsilon
(10092021 03:13 AM)lrdheat Wrote: I’m on 2.1.14588 What function are you using to calculate the integral? I use int. Tom L Cui bono? 

10092021, 02:01 PM
Post: #15




RE: Accuracy of Integral with epsilon
(10092021 01:41 AM)robve Wrote:(10082021 11:28 PM)lrdheat Wrote: Did not think the integral was remarkable, but was fascinated by how different calculators handled it as written. My Prime in home shows 40.0000000143 (!). Impressed with the WP 34S. WP 34S results with a IEEE 754 based implementation matching the Python code given by the author: eps=10^4: abs error = 0.0051939; est. error = 0.0155524; evals = 7 eps=10^5: abs error = 0.000130107; est. error = 0.0156835; evals = 9 eps=10^6: abs error = 0.000121986; est. error = 6.30233e06; evals = 17 eps=10^7: abs error = 2.14894e06; est. error = 3.30639e06; evals = 19 eps=10^8: abs error = 2.14894e06; est. error = 3.30639e06; evals = 19 eps=10^9: abs error = 5.81537e09; est. error = 3.35997e06; evals = 21 Note that the given eps is not necessarily smaller than the abs error obtained ("under delivers"). Only the last result delivers a decent estimate. Whether or not the convergence check is satisfied with the estimated error is a detail that may differ in the implementation(s) of TanhSinh, thus eps=10^4 may produce a better answer for some, but to get the better answer 19, 21 or more evaluations are needed to get close enough. The qthsh implementation (discussed on the HP forum thread) of the same algorithm with some differences has 27 evaluations for eps=10^4 with small abs error 5.63299e09. Qthsh approaches integration boundaries more consistently, evaluating additional boundary points by handling singularities instead of terminating early. Otherwise the algorithms are the same. This explains why qthsh has 27 versus WP 34S 21 evaluations. Like the WP 34S, qthsh was verified and compared to WP 34S and other algorithms, including a verification of accuracy with over 800 integrals.  Rob "I count on old friends to remain rational" 

10092021, 03:25 PM
Post: #16




RE: Accuracy of Integral with epsilon
The Prime and TI Nspire CAS emulators both (with integrand not simplified) produce an exact result of 40 (not 40.). Surprisingly, the CASIO Classpad, while readily simplifying the integrand by itself to “2”, it’s CAS, when presented with the integral with the integrand not simplified will not give an exact result.


10092021, 05:40 PM
Post: #17




RE: Accuracy of Integral with epsilon
(10092021 08:47 AM)Thomas Okken Wrote: It seems what's really being tested here is how the integrator handles the slight noisiness of the function at the low end of the range, when sinh(x) and cosh(x) become large but with opposite signs, so numerical cancellation becomes an issue. Indeed. Evaluating the integrand numerically, the HP 50g returns 2. for positive values of x but 1.99999690584 for x = 10. 

10092021, 08:37 PM
Post: #18




RE: Accuracy of Integral with epsilon
(10092021 02:27 AM)Paul Dale Wrote:(10092021 02:22 AM)Valentin Albillo Wrote: Are you saying that a complex stateoftheart TanhSinh method requires 27 evaluations of the integrand (the constant 2) to integrate the constant 2 from 10 to 10 ? Really !? There is never a guarantee that a given function is constant by probing a few points. That should be obvious, no? Naive implementations of GaussLegendre quadrature methods may stop at two or three points to fit a constant (polynomial). However, it would be detrimental to give up that early. Evidently, calculators run quite a bit longer on this function for a good reason. My point about TanhSinh was simply that this method doesn't require a huge number of points to converge accurately and it is used in the WP 34S... so make your own conclusions. TanhSinh simply works well for this type of function based on its syntactical structure (duh). More importantly, the function is NOT numerically constant in its present nonsimplified form, which is to be expected: noise increases quickly towards 10 and beyond. This noise is already in the order of 10^2 at x=18. Tabulating the error in Excel from 10 to 8.9 gives: 5.31248E08 6.01538E08 5.39791E08 1.77564E08 1.86315E08 1.64802E08 5.18465E10 8.54069E09 1.47008E08 1.40072E08 6.26628E09 5.77673E09 This amount of noise is sufficient to trip up quadrature methods. Of course, the noise may differ with nonIEEE 754 double precision and different implementations of the constituent functions EXP, SINH and COSH. However, a general consequence of the noise and floating point limitations, you can either get lucky to get to the exact result 40 with a few points or unlucky, which can cost you a great deal of time wasted to evaluate points in the presence of noise. The WP 34S and qthsh points reported show exactly what I mean. Also, attempting to push the accuracy beyond 10^8 (or about) is pointless. There is too much noise to make definitive conclusions.  Rob "I count on old friends to remain rational" 

10092021, 09:50 PM
Post: #19




RE: Accuracy of Integral with epsilon
(10092021 02:22 AM)Valentin Albillo Wrote: Someone here is not paying attention. You do realize that it is a fallacy to assume exp(x) and cosh(x)+sinh(x) equate numerically, in addition to algebraically? For one, sinh(x) typically uses expm1(x)=exp(x)1 for higher accuracy sinh(x) = (z+z/(z+1))/2 with z=expm1(x). This asymmetrically affects cosh(x)+sinh(x) for x<0. The difference grows with smaller x<0, which is easy to see when we plug in expm1(x).  Rob "I count on old friends to remain rational" 

10102021, 11:51 AM
(This post was last modified: 10142021 12:30 PM by Albert Chan.)
Post: #20




RE: Accuracy of Integral with epsilon
(10092021 08:37 PM)robve Wrote: There is never a guarantee that a given function is constant by probing a few points. That should be obvious, no? TanhSinh (and other DE variants) does not see a constant function ... it see a bellshaped curve. CAS> 'int(1,x,1,1)' (x = tanh(sinh(t))) \(\int _{\infty }^{+\infty }\mathrm{cosh}\left(t\right)\cdot (1\left(\mathrm{tanh}\left(\mathrm{sinh}\left(t\right)\right)\right)^{2})\, dt\) (10092021 09:50 PM)robve Wrote: You do realize that it is a fallacy to assume exp(x) and cosh(x)+sinh(x) equate numerically, in addition to algebraically? It is interesting that my MAPM implementation of raw_exp() actually use above "identity" Of course, x is limited to tiny value, below 1E4, so catastrophic cancellation is not an issue. sinh(x) converge with half as many terms, compared with exp(x) exp(x) = x + x^2/2! + x^3/3! + ... sinh(x) = x + x^3/3! + x^5/5! + ... 

« Next Oldest  Next Newest »

User(s) browsing this thread: