|
1983
|
1 /* Numeric types for XEmacs using the GNU MP library.
|
|
|
2 Copyright (C) 2004 Jerry James.
|
|
|
3
|
|
|
4 This file is part of XEmacs.
|
|
|
5
|
|
|
6 XEmacs is free software; you can redistribute it and/or modify it
|
|
|
7 under the terms of the GNU General Public License as published by the
|
|
|
8 Free Software Foundation; either version 2, or (at your option) any
|
|
|
9 later version.
|
|
|
10
|
|
|
11 XEmacs is distributed in the hope that it will be useful, but WITHOUT
|
|
|
12 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
|
|
|
13 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
|
|
|
14 for more details.
|
|
|
15
|
|
|
16 You should have received a copy of the GNU General Public License
|
|
|
17 along with XEmacs; see the file COPYING. If not, write to
|
|
|
18 the Free Software Foundation, Inc., 59 Temple Place - Suite 330,
|
|
|
19 Boston, MA 02111-1307, USA. */
|
|
|
20
|
|
|
21 /* Synched up with: Not in FSF. */
|
|
|
22
|
|
|
23 #include <config.h>
|
|
|
24 #include <limits.h>
|
|
|
25 #include <math.h>
|
|
|
26 #include "lisp.h"
|
|
1995
|
27 #include "sysproc.h" /* For qxe_getpid */
|
|
1983
|
28
|
|
|
29 static mpf_t float_print_min, float_print_max;
|
|
|
30 gmp_randstate_t random_state;
|
|
|
31
|
|
|
32 CIbyte *
|
|
|
33 bigfloat_to_string(mpf_t f, int base)
|
|
|
34 {
|
|
|
35 mp_exp_t expt;
|
|
|
36 CIbyte *str = mpf_get_str (NULL, &expt, base, 0, f);
|
|
|
37 const int sign = mpf_sgn (f);
|
|
|
38 const int neg = (sign < 0) ? 1 : 0;
|
|
|
39 int len = strlen (str) + 1; /* Count the null terminator */
|
|
|
40
|
|
|
41 if (sign == 0 || (mpf_cmp (float_print_min, f) <= 0 &&
|
|
|
42 mpf_cmp (f, float_print_max) <= 0))
|
|
|
43 {
|
|
|
44 /* Move digits down to insert a radix point */
|
|
|
45 if (expt <= 0)
|
|
|
46 {
|
|
|
47 /* We need room for a radix point and leading zeroes */
|
|
|
48 const int space = -expt + 2;
|
|
|
49 XREALLOC_ARRAY (str, CIbyte, len + space);
|
|
|
50 memmove (&str[space + neg], &str[neg], len - neg);
|
|
|
51 memset (&str[neg], '0', space);
|
|
|
52 str[neg + 1] = '.';
|
|
|
53 len += space;
|
|
|
54 }
|
|
|
55 else
|
|
|
56 {
|
|
|
57 /* We just need room for a radix point */
|
|
|
58 XREALLOC_ARRAY (str, CIbyte, len + 1);
|
|
|
59 memmove (&str[expt + neg + 1], &str[expt + neg], len - (expt + neg));
|
|
|
60 str[expt + neg] = '.';
|
|
|
61 len++;
|
|
|
62 }
|
|
|
63 }
|
|
|
64 else
|
|
|
65 {
|
|
|
66 /* Computerized scientific notation */
|
|
|
67 /* We need room for a radix point, format identifier, and exponent */
|
|
|
68 const int space = (expt < 0)
|
|
2956
|
69 ? (int)(log ((double) (-expt)) / log ((double) base)) + 3
|
|
|
70 : (int)(log ((double) expt) / log ((double) base)) + 2;
|
|
1983
|
71 XREALLOC_ARRAY (str, CIbyte, len + space);
|
|
|
72 memmove (&str[neg + 2], &str[neg + 1], len - neg);
|
|
|
73 str[len + 1] = 'l';
|
|
|
74 sprintf (&str[len + 2], "%ld", expt);
|
|
|
75 }
|
|
|
76 return str;
|
|
|
77 }
|
|
|
78
|
|
|
79 /* We need the next two functions since GNU MP insists on giving us an extra
|
|
|
80 parameter. */
|
|
2286
|
81 static void *gmp_realloc (void *ptr, size_t UNUSED (old_size), size_t new_size)
|
|
1983
|
82 {
|
|
|
83 return xrealloc (ptr, new_size);
|
|
|
84 }
|
|
|
85
|
|
2286
|
86 static void gmp_free (void *ptr, size_t UNUSED (size))
|
|
1983
|
87 {
|
|
|
88 xfree (ptr, void *);
|
|
|
89 }
|
|
|
90
|
|
|
91 void
|
|
|
92 init_number_gmp ()
|
|
|
93 {
|
|
2367
|
94 mp_set_memory_functions ((void *(*) (size_t)) xmalloc, gmp_realloc,
|
|
|
95 gmp_free);
|
|
1983
|
96
|
|
|
97 /* The smallest number that is printed without exponents */
|
|
|
98 mpf_init_set_d (float_print_min, 0.001);
|
|
|
99
|
|
|
100 /* The largest number that is printed without exponents */
|
|
|
101 mpf_init_set_ui (float_print_max, 10000000UL);
|
|
|
102
|
|
|
103 /* Prepare the bignum/bigfloat random number generator */
|
|
|
104 gmp_randinit_default (random_state);
|
|
|
105 gmp_randseed_ui (random_state, qxe_getpid () + time (NULL));
|
|
|
106 }
|
