16ded16be15dd865a9b21ea304d5273c8be299c87Steve Block;; Copyright 2010 the V8 project authors. All rights reserved. 26ded16be15dd865a9b21ea304d5273c8be299c87Steve Block;; Redistribution and use in source and binary forms, with or without 36ded16be15dd865a9b21ea304d5273c8be299c87Steve Block;; modification, are permitted provided that the following conditions are 46ded16be15dd865a9b21ea304d5273c8be299c87Steve Block;; met: 56ded16be15dd865a9b21ea304d5273c8be299c87Steve Block;; 66ded16be15dd865a9b21ea304d5273c8be299c87Steve Block;; * Redistributions of source code must retain the above copyright 76ded16be15dd865a9b21ea304d5273c8be299c87Steve Block;; notice, this list of conditions and the following disclaimer. 86ded16be15dd865a9b21ea304d5273c8be299c87Steve Block;; * Redistributions in binary form must reproduce the above 96ded16be15dd865a9b21ea304d5273c8be299c87Steve Block;; copyright notice, this list of conditions and the following 106ded16be15dd865a9b21ea304d5273c8be299c87Steve Block;; disclaimer in the documentation and/or other materials provided 116ded16be15dd865a9b21ea304d5273c8be299c87Steve Block;; with the distribution. 126ded16be15dd865a9b21ea304d5273c8be299c87Steve Block;; * Neither the name of Google Inc. nor the names of its 136ded16be15dd865a9b21ea304d5273c8be299c87Steve Block;; contributors may be used to endorse or promote products derived 146ded16be15dd865a9b21ea304d5273c8be299c87Steve Block;; from this software without specific prior written permission. 156ded16be15dd865a9b21ea304d5273c8be299c87Steve Block;; 166ded16be15dd865a9b21ea304d5273c8be299c87Steve Block;; THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS 176ded16be15dd865a9b21ea304d5273c8be299c87Steve Block;; "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT 186ded16be15dd865a9b21ea304d5273c8be299c87Steve Block;; LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR 196ded16be15dd865a9b21ea304d5273c8be299c87Steve Block;; A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT 206ded16be15dd865a9b21ea304d5273c8be299c87Steve Block;; OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, 216ded16be15dd865a9b21ea304d5273c8be299c87Steve Block;; SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT 226ded16be15dd865a9b21ea304d5273c8be299c87Steve Block;; LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, 236ded16be15dd865a9b21ea304d5273c8be299c87Steve Block;; DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY 246ded16be15dd865a9b21ea304d5273c8be299c87Steve Block;; THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT 256ded16be15dd865a9b21ea304d5273c8be299c87Steve Block;; (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE 266ded16be15dd865a9b21ea304d5273c8be299c87Steve Block;; OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. 276ded16be15dd865a9b21ea304d5273c8be299c87Steve Block 286ded16be15dd865a9b21ea304d5273c8be299c87Steve Block;; This is a Scheme script for the Bigloo compiler. Bigloo must be compiled with 296ded16be15dd865a9b21ea304d5273c8be299c87Steve Block;; support for bignums. The compilation of the script can be done as follows: 306ded16be15dd865a9b21ea304d5273c8be299c87Steve Block;; bigloo -static-bigloo -o generate-ten-powers generate-ten-powers.scm 316ded16be15dd865a9b21ea304d5273c8be299c87Steve Block;; 326ded16be15dd865a9b21ea304d5273c8be299c87Steve Block;; Generate approximations of 10^k. 336ded16be15dd865a9b21ea304d5273c8be299c87Steve Block 346ded16be15dd865a9b21ea304d5273c8be299c87Steve Block(module gen-ten-powers 356ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (static (class Cached-Fast 366ded16be15dd865a9b21ea304d5273c8be299c87Steve Block v::bignum 376ded16be15dd865a9b21ea304d5273c8be299c87Steve Block e::bint 386ded16be15dd865a9b21ea304d5273c8be299c87Steve Block exact?::bool)) 396ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (main my-main)) 406ded16be15dd865a9b21ea304d5273c8be299c87Steve Block 416ded16be15dd865a9b21ea304d5273c8be299c87Steve Block 426ded16be15dd865a9b21ea304d5273c8be299c87Steve Block;;----------------bignum shifts ----------------------------------------------- 436ded16be15dd865a9b21ea304d5273c8be299c87Steve Block(define (bit-lshbx::bignum x::bignum by::bint) 446ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (if (by 0) 456ded16be15dd865a9b21ea304d5273c8be299c87Steve Block #z0 466ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (*bx x (exptbx #z2 (fixnum->bignum by))))) 476ded16be15dd865a9b21ea304d5273c8be299c87Steve Block 486ded16be15dd865a9b21ea304d5273c8be299c87Steve Block(define (bit-rshbx::bignum x::bignum by::bint) 496ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (if ( by 0) 506ded16be15dd865a9b21ea304d5273c8be299c87Steve Block #z0 516ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (/bx x (exptbx #z2 (fixnum->bignum by))))) 526ded16be15dd865a9b21ea304d5273c8be299c87Steve Block 536ded16be15dd865a9b21ea304d5273c8be299c87Steve Block;;----------------the actual power generation ------------------------------- 546ded16be15dd865a9b21ea304d5273c8be299c87Steve Block 556ded16be15dd865a9b21ea304d5273c8be299c87Steve Block;; e should be an indication. it might be too small. 566ded16be15dd865a9b21ea304d5273c8be299c87Steve Block(define (round-n-cut n e nb-bits) 576ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (define max-container (- (bit-lshbx #z1 nb-bits) 1)) 586ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (define (round n) 596ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (case *round* 606ded16be15dd865a9b21ea304d5273c8be299c87Steve Block ((down) n) 616ded16be15dd865a9b21ea304d5273c8be299c87Steve Block ((up) 626ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (+bx n 636ded16be15dd865a9b21ea304d5273c8be299c87Steve Block ;; with the -1 it will only round up if the cut off part is 646ded16be15dd865a9b21ea304d5273c8be299c87Steve Block ;; non-zero 656ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (-bx (bit-lshbx #z1 666ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (-fx (+fx e nb-bits) 1)) 676ded16be15dd865a9b21ea304d5273c8be299c87Steve Block #z1))) 686ded16be15dd865a9b21ea304d5273c8be299c87Steve Block ((round) 696ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (+bx n 706ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (bit-lshbx #z1 716ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (-fx (+fx e nb-bits) 2)))))) 726ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (let* ((shift (-fx (+fx e nb-bits) 1)) 736ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (cut (bit-rshbx (round n) shift)) 746ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (exact? (=bx n (bit-lshbx cut shift)))) 756ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (if (<=bx cut max-container) 766ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (values cut e exact?) 776ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (round-n-cut n (+fx e 1) nb-bits)))) 786ded16be15dd865a9b21ea304d5273c8be299c87Steve Block 796ded16be15dd865a9b21ea304d5273c8be299c87Steve Block(define (rounded-/bx x y) 806ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (case *round* 816ded16be15dd865a9b21ea304d5273c8be299c87Steve Block ((down) (/bx x y)) 826ded16be15dd865a9b21ea304d5273c8be299c87Steve Block ((up) (+bx (/bx x y) #z1)) 836ded16be15dd865a9b21ea304d5273c8be299c87Steve Block ((round) (let ((tmp (/bx (*bx #z2 x) y))) 846ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (if (zerobx? (remainderbx tmp #z2)) 856ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (/bx tmp #z2) 866ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (+bx (/bx tmp #z2) #z1)))))) 876ded16be15dd865a9b21ea304d5273c8be299c87Steve Block 886ded16be15dd865a9b21ea304d5273c8be299c87Steve Block(define (generate-powers from to mantissa-size) 896ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (let* ((nb-bits mantissa-size) 906ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (offset (- from)) 916ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (nb-elements (+ (- from) to 1)) 926ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (vec (make-vector nb-elements)) 936ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (max-container (- (bit-lshbx #z1 nb-bits) 1))) 946ded16be15dd865a9b21ea304d5273c8be299c87Steve Block ;; the negative ones. 10^-1, 10^-2, etc. 956ded16be15dd865a9b21ea304d5273c8be299c87Steve Block ;; We already know, that we can't be exact, so exact? will always be #f. 966ded16be15dd865a9b21ea304d5273c8be299c87Steve Block ;; Basically we will have a ten^i that we will *10 at each iteration. We 976ded16be15dd865a9b21ea304d5273c8be299c87Steve Block ;; want to create the matissa of 1/ten^i. However the mantissa must be 986ded16be15dd865a9b21ea304d5273c8be299c87Steve Block ;; normalized (start with a 1). -> we have to shift the number. 996ded16be15dd865a9b21ea304d5273c8be299c87Steve Block ;; We shift by multiplying with two^e. -> We encode two^e*(1/ten^i) == 1006ded16be15dd865a9b21ea304d5273c8be299c87Steve Block ;; two^e/ten^i. 1016ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (let loop ((i 1) 1026ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (ten^i #z10) 1036ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (two^e #z1) 1046ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (e 0)) 1056ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (unless (< (- i) from) 1066ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (if (>bx (/bx (*bx #z2 two^e) ten^i) max-container) 1076ded16be15dd865a9b21ea304d5273c8be299c87Steve Block ;; another shift would make the number too big. We are 1086ded16be15dd865a9b21ea304d5273c8be299c87Steve Block ;; hence normalized now. 1096ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (begin 1106ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (vector-set! vec (-fx offset i) 1116ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (instantiate::Cached-Fast 1126ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (v (rounded-/bx two^e ten^i)) 1136ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (e (negfx e)) 1146ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (exact? #f))) 1156ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (loop (+fx i 1) (*bx ten^i #z10) two^e e)) 1166ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (loop i ten^i (bit-lshbx two^e 1) (+fx e 1))))) 1176ded16be15dd865a9b21ea304d5273c8be299c87Steve Block ;; the positive ones 10^0, 10^1, etc. 1186ded16be15dd865a9b21ea304d5273c8be299c87Steve Block ;; start with 1.0. mantissa: 10...0 (1 followed by nb-bits-1 bits) 1196ded16be15dd865a9b21ea304d5273c8be299c87Steve Block ;; -> e = -(nb-bits-1) 1206ded16be15dd865a9b21ea304d5273c8be299c87Steve Block ;; exact? is true when the container can still hold the complete 10^i 1216ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (let loop ((i 0) 1226ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (n (bit-lshbx #z1 (-fx nb-bits 1))) 1236ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (e (-fx 1 nb-bits))) 1246ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (when (<= i to) 1256ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (receive (cut e exact?) 1266ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (round-n-cut n e nb-bits) 1276ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (vector-set! vec (+fx i offset) 1286ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (instantiate::Cached-Fast 1296ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (v cut) 1306ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (e e) 1316ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (exact? exact?))) 1326ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (loop (+fx i 1) (*bx n #z10) e)))) 1336ded16be15dd865a9b21ea304d5273c8be299c87Steve Block vec)) 1346ded16be15dd865a9b21ea304d5273c8be299c87Steve Block 1356ded16be15dd865a9b21ea304d5273c8be299c87Steve Block(define (print-c powers from to struct-type 1366ded16be15dd865a9b21ea304d5273c8be299c87Steve Block cache-name max-distance-name offset-name macro64) 1376ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (define (display-power power k) 1386ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (with-access::Cached-Fast power (v e exact?) 1396ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (let ((tmp-p (open-output-string))) 1406ded16be15dd865a9b21ea304d5273c8be299c87Steve Block ;; really hackish way of getting the digits 1416ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (display (format "~x" v) tmp-p) 1426ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (let ((str (close-output-port tmp-p))) 1436ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (printf " {~a(0x~a, ~a), ~a, ~a},\n" 1446ded16be15dd865a9b21ea304d5273c8be299c87Steve Block macro64 1456ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (substring str 0 8) 1466ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (substring str 8 16) 1476ded16be15dd865a9b21ea304d5273c8be299c87Steve Block e 1486ded16be15dd865a9b21ea304d5273c8be299c87Steve Block k))))) 1496ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (define (print-powers-reduced n) 1506ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (print "static const " struct-type " " cache-name 1516ded16be15dd865a9b21ea304d5273c8be299c87Steve Block "(" n ")" 1526ded16be15dd865a9b21ea304d5273c8be299c87Steve Block "[] = {") 1536ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (let loop ((i 0) 1546ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (nb-elements 0) 1556ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (last-e 0) 1566ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (max-distance 0)) 1576ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (cond 1586ded16be15dd865a9b21ea304d5273c8be299c87Steve Block ((>= i (vector-length powers)) 1596ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (print " };") 1606ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (print "static const int " max-distance-name "(" n ") = " 1616ded16be15dd865a9b21ea304d5273c8be299c87Steve Block max-distance ";") 1626ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (print "// nb elements (" n "): " nb-elements)) 1636ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (else 1646ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (let* ((power (vector-ref powers i)) 1656ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (e (Cached-Fast-e power))) 1666ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (display-power power (+ i from)) 1676ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (loop (+ i n) 1686ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (+ nb-elements 1) 1696ded16be15dd865a9b21ea304d5273c8be299c87Steve Block e 1706ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (cond 1716ded16be15dd865a9b21ea304d5273c8be299c87Steve Block ((=fx i 0) max-distance) 1726ded16be15dd865a9b21ea304d5273c8be299c87Steve Block ((> (- e last-e) max-distance) (- e last-e)) 1736ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (else max-distance)))))))) 1746ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (print "// Copyright 2010 the V8 project authors. All rights reserved.") 1756ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (print "// ------------ GENERATED FILE ----------------") 1766ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (print "// command used:") 1776ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (print "// " 1786ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (apply string-append (map (lambda (str) 1796ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (string-append " " str)) 1806ded16be15dd865a9b21ea304d5273c8be299c87Steve Block *main-args*)) 1816ded16be15dd865a9b21ea304d5273c8be299c87Steve Block " // NOLINT") 1826ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (print) 1836ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (print 1846ded16be15dd865a9b21ea304d5273c8be299c87Steve Block "// This file is intended to be included inside another .h or .cc files\n" 1856ded16be15dd865a9b21ea304d5273c8be299c87Steve Block "// with the following defines set:\n" 1866ded16be15dd865a9b21ea304d5273c8be299c87Steve Block "// GRISU_CACHE_STRUCT: should expand to the name of a struct that will\n" 1876ded16be15dd865a9b21ea304d5273c8be299c87Steve Block "// hold the cached powers of ten. Each entry will hold a 64-bit\n" 1886ded16be15dd865a9b21ea304d5273c8be299c87Steve Block "// significand, a 16-bit signed binary exponent, and a 16-bit\n" 1896ded16be15dd865a9b21ea304d5273c8be299c87Steve Block "// signed decimal exponent. Each entry will be constructed as follows:\n" 1906ded16be15dd865a9b21ea304d5273c8be299c87Steve Block "// { significand, binary_exponent, decimal_exponent }.\n" 1916ded16be15dd865a9b21ea304d5273c8be299c87Steve Block "// GRISU_CACHE_NAME(i): generates the name for the different caches.\n" 1926ded16be15dd865a9b21ea304d5273c8be299c87Steve Block "// The parameter i will be a number in the range 1-20. A cache will\n" 1936ded16be15dd865a9b21ea304d5273c8be299c87Steve Block "// hold every i'th element of a full cache. GRISU_CACHE_NAME(1) will\n" 1946ded16be15dd865a9b21ea304d5273c8be299c87Steve Block "// thus hold all elements. The higher i the fewer elements it has.\n" 1956ded16be15dd865a9b21ea304d5273c8be299c87Steve Block "// Ideally the user should only reference one cache and let the\n" 1966ded16be15dd865a9b21ea304d5273c8be299c87Steve Block "// compiler remove the unused ones.\n" 1976ded16be15dd865a9b21ea304d5273c8be299c87Steve Block "// GRISU_CACHE_MAX_DISTANCE(i): generates the name for the maximum\n" 1986ded16be15dd865a9b21ea304d5273c8be299c87Steve Block "// binary exponent distance between all elements of a given cache.\n" 1996ded16be15dd865a9b21ea304d5273c8be299c87Steve Block "// GRISU_CACHE_OFFSET: is used as variable name for the decimal\n" 2006ded16be15dd865a9b21ea304d5273c8be299c87Steve Block "// exponent offset. It is equal to -cache[0].decimal_exponent.\n" 2016ded16be15dd865a9b21ea304d5273c8be299c87Steve Block "// GRISU_UINT64_C: used to construct 64-bit values in a platform\n" 2026ded16be15dd865a9b21ea304d5273c8be299c87Steve Block "// independent way. In order to encode 0x123456789ABCDEF0 the macro\n" 2036ded16be15dd865a9b21ea304d5273c8be299c87Steve Block "// will be invoked as follows: GRISU_UINT64_C(0x12345678,9ABCDEF0).\n") 2046ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (print) 2056ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (print-powers-reduced 1) 2066ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (print-powers-reduced 2) 2076ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (print-powers-reduced 3) 2086ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (print-powers-reduced 4) 2096ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (print-powers-reduced 5) 2106ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (print-powers-reduced 6) 2116ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (print-powers-reduced 7) 2126ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (print-powers-reduced 8) 2136ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (print-powers-reduced 9) 2146ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (print-powers-reduced 10) 2156ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (print-powers-reduced 11) 2166ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (print-powers-reduced 12) 2176ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (print-powers-reduced 13) 2186ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (print-powers-reduced 14) 2196ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (print-powers-reduced 15) 2206ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (print-powers-reduced 16) 2216ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (print-powers-reduced 17) 2226ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (print-powers-reduced 18) 2236ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (print-powers-reduced 19) 2246ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (print-powers-reduced 20) 2256ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (print "static const int GRISU_CACHE_OFFSET = " (- from) ";")) 2266ded16be15dd865a9b21ea304d5273c8be299c87Steve Block 2276ded16be15dd865a9b21ea304d5273c8be299c87Steve Block;;----------------main -------------------------------------------------------- 2286ded16be15dd865a9b21ea304d5273c8be299c87Steve Block(define *main-args* #f) 2296ded16be15dd865a9b21ea304d5273c8be299c87Steve Block(define *mantissa-size* #f) 2306ded16be15dd865a9b21ea304d5273c8be299c87Steve Block(define *dest* #f) 2316ded16be15dd865a9b21ea304d5273c8be299c87Steve Block(define *round* #f) 2326ded16be15dd865a9b21ea304d5273c8be299c87Steve Block(define *from* #f) 2336ded16be15dd865a9b21ea304d5273c8be299c87Steve Block(define *to* #f) 2346ded16be15dd865a9b21ea304d5273c8be299c87Steve Block 2356ded16be15dd865a9b21ea304d5273c8be299c87Steve Block(define (my-main args) 2366ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (set! *main-args* args) 2376ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (args-parse (cdr args) 2386ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (section "Help") 2396ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (("?") (args-parse-usage #f)) 2406ded16be15dd865a9b21ea304d5273c8be299c87Steve Block ((("-h" "--help") (help "?, -h, --help" "This help message")) 2416ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (args-parse-usage #f)) 2426ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (section "Misc") 2436ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (("-o" ?file (help "The output file")) 2446ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (set! *dest* file)) 2456ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (("--mantissa-size" ?size (help "Container-size in bits")) 2466ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (set! *mantissa-size* (string->number size))) 2476ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (("--round" ?direction (help "Round bignums (down, round or up)")) 2486ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (set! *round* (string->symbol direction))) 2496ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (("--from" ?from (help "start at 10^from")) 2506ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (set! *from* (string->number from))) 2516ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (("--to" ?to (help "go up to 10^to")) 2526ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (set! *to* (string->number to))) 2536ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (else 2546ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (print "Illegal argument `" else "'. Usage:") 2556ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (args-parse-usage #f))) 2566ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (when (not *from*) 2576ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (error "generate-ten-powers" 2586ded16be15dd865a9b21ea304d5273c8be299c87Steve Block "Missing from" 2596ded16be15dd865a9b21ea304d5273c8be299c87Steve Block #f)) 2606ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (when (not *to*) 2616ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (error "generate-ten-powers" 2626ded16be15dd865a9b21ea304d5273c8be299c87Steve Block "Missing to" 2636ded16be15dd865a9b21ea304d5273c8be299c87Steve Block #f)) 2646ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (when (not *mantissa-size*) 2656ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (error "generate-ten-powers" 2666ded16be15dd865a9b21ea304d5273c8be299c87Steve Block "Missing mantissa size" 2676ded16be15dd865a9b21ea304d5273c8be299c87Steve Block #f)) 2686ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (when (not (memv *round* '(up down round))) 2696ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (error "generate-ten-powers" 2706ded16be15dd865a9b21ea304d5273c8be299c87Steve Block "Missing round-method" 2716ded16be15dd865a9b21ea304d5273c8be299c87Steve Block *round*)) 2726ded16be15dd865a9b21ea304d5273c8be299c87Steve Block 2736ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (let ((dividers (generate-powers *from* *to* *mantissa-size*)) 2746ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (p (if (not *dest*) 2756ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (current-output-port) 2766ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (open-output-file *dest*)))) 2776ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (unwind-protect 2786ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (with-output-to-port p 2796ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (lambda () 2806ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (print-c dividers *from* *to* 2816ded16be15dd865a9b21ea304d5273c8be299c87Steve Block "GRISU_CACHE_STRUCT" "GRISU_CACHE_NAME" 2826ded16be15dd865a9b21ea304d5273c8be299c87Steve Block "GRISU_CACHE_MAX_DISTANCE" "GRISU_CACHE_OFFSET" 2836ded16be15dd865a9b21ea304d5273c8be299c87Steve Block "GRISU_UINT64_C" 2846ded16be15dd865a9b21ea304d5273c8be299c87Steve Block ))) 2856ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (if *dest* 2866ded16be15dd865a9b21ea304d5273c8be299c87Steve Block (close-output-port p))))) 287