annotate lisp/obsolete/float.el @ 86769:fc4b337f4d87

Move to etc/nxml
author Glenn Morris <rgm@gnu.org>
date Wed, 28 Nov 2007 05:12:01 +0000
parents 1e016d63f292
children e68f97372143 f55f9811f5d7
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
64758
1b36b4662ec6 Remove spurious character.
Juanma Barranquero <lekktu@gmail.com>
parents: 64751
diff changeset
1 ;;; float.el --- obsolete floating point arithmetic package
51342
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
2
74509
f7702c5f335d Update copyright years.
Glenn Morris <rgm@gnu.org>
parents: 68640
diff changeset
3 ;; Copyright (C) 1986, 2001, 2002, 2003, 2004, 2005,
75347
e3694f1cb928 Add 2007 to copyright years.
Glenn Morris <rgm@gnu.org>
parents: 74509
diff changeset
4 ;; 2006, 2007 Free Software Foundation, Inc.
51342
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
5
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
6 ;; Author: Bill Rosenblatt
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
7 ;; Maintainer: FSF
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
8 ;; Keywords: extensions
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
9
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
10 ;; This file is part of GNU Emacs.
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
11
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
12 ;; GNU Emacs is free software; you can redistribute it and/or modify
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
13 ;; it under the terms of the GNU General Public License as published by
78228
1e016d63f292 Switch license to GPLv3 or later.
Glenn Morris <rgm@gnu.org>
parents: 75347
diff changeset
14 ;; the Free Software Foundation; either version 3, or (at your option)
51342
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
15 ;; any later version.
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
16
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
17 ;; GNU Emacs is distributed in the hope that it will be useful,
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
18 ;; but WITHOUT ANY WARRANTY; without even the implied warranty of
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
19 ;; MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
20 ;; GNU General Public License for more details.
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
21
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
22 ;; You should have received a copy of the GNU General Public License
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
23 ;; along with GNU Emacs; see the file COPYING. If not, write to the
64085
18a818a2ee7c Update FSF's address.
Lute Kamstra <lute@gnu.org>
parents: 62194
diff changeset
24 ;; Free Software Foundation, Inc., 51 Franklin Street, Fifth Floor,
18a818a2ee7c Update FSF's address.
Lute Kamstra <lute@gnu.org>
parents: 62194
diff changeset
25 ;; Boston, MA 02110-1301, USA.
51342
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
26
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
27 ;;; Commentary:
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
28
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
29 ;; Floating point numbers are represented by dot-pairs (mant . exp)
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
30 ;; where mant is the 24-bit signed integral mantissa and exp is the
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
31 ;; base 2 exponent.
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
32 ;;
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
33 ;; Emacs LISP supports a 24-bit signed integer data type, which has a
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
34 ;; range of -(2**23) to +(2**23)-1, or -8388608 to 8388607 decimal.
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
35 ;; This gives six significant decimal digit accuracy. Exponents can
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
36 ;; be anything in the range -(2**23) to +(2**23)-1.
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
37 ;;
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
38 ;; User interface:
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
39 ;; function f converts from integer to floating point
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
40 ;; function string-to-float converts from string to floating point
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
41 ;; function fint converts a floating point to integer (with truncation)
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
42 ;; function float-to-string converts from floating point to string
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
43 ;;
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
44 ;; Caveats:
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
45 ;; - Exponents outside of the range of +/-100 or so will cause certain
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
46 ;; functions (especially conversion routines) to take forever.
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
47 ;; - Very little checking is done for fixed point overflow/underflow.
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
48 ;; - No checking is done for over/underflow of the exponent
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
49 ;; (hardly necessary when exponent can be 2**23).
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
50 ;;
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
51 ;;
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
52 ;; Bill Rosenblatt
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
53 ;; June 20, 1986
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
54 ;;
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
55
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
56 ;;; Code:
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
57
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
58 ;; fundamental implementation constants
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
59 (defconst exp-base 2
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
60 "Base of exponent in this floating point representation.")
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
61
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
62 (defconst mantissa-bits 24
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
63 "Number of significant bits in this floating point representation.")
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
64
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
65 (defconst decimal-digits 6
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
66 "Number of decimal digits expected to be accurate.")
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
67
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
68 (defconst expt-digits 2
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
69 "Maximum permitted digits in a scientific notation exponent.")
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
70
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
71 ;; other constants
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
72 (defconst maxbit (1- mantissa-bits)
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
73 "Number of highest bit")
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
74
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
75 (defconst mantissa-maxval (1- (ash 1 maxbit))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
76 "Maximum permissible value of mantissa")
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
77
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
78 (defconst mantissa-minval (ash 1 maxbit)
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
79 "Minimum permissible value of mantissa")
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
80
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
81 (defconst floating-point-regexp
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
82 "^[ \t]*\\(-?\\)\\([0-9]*\\)\
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
83 \\(\\.\\([0-9]*\\)\\|\\)\
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
84 \\(\\(\\([Ee]\\)\\(-?\\)\\([0-9][0-9]*\\)\\)\\|\\)[ \t]*$"
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
85 "Regular expression to match floating point numbers. Extract matches:
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
86 1 - minus sign
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
87 2 - integer part
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
88 4 - fractional part
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
89 8 - minus sign for power of ten
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
90 9 - power of ten
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
91 ")
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
92
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
93 (defconst high-bit-mask (ash 1 maxbit)
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
94 "Masks all bits except the high-order (sign) bit.")
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
95
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
96 (defconst second-bit-mask (ash 1 (1- maxbit))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
97 "Masks all bits except the highest-order magnitude bit")
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
98
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
99 ;; various useful floating point constants
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
100 (defconst _f0 '(0 . 1))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
101
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
102 (defconst _f1/2 '(4194304 . -23))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
103
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
104 (defconst _f1 '(4194304 . -22))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
105
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
106 (defconst _f10 '(5242880 . -19))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
107
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
108 ;; support for decimal conversion routines
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
109 (defvar powers-of-10 (make-vector (1+ decimal-digits) _f1))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
110 (aset powers-of-10 1 _f10)
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
111 (aset powers-of-10 2 '(6553600 . -16))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
112 (aset powers-of-10 3 '(8192000 . -13))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
113 (aset powers-of-10 4 '(5120000 . -9))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
114 (aset powers-of-10 5 '(6400000 . -6))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
115 (aset powers-of-10 6 '(8000000 . -3))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
116
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
117 (defconst all-decimal-digs-minval (aref powers-of-10 (1- decimal-digits)))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
118 (defconst highest-power-of-10 (aref powers-of-10 decimal-digits))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
119
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
120 (defun fashl (fnum) ; floating-point arithmetic shift left
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
121 (cons (ash (car fnum) 1) (1- (cdr fnum))))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
122
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
123 (defun fashr (fnum) ; floating point arithmetic shift right
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
124 (cons (ash (car fnum) -1) (1+ (cdr fnum))))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
125
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
126 (defun normalize (fnum)
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
127 (if (> (car fnum) 0) ; make sure next-to-highest bit is set
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
128 (while (zerop (logand (car fnum) second-bit-mask))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
129 (setq fnum (fashl fnum)))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
130 (if (< (car fnum) 0) ; make sure highest bit is set
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
131 (while (zerop (logand (car fnum) high-bit-mask))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
132 (setq fnum (fashl fnum)))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
133 (setq fnum _f0))) ; "standard 0"
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
134 fnum)
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
135
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
136 (defun abs (n) ; integer absolute value
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
137 (if (>= n 0) n (- n)))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
138
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
139 (defun fabs (fnum) ; re-normalize after taking abs value
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
140 (normalize (cons (abs (car fnum)) (cdr fnum))))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
141
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
142 (defun xor (a b) ; logical exclusive or
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
143 (and (or a b) (not (and a b))))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
144
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
145 (defun same-sign (a b) ; two f-p numbers have same sign?
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
146 (not (xor (natnump (car a)) (natnump (car b)))))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
147
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
148 (defun extract-match (str i) ; used after string-match
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
149 (condition-case ()
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
150 (substring str (match-beginning i) (match-end i))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
151 (error "")))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
152
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
153 ;; support for the multiplication function
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
154 (defconst halfword-bits (/ mantissa-bits 2)) ; bits in a halfword
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
155 (defconst masklo (1- (ash 1 halfword-bits))) ; isolate the lower halfword
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
156 (defconst maskhi (lognot masklo)) ; isolate the upper halfword
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
157 (defconst round-limit (ash 1 (/ halfword-bits 2)))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
158
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
159 (defun hihalf (n) ; return high halfword, shifted down
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
160 (ash (logand n maskhi) (- halfword-bits)))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
161
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
162 (defun lohalf (n) ; return low halfword
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
163 (logand n masklo))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
164
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
165 ;; Visible functions
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
166
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
167 ;; Arithmetic functions
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
168 (defun f+ (a1 a2)
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
169 "Returns the sum of two floating point numbers."
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
170 (let ((f1 (fmax a1 a2))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
171 (f2 (fmin a1 a2)))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
172 (if (same-sign a1 a2)
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
173 (setq f1 (fashr f1) ; shift right to avoid overflow
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
174 f2 (fashr f2)))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
175 (normalize
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
176 (cons (+ (car f1) (ash (car f2) (- (cdr f2) (cdr f1))))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
177 (cdr f1)))))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
178
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
179 (defun f- (a1 &optional a2) ; unary or binary minus
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
180 "Returns the difference of two floating point numbers."
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
181 (if a2
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
182 (f+ a1 (f- a2))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
183 (normalize (cons (- (car a1)) (cdr a1)))))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
184
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
185 (defun f* (a1 a2) ; multiply in halfword chunks
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
186 "Returns the product of two floating point numbers."
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
187 (let* ((i1 (car (fabs a1)))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
188 (i2 (car (fabs a2)))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
189 (sign (not (same-sign a1 a2)))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
190 (prodlo (+ (hihalf (* (lohalf i1) (lohalf i2)))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
191 (lohalf (* (hihalf i1) (lohalf i2)))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
192 (lohalf (* (lohalf i1) (hihalf i2)))))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
193 (prodhi (+ (* (hihalf i1) (hihalf i2))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
194 (hihalf (* (hihalf i1) (lohalf i2)))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
195 (hihalf (* (lohalf i1) (hihalf i2)))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
196 (hihalf prodlo))))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
197 (if (> (lohalf prodlo) round-limit)
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
198 (setq prodhi (1+ prodhi))) ; round off truncated bits
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
199 (normalize
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
200 (cons (if sign (- prodhi) prodhi)
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
201 (+ (cdr (fabs a1)) (cdr (fabs a2)) mantissa-bits)))))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
202
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
203 (defun f/ (a1 a2) ; SLOW subtract-and-shift algorithm
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
204 "Returns the quotient of two floating point numbers."
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
205 (if (zerop (car a2)) ; if divide by 0
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
206 (signal 'arith-error (list "attempt to divide by zero" a1 a2))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
207 (let ((bits (1- maxbit))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
208 (quotient 0)
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
209 (dividend (car (fabs a1)))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
210 (divisor (car (fabs a2)))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
211 (sign (not (same-sign a1 a2))))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
212 (while (natnump bits)
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
213 (if (< (- dividend divisor) 0)
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
214 (setq quotient (ash quotient 1))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
215 (setq quotient (1+ (ash quotient 1))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
216 dividend (- dividend divisor)))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
217 (setq dividend (ash dividend 1)
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
218 bits (1- bits)))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
219 (normalize
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
220 (cons (if sign (- quotient) quotient)
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
221 (- (cdr (fabs a1)) (cdr (fabs a2)) (1- maxbit)))))))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
222
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
223 (defun f% (a1 a2)
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
224 "Returns the remainder of first floating point number divided by second."
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
225 (f- a1 (f* (ftrunc (f/ a1 a2)) a2)))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
226
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
227
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
228 ;; Comparison functions
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
229 (defun f= (a1 a2)
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
230 "Returns t if two floating point numbers are equal, nil otherwise."
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
231 (equal a1 a2))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
232
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
233 (defun f> (a1 a2)
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
234 "Returns t if first floating point number is greater than second,
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
235 nil otherwise."
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
236 (cond ((and (natnump (car a1)) (< (car a2) 0))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
237 t) ; a1 nonnegative, a2 negative
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
238 ((and (> (car a1) 0) (<= (car a2) 0))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
239 t) ; a1 positive, a2 nonpositive
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
240 ((and (<= (car a1) 0) (natnump (car a2)))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
241 nil) ; a1 nonpos, a2 nonneg
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
242 ((/= (cdr a1) (cdr a2)) ; same signs. exponents differ
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
243 (> (cdr a1) (cdr a2))) ; compare the mantissas.
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
244 (t
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
245 (> (car a1) (car a2))))) ; same exponents.
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
246
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
247 (defun f>= (a1 a2)
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
248 "Returns t if first floating point number is greater than or equal to
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
249 second, nil otherwise."
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
250 (or (f> a1 a2) (f= a1 a2)))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
251
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
252 (defun f< (a1 a2)
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
253 "Returns t if first floating point number is less than second,
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
254 nil otherwise."
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
255 (not (f>= a1 a2)))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
256
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
257 (defun f<= (a1 a2)
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
258 "Returns t if first floating point number is less than or equal to
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
259 second, nil otherwise."
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
260 (not (f> a1 a2)))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
261
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
262 (defun f/= (a1 a2)
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
263 "Returns t if first floating point number is not equal to second,
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
264 nil otherwise."
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
265 (not (f= a1 a2)))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
266
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
267 (defun fmin (a1 a2)
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
268 "Returns the minimum of two floating point numbers."
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
269 (if (f< a1 a2) a1 a2))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
270
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
271 (defun fmax (a1 a2)
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
272 "Returns the maximum of two floating point numbers."
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
273 (if (f> a1 a2) a1 a2))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
274
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
275 (defun fzerop (fnum)
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
276 "Returns t if the floating point number is zero, nil otherwise."
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
277 (= (car fnum) 0))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
278
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
279 (defun floatp (fnum)
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
280 "Returns t if the arg is a floating point number, nil otherwise."
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
281 (and (consp fnum) (integerp (car fnum)) (integerp (cdr fnum))))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
282
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
283 ;; Conversion routines
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
284 (defun f (int)
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
285 "Convert the integer argument to floating point, like a C cast operator."
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
286 (normalize (cons int '0)))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
287
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
288 (defun int-to-hex-string (int)
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
289 "Convert the integer argument to a C-style hexadecimal string."
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
290 (let ((shiftval -20)
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
291 (str "0x")
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
292 (hex-chars "0123456789ABCDEF"))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
293 (while (<= shiftval 0)
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
294 (setq str (concat str (char-to-string
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
295 (aref hex-chars
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
296 (logand (lsh int shiftval) 15))))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
297 shiftval (+ shiftval 4)))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
298 str))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
299
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
300 (defun ftrunc (fnum) ; truncate fractional part
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
301 "Truncate the fractional part of a floating point number."
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
302 (cond ((natnump (cdr fnum)) ; it's all integer, return number as is
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
303 fnum)
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
304 ((<= (cdr fnum) (- maxbit)) ; it's all fractional, return 0
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
305 '(0 . 1))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
306 (t ; otherwise mask out fractional bits
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
307 (let ((mant (car fnum)) (exp (cdr fnum)))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
308 (normalize
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
309 (cons (if (natnump mant) ; if negative, use absolute value
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
310 (ash (ash mant exp) (- exp))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
311 (- (ash (ash (- mant) exp) (- exp))))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
312 exp))))))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
313
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
314 (defun fint (fnum) ; truncate and convert to integer
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
315 "Convert the floating point number to integer, with truncation,
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
316 like a C cast operator."
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
317 (let* ((tf (ftrunc fnum)) (tint (car tf)) (texp (cdr tf)))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
318 (cond ((>= texp mantissa-bits) ; too high, return "maxint"
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
319 mantissa-maxval)
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
320 ((<= texp (- mantissa-bits)) ; too low, return "minint"
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
321 mantissa-minval)
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
322 (t ; in range
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
323 (ash tint texp))))) ; shift so that exponent is 0
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
324
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
325 (defun float-to-string (fnum &optional sci)
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
326 "Convert the floating point number to a decimal string.
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
327 Optional second argument non-nil means use scientific notation."
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
328 (let* ((value (fabs fnum)) (sign (< (car fnum) 0))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
329 (power 0) (result 0) (str "")
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
330 (temp 0) (pow10 _f1))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
331
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
332 (if (f= fnum _f0)
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
333 "0"
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
334 (if (f>= value _f1) ; find largest power of 10 <= value
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
335 (progn ; value >= 1, power is positive
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
336 (while (f<= (setq temp (f* pow10 highest-power-of-10)) value)
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
337 (setq pow10 temp
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
338 power (+ power decimal-digits)))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
339 (while (f<= (setq temp (f* pow10 _f10)) value)
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
340 (setq pow10 temp
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
341 power (1+ power))))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
342 (progn ; value < 1, power is negative
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
343 (while (f> (setq temp (f/ pow10 highest-power-of-10)) value)
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
344 (setq pow10 temp
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
345 power (- power decimal-digits)))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
346 (while (f> pow10 value)
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
347 (setq pow10 (f/ pow10 _f10)
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
348 power (1- power)))))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
349 ; get value in range 100000 to 999999
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
350 (setq value (f* (f/ value pow10) all-decimal-digs-minval)
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
351 result (ftrunc value))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
352 (let (int)
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
353 (if (f> (f- value result) _f1/2) ; round up if remainder > 0.5
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
354 (setq int (1+ (fint result)))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
355 (setq int (fint result)))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
356 (setq str (int-to-string int))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
357 (if (>= int 1000000)
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
358 (setq power (1+ power))))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
359
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
360 (if sci ; scientific notation
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
361 (setq str (concat (substring str 0 1) "." (substring str 1)
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
362 "E" (int-to-string power)))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
363
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
364 ; regular decimal string
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
365 (cond ((>= power (1- decimal-digits))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
366 ; large power, append zeroes
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
367 (let ((zeroes (- power decimal-digits)))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
368 (while (natnump zeroes)
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
369 (setq str (concat str "0")
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
370 zeroes (1- zeroes)))))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
371
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
372 ; negative power, prepend decimal
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
373 ((< power 0) ; point and zeroes
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
374 (let ((zeroes (- (- power) 2)))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
375 (while (natnump zeroes)
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
376 (setq str (concat "0" str)
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
377 zeroes (1- zeroes)))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
378 (setq str (concat "0." str))))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
379
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
380 (t ; in range, insert decimal point
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
381 (setq str (concat
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
382 (substring str 0 (1+ power))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
383 "."
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
384 (substring str (1+ power)))))))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
385
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
386 (if sign ; if negative, prepend minus sign
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
387 (concat "-" str)
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
388 str))))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
389
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
390
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
391 ;; string to float conversion.
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
392 ;; accepts scientific notation, but ignores anything after the first two
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
393 ;; digits of the exponent.
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
394 (defun string-to-float (str)
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
395 "Convert the string to a floating point number.
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
396 Accepts a decimal string in scientific notation, with exponent preceded
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
397 by either E or e. Only the six most significant digits of the integer
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
398 and fractional parts are used; only the first two digits of the exponent
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
399 are used. Negative signs preceding both the decimal number and the exponent
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
400 are recognized."
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
401
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
402 (if (string-match floating-point-regexp str 0)
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
403 (let (power)
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
404 (f*
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
405 ; calculate the mantissa
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
406 (let* ((int-subst (extract-match str 2))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
407 (fract-subst (extract-match str 4))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
408 (digit-string (concat int-subst fract-subst))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
409 (mant-sign (equal (extract-match str 1) "-"))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
410 (leading-0s 0) (round-up nil))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
411
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
412 ; get rid of leading 0's
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
413 (setq power (- (length int-subst) decimal-digits))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
414 (while (and (< leading-0s (length digit-string))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
415 (= (aref digit-string leading-0s) ?0))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
416 (setq leading-0s (1+ leading-0s)))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
417 (setq power (- power leading-0s)
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
418 digit-string (substring digit-string leading-0s))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
419
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
420 ; if more than 6 digits, round off
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
421 (if (> (length digit-string) decimal-digits)
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
422 (setq round-up (>= (aref digit-string decimal-digits) ?5)
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
423 digit-string (substring digit-string 0 decimal-digits))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
424 (setq power (+ power (- decimal-digits (length digit-string)))))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
425
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
426 ; round up and add minus sign, if necessary
62194
dda44460de33 (string-to-float): Replace `string-to-int' by `string-to-number'.
Juanma Barranquero <lekktu@gmail.com>
parents: 52401
diff changeset
427 (f (* (+ (string-to-number digit-string)
51342
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
428 (if round-up 1 0))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
429 (if mant-sign -1 1))))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
430
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
431 ; calculate the exponent (power of ten)
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
432 (let* ((expt-subst (extract-match str 9))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
433 (expt-sign (equal (extract-match str 8) "-"))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
434 (expt 0) (chunks 0) (tens 0) (exponent _f1)
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
435 (func 'f*))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
436
62194
dda44460de33 (string-to-float): Replace `string-to-int' by `string-to-number'.
Juanma Barranquero <lekktu@gmail.com>
parents: 52401
diff changeset
437 (setq expt (+ (* (string-to-number
51342
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
438 (substring expt-subst 0
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
439 (min expt-digits (length expt-subst))))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
440 (if expt-sign -1 1))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
441 power))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
442 (if (< expt 0) ; if power of 10 negative
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
443 (setq expt (- expt) ; take abs val of exponent
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
444 func 'f/)) ; and set up to divide, not multiply
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
445
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
446 (setq chunks (/ expt decimal-digits)
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
447 tens (% expt decimal-digits))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
448 ; divide or multiply by "chunks" of 10**6
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
449 (while (> chunks 0)
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
450 (setq exponent (funcall func exponent highest-power-of-10)
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
451 chunks (1- chunks)))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
452 ; divide or multiply by remaining power of ten
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
453 (funcall func exponent (aref powers-of-10 tens)))))
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
454
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
455 _f0)) ; if invalid, return 0
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
456
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
457 (provide 'float)
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
458
52401
695cf19ef79e Add arch taglines
Miles Bader <miles@gnu.org>
parents: 51342
diff changeset
459 ;;; arch-tag: cc0c89c6-5718-49af-978e-585f6b14e347
51342
Juanma Barranquero <lekktu@gmail.com>
parents:
diff changeset
460 ;;; float.el ends here