User Tools

Site Tools


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