Casio fx991EX vs Hp 50g speed difference

07132018, 08:50 PM
(This post was last modified: 07132018 08:52 PM by Tugdual.)
Post: #1




Casio fx991EX vs Hp 50g speed difference
Found something weird.
Was speed testing the fx991ex vs 50g on \(\int_0^{500}e^{x} dx \) Results: fx991EX: result:1 time:6.64s 50g using >NUM : result:1 time:46.69s Huge difference... Pressing EVAL on the 50g returns instantly \(  \frac {1e^{500}}{e^{500}} \) which seems to indicate that CAS was used since this is the exact result. I'm a bit surprised the 50g performs so poorly since I doubt the fx991ex uses any smart CAS behind the scene. Any thoughts on that? 

07132018, 10:36 PM
(This post was last modified: 07132018 10:44 PM by StephenG1CMZ.)
Post: #2




RE: Casio fx991EX vs Hp 50g speed difference
On HP Prime I get result 1(1/e^500), in 20 ms (Android), using CAS.
Stephen Lewkowicz (G1CMZ) 

07132018, 11:13 PM
Post: #3




RE: Casio fx991EX vs Hp 50g speed difference
(07132018 08:50 PM)Tugdual Wrote: I'm a bit surprised the 50g performs so poorly since I doubt the fx991ex uses any smart CAS behind the scene. Any thoughts on that? I can't really say why the HP 50g takes awhile to get an approximate answer for the given integral. I do know that if you take it out of Standard format and put it into a smaller number format like Fix 4, it would get a numeric answer quicker with an accuracy up to the given number format. But you probably knew that already. (07132018 08:50 PM)Tugdual Wrote: Pressing EVAL on the 50g returns instantly \(  \frac {1e^{500}}{e^{500}} \) which seems to indicate that CAS was used since this is the exact result. Weird right? The difference between Exact and Approx. speeds in this case is very interesting and I don't know why that is. I get the answer in 0.70 milliseconds if I want an exact result. And when I press >NUM on the exact result, it solves the exact answer in .06 milliseconds. This means I can get the approximate result of the integral in .76 milliseconds, beating the Casio. But every attempt to solve the integral approximately, I get a 47 second wait. I'm sorry, I don't have a guess on why it does this. (07132018 10:36 PM)StephenG1CMZ Wrote: On HP Prime I get result 1(1/e^500), in 20 ms (Android), using CAS. The Prime is great isn't it? But the 50g getting the same answer in .70 milliseconds is also impressive to me. 

07142018, 01:45 AM
Post: #4




RE: Casio fx991EX vs Hp 50g speed difference
(07132018 08:50 PM)Tugdual Wrote: Found something weird. The function e^x is not very wellbehaved for integration over the interval [0, 500], being close to zero for most of it. Everything happens in a small interval near 0, and that kind of thing makes some numerical integrators converge rather slowly. On my HP42S, which I'm guessing uses the same Romberg integration algorithm as the RPL numerical integrator, the integral in question takes 2047 evaluations of the function (with ACC=0, which is full accuracy), taking 4 minutes and 4 seconds. 

07142018, 03:17 AM
(This post was last modified: 07142018 03:18 AM by Albert Chan.)
Post: #5




RE: Casio fx991EX vs Hp 50g speed difference
I wonder if Casio flush values below 1e99 to zero, which might save time
for integration (flushtozero section "converge" to zero area right the way) For a fair test, try limit of 0 to 200 Just curious, is the final result same for both calculator ? 

07142018, 10:27 AM
Post: #6




RE: Casio fx991EX vs Hp 50g speed difference
another result that may be interesting.
http://www.hpmuseum.org/forum/thread9750.html (sum, trig, power/root , exp/log) Wikis are great, Contribute :) 

07142018, 03:10 PM
(This post was last modified: 07142018 03:11 PM by ijabbott.)
Post: #7




RE: Casio fx991EX vs Hp 50g speed difference
(07142018 10:27 AM)pier4r Wrote: another result that may be interesting. I appear to have one of the "fast" batch. The diagnostics show "CY235 VerB". To enter diagnostics mode, while holding "SHIFT" and "7" together, press "ON", then release all buttons. You should see: Quote:DIAGNOSTIC (If it doesn't work, press and release "ON" first, then try again.) Pressing "9" in diagnostic mode shows: Quote:8888888888888888 Then pressing "SHIFT" repeatedly goes through some display tests, ending with: Quote:CY235 VerB (The version may be different on the "slow" batch.) Pressing "AC" shows "00" in the top corner of the display. If you press the keys in the correct order, it increments to "01", "02" etc. On the initial diagnostic screen, pressing "8" goes into a more brief keyboard test where various labels get highlighted when you press the corresponding key (although the \(x\) key is labelled as "[logab]" for some reason!). After pressing all the keys, the display changes to: Quote:Solar MODEL OK! — Ian Abbott 

07162018, 07:07 AM
(This post was last modified: 07162018 07:08 AM by Tugdual.)
Post: #8




RE: Casio fx991EX vs Hp 50g speed difference
Just for the records mine says
CY235 VerF And seems to belong to the fast batch 

07182018, 04:28 AM
(This post was last modified: 07182018 05:50 AM by brickviking.)
Post: #9




RE: Casio fx991EX vs Hp 50g speed difference
(07132018 08:50 PM)Tugdual Wrote: Found something weird. My upgraded 9750GII gave me 1 for that formula: \( \int \)(e^X,0,500), almost immediately. If I add a tolerance value: \( \int \)(e^X,0,500,0.01) I get 1.000000014 in about 0.5 seconds for tolerance values of 0.009 to 0.1724, and 1.000050478 (for a tolerance value of >0.1724 to 0.66558). Larger tolerance values affect the result more. It'd be interesting to see whether the 991EX has this tolerance value in the formula it uses. (07132018 08:50 PM)Tugdual Wrote: Pressing EVAL on the 50g returns instantly \(  \frac {1e^{500}}{e^{500}} \) which seems to indicate that CAS was used since this is the exact result. My HP50G gave me \(\frac {e^{500}1}{e^{500}} \) instantly instead of \(  \frac {1e^{500}}{e^{500}} \) when I hit EVAL from the equation writer. I presume those two are equivalent. I also observed the Approx/Exact discrepancy. (Post 258) Regards, BrickViking HP50g Casio fx9750G+ Casio fx9750GII (SH4a) 

07252018, 01:50 AM
Post: #10




RE: Casio fx991EX vs Hp 50g speed difference
(07132018 08:50 PM)Tugdual Wrote: Found something weird. Try this: \(\int_0^{200}sin(x) dx \) = 1  cos(200) ~ 0.51281 If 50g returned 96.581, that means: 1. No symbolic processing used for integration 2. Integration by Simpson rule, thus poor performance (and wrong!) 

07252018, 02:44 AM
Post: #11




RE: Casio fx991EX vs Hp 50g speed difference
When I enter the integral in exact mode, I get (COS(200)1) in .51 seconds. Then I proceed to get the numerical answer of 0.512812324993 in .04 seconds by executing >NUM on the equation. Thus, computing it in .55 seconds.
When I enter the integral in approx. mode, I get the answer of 0.512812324889 in 111.36 seconds with an error of 1.3E9. 

07252018, 12:55 PM
Post: #12




RE: Casio fx991EX vs Hp 50g speed difference
(07252018 01:50 AM)Albert Chan Wrote: Try this: \(\int_0^{200}sin(x) dx \) = 1  cos(200) ~ 0.51281 My mistake, Simpson Rule + Romberg's Extrapolation led to a wrong answer Code: Intervals Simpson Rule Romberg Extrapolation Simpson Rule alone will eventually converge: 0.51281 (2048 intervals), 0.5128123250 (32768 intervals) 

07252018, 06:32 PM
Post: #13




RE: Casio fx991EX vs Hp 50g speed difference
Just learned there is a trick to use Romberg's method for periodic function:
Casio keisan online Romberg Integration calculator, option nonlinear substitution Code: def nonlinear(f, a, b): Just integrate above transformed f, from 1 to 1: \(\int_0^{200}sin(x) dx \) = \(\int_{1}^{1}150(1u^2) sin[50u (3  u^2) + 100] du \) On my Casio FX115MS, it had gotten the sin(x) integral wrong: 96.581 With transformed integral, after 3 minutes 45 seconds, I got 0.51281 :) Amazingly, Mathematica can solve the messy transformed integral symbolically: 2 (sin(100))^2 

07252018, 08:07 PM
Post: #14




RE: Casio fx991EX vs Hp 50g speed difference
(07252018 06:32 PM)Albert Chan Wrote: Just integrate above transformed f, from 1 to 1: My 50g reduces the transformed integral to (COS(200)1) in 1.7 seconds. If I evaluate the same transformed integral in approx. mode, I get the answer of .512812325001 in 61.45 seconds. 

07252018, 10:13 PM
Post: #15




RE: Casio fx991EX vs Hp 50g speed difference
(07252018 06:32 PM)Albert Chan Wrote: Just integrate above transformed f, from 1 to 1: (07252018 02:44 AM)Carsen Wrote: When I enter the integral in approx. mode, I get the answer of 0.512812324889 in 111.36 seconds with an error of 1.3E9. (07252018 08:07 PM)Carsen Wrote: If I evaluate the transformed integral in approx. mode, I get the answer of .512812325001 in 61.45 seconds. For transformed integral, it only take half the time ? And, compare against 1  cos(200), more accurate ? Is it tested on the same calculator ? What if it were applied not just to periodic function, but others ? Say, this thread original post ? \(\int_0^{500}e^{x}dx \) = \(\int_{1}^{1}375(1u^2) e^{(125u (3  u^2) + 250)} du \) 

07262018, 03:01 AM
Post: #16




RE: Casio fx991EX vs Hp 50g speed difference
(07252018 10:13 PM)Albert Chan Wrote: For transformed integral, it only take half the time ? No. It took 0.51 seconds to solve the SIN(X) integral symbolically (This time is found in post #11 in this thread. I apologize for not defining what integral I evaluated in that post). On the other hand, it took 1.7 seconds to solve the transformed SIN(X) integral symbolically (This time is found in post #14 in this thread). Both of them returned the same answer of (COS(200)1). (07252018 10:13 PM)Albert Chan Wrote: And, compare against 1  cos(200), more accurate ? I would say that 1COS(200) has the same accuracy as (COS(200)1). They both are the same. (07252018 10:13 PM)Albert Chan Wrote: Is it tested on the same calculator ? Yes. All is tested on the same HP 50g. (07252018 10:13 PM)Albert Chan Wrote: What if it were applied not just to periodic function, but others ? Same as the first integral of SIN(X). Symbolically, the 50g solves the integral e^x in .76 seconds, while it solves its transformed counterpart in 2.84 seconds. Again, both get the same answer of (e^5001)/e^500. Here are some pictures of my results. 

07262018, 03:28 AM
Post: #17




RE: Casio fx991EX vs Hp 50g speed difference
Hi, Carsen
I might have phrase the question poorly. I am not interested in the performance of symbolic solver. Taking longer to solve the transformed integral is not surprising. I am more interested in numerical integration. I noticed your timings in the approx mode, transformed integral solving time is halved, with better accuracy. What is the timings for the transformed exponential integral (in approx mode) ? Tried the transformed exponential integral on my old computer. It speed up both Romberg and Gaussian Quadrature code. 

07262018, 04:53 PM
Post: #18




RE: Casio fx991EX vs Hp 50g speed difference
(07262018 03:28 AM)Albert Chan Wrote: I might have phrase the question poorly. Possibly but this was the main reason. (07252018 10:13 PM)Albert Chan Wrote: And, compare against 1  cos(200), more accurate ? I saw the 1COS(200) and immediately thought, "that's a symbolic answer so that must be what we are talking about". Of course, I didn't really understand how you could mention accuracy with an symbolic answer but I didn't give that much thought. I should think a little more to see if what I'm am doing (posting, solving, fixing, anything in general). But now that I understand the miscommunication, here is the result. The time it took to evaluate the approx. answer of the transformed integral of e^x is 234.53 seconds. That's 3 minutes and 54.53 seconds. 

07262018, 07:08 PM
Post: #19




RE: Casio fx991EX vs Hp 50g speed difference
Thanks, Carsen for the test.
I think the transformed sine integral is faster because it converge slightly faster, skipping 1 Simpson iteration (that is why the halved timing). But, that is just lucky. It could be the other way around, like the transformed exponential integral. 5X slow down possibly caused by 2 extra Simpson iteration for convergence. Transformed integral is not really for performance, but to avoid false convergence. 

09032018, 07:40 PM
Post: #20




RE: Casio fx991EX vs Hp 50g speed difference
I know for me personally when I solve the exponential integral on the TInspire it takes a really long time depending on the value chosen for x:
https://en.wikipedia.org/wiki/Exponential_integral When I program it into my WP34s (real hardware) it takes probably 1/20th the time. This is because I use an algorithm specifically designed to solve the exponential integral quickly. I found in my personal research in regards to integrals of exponential functions that algorithms that are really good at them are generally worse at other types of integrals and vice versa. The HP 50g probably just uses an algorithm that is not very efficient at solving integrals of exponential functions. You could test this by programming the 50g with your own custom solver that is good with exponentials, I'm sure it would be must faster. 

« Next Oldest  Next Newest »

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