home *** CD-ROM | disk | FTP | other *** search
/ The CDPD Public Domain Collection for CDTV 3 / CDPDIII.bin / pd / programming / gnusmalltalk / float.st < prev    next >
Text File  |  1992-02-15  |  6KB  |  233 lines
  1. "======================================================================
  2. |
  3. |   Float Method Definitions
  4. |
  5.  ======================================================================"
  6.  
  7.  
  8. "======================================================================
  9. |
  10. | Copyright (C) 1990, 1991, 1992 Free Software Foundation, Inc.
  11. | Written by Steve Byrne.
  12. |
  13. | This file is part of GNU Smalltalk.
  14. |
  15. | GNU Smalltalk is free software; you can redistribute it and/or modify it
  16. | under the terms of the GNU General Public License as published by the Free
  17. | Software Foundation; either version 1, or (at your option) any later version.
  18. | 
  19. | GNU Smalltalk is distributed in the hope that it will be useful, but WITHOUT
  20. | ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
  21. | FOR A PARTICULAR PURPOSE.  See the GNU General Public License for more
  22. | details.
  23. | 
  24. | You should have received a copy of the GNU General Public License along with
  25. | GNU Smalltalk; see the file COPYING.  If not, write to the Free Software
  26. | Foundation, 675 Mass Ave, Cambridge, MA 02139, USA.  
  27. |
  28.  ======================================================================"
  29.  
  30.  
  31. "
  32. |     Change Log
  33. | ============================================================================
  34. | Author       Date       Change 
  35. | sbb         24 Mar 91      Fixed up floating point printing to take account of
  36. |              the number of digits printed in the integer part when
  37. |              printing the fractional part.
  38. |
  39. | sbb         16 Mar 91      Class creation now separate statement.
  40. |
  41. | sbb         20 Sep 90      Fixed storeOn: to not be recursive.
  42. |
  43. | sbyrne     15 Apr 90      Added asFloat...I could have sworn this was already
  44. |                         here...must have been lost in an edit.
  45. |
  46. | sbyrne     25 Apr 89      created.
  47. |
  48. "
  49.  
  50. Number variableWordSubclass: #Float "not really a variable class"
  51.        instanceVariableNames: ''
  52.        classVariableNames: ''
  53.        poolDictionaries: ''
  54.        category: nil
  55. !
  56.  
  57. Float comment: 
  58. 'My instances represent floating point numbers that have 64 bits of 
  59. precision (well, less than that in precision; they are precisely the same
  60. as C''s "double" datatype).  Besides the standard numerical operations,
  61. I provide transcendental operations too.'  !
  62.  
  63. !Float class methodsFor: 'basic'!
  64.  
  65. pi
  66.     "Returns the value of pi"
  67.     ^3.14159265358979323846264338327950288
  68. !!
  69.  
  70.  
  71.  
  72. !Float methodsFor: 'arithmetic'!
  73.  
  74. // aNumber
  75.     "Return the integer quotient of dividing the receiver by aNumber with
  76.     truncation towards negative infinity."
  77.     ^(self / aNumber) floor
  78. !
  79.  
  80. \\ aNumber
  81.     "Return the remainder of dividing the receiver by aNumber with
  82.     truncation towards negative infinity."
  83.     ^(self - (self // aNumber * aNumber)) floor
  84. !
  85.  
  86. integerPart
  87.     ^self - self fractionPart
  88. !!
  89.  
  90.  
  91.  
  92. !Float methodsFor: 'coercing methods'!
  93.  
  94. coerce: aNumber
  95.     ^aNumber asFloat
  96. !
  97.  
  98. generality
  99.     ^4
  100. !
  101.  
  102. asFloat
  103.     "Just defined for completeness."
  104.     ^self
  105. !
  106.  
  107. hash
  108.     "Returns the hash value for the receiver (a Float).  If it's representable
  109.     as an integer, we use that value, since it's likely to be much more random,
  110.     but if not, we use the exponent.  Another approach that may be used in 
  111.     the future is to scale the number so that it's within the integer range if
  112.     it's not normally and then return the integer part"
  113.     | expt |
  114.     expt _ self exponent.
  115.     expt >= 0 & (expt < 31)
  116.         ifTrue: [ ^self truncated ]
  117.     ifFalse: [ ^expt ]
  118. !!
  119.  
  120.  
  121.  
  122. !Float methodsFor: 'copying'!
  123.  
  124. shallowCopy
  125.     ^self
  126. !
  127.  
  128. deepCopy
  129.     ^self
  130. !!
  131.  
  132.  
  133.  
  134.  
  135. !Float methodsFor: 'printing'!
  136.  
  137. printOn: aStream
  138.     | num exp |
  139.     (self exponent between: -50 and: 49)
  140.         ifTrue: [ self printFloatOn: aStream ]
  141.     ifFalse: [ exp _ 0.
  142.                num _ self.
  143.            num abs < 1.0
  144.               ifTrue: [ [ num abs < 1.0 ] whileTrue:
  145.                     [ num _ num * 10.0.
  146.                       exp _ exp - 1 ] ]
  147.               ifFalse: [ [num abs >= 10.0] whileTrue:
  148.                               [ num _ num / 10.0.
  149.                       exp _ exp + 1 ] ].
  150.                    num printFloatOn: aStream.
  151.            aStream nextPut: $e.
  152.            exp printOn: aStream ].
  153. !!
  154.  
  155.  
  156.  
  157. !Float methodsFor: 'storing'!
  158.  
  159. storeOn: aStream
  160.     self printOn: aStream
  161. !!
  162.  
  163.  
  164.  
  165. !Float methodsFor: 'private'!
  166.  
  167. printFloatOn: aStream
  168.     | num str intDigits |
  169.     num _ self.
  170.     num < 0.0 ifTrue:
  171.         [ aStream nextPut: $-.
  172.       num _ num negated. ].
  173.     intDigits _ num printIntegerPartOn: aStream.
  174.     aStream nextPut: $..
  175.     num printFracPartOn: aStream digitsPrinted: intDigits.
  176. !
  177.  
  178. printIntegerPartOn: aStream
  179.     | num digitString numDigits |
  180.     num _ self integerPart.
  181.     (digitString _ num revDigitsBase: 10.0) reverseDo:
  182.         [ :char | aStream nextPut: char ].
  183.     numDigits _ digitString size.
  184.     $0 = (digitString at: 1)    "we've been forced to emit a 0; doesn't count"
  185.     ifTrue: [ numDigits _ numDigits - 1 ].
  186.     ^numDigits
  187. !
  188.  
  189. printFracPartOn: aStream digitsPrinted: numDigits
  190.     | num count zeroCount digit digitLimit |
  191.     num _ self fractionPart.
  192.     num = 0.0 ifTrue: [ ^aStream nextPut: $0 ].
  193.     zeroCount _ count _ 0.
  194.  
  195.     " produce the digits, up to a maximum of 15, suppressing trailing
  196.       zeros (ln(2^52)/ln(10)) "
  197.     digitLimit _ 15 - numDigits.
  198.     [ (num ~= 0.0) & (count < digitLimit) ] whileTrue:
  199.         [ num _ num * 10.0.
  200.       digit _ num floor.
  201.       digit = 0        "for trailing zero suppression"
  202.         ifTrue: [ zeroCount _ zeroCount + 1 ]
  203.         ifFalse: [ aStream next: zeroCount put: $0.
  204.                        aStream nextPut: (Character digitValue: digit).
  205.                    zeroCount _ 0 ].
  206.       num _ num fractionPart.
  207.       count _ count + 1 ].
  208.  
  209.     zeroCount = digitLimit ifTrue: [ aStream nextPut: $0 ]
  210. !
  211.  
  212. "AARRGGHH!  Why can't I share this code with integer?  Grrr...."
  213. revDigitsBase: b    
  214.     | str num digit |
  215.     str _ WriteStream on: (String new: 10).
  216.     self = 0.0
  217.         ifTrue: [ str nextPut: $0 ]
  218.     ifFalse: [
  219.             num _ self.
  220.         [ num = 0.0 ] whileFalse:
  221.             [ "stdout nextPutAll: '------------'; nl.
  222.                   num basicPrint.
  223.           (num/b) basicPrint.
  224.           (num/b) fractionPart basicPrint.
  225.           ((num/b) fractionPart * b) basicPrint.
  226.           ((num/b) fractionPart * b) rounded basicPrint."
  227.           digit _ ((num / b) fractionPart * b) rounded.
  228.           str nextPut: (Character digitValue: digit).
  229.           num _ (num / b) integerPart ] ].
  230.     ^str contents
  231. ! !
  232.  
  233.