HP Calculator Wiki

User Tools

Site Tools

You are here: start » benchmarks » savage » accuracy
Trace:
benchmarks:savage:accuracy

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

— benchmarks:savage:accuracy [2013/09/13 07:36] (current)
pier4r created
Line 1: Line 1:
 +See the remarks on [[benchmarks:savage|savage benchmark]].  For that reason we "partially" dismiss the importance of the accuracy in the savage test.
  
 +===== Accuracy section =====
 +
 +==== Physical calculators ====
 +  - HP 50g, 2.15, HPGCC
 +    * Source: [[http://www.hpmuseum.org/cgi-sys/cgiwrap/hpmuseum/archv016.cgi?read=107898|New Bench ? :-)]]
 +    * ? digits
 +    * Time: around 0.52 secs @12mhz
 +    * Result: 2500 ((is it just a lucky rounding or there is an error but we don't check entirely the result? One should do 2500-result))
 +    * Absolute error: 0
 +    * Code: ?
 +  - Hp 9g
 +    * Source: [[http://www.hpmuseum.org/cgi-sys/cgiwrap/hpmuseum/archv016.cgi?read=103356|What should we get, Part 2 (MoHPC)]]
 +    * 22 digits (confirmed?)
 +    * Time: ?
 +    * Result: 2500 ((is it just a lucky rounding or there is an error but we don't check entirely the result? One should do 2500-result))
 +    * Absolute error: 0
 +    * Code: <code>
 +INPUT N;
 +A=1;
 +FOR(I=1;I<N;++I){
 +A=tan(tan-1(e^(ln(|/(A*A)))))+1;
 +}
 +PRINT A;
 +END; </code>
 +  - HP 34s, 3.2r3371
 +    * Source: [[http://www.hpmuseum.org/|Search for Posted by Paul Dale on 13 Sept 2013, 7:06 a.m.]]
 +    * ? digits
 +    * Time: ?
 +    * Result:
 +      * double 2499.999999999999999999999999999661
 +      * single 2500.000000001578
 +    * Absolute error: 
 +      * double 3.39 e-28
 +      * single 1.578 e-9
 +    * Code: ((<code>
 +LBL A
 +TICKS
 +2
 +4
 +9
 +9
 +# 001
 +LBL 00
 + x[^2]
 + [sqrt]
 + LN
 + e[^x]
 + ATAN
 + TAN
 + INC X
 + DSE Y
 + GTO 00
 +TICKS
 +RCL- T
 +SDR 001
 +x[<->] Y
 +END </code>))
 +  - sharp PC-E500S
 +    * Source: [[http://www.hpmuseum.org/cgi-sys/cgiwrap/hpmuseum/archv014.cgi?read=61403|Floating point didacticism (MoHPC)]]
 +    * ? digits
 +    * Time: ?
 +    * Result: 
 +      * single precision 2500.000056
 +      * double precision 2499.9999999999999999
 +    * Absolute error:
 +      * single precision 5.6 e-5
 +      * double precision 1 e-16
 +    * Code: ?
 +  - sharp PC-1475
 +    * Source: [[http://www.hpmuseum.org/cgi-sys/cgiwrap/hpmuseum/archv014.cgi?read=61403|Floating point didacticism (MoHPC)]]
 +    * ? digits
 +    * Time: ?
 +    * Result: 
 +      * single precision 2499.987313
 +      * double precision 2499.9999999999999925
 +    * Absolute error:
 +      * single precision 1.2687 e-2
 +      * double precision 7.5 e-15
 +    * Code: ?
 +  - HP 200 LX
 +    * Source: [[http://www.hpmuseum.org/cgi-sys/cgiwrap/hpmuseum/archv014.cgi?read=61403|Floating point didacticism (MoHPC)]]
 +    * ? digits
 +    * Time: around 210 secs
 +    * Result: 2499.999999998308
 +    * Absolute error: 1.692 e-9
 +    * Code: <code>
 +0*L(A,1)+SIGMA(I,1,N,1,
 +0*L(A,1+TAN(ATAN(EXP(LN(SQRT(g(A)*g(A))))))))
 ++L(B,g(A))-B
 +
 +2499 N and solving for B
 +</code>
 +  - TI 83+ SE
 +    * Source: [[http://www.hpmuseum.org/cgi-sys/cgiwrap/hpmuseum/archv016.cgi?read=103356|What should we get, Part 2 (MoHPC)]]
 +    * ? digits
 +    * Time: ?
 +    * Result: 2499.999992
 +    * Absolute error: 8 e-6
 +    * Code: <code>
 +1->A
 +For(I,1,2499)
 +tan(tan-1(e^(ln(\|(A^2)))))+1->A
 +End
 +Disp A
 +</code>
 +  - HP 50g, 2.15, userRPL
 +    * Source: [[http://www.hpmuseum.org/cgi-sys/cgiwrap/hpmuseum/archv017.cgi?read=116285|UserRPL Savage benchmarks]]
 +    * ? digits
 +    * Time: around 65 secs
 +    * Result: 2499.99948647
 +    * Absolute error: 5.1353 e-4
 +    * Code: ((<code>
 +@ "Savage benchmark" for 48 and 49 series.
 +@ 48 series checksum: # ECAh
 +@ 48 series size:        159
 +@ 49 series checksum: # B0C9h
 +@ 49 series size:        159.
 +\<<             @
 +  STD           @ Force standard display mode.
 +  RCLF          @ Get original flags.
 +  -55. SF       @ Force last arguments disabled.
 +  64. STWS      @ Force wordsize.
 +  RAD           @ Force radians mode.
 +  MEM DROP      @ Force a GC.
 +  TICKS         @ Initial system time.
 +  1.            @ Initial value.
 +  1. 2499.      @ Loop start/stop values.
 +  START         @
 +    DUP *       @ Square.
 +    \v/         @ Square root command.
 +    LN          @
 +    EXP         @
 +    ATAN        @
 +    TAN         @
 +    1. +        @
 +  NEXT          @
 +  TICKS         @ Ending system time.
 +  ROT           @ Move initial time to level 1.
 +  -             @ Elapsed time.
 +  B\->R         @ Convert binary to real.
 +  "Ticks"       @
 +  \->TAG        @
 +  DUP           @
 +  8192. /       @ Convert ticks to seconds.
 +  3. RND        @ Round to 3 decimal places.
 +  "Seconds"     @
 +  \->TAG        @
 +  4. ROLL STOF  @ Restore original flags.
 +\>></code>))
 +    * Source: [[http://www.hpmuseum.org/cgi-sys/cgiwrap/hpmuseum/archv021.cgi?read=239632|Good (simple?) calculation for benchmarking?]]
 +    * ? digits
 +    * Time: around 65 secs
 +    * Result: 2499.99946106
 +    * Absolute error: 5.3894 e-4
 +    * Code: ((<code>
 +<<
 +RAD TICKS
 +0.
 +1. 2500. START
 +1. + sqrt SQ LN EXP ATAN TAN
 +NEXT
 +SWAP TICKS SWAP - B->R 8192. /
 +DEG
 +800. 2. BEEP
 +>> </code>))
 +  - HP 49g
 +    * Source: [[http://www.hpmuseum.org/cgi-sys/cgiwrap/hpmuseum/archv017.cgi?read=116285|UserRPL Savage benchmarks]]
 +    * ? digits
 +    * Time: around 111 secs
 +    * Result: 2499.99948647
 +    * Absolute error: 5.1353 e-4
 +    * Code: See 50g
 +  - HP 49gx
 +    * Source: [[http://www.hpmuseum.org/cgi-sys/cgiwrap/hpmuseum/archv017.cgi?read=116285|UserRPL Savage benchmarks]]
 +    * ? digits
 +    * Time: around 118 secs
 +    * Result: 2499.99948647
 +    * Absolute error: 5.1353 e-4
 +    * Code: See 50g
 +  - HP 49gx
 +    * Source: [[http://www.hpmuseum.org/cgi-sys/cgiwrap/hpmuseum/archv017.cgi?read=116285|UserRPL Savage benchmarks]]
 +    * ? digits
 +    * Time: around 193 secs
 +    * Result: 2499.99948647
 +    * Absolute error: 5.1353 e-4
 +    * Code: See 50g
 +  - HP 28s
 +    * Source: [[http://www.hpmuseum.org/cgi-sys/cgiwrap/hpmuseum/archv017.cgi?read=116285|UserRPL Savage benchmarks]]
 +    * ? digits
 +    * Time: around 251 secs
 +    * Result: 2499.99948647
 +    * Absolute error: 5.1353 e-4
 +    * Code: ((<code>
 +@ "Savage benchmark" for 28 series.
 +@ A 28C must be in HEX mode when this program is entered to ensure
 +@ that the correct address is supplied to SYSEVAL!
 +\<<             @
 +  STD           @
 +  RCLF          @ Get original flags.
 +  31 CF         @ Force last arguments disabled.
 +  64 STWS       @ Force wordsize.
 +  RAD           @ Force radians mode.
 +  MEM DROP      @ Force a GC.
 +
 +@ Uncomment the binary integer in one of the following 3 lines.
 +  @ #123E @     @ For 28C ROM version 1BB.
 +  @ #1266 @     @ For 28C ROM version 1CC.
 +  @ #11CAh @    @ For 28S ROM version 2BB.
 +
 +
 +  SYSEVAL       @ Initial system time.
 +  1             @ Initial value.
 +  1 2499        @ Loop start/stop values.
 +  START         @
 +    DUP *       @ Square.
 +    \v/         @ Square root command.
 +    LN          @
 +    EXP         @
 +    ATAN        @
 +    TAN         @
 +    1  +        @
 +  NEXT          @
 +
 +
 +@ Uncomment the binary integer in one of the following 3 lines.
 +  @ #123E @     @ For 28C ROM version 1BB.
 +  @ #1266 @     @ For 28C ROM version 1CC.
 +  @ #11CAh @    @ For 28S ROM version 2BB.
 +
 +
 +  SYSEVAL       @ Ending system time.
 +  ROT           @ Move initial time to level 1.
 +  -             @ Elapsed time.
 +  B\->R         @ Convert binary to real.
 +  DUP           @
 +  \->STR        @ Convert real to character string.
 +  "Ticks="      @
 +  SWAP +        @
 +  SWAP          @
 +  8192 /        @ Convert ticks to seconds.
 +  3 FIX RND     @ Round to 3 decimal places.
 +  \->STR        @ Convert real to character string.
 +  "Seconds="    @
 +  SWAP +        @
 +  4 ROLL STOF   @ Restore original flags.
 +\>>
 +</code>))
 +  - HP71 Basic
 +    * Source: [[http://www.hpmuseum.org/cgi-sys/cgiwrap/hpmuseum/archv016.cgi?read=103356|What should we get, Part 2 (MoHPC)]]
 +    * 12 digits
 +    * Time: 365 [[http://www.hpmuseum.org/cgi-sys/cgiwrap/hpmuseum/archv021.cgi?read=239632|Good (simple?) calculation for benchmarking?]]
 +    * Result: 2499.99948647
 +    * Absolute error: 5.1353 e-4
 +    * Note: All the other Saturn based machines, such as the HP48, should get this result also (but maybe they have more digits!)
 +    * <code>
 +>RADIANS
 +>1 T=TIME @ A=0 @ FOR I=1 TO 2500 @ B=SQR(A+1) @ A=TAN(ATAN(EXP(LN(B*B)))) @ NEXT I
 +>2 DISP A;TIME-T
 +</code>
 +  - HP 28s
 +    * Source: [[http://www.hpmuseum.org/cgi-sys/cgiwrap/hpmuseum/archv017.cgi?read=116285|UserRPL Savage benchmarks]]
 +    * ? digits
 +    * Time: around 398 secs
 +    * Result: 2499.99948647
 +    * Absolute error: 5.1353 e-4
 +    * Code: see 28s
 +  - HP 19BII
 +    * Source: [[http://www.hpmuseum.org/cgi-sys/cgiwrap/hpmuseum/archv014.cgi?read=61403|Floating point didacticism (MoHPC)]]
 +    * ? digits
 +    * Time: around 360 secs
 +    * Result: 2500.00005173
 +    * Absolute error: 5.173 e-5
 +    * Code: <code>
 +0*L(A,1)+SIGMA(I,1,N,1,
 +0*L(A,1+TAN(ATAN(EXP(LN(SQRT(g(A)*g(A))))))))
 ++L(B,g(A))-B
 +
 +2499 N and solving for B
 +</code>
 +  - HP 39gII
 +    * Source: [[http://www.hpmuseum.org/cgi-sys/cgiwrap/hpmuseum/archv021.cgi?read=239632|Good (simple?) calculation for benchmarking?]]
 +    * ? digits
 +    * Time: 3.6 s
 +    * Result: 2499.99946106
 +    * Absolute error: 5.3894 e-4
 +    * Code: <code>
 +EXPORT SAVAGE()
 +BEGIN
 +A:=0
 +  FOR I FROM 1 TO 2500 DO
 +   A:=TAN(ATAN(EXP(LN(sqrt(A+1)²)))); 
 +  END:
 +END; </code>
 +  - HP 30b 
 +    * Source: [[http://www.hpmuseum.org/cgi-sys/cgiwrap/hpmuseum/archv020.cgi?read=193816|Savage Benchmark Results]]
 +    * ? digits
 +    * Time: around 6 secs
 +    * Result: 2499.99946106
 +    * Absolute error: 5.3894 e-4
 +    * Code: ?
 +   - HP 32S
 +    * Source: [[http://www.hpmuseum.org/cgi-sys/cgiwrap/hpmuseum/archv019.cgi?read=154497|HP 32SII vs 35S]]
 +    * ? digits
 +    * Time: around 355 s
 +    * Result: 2499.99946106
 +    * Absolute error: 5.3894 e-4
 +    * Code: ?
 +   - HP 32SII
 +    * Source: [[http://www.hpmuseum.org/cgi-sys/cgiwrap/hpmuseum/archv019.cgi?read=154497|HP 32SII vs 35S]]
 +    * ? digits
 +    * Time: around 390
 +    * Result: 2499.99946106
 +    * Absolute error: 5.3894 e-4
 +    * Code: ?
 +  - HP 42s
 +    * Source: [[http://www.hpmuseum.org/cgi-sys/cgiwrap/hpmuseum/archv020.cgi?read=193816|Savage Benchmark Results]]
 +    * ? digits
 +    * Time: around 600 secs
 +    * Result: 2499.99946106
 +    * Absolute error: 5.3894 e-4
 +    * Code: ((<code>
 +01 LBL "SB"
 +02 RAD
 +03 0
 +04 LBL 01
 +05 1
 +06 +
 +07 SQRT
 +08 X^2
 +09 LN
 +10 E^X
 +11 ATAN
 +12 TAN
 +13 DSE ST Y
 +14 GTO 01
 +15 BEEP
 +16 DEG
 +17 END
 +</code>))
 +  - HP 86b
 +    * Source: [[http://www.hpmuseum.org/cgi-sys/cgiwrap/hpmuseum/archv014.cgi?read=61403|Floating point didacticism (MoHPC)]]
 +    * ? digits
 +    * Time: around 82 secs
 +    * Result: 2499.99942403
 +    * Absolute error: 5.7597 e-4
 +    * Code: <code>
 +5 RAD
 +7 t=TIME 
 +10 A=1
 +20 FOR I=1 TO 2499
 +30 A=TAN (ATN (EXP (LOG (SQR (A*A)))))+1
 +40 NEXT I
 +50 PRINT A
 +55 PRINT HMS$ (TIME -t)
 +60 END
 +</code>
 +  - HP 33s
 +    * Source: [[http://www.hpmuseum.org/cgi-sys/cgiwrap/hpmuseum/archv016.cgi?read=103356|What should we get, Part 2 (MoHPC)]]
 +    * ? digits
 +    * Time: 475 [[http://www.hpmuseum.org/cgi-sys/cgiwrap/hpmuseum/archv019.cgi?read=154497|HP 32SII vs 35S]]
 +    * Result: 2499.99942403
 +    * Absolute error: 5.7597 e-4
 +    * Code: see the note ((<code>
 +LBL A
 +1
 +STO A
 +STO I
 +LBL B
 +RCL A
 +x^2
 +\|x (Square root of x)
 +LN
 +e^x
 +ATAN
 +TAN
 +1
 ++
 +STO A
 +1
 +STO+ I
 +RCL I
 +2,499
 +x>=y?
 +GTO B
 +RCL A
 +RTN
 +</code>))
 +  - HP 35s
 +    * Source: [[http://www.hpmuseum.org/cgi-sys/cgiwrap/hpmuseum/archv019.cgi?read=154497|HP 32SII]]
 +    * ? digits
 +    * Time: around 630 secs
 +    * Result: 2499.99939863
 +    * Absolute error: 6.0137 e-4
 +    * Code: ?
 +  - sharp PC-1262
 +    * Source: [[http://www.hpmuseum.org/cgi-sys/cgiwrap/hpmuseum/archv014.cgi?read=61403|Floating point didacticism (MoHPC)]]
 +    * ? digits
 +    * Time: ?
 +    * Result: 2499.98731
 +    * Absolute error: 1.269 e-2
 +    * Code: ?
 +  - sharp PC-1403
 +    * Source: [[http://www.hpmuseum.org/cgi-sys/cgiwrap/hpmuseum/archv014.cgi?read=61403|Floating point didacticism (MoHPC)]]
 +    * ? digits
 +    * Time: ?
 +    * Result: 2499.98731
 +    * Absolute error: 1.269 e-2
 +    * Code: ?
 +  - HP 67
 +    * Source: [[http://www.hpmuseum.org/cgi-sys/cgiwrap/hpmuseum/archv019.cgi?read=15670|which calculators use the improved log/exp algorithms]]
 +    * ? digits
 +    * Time: around 5'840 s
 +    * Result: 2499.976148
 +    * Absolute error: 2.3852 e-2
 +    * Code: ((<code>
 +01 LBL B
 +02 RAD
 +03 ST I
 +04 0
 +05 LBL 1
 +06 1
 +07 +
 +08 SQRT
 +09 X^2
 +10 LN
 +11 e^X
 +12 ATAN
 +13 TAN
 +14 DSZ I
 +15 GTO 1
 +16 DEG
 +17 RTN
 +</code>))
 +  - HP 19c
 +    * Source: [[http://www.hpmuseum.org/cgi-sys/cgiwrap/hpmuseum/archv019.cgi?read=15670|which calculators use the improved log/exp algorithms]]
 +    * ? digits
 +    * Time: around 6'600  s
 +    * Result: 2499.976148
 +    * Absolute error: 2.3852 e-2
 +    * Code: ?
 +  - HP 11C
 +    * Source: [[http://www.hpmuseum.org/cgi-sys/cgiwrap/hpmuseum/archv014.cgi?read=61403|Floating point didacticism (MoHPC)]]
 +    * ? digits
 +    * Time: ?
 +    * Result: 2499.970322
 +    * Absolute error: 2.9678 e-2
 +    * Code: ?
 +  - HP 41C
 +    * Source: [[http://www.hpmuseum.org/cgi-sys/cgiwrap/hpmuseum/archv014.cgi?read=61403|Floating point didacticism (MoHPC)]]
 +    * ? digits
 +    * Time: ?
 +    * Result: 2499.970322
 +    * Absolute error: 2.9678 e-2
 +    * Code: ?
 +  - HP 15c LE
 +    * Source: [[http://www.hpmuseum.org/cgi-sys/cgiwrap/hpmuseum/archv020.cgi?read=193816|Savage Benchmark Results]]
 +    * ? digits
 +    * Time: around 48 secs
 +    * Result: 2499.969898
 +    * Absolute error: 3.0102 e-2
 +    * Code: ((<code>
 +001 LBL B
 +002 STO 0
 +003 RAD
 +004 FIX 9
 +005 0
 +006 LBL 0
 +007 1
 +008 +
 +009 SQRT
 +010 X^2
 +011 LN
 +012 e^X
 +013 ARCTAN
 +014 TAN
 +015 DSE 0
 +016 GTO 0
 +017 DEG
 +018 RTN
 +</code>))
 +  - HP 41cx
 +    * Source: [[http://www.hpmuseum.org/cgi-sys/cgiwrap/hpmuseum/archv019.cgi?read=155670|which calculators use the improved log/exp algorithms]]
 +    * ? digits
 +    * Time: around 3'300 secs
 +    * Result: 2499.969898
 +    * Absolute error: 3.0102 e-2
 +    * Code: ((<code>
 +01 LBL B
 +02 RAD
 +03 STO I
 +04 0
 +05 LBL 1
 +06 1
 +07 +
 +08 SQRT
 +09 X^2
 +10 LN
 +11 e^X
 +12 ATAN
 +13 TAN
 +14 DSE I
 +15 GTO 1
 +16 DEG
 +17 RTN
 +</code>))
 +  - HP 15c 
 +    * Source: [[http://www.hpmuseum.org/cgi-sys/cgiwrap/hpmuseum/archv020.cgi?read=193816|Savage Benchmark Results]]
 +    * ? digits
 +    * Time: around 5'840 secs
 +    * Result: 2499.969898
 +    * Absolute error: 3.0102 e-2
 +    * Code: see 15c le
 +  - HP 34c
 +    * Source: [[http://www.hpmuseum.org/cgi-sys/cgiwrap/hpmuseum/archv019.cgi?read=155670|which calculators use the improved log/exp algorithms]]
 +    * ? digits
 +    * Time: around 7'530 secs
 +    * Result: 2499.969898
 +    * Absolute error: 3.0102 e-2
 +    * Code: ((<code>
 +01 LBL B
 +02 RAD
 +03 STO I
 +04 0
 +05 LBL 1
 +06 1
 +07 +
 +08 SQRT
 +09 X^2
 +10 LN
 +11 e^X
 +12 ATAN
 +13 TAN
 +14 DSE I
 +15 GTO 1
 +16 DEG
 +17 RTN
 +</code>))
 +  - HP 25
 +    * Source: [[http://www.hpmuseum.org/cgi-sys/cgiwrap/hpmuseum/archv019.cgi?read=155670|which calculators use the improved log/exp algorithms]]
 +    * ? digits
 +    * Time: around 5'940  s
 +    * Result: 2499.079996
 +    * Absolute error: 9.20004 e-1
 +    * Code: ((<code>
 +01 RAD
 +02 STO 0
 +03 0
 +04 STO 1
 +05 RCL 1
 +06 1
 +07 +
 +08 SQRT
 +09 X^2
 +10 LN
 +11 e^X
 +12 ATAN
 +13 TAN
 +14 STO 1
 +15 RCL 0
 +16 1
 +17 -
 +18 STO 0
 +19 X=0?
 +20 GTO 22
 +21 GTO 05
 +22 RCL 1
 +23 DEG
 +24 GTO 00
 +</code>))
 +
 +==== Emulators on mobile/handheld devices smaller than 7'' (7'' included) ====
 +  - Free 42 (other info?)
 +    * Source: [[http://www.hpmuseum.org/cgi-sys/cgiwrap/hpmuseum/archv016.cgi?read=103356|What should we get, Part 2 (MoHPC)]]
 +    * ? digits
 +    * Time: ?
 +    * Result: 2500.000000000000729718363
 +    * Absolute error: 7.29718363 e-15
 +    * Code ?
benchmarks/savage/accuracy.txt · Last modified: 2013/09/13 07:36 by pier4r

Page Tools

Except where otherwise noted, content on this wiki is licensed under the following license: CC Attribution-Noncommercial-Share Alike 4.0 International
CC Attribution-Noncommercial-Share Alike 4.0 International Donate Powered by PHP Valid HTML5 Valid CSS Driven by DokuWiki