st/pidigits.st

   1
   2 "
   3 * The Computer Language Shootout
   4     http://shootout.alioth.debian.org/
   5     contributed by Isaac Gouy
   6     modified by Eliot Miranda *
   7
   8 Redistribution and use in source and binary forms, with or without modification, are permitted provided that the following conditions are met:
   9
  10     * Redistributions of source code must retain the above copyright notice, this list of conditions and the following disclaimer.
  11     * Redistributions in binary form must reproduce the above copyright notice, this list of conditions and the following disclaimer in the documentation and/or other materials provided with the distribution.
  12     * Neither the name of 'The Computer Language Benchmarks Game' nor the name of 'The Computer Language Shootout Benchmarks' nor the names of its contributors may be used to endorse or promote products derived from this software without specific prior written permission.
  13
  14 THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS 'AS IS' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE."
  15
  16 PiDigitSpigot method!
  17 consume: aTransformation
  18    ^z * aTransformation!
  19
  20 PiDigitSpigot method!
  21 digit
  22    ^(z extract: 3) floor!
  23
  24 PiDigitSpigot method!
  25 isSafe: aDigit
  26    ^ aDigit = (z extract: 4) floor!
  27
  28 PiDigitSpigot method!
  29 produce: anInteger
  30    inverse q: 10 r: -10 * anInteger s: 0 t: 1.
  31    ^inverse * z!
  32
  33 PiDigitSpigot method!
  34 next
  35    | y |
  36    ^(self isSafe: (y := self digit))
  37       ifTrue: [z := self produce: y. self assert: [ y ~= 0 ]. y]
  38       ifFalse: [z := self consume: x next. y := self next. self assert: [ y ~= 0 ]. y]!
  39
  40 PiDigitSpigot method!
  41 initialize
  42    z := Transformation unity.
  43    x := Transformation new.
  44    inverse := Transformation new.!
  45
  46 Number class method!
  47 pidigitsTo: v width: width to: output
  48    | n i pidigits |
  49    n := v.
  50    i := 0.
  51    pidigits := PiDigitSpigot new.
  52    [n > 0] whileTrue:
  53       [n < width
  54          ifTrue:
  55             [n timesRepeat: [output nextPutAll: pidigits next printString].
  56              n to: width do: [:each | output nextPut: Character space].
  57              i := i + n]
  58          ifFalse:
  59             [width timesRepeat: [output nextPutAll: pidigits next printString].
  60              i := i + width].
  61
  62       output nextPut: Character tab;
  63                          nextPut: $:;
  64                          nextPutAll: i printString;
  65              nextPut: Character cr.
  66
  67       n := n - width]!
  68
  69 Transformation method!
  70 * aTransformation
  71    ^ Transformation
  72       q: q * aTransformation q
  73       r: q * aTransformation r + (r * aTransformation t)
  74       s: s * aTransformation q + (t * aTransformation s)
  75       t: s * aTransformation r + (t * aTransformation t)!
  76
  77 Transformation method!
  78 extract: anInteger
  79    ^(q * anInteger + r) // (s * anInteger + t)!
  80
  81 Transformation method!
  82 next
  83    k := k + 1.
  84    q := k.
  85    r := 4 * k + 2.
  86    s := 0.
  87    t := 2 * k + 1.!
  88
  89 Transformation method!
  90 q
  91    ^q!
  92
  93 Transformation method!
  94 q: anInteger1 r: anInteger2 s: anInteger3 t: anInteger4
  95    q := anInteger1.
  96    r := anInteger2.
  97    s := anInteger3.
  98    t := anInteger4.
  99    k := 0.!
 100
 101 Transformation method!
 102 r
 103    ^r!
 104
 105 Transformation method!
 106 s
 107    ^s!
 108
 109 Transformation method!
 110 t
 111    ^t!
 112
 113 Transformation method!
 114 initialize
 115    q := 0.
 116    r := 0.
 117    s := 0.
 118    t := 0.
 119    k := 0.!
 120
 121 Transformation class method!
 122 q: anInteger1 r: anInteger2 s: anInteger3 t: anInteger4
 123    ^(self basicNew) q: anInteger1 r: anInteger2 s: anInteger3 t: anInteger4!
 124
 125 Transformation class method!
 126 unity
 127    ^self q: 1 r: 0 s: 0 t: 1!