annotate src/floatfns.c @ 95080:5eb7d7d544f6

(calc-register-alist): New variable. (calc-set-register,calc-get-register,calc-copy-to-register) (calc-insert-register,calc-add-to-register,calc-append-to-register) (calc-prepend-to-register): New functions.
author Jay Belanger <jay.p.belanger@gmail.com>
date Sun, 18 May 2008 20:34:02 +0000
parents 8971ddf55736
children ddedcecb18ef
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
1 /* Primitive operations on floating point for GNU Emacs Lisp interpreter.
75227
e90d04cd455a Update copyright for years from Emacs 21 to present (mainly adding
Glenn Morris <rgm@gnu.org>
parents: 71978
diff changeset
2 Copyright (C) 1988, 1993, 1994, 1999, 2001, 2002, 2003, 2004,
79759
fc2bcd2a8aad Add 2008 to copyright years.
Glenn Morris <rgm@gnu.org>
parents: 78260
diff changeset
3 2005, 2006, 2007, 2008 Free Software Foundation, Inc.
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
4
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
5 This file is part of GNU Emacs.
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
6
94963
8971ddf55736 Switch to recommended form of GPLv3 permissions notice.
Glenn Morris <rgm@gnu.org>
parents: 87649
diff changeset
7 GNU Emacs is free software: you can redistribute it and/or modify
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
8 it under the terms of the GNU General Public License as published by
94963
8971ddf55736 Switch to recommended form of GPLv3 permissions notice.
Glenn Morris <rgm@gnu.org>
parents: 87649
diff changeset
9 the Free Software Foundation, either version 3 of the License, or
8971ddf55736 Switch to recommended form of GPLv3 permissions notice.
Glenn Morris <rgm@gnu.org>
parents: 87649
diff changeset
10 (at your option) any later version.
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
11
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
12 GNU Emacs is distributed in the hope that it will be useful,
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
13 but WITHOUT ANY WARRANTY; without even the implied warranty of
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
14 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
15 GNU General Public License for more details.
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
16
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
17 You should have received a copy of the GNU General Public License
94963
8971ddf55736 Switch to recommended form of GPLv3 permissions notice.
Glenn Morris <rgm@gnu.org>
parents: 87649
diff changeset
18 along with GNU Emacs. If not, see <http://www.gnu.org/licenses/>. */
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
19
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
20
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
21 /* ANSI C requires only these float functions:
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
22 acos, asin, atan, atan2, ceil, cos, cosh, exp, fabs, floor, fmod,
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
23 frexp, ldexp, log, log10, modf, pow, sin, sinh, sqrt, tan, tanh.
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
24
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
25 Define HAVE_INVERSE_HYPERBOLIC if you have acosh, asinh, and atanh.
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
26 Define HAVE_CBRT if you have cbrt.
19737
cad2a27d4451 (emacs_rint): Define this,
Richard M. Stallman <rms@gnu.org>
parents: 18627
diff changeset
27 Define HAVE_RINT if you have a working rint.
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
28 If you don't define these, then the appropriate routines will be simulated.
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
29
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
30 Define HAVE_MATHERR if on a system supporting the SysV matherr callback.
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
31 (This should happen automatically.)
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
32
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
33 Define FLOAT_CHECK_ERRNO if the float library routines set errno.
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
34 This has no effect if HAVE_MATHERR is defined.
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
35
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
36 Define FLOAT_CATCH_SIGILL if the float library routines signal SIGILL.
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
37 (What systems actually do this? Please let us know.)
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
38
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
39 Define FLOAT_CHECK_DOMAIN if the float library doesn't handle errors by
14036
621a575db6f7 Comment fixes.
Karl Heuer <kwzh@gnu.org>
parents: 13442
diff changeset
40 either setting errno, or signaling SIGFPE/SIGILL. Otherwise, domain and
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
41 range checking will happen before calling the float routines. This has
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
42 no effect if HAVE_MATHERR is defined (since matherr will be called when
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
43 a domain error occurs.)
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
44 */
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
45
26088
b7aa6ac26872 Add support for large files, 64-bit Solaris, system locale codings.
Paul Eggert <eggert@twinsun.com>
parents: 25645
diff changeset
46 #include <config.h>
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
47 #include <signal.h>
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
48 #include "lisp.h"
638
40b255f55df3 *** empty log message ***
Jim Blandy <jimb@redhat.com>
parents: 624
diff changeset
49 #include "syssignal.h"
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
50
20122
923e1f635ace No need to include <float.h> before "lisp.h",
Paul Eggert <eggert@twinsun.com>
parents: 19737
diff changeset
51 #if STDC_HEADERS
923e1f635ace No need to include <float.h> before "lisp.h",
Paul Eggert <eggert@twinsun.com>
parents: 19737
diff changeset
52 #include <float.h>
923e1f635ace No need to include <float.h> before "lisp.h",
Paul Eggert <eggert@twinsun.com>
parents: 19737
diff changeset
53 #endif
923e1f635ace No need to include <float.h> before "lisp.h",
Paul Eggert <eggert@twinsun.com>
parents: 19737
diff changeset
54
16786
8907c00c0cc6 <float.h>: Include if STDC_HEADERS.
Paul Eggert <eggert@twinsun.com>
parents: 16220
diff changeset
55 /* If IEEE_FLOATING_POINT isn't defined, default it from FLT_*. */
8907c00c0cc6 <float.h>: Include if STDC_HEADERS.
Paul Eggert <eggert@twinsun.com>
parents: 16220
diff changeset
56 #ifndef IEEE_FLOATING_POINT
8907c00c0cc6 <float.h>: Include if STDC_HEADERS.
Paul Eggert <eggert@twinsun.com>
parents: 16220
diff changeset
57 #if (FLT_RADIX == 2 && FLT_MANT_DIG == 24 \
8907c00c0cc6 <float.h>: Include if STDC_HEADERS.
Paul Eggert <eggert@twinsun.com>
parents: 16220
diff changeset
58 && FLT_MIN_EXP == -125 && FLT_MAX_EXP == 128)
8907c00c0cc6 <float.h>: Include if STDC_HEADERS.
Paul Eggert <eggert@twinsun.com>
parents: 16220
diff changeset
59 #define IEEE_FLOATING_POINT 1
8907c00c0cc6 <float.h>: Include if STDC_HEADERS.
Paul Eggert <eggert@twinsun.com>
parents: 16220
diff changeset
60 #else
8907c00c0cc6 <float.h>: Include if STDC_HEADERS.
Paul Eggert <eggert@twinsun.com>
parents: 16220
diff changeset
61 #define IEEE_FLOATING_POINT 0
8907c00c0cc6 <float.h>: Include if STDC_HEADERS.
Paul Eggert <eggert@twinsun.com>
parents: 16220
diff changeset
62 #endif
8907c00c0cc6 <float.h>: Include if STDC_HEADERS.
Paul Eggert <eggert@twinsun.com>
parents: 16220
diff changeset
63 #endif
8907c00c0cc6 <float.h>: Include if STDC_HEADERS.
Paul Eggert <eggert@twinsun.com>
parents: 16220
diff changeset
64
4843
136e32c763e3 [hpux 7] (_MAXLDBL, _NMAXLDBL): New macro definitions.
Richard M. Stallman <rms@gnu.org>
parents: 4696
diff changeset
65 /* Work around a problem that happens because math.h on hpux 7
136e32c763e3 [hpux 7] (_MAXLDBL, _NMAXLDBL): New macro definitions.
Richard M. Stallman <rms@gnu.org>
parents: 4696
diff changeset
66 defines two static variables--which, in Emacs, are not really static,
136e32c763e3 [hpux 7] (_MAXLDBL, _NMAXLDBL): New macro definitions.
Richard M. Stallman <rms@gnu.org>
parents: 4696
diff changeset
67 because `static' is defined as nothing. The problem is that they are
136e32c763e3 [hpux 7] (_MAXLDBL, _NMAXLDBL): New macro definitions.
Richard M. Stallman <rms@gnu.org>
parents: 4696
diff changeset
68 defined both here and in lread.c.
136e32c763e3 [hpux 7] (_MAXLDBL, _NMAXLDBL): New macro definitions.
Richard M. Stallman <rms@gnu.org>
parents: 4696
diff changeset
69 These macros prevent the name conflict. */
136e32c763e3 [hpux 7] (_MAXLDBL, _NMAXLDBL): New macro definitions.
Richard M. Stallman <rms@gnu.org>
parents: 4696
diff changeset
70 #if defined (HPUX) && !defined (HPUX8)
136e32c763e3 [hpux 7] (_MAXLDBL, _NMAXLDBL): New macro definitions.
Richard M. Stallman <rms@gnu.org>
parents: 4696
diff changeset
71 #define _MAXLDBL floatfns_maxldbl
136e32c763e3 [hpux 7] (_MAXLDBL, _NMAXLDBL): New macro definitions.
Richard M. Stallman <rms@gnu.org>
parents: 4696
diff changeset
72 #define _NMAXLDBL floatfns_nmaxldbl
136e32c763e3 [hpux 7] (_MAXLDBL, _NMAXLDBL): New macro definitions.
Richard M. Stallman <rms@gnu.org>
parents: 4696
diff changeset
73 #endif
136e32c763e3 [hpux 7] (_MAXLDBL, _NMAXLDBL): New macro definitions.
Richard M. Stallman <rms@gnu.org>
parents: 4696
diff changeset
74
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
75 #include <math.h>
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
76
4881
e53d4ae93675 Declare `logb' only if HAVE_LOGB is defined.
Brian Fox <bfox@gnu.org>
parents: 4843
diff changeset
77 /* This declaration is omitted on some systems, like Ultrix. */
7448
bf93ac2d8409 Don't declare logb if it is a macro.
Richard M. Stallman <rms@gnu.org>
parents: 7361
diff changeset
78 #if !defined (HPUX) && defined (HAVE_LOGB) && !defined (logb)
2205
c021f53fe7e5 * floatfns.c (logb): Add extern declaration for this.
Jim Blandy <jimb@redhat.com>
parents: 2129
diff changeset
79 extern double logb ();
7448
bf93ac2d8409 Don't declare logb if it is a macro.
Richard M. Stallman <rms@gnu.org>
parents: 7361
diff changeset
80 #endif /* not HPUX and HAVE_LOGB and no logb macro */
2205
c021f53fe7e5 * floatfns.c (logb): Add extern declaration for this.
Jim Blandy <jimb@redhat.com>
parents: 2129
diff changeset
81
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
82 #if defined(DOMAIN) && defined(SING) && defined(OVERFLOW)
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
83 /* If those are defined, then this is probably a `matherr' machine. */
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
84 # ifndef HAVE_MATHERR
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
85 # define HAVE_MATHERR
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
86 # endif
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
87 #endif
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
88
3027
7ed290bef028 Fix typo in last change.
Richard M. Stallman <rms@gnu.org>
parents: 3025
diff changeset
89 #ifdef NO_MATHERR
3025
d09f68ccd145 [NO_MATHERR]: Undef HAVE_MATHERR.
Richard M. Stallman <rms@gnu.org>
parents: 2961
diff changeset
90 #undef HAVE_MATHERR
d09f68ccd145 [NO_MATHERR]: Undef HAVE_MATHERR.
Richard M. Stallman <rms@gnu.org>
parents: 2961
diff changeset
91 #endif
d09f68ccd145 [NO_MATHERR]: Undef HAVE_MATHERR.
Richard M. Stallman <rms@gnu.org>
parents: 2961
diff changeset
92
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
93 #ifdef HAVE_MATHERR
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
94 # ifdef FLOAT_CHECK_ERRNO
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
95 # undef FLOAT_CHECK_ERRNO
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
96 # endif
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
97 # ifdef FLOAT_CHECK_DOMAIN
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
98 # undef FLOAT_CHECK_DOMAIN
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
99 # endif
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
100 #endif
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
101
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
102 #ifndef NO_FLOAT_CHECK_ERRNO
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
103 #define FLOAT_CHECK_ERRNO
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
104 #endif
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
105
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
106 #ifdef FLOAT_CHECK_ERRNO
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
107 # include <errno.h>
485
8c615e453683 *** empty log message ***
Jim Blandy <jimb@redhat.com>
parents: 102
diff changeset
108
31095
e19d38e14720 [USE_CRT_DLL]: Remove unnecessary extern, which screws
Andrew Innes <andrewi@gnu.org>
parents: 27727
diff changeset
109 #ifndef USE_CRT_DLL
485
8c615e453683 *** empty log message ***
Jim Blandy <jimb@redhat.com>
parents: 102
diff changeset
110 extern int errno;
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
111 #endif
31095
e19d38e14720 [USE_CRT_DLL]: Remove unnecessary extern, which screws
Andrew Innes <andrewi@gnu.org>
parents: 27727
diff changeset
112 #endif
485
8c615e453683 *** empty log message ***
Jim Blandy <jimb@redhat.com>
parents: 102
diff changeset
113
8c615e453683 *** empty log message ***
Jim Blandy <jimb@redhat.com>
parents: 102
diff changeset
114 /* Avoid traps on VMS from sinh and cosh.
8c615e453683 *** empty log message ***
Jim Blandy <jimb@redhat.com>
parents: 102
diff changeset
115 All the other functions set errno instead. */
8c615e453683 *** empty log message ***
Jim Blandy <jimb@redhat.com>
parents: 102
diff changeset
116
8c615e453683 *** empty log message ***
Jim Blandy <jimb@redhat.com>
parents: 102
diff changeset
117 #ifdef VMS
8c615e453683 *** empty log message ***
Jim Blandy <jimb@redhat.com>
parents: 102
diff changeset
118 #undef cosh
8c615e453683 *** empty log message ***
Jim Blandy <jimb@redhat.com>
parents: 102
diff changeset
119 #undef sinh
8c615e453683 *** empty log message ***
Jim Blandy <jimb@redhat.com>
parents: 102
diff changeset
120 #define cosh(x) ((exp(x)+exp(-x))*0.5)
8c615e453683 *** empty log message ***
Jim Blandy <jimb@redhat.com>
parents: 102
diff changeset
121 #define sinh(x) ((exp(x)-exp(-x))*0.5)
8c615e453683 *** empty log message ***
Jim Blandy <jimb@redhat.com>
parents: 102
diff changeset
122 #endif /* VMS */
8c615e453683 *** empty log message ***
Jim Blandy <jimb@redhat.com>
parents: 102
diff changeset
123
39814
a9c8422074c9 float_error() is defined iff FLOAT_CATCH_SIGILL is defined.
Pavel Janík <Pavel@Janik.cz>
parents: 39638
diff changeset
124 #ifdef FLOAT_CATCH_SIGILL
621
eca8812e61cd *** empty log message ***
Jim Blandy <jimb@redhat.com>
parents: 485
diff changeset
125 static SIGTYPE float_error ();
39814
a9c8422074c9 float_error() is defined iff FLOAT_CATCH_SIGILL is defined.
Pavel Janík <Pavel@Janik.cz>
parents: 39638
diff changeset
126 #endif
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
127
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
128 /* Nonzero while executing in floating point.
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
129 This tells float_error what to do. */
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
130
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
131 static int in_float;
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
132
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
133 /* If an argument is out of range for a mathematical function,
16207
f414d3a2fd58 Comment change.
Richard M. Stallman <rms@gnu.org>
parents: 15694
diff changeset
134 here is the actual argument value to use in the error message.
f414d3a2fd58 Comment change.
Richard M. Stallman <rms@gnu.org>
parents: 15694
diff changeset
135 These variables are used only across the floating point library call
f414d3a2fd58 Comment change.
Richard M. Stallman <rms@gnu.org>
parents: 15694
diff changeset
136 so there is no need to staticpro them. */
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
137
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
138 static Lisp_Object float_error_arg, float_error_arg2;
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
139
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
140 static char *float_error_fn_name;
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
141
485
8c615e453683 *** empty log message ***
Jim Blandy <jimb@redhat.com>
parents: 102
diff changeset
142 /* Evaluate the floating point expression D, recording NUM
8c615e453683 *** empty log message ***
Jim Blandy <jimb@redhat.com>
parents: 102
diff changeset
143 as the original argument for error messages.
8c615e453683 *** empty log message ***
Jim Blandy <jimb@redhat.com>
parents: 102
diff changeset
144 D is normally an assignment expression.
1918
699ae3079c09 * floatfns.c (Flogb): Always implement this by calling Flog, even
Jim Blandy <jimb@redhat.com>
parents: 1715
diff changeset
145 Handle errors which may result in signals or may set errno.
699ae3079c09 * floatfns.c (Flogb): Always implement this by calling Flog, even
Jim Blandy <jimb@redhat.com>
parents: 1715
diff changeset
146
699ae3079c09 * floatfns.c (Flogb): Always implement this by calling Flog, even
Jim Blandy <jimb@redhat.com>
parents: 1715
diff changeset
147 Note that float_error may be declared to return void, so you can't
699ae3079c09 * floatfns.c (Flogb): Always implement this by calling Flog, even
Jim Blandy <jimb@redhat.com>
parents: 1715
diff changeset
148 just cast the zero after the colon to (SIGTYPE) to make the types
699ae3079c09 * floatfns.c (Flogb): Always implement this by calling Flog, even
Jim Blandy <jimb@redhat.com>
parents: 1715
diff changeset
149 check properly. */
485
8c615e453683 *** empty log message ***
Jim Blandy <jimb@redhat.com>
parents: 102
diff changeset
150
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
151 #ifdef FLOAT_CHECK_ERRNO
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
152 #define IN_FLOAT(d, name, num) \
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
153 do { \
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
154 float_error_arg = num; \
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
155 float_error_fn_name = name; \
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
156 in_float = 1; errno = 0; (d); in_float = 0; \
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
157 switch (errno) { \
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
158 case 0: break; \
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
159 case EDOM: domain_error (float_error_fn_name, float_error_arg); \
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
160 case ERANGE: range_error (float_error_fn_name, float_error_arg); \
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
161 default: arith_error (float_error_fn_name, float_error_arg); \
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
162 } \
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
163 } while (0)
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
164 #define IN_FLOAT2(d, name, num, num2) \
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
165 do { \
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
166 float_error_arg = num; \
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
167 float_error_arg2 = num2; \
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
168 float_error_fn_name = name; \
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
169 in_float = 1; errno = 0; (d); in_float = 0; \
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
170 switch (errno) { \
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
171 case 0: break; \
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
172 case EDOM: domain_error (float_error_fn_name, float_error_arg); \
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
173 case ERANGE: range_error (float_error_fn_name, float_error_arg); \
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
174 default: arith_error (float_error_fn_name, float_error_arg); \
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
175 } \
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
176 } while (0)
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
177 #else
3407
d00aaf536bfd [!FLOAT_CHECK_ERRNO] (IN_FLOAT): New definition.
Richard M. Stallman <rms@gnu.org>
parents: 3094
diff changeset
178 #define IN_FLOAT(d, name, num) (in_float = 1, (d), in_float = 0)
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
179 #define IN_FLOAT2(d, name, num, num2) (in_float = 1, (d), in_float = 0)
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
180 #endif
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
181
6375
212dcd2c06e4 (FLOAT_TO_INT, FLOAT_TO_INT2, range_error2): New macros.
Karl Heuer <kwzh@gnu.org>
parents: 6359
diff changeset
182 /* Convert float to Lisp_Int if it fits, else signal a range error
212dcd2c06e4 (FLOAT_TO_INT, FLOAT_TO_INT2, range_error2): New macros.
Karl Heuer <kwzh@gnu.org>
parents: 6359
diff changeset
183 using the given arguments. */
212dcd2c06e4 (FLOAT_TO_INT, FLOAT_TO_INT2, range_error2): New macros.
Karl Heuer <kwzh@gnu.org>
parents: 6359
diff changeset
184 #define FLOAT_TO_INT(x, i, name, num) \
212dcd2c06e4 (FLOAT_TO_INT, FLOAT_TO_INT2, range_error2): New macros.
Karl Heuer <kwzh@gnu.org>
parents: 6359
diff changeset
185 do \
212dcd2c06e4 (FLOAT_TO_INT, FLOAT_TO_INT2, range_error2): New macros.
Karl Heuer <kwzh@gnu.org>
parents: 6359
diff changeset
186 { \
39638
911382113c16 (FLOAT_TO_INT2, FLOAT_TO_INT): Use FIXNUM_OVERFLOW_P.
Gerd Moellmann <gerd@gnu.org>
parents: 31095
diff changeset
187 if (FIXNUM_OVERFLOW_P (x)) \
6375
212dcd2c06e4 (FLOAT_TO_INT, FLOAT_TO_INT2, range_error2): New macros.
Karl Heuer <kwzh@gnu.org>
parents: 6359
diff changeset
188 range_error (name, num); \
9267
1b685d477c10 (FLOAT_TO_INT, FLOAT_TO_INT2, Fexpt, Flogb, Ffloor): Use new accessor macros
Karl Heuer <kwzh@gnu.org>
parents: 9129
diff changeset
189 XSETINT (i, (EMACS_INT)(x)); \
6375
212dcd2c06e4 (FLOAT_TO_INT, FLOAT_TO_INT2, range_error2): New macros.
Karl Heuer <kwzh@gnu.org>
parents: 6359
diff changeset
190 } \
212dcd2c06e4 (FLOAT_TO_INT, FLOAT_TO_INT2, range_error2): New macros.
Karl Heuer <kwzh@gnu.org>
parents: 6359
diff changeset
191 while (0)
212dcd2c06e4 (FLOAT_TO_INT, FLOAT_TO_INT2, range_error2): New macros.
Karl Heuer <kwzh@gnu.org>
parents: 6359
diff changeset
192 #define FLOAT_TO_INT2(x, i, name, num1, num2) \
212dcd2c06e4 (FLOAT_TO_INT, FLOAT_TO_INT2, range_error2): New macros.
Karl Heuer <kwzh@gnu.org>
parents: 6359
diff changeset
193 do \
212dcd2c06e4 (FLOAT_TO_INT, FLOAT_TO_INT2, range_error2): New macros.
Karl Heuer <kwzh@gnu.org>
parents: 6359
diff changeset
194 { \
39638
911382113c16 (FLOAT_TO_INT2, FLOAT_TO_INT): Use FIXNUM_OVERFLOW_P.
Gerd Moellmann <gerd@gnu.org>
parents: 31095
diff changeset
195 if (FIXNUM_OVERFLOW_P (x)) \
6375
212dcd2c06e4 (FLOAT_TO_INT, FLOAT_TO_INT2, range_error2): New macros.
Karl Heuer <kwzh@gnu.org>
parents: 6359
diff changeset
196 range_error2 (name, num1, num2); \
9267
1b685d477c10 (FLOAT_TO_INT, FLOAT_TO_INT2, Fexpt, Flogb, Ffloor): Use new accessor macros
Karl Heuer <kwzh@gnu.org>
parents: 9129
diff changeset
197 XSETINT (i, (EMACS_INT)(x)); \
6375
212dcd2c06e4 (FLOAT_TO_INT, FLOAT_TO_INT2, range_error2): New macros.
Karl Heuer <kwzh@gnu.org>
parents: 6359
diff changeset
198 } \
212dcd2c06e4 (FLOAT_TO_INT, FLOAT_TO_INT2, range_error2): New macros.
Karl Heuer <kwzh@gnu.org>
parents: 6359
diff changeset
199 while (0)
212dcd2c06e4 (FLOAT_TO_INT, FLOAT_TO_INT2, range_error2): New macros.
Karl Heuer <kwzh@gnu.org>
parents: 6359
diff changeset
200
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
201 #define arith_error(op,arg) \
71978
6b2f430f5341 (arith_error, range_error, domain_error): Use xsignal2.
Kim F. Storm <storm@cua.dk>
parents: 68651
diff changeset
202 xsignal2 (Qarith_error, build_string ((op)), (arg))
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
203 #define range_error(op,arg) \
71978
6b2f430f5341 (arith_error, range_error, domain_error): Use xsignal2.
Kim F. Storm <storm@cua.dk>
parents: 68651
diff changeset
204 xsignal2 (Qrange_error, build_string ((op)), (arg))
6375
212dcd2c06e4 (FLOAT_TO_INT, FLOAT_TO_INT2, range_error2): New macros.
Karl Heuer <kwzh@gnu.org>
parents: 6359
diff changeset
205 #define range_error2(op,a1,a2) \
71978
6b2f430f5341 (arith_error, range_error, domain_error): Use xsignal2.
Kim F. Storm <storm@cua.dk>
parents: 68651
diff changeset
206 xsignal3 (Qrange_error, build_string ((op)), (a1), (a2))
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
207 #define domain_error(op,arg) \
71978
6b2f430f5341 (arith_error, range_error, domain_error): Use xsignal2.
Kim F. Storm <storm@cua.dk>
parents: 68651
diff changeset
208 xsignal2 (Qdomain_error, build_string ((op)), (arg))
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
209 #define domain_error2(op,a1,a2) \
71978
6b2f430f5341 (arith_error, range_error, domain_error): Use xsignal2.
Kim F. Storm <storm@cua.dk>
parents: 68651
diff changeset
210 xsignal3 (Qdomain_error, build_string ((op)), (a1), (a2))
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
211
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
212 /* Extract a Lisp number as a `double', or signal an error. */
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
213
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
214 double
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
215 extract_float (num)
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
216 Lisp_Object num;
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
217 {
40656
cdfd4d09b79a Update usage of CHECK_ macros (remove unused second argument).
Pavel Janík <Pavel@Janik.cz>
parents: 39814
diff changeset
218 CHECK_NUMBER_OR_FLOAT (num);
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
219
9129
b2d1d925d5cc (extract_float, Fexpt, Fabs, Ffloat, Fceiling, Ffloor, Fround, Ftruncate): Use
Karl Heuer <kwzh@gnu.org>
parents: 8825
diff changeset
220 if (FLOATP (num))
25645
a14111a2a100 Use XCAR, XCDR, XFLOAT_DATA instead of explicit member access.
Ken Raeburn <raeburn@raeburn.org>
parents: 21514
diff changeset
221 return XFLOAT_DATA (num);
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
222 return (double) XINT (num);
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
223 }
683
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
224
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
225 /* Trig functions. */
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
226
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
227 DEFUN ("acos", Facos, Sacos, 1, 1, 0,
41001
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
228 doc: /* Return the inverse cosine of ARG. */)
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
229 (arg)
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
230 register Lisp_Object arg;
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
231 {
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
232 double d = extract_float (arg);
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
233 #ifdef FLOAT_CHECK_DOMAIN
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
234 if (d > 1.0 || d < -1.0)
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
235 domain_error ("acos", arg);
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
236 #endif
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
237 IN_FLOAT (d = acos (d), "acos", arg);
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
238 return make_float (d);
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
239 }
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
240
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
241 DEFUN ("asin", Fasin, Sasin, 1, 1, 0,
41001
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
242 doc: /* Return the inverse sine of ARG. */)
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
243 (arg)
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
244 register Lisp_Object arg;
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
245 {
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
246 double d = extract_float (arg);
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
247 #ifdef FLOAT_CHECK_DOMAIN
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
248 if (d > 1.0 || d < -1.0)
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
249 domain_error ("asin", arg);
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
250 #endif
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
251 IN_FLOAT (d = asin (d), "asin", arg);
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
252 return make_float (d);
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
253 }
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
254
43413
0f448bd1bf9a (Fatan): Accept an optional second arg and call atan2 if passed 2 args.
Eli Zaretskii <eliz@gnu.org>
parents: 43189
diff changeset
255 DEFUN ("atan", Fatan, Satan, 1, 2, 0,
0f448bd1bf9a (Fatan): Accept an optional second arg and call atan2 if passed 2 args.
Eli Zaretskii <eliz@gnu.org>
parents: 43189
diff changeset
256 doc: /* Return the inverse tangent of the arguments.
0f448bd1bf9a (Fatan): Accept an optional second arg and call atan2 if passed 2 args.
Eli Zaretskii <eliz@gnu.org>
parents: 43189
diff changeset
257 If only one argument Y is given, return the inverse tangent of Y.
0f448bd1bf9a (Fatan): Accept an optional second arg and call atan2 if passed 2 args.
Eli Zaretskii <eliz@gnu.org>
parents: 43189
diff changeset
258 If two arguments Y and X are given, return the inverse tangent of Y
0f448bd1bf9a (Fatan): Accept an optional second arg and call atan2 if passed 2 args.
Eli Zaretskii <eliz@gnu.org>
parents: 43189
diff changeset
259 divided by X, i.e. the angle in radians between the vector (X, Y)
0f448bd1bf9a (Fatan): Accept an optional second arg and call atan2 if passed 2 args.
Eli Zaretskii <eliz@gnu.org>
parents: 43189
diff changeset
260 and the x-axis. */)
0f448bd1bf9a (Fatan): Accept an optional second arg and call atan2 if passed 2 args.
Eli Zaretskii <eliz@gnu.org>
parents: 43189
diff changeset
261 (y, x)
0f448bd1bf9a (Fatan): Accept an optional second arg and call atan2 if passed 2 args.
Eli Zaretskii <eliz@gnu.org>
parents: 43189
diff changeset
262 register Lisp_Object y, x;
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
263 {
43413
0f448bd1bf9a (Fatan): Accept an optional second arg and call atan2 if passed 2 args.
Eli Zaretskii <eliz@gnu.org>
parents: 43189
diff changeset
264 double d = extract_float (y);
0f448bd1bf9a (Fatan): Accept an optional second arg and call atan2 if passed 2 args.
Eli Zaretskii <eliz@gnu.org>
parents: 43189
diff changeset
265
0f448bd1bf9a (Fatan): Accept an optional second arg and call atan2 if passed 2 args.
Eli Zaretskii <eliz@gnu.org>
parents: 43189
diff changeset
266 if (NILP (x))
0f448bd1bf9a (Fatan): Accept an optional second arg and call atan2 if passed 2 args.
Eli Zaretskii <eliz@gnu.org>
parents: 43189
diff changeset
267 IN_FLOAT (d = atan (d), "atan", y);
0f448bd1bf9a (Fatan): Accept an optional second arg and call atan2 if passed 2 args.
Eli Zaretskii <eliz@gnu.org>
parents: 43189
diff changeset
268 else
0f448bd1bf9a (Fatan): Accept an optional second arg and call atan2 if passed 2 args.
Eli Zaretskii <eliz@gnu.org>
parents: 43189
diff changeset
269 {
0f448bd1bf9a (Fatan): Accept an optional second arg and call atan2 if passed 2 args.
Eli Zaretskii <eliz@gnu.org>
parents: 43189
diff changeset
270 double d2 = extract_float (x);
0f448bd1bf9a (Fatan): Accept an optional second arg and call atan2 if passed 2 args.
Eli Zaretskii <eliz@gnu.org>
parents: 43189
diff changeset
271
0f448bd1bf9a (Fatan): Accept an optional second arg and call atan2 if passed 2 args.
Eli Zaretskii <eliz@gnu.org>
parents: 43189
diff changeset
272 IN_FLOAT2 (d = atan2 (d, d2), "atan", y, x);
0f448bd1bf9a (Fatan): Accept an optional second arg and call atan2 if passed 2 args.
Eli Zaretskii <eliz@gnu.org>
parents: 43189
diff changeset
273 }
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
274 return make_float (d);
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
275 }
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
276
683
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
277 DEFUN ("cos", Fcos, Scos, 1, 1, 0,
41001
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
278 doc: /* Return the cosine of ARG. */)
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
279 (arg)
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
280 register Lisp_Object arg;
683
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
281 {
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
282 double d = extract_float (arg);
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
283 IN_FLOAT (d = cos (d), "cos", arg);
683
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
284 return make_float (d);
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
285 }
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
286
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
287 DEFUN ("sin", Fsin, Ssin, 1, 1, 0,
41001
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
288 doc: /* Return the sine of ARG. */)
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
289 (arg)
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
290 register Lisp_Object arg;
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
291 {
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
292 double d = extract_float (arg);
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
293 IN_FLOAT (d = sin (d), "sin", arg);
683
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
294 return make_float (d);
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
295 }
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
296
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
297 DEFUN ("tan", Ftan, Stan, 1, 1, 0,
41001
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
298 doc: /* Return the tangent of ARG. */)
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
299 (arg)
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
300 register Lisp_Object arg;
683
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
301 {
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
302 double d = extract_float (arg);
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
303 double c = cos (d);
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
304 #ifdef FLOAT_CHECK_DOMAIN
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
305 if (c == 0.0)
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
306 domain_error ("tan", arg);
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
307 #endif
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
308 IN_FLOAT (d = sin (d) / c, "tan", arg);
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
309 return make_float (d);
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
310 }
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
311
683
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
312 #if 0 /* Leave these out unless we find there's a reason for them. */
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
313
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
314 DEFUN ("bessel-j0", Fbessel_j0, Sbessel_j0, 1, 1, 0,
41001
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
315 doc: /* Return the bessel function j0 of ARG. */)
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
316 (arg)
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
317 register Lisp_Object arg;
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
318 {
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
319 double d = extract_float (arg);
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
320 IN_FLOAT (d = j0 (d), "bessel-j0", arg);
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
321 return make_float (d);
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
322 }
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
323
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
324 DEFUN ("bessel-j1", Fbessel_j1, Sbessel_j1, 1, 1, 0,
41001
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
325 doc: /* Return the bessel function j1 of ARG. */)
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
326 (arg)
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
327 register Lisp_Object arg;
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
328 {
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
329 double d = extract_float (arg);
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
330 IN_FLOAT (d = j1 (d), "bessel-j1", arg);
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
331 return make_float (d);
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
332 }
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
333
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
334 DEFUN ("bessel-jn", Fbessel_jn, Sbessel_jn, 2, 2, 0,
41001
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
335 doc: /* Return the order N bessel function output jn of ARG.
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
336 The first arg (the order) is truncated to an integer. */)
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
337 (n, arg)
14076
b339e9aaea4b (Fbessel_jn, Fbessel_yn): Harmonize arguments with documentation.
Erik Naggum <erik@naggum.no>
parents: 14036
diff changeset
338 register Lisp_Object n, arg;
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
339 {
14076
b339e9aaea4b (Fbessel_jn, Fbessel_yn): Harmonize arguments with documentation.
Erik Naggum <erik@naggum.no>
parents: 14036
diff changeset
340 int i1 = extract_float (n);
b339e9aaea4b (Fbessel_jn, Fbessel_yn): Harmonize arguments with documentation.
Erik Naggum <erik@naggum.no>
parents: 14036
diff changeset
341 double f2 = extract_float (arg);
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
342
14076
b339e9aaea4b (Fbessel_jn, Fbessel_yn): Harmonize arguments with documentation.
Erik Naggum <erik@naggum.no>
parents: 14036
diff changeset
343 IN_FLOAT (f2 = jn (i1, f2), "bessel-jn", n);
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
344 return make_float (f2);
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
345 }
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
346
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
347 DEFUN ("bessel-y0", Fbessel_y0, Sbessel_y0, 1, 1, 0,
41001
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
348 doc: /* Return the bessel function y0 of ARG. */)
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
349 (arg)
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
350 register Lisp_Object arg;
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
351 {
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
352 double d = extract_float (arg);
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
353 IN_FLOAT (d = y0 (d), "bessel-y0", arg);
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
354 return make_float (d);
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
355 }
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
356
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
357 DEFUN ("bessel-y1", Fbessel_y1, Sbessel_y1, 1, 1, 0,
41001
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
358 doc: /* Return the bessel function y1 of ARG. */)
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
359 (arg)
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
360 register Lisp_Object arg;
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
361 {
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
362 double d = extract_float (arg);
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
363 IN_FLOAT (d = y1 (d), "bessel-y0", arg);
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
364 return make_float (d);
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
365 }
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
366
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
367 DEFUN ("bessel-yn", Fbessel_yn, Sbessel_yn, 2, 2, 0,
41001
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
368 doc: /* Return the order N bessel function output yn of ARG.
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
369 The first arg (the order) is truncated to an integer. */)
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
370 (n, arg)
14076
b339e9aaea4b (Fbessel_jn, Fbessel_yn): Harmonize arguments with documentation.
Erik Naggum <erik@naggum.no>
parents: 14036
diff changeset
371 register Lisp_Object n, arg;
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
372 {
14076
b339e9aaea4b (Fbessel_jn, Fbessel_yn): Harmonize arguments with documentation.
Erik Naggum <erik@naggum.no>
parents: 14036
diff changeset
373 int i1 = extract_float (n);
b339e9aaea4b (Fbessel_jn, Fbessel_yn): Harmonize arguments with documentation.
Erik Naggum <erik@naggum.no>
parents: 14036
diff changeset
374 double f2 = extract_float (arg);
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
375
14076
b339e9aaea4b (Fbessel_jn, Fbessel_yn): Harmonize arguments with documentation.
Erik Naggum <erik@naggum.no>
parents: 14036
diff changeset
376 IN_FLOAT (f2 = yn (i1, f2), "bessel-yn", n);
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
377 return make_float (f2);
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
378 }
683
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
379
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
380 #endif
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
381
683
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
382 #if 0 /* Leave these out unless we see they are worth having. */
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
383
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
384 DEFUN ("erf", Ferf, Serf, 1, 1, 0,
41001
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
385 doc: /* Return the mathematical error function of ARG. */)
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
386 (arg)
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
387 register Lisp_Object arg;
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
388 {
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
389 double d = extract_float (arg);
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
390 IN_FLOAT (d = erf (d), "erf", arg);
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
391 return make_float (d);
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
392 }
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
393
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
394 DEFUN ("erfc", Ferfc, Serfc, 1, 1, 0,
41001
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
395 doc: /* Return the complementary error function of ARG. */)
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
396 (arg)
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
397 register Lisp_Object arg;
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
398 {
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
399 double d = extract_float (arg);
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
400 IN_FLOAT (d = erfc (d), "erfc", arg);
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
401 return make_float (d);
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
402 }
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
403
683
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
404 DEFUN ("log-gamma", Flog_gamma, Slog_gamma, 1, 1, 0,
41001
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
405 doc: /* Return the log gamma of ARG. */)
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
406 (arg)
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
407 register Lisp_Object arg;
683
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
408 {
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
409 double d = extract_float (arg);
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
410 IN_FLOAT (d = lgamma (d), "log-gamma", arg);
683
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
411 return make_float (d);
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
412 }
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
413
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
414 DEFUN ("cube-root", Fcube_root, Scube_root, 1, 1, 0,
41001
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
415 doc: /* Return the cube root of ARG. */)
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
416 (arg)
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
417 register Lisp_Object arg;
683
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
418 {
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
419 double d = extract_float (arg);
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
420 #ifdef HAVE_CBRT
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
421 IN_FLOAT (d = cbrt (d), "cube-root", arg);
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
422 #else
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
423 if (d >= 0.0)
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
424 IN_FLOAT (d = pow (d, 1.0/3.0), "cube-root", arg);
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
425 else
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
426 IN_FLOAT (d = -pow (-d, 1.0/3.0), "cube-root", arg);
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
427 #endif
683
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
428 return make_float (d);
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
429 }
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
430
694
714b8017cc6b entered into RCS
Richard M. Stallman <rms@gnu.org>
parents: 683
diff changeset
431 #endif
714b8017cc6b entered into RCS
Richard M. Stallman <rms@gnu.org>
parents: 683
diff changeset
432
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
433 DEFUN ("exp", Fexp, Sexp, 1, 1, 0,
41001
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
434 doc: /* Return the exponential base e of ARG. */)
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
435 (arg)
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
436 register Lisp_Object arg;
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
437 {
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
438 double d = extract_float (arg);
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
439 #ifdef FLOAT_CHECK_DOMAIN
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
440 if (d > 709.7827) /* Assume IEEE doubles here */
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
441 range_error ("exp", arg);
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
442 else if (d < -709.0)
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
443 return make_float (0.0);
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
444 else
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
445 #endif
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
446 IN_FLOAT (d = exp (d), "exp", arg);
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
447 return make_float (d);
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
448 }
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
449
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
450 DEFUN ("expt", Fexpt, Sexpt, 2, 2, 0,
41001
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
451 doc: /* Return the exponential ARG1 ** ARG2. */)
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
452 (arg1, arg2)
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
453 register Lisp_Object arg1, arg2;
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
454 {
84998
259a9fdd1450 (Fexpt): Manually check for overflows, so that a power
Michaël Cadilhac <michael.cadilhac@lrde.org>
parents: 78260
diff changeset
455 double f1, f2, f3;
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
456
40656
cdfd4d09b79a Update usage of CHECK_ macros (remove unused second argument).
Pavel Janík <Pavel@Janik.cz>
parents: 39814
diff changeset
457 CHECK_NUMBER_OR_FLOAT (arg1);
cdfd4d09b79a Update usage of CHECK_ macros (remove unused second argument).
Pavel Janík <Pavel@Janik.cz>
parents: 39814
diff changeset
458 CHECK_NUMBER_OR_FLOAT (arg2);
9129
b2d1d925d5cc (extract_float, Fexpt, Fabs, Ffloat, Fceiling, Ffloor, Fround, Ftruncate): Use
Karl Heuer <kwzh@gnu.org>
parents: 8825
diff changeset
459 if (INTEGERP (arg1) /* common lisp spec */
64810
3529fe766935 (Fexpt): Use floats for negative exponent.
Thien-Thi Nguyen <ttn@gnuvola.org>
parents: 64778
diff changeset
460 && INTEGERP (arg2) /* don't promote, if both are ints, and */
3529fe766935 (Fexpt): Use floats for negative exponent.
Thien-Thi Nguyen <ttn@gnuvola.org>
parents: 64778
diff changeset
461 && 0 <= XINT (arg2)) /* we are sure the result is not fractional */
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
462 { /* this can be improved by pre-calculating */
12537
476296adb950 (Fexpt): Use EMACS_INT for integer calculation.
Karl Heuer <kwzh@gnu.org>
parents: 11859
diff changeset
463 EMACS_INT acc, x, y; /* some binary powers of x then accumulating */
3673
8ab0a7453577 (Fexpt): New local `val' for making integer to return.
Richard M. Stallman <rms@gnu.org>
parents: 3591
diff changeset
464 Lisp_Object val;
8ab0a7453577 (Fexpt): New local `val' for making integer to return.
Richard M. Stallman <rms@gnu.org>
parents: 3591
diff changeset
465
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
466 x = XINT (arg1);
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
467 y = XINT (arg2);
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
468 acc = 1;
49600
23a1cea22d13 Trailing whitespace deleted.
Juanma Barranquero <lekktu@gmail.com>
parents: 43413
diff changeset
469
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
470 if (y < 0)
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
471 {
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
472 if (x == 1)
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
473 acc = 1;
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
474 else if (x == -1)
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
475 acc = (y & 1) ? -1 : 1;
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
476 else
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
477 acc = 0;
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
478 }
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
479 else
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
480 {
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
481 while (y > 0)
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
482 {
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
483 if (y & 1)
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
484 acc *= x;
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
485 x *= x;
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
486 y = (unsigned)y >> 1;
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
487 }
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
488 }
9267
1b685d477c10 (FLOAT_TO_INT, FLOAT_TO_INT2, Fexpt, Flogb, Ffloor): Use new accessor macros
Karl Heuer <kwzh@gnu.org>
parents: 9129
diff changeset
489 XSETINT (val, acc);
3673
8ab0a7453577 (Fexpt): New local `val' for making integer to return.
Richard M. Stallman <rms@gnu.org>
parents: 3591
diff changeset
490 return val;
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
491 }
25645
a14111a2a100 Use XCAR, XCDR, XFLOAT_DATA instead of explicit member access.
Ken Raeburn <raeburn@raeburn.org>
parents: 21514
diff changeset
492 f1 = FLOATP (arg1) ? XFLOAT_DATA (arg1) : XINT (arg1);
a14111a2a100 Use XCAR, XCDR, XFLOAT_DATA instead of explicit member access.
Ken Raeburn <raeburn@raeburn.org>
parents: 21514
diff changeset
493 f2 = FLOATP (arg2) ? XFLOAT_DATA (arg2) : XINT (arg2);
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
494 /* Really should check for overflow, too */
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
495 if (f1 == 0.0 && f2 == 0.0)
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
496 f1 = 1.0;
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
497 #ifdef FLOAT_CHECK_DOMAIN
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
498 else if ((f1 == 0.0 && f2 < 0.0) || (f1 < 0 && f2 != floor(f2)))
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
499 domain_error2 ("expt", arg1, arg2);
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
500 #endif
84998
259a9fdd1450 (Fexpt): Manually check for overflows, so that a power
Michaël Cadilhac <michael.cadilhac@lrde.org>
parents: 78260
diff changeset
501 IN_FLOAT2 (f3 = pow (f1, f2), "expt", arg1, arg2);
259a9fdd1450 (Fexpt): Manually check for overflows, so that a power
Michaël Cadilhac <michael.cadilhac@lrde.org>
parents: 78260
diff changeset
502 /* Check for overflow in the result. */
259a9fdd1450 (Fexpt): Manually check for overflows, so that a power
Michaël Cadilhac <michael.cadilhac@lrde.org>
parents: 78260
diff changeset
503 if (f1 != 0.0 && f3 == 0.0)
259a9fdd1450 (Fexpt): Manually check for overflows, so that a power
Michaël Cadilhac <michael.cadilhac@lrde.org>
parents: 78260
diff changeset
504 range_error ("expt", arg1);
259a9fdd1450 (Fexpt): Manually check for overflows, so that a power
Michaël Cadilhac <michael.cadilhac@lrde.org>
parents: 78260
diff changeset
505 return make_float (f3);
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
506 }
683
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
507
1005
70ed307d9047 * floatfns.c (Fexpm1, Flog1p): Function removed; it's not widely
Jim Blandy <jimb@redhat.com>
parents: 694
diff changeset
508 DEFUN ("log", Flog, Slog, 1, 2, 0,
41001
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
509 doc: /* Return the natural logarithm of ARG.
67496
199bdb85a35f (Flog): Doc fix.
Richard M. Stallman <rms@gnu.org>
parents: 64810
diff changeset
510 If the optional argument BASE is given, return log ARG using that base. */)
41001
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
511 (arg, base)
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
512 register Lisp_Object arg, base;
683
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
513 {
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
514 double d = extract_float (arg);
1005
70ed307d9047 * floatfns.c (Fexpm1, Flog1p): Function removed; it's not widely
Jim Blandy <jimb@redhat.com>
parents: 694
diff changeset
515
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
516 #ifdef FLOAT_CHECK_DOMAIN
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
517 if (d <= 0.0)
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
518 domain_error2 ("log", arg, base);
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
519 #endif
1005
70ed307d9047 * floatfns.c (Fexpm1, Flog1p): Function removed; it's not widely
Jim Blandy <jimb@redhat.com>
parents: 694
diff changeset
520 if (NILP (base))
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
521 IN_FLOAT (d = log (d), "log", arg);
1005
70ed307d9047 * floatfns.c (Fexpm1, Flog1p): Function removed; it's not widely
Jim Blandy <jimb@redhat.com>
parents: 694
diff changeset
522 else
70ed307d9047 * floatfns.c (Fexpm1, Flog1p): Function removed; it's not widely
Jim Blandy <jimb@redhat.com>
parents: 694
diff changeset
523 {
70ed307d9047 * floatfns.c (Fexpm1, Flog1p): Function removed; it's not widely
Jim Blandy <jimb@redhat.com>
parents: 694
diff changeset
524 double b = extract_float (base);
70ed307d9047 * floatfns.c (Fexpm1, Flog1p): Function removed; it's not widely
Jim Blandy <jimb@redhat.com>
parents: 694
diff changeset
525
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
526 #ifdef FLOAT_CHECK_DOMAIN
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
527 if (b <= 0.0 || b == 1.0)
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
528 domain_error2 ("log", arg, base);
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
529 #endif
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
530 if (b == 10.0)
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
531 IN_FLOAT2 (d = log10 (d), "log", arg, base);
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
532 else
3407
d00aaf536bfd [!FLOAT_CHECK_ERRNO] (IN_FLOAT): New definition.
Richard M. Stallman <rms@gnu.org>
parents: 3094
diff changeset
533 IN_FLOAT2 (d = log (d) / log (b), "log", arg, base);
1005
70ed307d9047 * floatfns.c (Fexpm1, Flog1p): Function removed; it's not widely
Jim Blandy <jimb@redhat.com>
parents: 694
diff changeset
534 }
683
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
535 return make_float (d);
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
536 }
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
537
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
538 DEFUN ("log10", Flog10, Slog10, 1, 1, 0,
41001
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
539 doc: /* Return the logarithm base 10 of ARG. */)
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
540 (arg)
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
541 register Lisp_Object arg;
683
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
542 {
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
543 double d = extract_float (arg);
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
544 #ifdef FLOAT_CHECK_DOMAIN
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
545 if (d <= 0.0)
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
546 domain_error ("log10", arg);
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
547 #endif
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
548 IN_FLOAT (d = log10 (d), "log10", arg);
683
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
549 return make_float (d);
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
550 }
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
551
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
552 DEFUN ("sqrt", Fsqrt, Ssqrt, 1, 1, 0,
41001
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
553 doc: /* Return the square root of ARG. */)
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
554 (arg)
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
555 register Lisp_Object arg;
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
556 {
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
557 double d = extract_float (arg);
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
558 #ifdef FLOAT_CHECK_DOMAIN
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
559 if (d < 0.0)
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
560 domain_error ("sqrt", arg);
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
561 #endif
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
562 IN_FLOAT (d = sqrt (d), "sqrt", arg);
683
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
563 return make_float (d);
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
564 }
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
565
694
714b8017cc6b entered into RCS
Richard M. Stallman <rms@gnu.org>
parents: 683
diff changeset
566 #if 0 /* Not clearly worth adding. */
683
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
567
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
568 DEFUN ("acosh", Facosh, Sacosh, 1, 1, 0,
41001
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
569 doc: /* Return the inverse hyperbolic cosine of ARG. */)
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
570 (arg)
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
571 register Lisp_Object arg;
683
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
572 {
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
573 double d = extract_float (arg);
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
574 #ifdef FLOAT_CHECK_DOMAIN
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
575 if (d < 1.0)
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
576 domain_error ("acosh", arg);
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
577 #endif
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
578 #ifdef HAVE_INVERSE_HYPERBOLIC
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
579 IN_FLOAT (d = acosh (d), "acosh", arg);
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
580 #else
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
581 IN_FLOAT (d = log (d + sqrt (d*d - 1.0)), "acosh", arg);
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
582 #endif
683
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
583 return make_float (d);
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
584 }
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
585
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
586 DEFUN ("asinh", Fasinh, Sasinh, 1, 1, 0,
41001
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
587 doc: /* Return the inverse hyperbolic sine of ARG. */)
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
588 (arg)
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
589 register Lisp_Object arg;
683
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
590 {
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
591 double d = extract_float (arg);
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
592 #ifdef HAVE_INVERSE_HYPERBOLIC
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
593 IN_FLOAT (d = asinh (d), "asinh", arg);
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
594 #else
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
595 IN_FLOAT (d = log (d + sqrt (d*d + 1.0)), "asinh", arg);
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
596 #endif
683
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
597 return make_float (d);
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
598 }
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
599
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
600 DEFUN ("atanh", Fatanh, Satanh, 1, 1, 0,
41001
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
601 doc: /* Return the inverse hyperbolic tangent of ARG. */)
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
602 (arg)
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
603 register Lisp_Object arg;
683
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
604 {
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
605 double d = extract_float (arg);
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
606 #ifdef FLOAT_CHECK_DOMAIN
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
607 if (d >= 1.0 || d <= -1.0)
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
608 domain_error ("atanh", arg);
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
609 #endif
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
610 #ifdef HAVE_INVERSE_HYPERBOLIC
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
611 IN_FLOAT (d = atanh (d), "atanh", arg);
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
612 #else
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
613 IN_FLOAT (d = 0.5 * log ((1.0 + d) / (1.0 - d)), "atanh", arg);
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
614 #endif
683
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
615 return make_float (d);
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
616 }
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
617
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
618 DEFUN ("cosh", Fcosh, Scosh, 1, 1, 0,
41001
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
619 doc: /* Return the hyperbolic cosine of ARG. */)
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
620 (arg)
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
621 register Lisp_Object arg;
683
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
622 {
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
623 double d = extract_float (arg);
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
624 #ifdef FLOAT_CHECK_DOMAIN
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
625 if (d > 710.0 || d < -710.0)
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
626 range_error ("cosh", arg);
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
627 #endif
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
628 IN_FLOAT (d = cosh (d), "cosh", arg);
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
629 return make_float (d);
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
630 }
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
631
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
632 DEFUN ("sinh", Fsinh, Ssinh, 1, 1, 0,
41001
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
633 doc: /* Return the hyperbolic sine of ARG. */)
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
634 (arg)
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
635 register Lisp_Object arg;
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
636 {
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
637 double d = extract_float (arg);
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
638 #ifdef FLOAT_CHECK_DOMAIN
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
639 if (d > 710.0 || d < -710.0)
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
640 range_error ("sinh", arg);
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
641 #endif
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
642 IN_FLOAT (d = sinh (d), "sinh", arg);
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
643 return make_float (d);
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
644 }
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
645
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
646 DEFUN ("tanh", Ftanh, Stanh, 1, 1, 0,
41001
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
647 doc: /* Return the hyperbolic tangent of ARG. */)
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
648 (arg)
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
649 register Lisp_Object arg;
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
650 {
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
651 double d = extract_float (arg);
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
652 IN_FLOAT (d = tanh (d), "tanh", arg);
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
653 return make_float (d);
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
654 }
683
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
655 #endif
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
656
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
657 DEFUN ("abs", Fabs, Sabs, 1, 1, 0,
41001
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
658 doc: /* Return the absolute value of ARG. */)
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
659 (arg)
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
660 register Lisp_Object arg;
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
661 {
40656
cdfd4d09b79a Update usage of CHECK_ macros (remove unused second argument).
Pavel Janík <Pavel@Janik.cz>
parents: 39814
diff changeset
662 CHECK_NUMBER_OR_FLOAT (arg);
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
663
9129
b2d1d925d5cc (extract_float, Fexpt, Fabs, Ffloat, Fceiling, Ffloor, Fround, Ftruncate): Use
Karl Heuer <kwzh@gnu.org>
parents: 8825
diff changeset
664 if (FLOATP (arg))
25645
a14111a2a100 Use XCAR, XCDR, XFLOAT_DATA instead of explicit member access.
Ken Raeburn <raeburn@raeburn.org>
parents: 21514
diff changeset
665 IN_FLOAT (arg = make_float (fabs (XFLOAT_DATA (arg))), "abs", arg);
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
666 else if (XINT (arg) < 0)
9365
0b431684f97b (Fabs): Don't use XFASTINT when negative.
Karl Heuer <kwzh@gnu.org>
parents: 9267
diff changeset
667 XSETINT (arg, - XINT (arg));
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
668
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
669 return arg;
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
670 }
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
671
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
672 DEFUN ("float", Ffloat, Sfloat, 1, 1, 0,
41001
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
673 doc: /* Return the floating point number equal to ARG. */)
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
674 (arg)
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
675 register Lisp_Object arg;
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
676 {
40656
cdfd4d09b79a Update usage of CHECK_ macros (remove unused second argument).
Pavel Janík <Pavel@Janik.cz>
parents: 39814
diff changeset
677 CHECK_NUMBER_OR_FLOAT (arg);
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
678
9129
b2d1d925d5cc (extract_float, Fexpt, Fabs, Ffloat, Fceiling, Ffloor, Fround, Ftruncate): Use
Karl Heuer <kwzh@gnu.org>
parents: 8825
diff changeset
679 if (INTEGERP (arg))
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
680 return make_float ((double) XINT (arg));
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
681 else /* give 'em the same float back */
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
682 return arg;
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
683 }
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
684
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
685 DEFUN ("logb", Flogb, Slogb, 1, 1, 0,
41001
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
686 doc: /* Returns largest integer <= the base 2 log of the magnitude of ARG.
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
687 This is the same as the exponent of a float. */)
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
688 (arg)
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
689 Lisp_Object arg;
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
690 {
2119
4077ef8ad483 * floatfns.c (Flogb): Undo the change of Feb 22.
Jim Blandy <jimb@redhat.com>
parents: 2094
diff changeset
691 Lisp_Object val;
11265
b2291296ec87 (Flogb): Use EMACS_INT for `value'.
Richard M. Stallman <rms@gnu.org>
parents: 9365
diff changeset
692 EMACS_INT value;
2129
6741f5f8ed54 (Flogb): Fix arg names. Don't confuse Lisp_Object with integer.
Richard M. Stallman <rms@gnu.org>
parents: 2128
diff changeset
693 double f = extract_float (arg);
2119
4077ef8ad483 * floatfns.c (Flogb): Undo the change of Feb 22.
Jim Blandy <jimb@redhat.com>
parents: 2094
diff changeset
694
6359
800c035273e9 (Flogb): Check for 0.0. Emulate logb if needed.
Karl Heuer <kwzh@gnu.org>
parents: 6314
diff changeset
695 if (f == 0.0)
53002
3e336113348d (Flogb): Don't use VALMASK.
Stefan Monnier <monnier@iro.umontreal.ca>
parents: 52401
diff changeset
696 value = MOST_NEGATIVE_FIXNUM;
6359
800c035273e9 (Flogb): Check for 0.0. Emulate logb if needed.
Karl Heuer <kwzh@gnu.org>
parents: 6314
diff changeset
697 else
800c035273e9 (Flogb): Check for 0.0. Emulate logb if needed.
Karl Heuer <kwzh@gnu.org>
parents: 6314
diff changeset
698 {
4501
9352d7d021c2 * config.h.in (HAVE_LOGB, HAVE_FREXP): Add #undefs for the
Jim Blandy <jimb@redhat.com>
parents: 3673
diff changeset
699 #ifdef HAVE_LOGB
6359
800c035273e9 (Flogb): Check for 0.0. Emulate logb if needed.
Karl Heuer <kwzh@gnu.org>
parents: 6314
diff changeset
700 IN_FLOAT (value = logb (f), "logb", arg);
4501
9352d7d021c2 * config.h.in (HAVE_LOGB, HAVE_FREXP): Add #undefs for the
Jim Blandy <jimb@redhat.com>
parents: 3673
diff changeset
701 #else
9352d7d021c2 * config.h.in (HAVE_LOGB, HAVE_FREXP): Add #undefs for the
Jim Blandy <jimb@redhat.com>
parents: 3673
diff changeset
702 #ifdef HAVE_FREXP
11859
8d57babf4b36 (Flogb): frexp needs a pointer to int, not EMACS_INT.
Karl Heuer <kwzh@gnu.org>
parents: 11265
diff changeset
703 int ivalue;
8d57babf4b36 (Flogb): frexp needs a pointer to int, not EMACS_INT.
Karl Heuer <kwzh@gnu.org>
parents: 11265
diff changeset
704 IN_FLOAT (frexp (f, &ivalue), "logb", arg);
8d57babf4b36 (Flogb): frexp needs a pointer to int, not EMACS_INT.
Karl Heuer <kwzh@gnu.org>
parents: 11265
diff changeset
705 value = ivalue - 1;
2205
c021f53fe7e5 * floatfns.c (logb): Add extern declaration for this.
Jim Blandy <jimb@redhat.com>
parents: 2129
diff changeset
706 #else
6359
800c035273e9 (Flogb): Check for 0.0. Emulate logb if needed.
Karl Heuer <kwzh@gnu.org>
parents: 6314
diff changeset
707 int i;
800c035273e9 (Flogb): Check for 0.0. Emulate logb if needed.
Karl Heuer <kwzh@gnu.org>
parents: 6314
diff changeset
708 double d;
800c035273e9 (Flogb): Check for 0.0. Emulate logb if needed.
Karl Heuer <kwzh@gnu.org>
parents: 6314
diff changeset
709 if (f < 0.0)
800c035273e9 (Flogb): Check for 0.0. Emulate logb if needed.
Karl Heuer <kwzh@gnu.org>
parents: 6314
diff changeset
710 f = -f;
800c035273e9 (Flogb): Check for 0.0. Emulate logb if needed.
Karl Heuer <kwzh@gnu.org>
parents: 6314
diff changeset
711 value = -1;
800c035273e9 (Flogb): Check for 0.0. Emulate logb if needed.
Karl Heuer <kwzh@gnu.org>
parents: 6314
diff changeset
712 while (f < 0.5)
800c035273e9 (Flogb): Check for 0.0. Emulate logb if needed.
Karl Heuer <kwzh@gnu.org>
parents: 6314
diff changeset
713 {
800c035273e9 (Flogb): Check for 0.0. Emulate logb if needed.
Karl Heuer <kwzh@gnu.org>
parents: 6314
diff changeset
714 for (i = 1, d = 0.5; d * d >= f; i += i)
800c035273e9 (Flogb): Check for 0.0. Emulate logb if needed.
Karl Heuer <kwzh@gnu.org>
parents: 6314
diff changeset
715 d *= d;
800c035273e9 (Flogb): Check for 0.0. Emulate logb if needed.
Karl Heuer <kwzh@gnu.org>
parents: 6314
diff changeset
716 f /= d;
800c035273e9 (Flogb): Check for 0.0. Emulate logb if needed.
Karl Heuer <kwzh@gnu.org>
parents: 6314
diff changeset
717 value -= i;
800c035273e9 (Flogb): Check for 0.0. Emulate logb if needed.
Karl Heuer <kwzh@gnu.org>
parents: 6314
diff changeset
718 }
800c035273e9 (Flogb): Check for 0.0. Emulate logb if needed.
Karl Heuer <kwzh@gnu.org>
parents: 6314
diff changeset
719 while (f >= 1.0)
800c035273e9 (Flogb): Check for 0.0. Emulate logb if needed.
Karl Heuer <kwzh@gnu.org>
parents: 6314
diff changeset
720 {
800c035273e9 (Flogb): Check for 0.0. Emulate logb if needed.
Karl Heuer <kwzh@gnu.org>
parents: 6314
diff changeset
721 for (i = 1, d = 2.0; d * d <= f; i += i)
800c035273e9 (Flogb): Check for 0.0. Emulate logb if needed.
Karl Heuer <kwzh@gnu.org>
parents: 6314
diff changeset
722 d *= d;
800c035273e9 (Flogb): Check for 0.0. Emulate logb if needed.
Karl Heuer <kwzh@gnu.org>
parents: 6314
diff changeset
723 f /= d;
800c035273e9 (Flogb): Check for 0.0. Emulate logb if needed.
Karl Heuer <kwzh@gnu.org>
parents: 6314
diff changeset
724 value += i;
800c035273e9 (Flogb): Check for 0.0. Emulate logb if needed.
Karl Heuer <kwzh@gnu.org>
parents: 6314
diff changeset
725 }
4501
9352d7d021c2 * config.h.in (HAVE_LOGB, HAVE_FREXP): Add #undefs for the
Jim Blandy <jimb@redhat.com>
parents: 3673
diff changeset
726 #endif
2205
c021f53fe7e5 * floatfns.c (logb): Add extern declaration for this.
Jim Blandy <jimb@redhat.com>
parents: 2129
diff changeset
727 #endif
6359
800c035273e9 (Flogb): Check for 0.0. Emulate logb if needed.
Karl Heuer <kwzh@gnu.org>
parents: 6314
diff changeset
728 }
9267
1b685d477c10 (FLOAT_TO_INT, FLOAT_TO_INT2, Fexpt, Flogb, Ffloor): Use new accessor macros
Karl Heuer <kwzh@gnu.org>
parents: 9129
diff changeset
729 XSETINT (val, value);
2119
4077ef8ad483 * floatfns.c (Flogb): Undo the change of Feb 22.
Jim Blandy <jimb@redhat.com>
parents: 2094
diff changeset
730 return val;
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
731 }
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
732
4506
6131dad14f6f (Ffloor): Optional second operand specifies divisor, as in Common Lisp.
Paul Eggert <eggert@twinsun.com>
parents: 4501
diff changeset
733
16857
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
734 /* the rounding functions */
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
735
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
736 static Lisp_Object
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
737 rounding_driver (arg, divisor, double_round, int_round2, name)
4506
6131dad14f6f (Ffloor): Optional second operand specifies divisor, as in Common Lisp.
Paul Eggert <eggert@twinsun.com>
parents: 4501
diff changeset
738 register Lisp_Object arg, divisor;
16857
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
739 double (*double_round) ();
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
740 EMACS_INT (*int_round2) ();
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
741 char *name;
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
742 {
40656
cdfd4d09b79a Update usage of CHECK_ macros (remove unused second argument).
Pavel Janík <Pavel@Janik.cz>
parents: 39814
diff changeset
743 CHECK_NUMBER_OR_FLOAT (arg);
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
744
4506
6131dad14f6f (Ffloor): Optional second operand specifies divisor, as in Common Lisp.
Paul Eggert <eggert@twinsun.com>
parents: 4501
diff changeset
745 if (! NILP (divisor))
6131dad14f6f (Ffloor): Optional second operand specifies divisor, as in Common Lisp.
Paul Eggert <eggert@twinsun.com>
parents: 4501
diff changeset
746 {
12537
476296adb950 (Fexpt): Use EMACS_INT for integer calculation.
Karl Heuer <kwzh@gnu.org>
parents: 11859
diff changeset
747 EMACS_INT i1, i2;
4506
6131dad14f6f (Ffloor): Optional second operand specifies divisor, as in Common Lisp.
Paul Eggert <eggert@twinsun.com>
parents: 4501
diff changeset
748
40656
cdfd4d09b79a Update usage of CHECK_ macros (remove unused second argument).
Pavel Janík <Pavel@Janik.cz>
parents: 39814
diff changeset
749 CHECK_NUMBER_OR_FLOAT (divisor);
4506
6131dad14f6f (Ffloor): Optional second operand specifies divisor, as in Common Lisp.
Paul Eggert <eggert@twinsun.com>
parents: 4501
diff changeset
750
9129
b2d1d925d5cc (extract_float, Fexpt, Fabs, Ffloat, Fceiling, Ffloor, Fround, Ftruncate): Use
Karl Heuer <kwzh@gnu.org>
parents: 8825
diff changeset
751 if (FLOATP (arg) || FLOATP (divisor))
4506
6131dad14f6f (Ffloor): Optional second operand specifies divisor, as in Common Lisp.
Paul Eggert <eggert@twinsun.com>
parents: 4501
diff changeset
752 {
6131dad14f6f (Ffloor): Optional second operand specifies divisor, as in Common Lisp.
Paul Eggert <eggert@twinsun.com>
parents: 4501
diff changeset
753 double f1, f2;
6131dad14f6f (Ffloor): Optional second operand specifies divisor, as in Common Lisp.
Paul Eggert <eggert@twinsun.com>
parents: 4501
diff changeset
754
25645
a14111a2a100 Use XCAR, XCDR, XFLOAT_DATA instead of explicit member access.
Ken Raeburn <raeburn@raeburn.org>
parents: 21514
diff changeset
755 f1 = FLOATP (arg) ? XFLOAT_DATA (arg) : XINT (arg);
a14111a2a100 Use XCAR, XCDR, XFLOAT_DATA instead of explicit member access.
Ken Raeburn <raeburn@raeburn.org>
parents: 21514
diff changeset
756 f2 = (FLOATP (divisor) ? XFLOAT_DATA (divisor) : XINT (divisor));
16786
8907c00c0cc6 <float.h>: Include if STDC_HEADERS.
Paul Eggert <eggert@twinsun.com>
parents: 16220
diff changeset
757 if (! IEEE_FLOATING_POINT && f2 == 0)
71978
6b2f430f5341 (arith_error, range_error, domain_error): Use xsignal2.
Kim F. Storm <storm@cua.dk>
parents: 68651
diff changeset
758 xsignal0 (Qarith_error);
4506
6131dad14f6f (Ffloor): Optional second operand specifies divisor, as in Common Lisp.
Paul Eggert <eggert@twinsun.com>
parents: 4501
diff changeset
759
16857
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
760 IN_FLOAT2 (f1 = (*double_round) (f1 / f2), name, arg, divisor);
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
761 FLOAT_TO_INT2 (f1, arg, name, arg, divisor);
4506
6131dad14f6f (Ffloor): Optional second operand specifies divisor, as in Common Lisp.
Paul Eggert <eggert@twinsun.com>
parents: 4501
diff changeset
762 return arg;
6131dad14f6f (Ffloor): Optional second operand specifies divisor, as in Common Lisp.
Paul Eggert <eggert@twinsun.com>
parents: 4501
diff changeset
763 }
6131dad14f6f (Ffloor): Optional second operand specifies divisor, as in Common Lisp.
Paul Eggert <eggert@twinsun.com>
parents: 4501
diff changeset
764
6131dad14f6f (Ffloor): Optional second operand specifies divisor, as in Common Lisp.
Paul Eggert <eggert@twinsun.com>
parents: 4501
diff changeset
765 i1 = XINT (arg);
6131dad14f6f (Ffloor): Optional second operand specifies divisor, as in Common Lisp.
Paul Eggert <eggert@twinsun.com>
parents: 4501
diff changeset
766 i2 = XINT (divisor);
6131dad14f6f (Ffloor): Optional second operand specifies divisor, as in Common Lisp.
Paul Eggert <eggert@twinsun.com>
parents: 4501
diff changeset
767
6131dad14f6f (Ffloor): Optional second operand specifies divisor, as in Common Lisp.
Paul Eggert <eggert@twinsun.com>
parents: 4501
diff changeset
768 if (i2 == 0)
71978
6b2f430f5341 (arith_error, range_error, domain_error): Use xsignal2.
Kim F. Storm <storm@cua.dk>
parents: 68651
diff changeset
769 xsignal0 (Qarith_error);
4506
6131dad14f6f (Ffloor): Optional second operand specifies divisor, as in Common Lisp.
Paul Eggert <eggert@twinsun.com>
parents: 4501
diff changeset
770
16857
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
771 XSETINT (arg, (*int_round2) (i1, i2));
4506
6131dad14f6f (Ffloor): Optional second operand specifies divisor, as in Common Lisp.
Paul Eggert <eggert@twinsun.com>
parents: 4501
diff changeset
772 return arg;
6131dad14f6f (Ffloor): Optional second operand specifies divisor, as in Common Lisp.
Paul Eggert <eggert@twinsun.com>
parents: 4501
diff changeset
773 }
6131dad14f6f (Ffloor): Optional second operand specifies divisor, as in Common Lisp.
Paul Eggert <eggert@twinsun.com>
parents: 4501
diff changeset
774
9129
b2d1d925d5cc (extract_float, Fexpt, Fabs, Ffloat, Fceiling, Ffloor, Fround, Ftruncate): Use
Karl Heuer <kwzh@gnu.org>
parents: 8825
diff changeset
775 if (FLOATP (arg))
6375
212dcd2c06e4 (FLOAT_TO_INT, FLOAT_TO_INT2, range_error2): New macros.
Karl Heuer <kwzh@gnu.org>
parents: 6359
diff changeset
776 {
212dcd2c06e4 (FLOAT_TO_INT, FLOAT_TO_INT2, range_error2): New macros.
Karl Heuer <kwzh@gnu.org>
parents: 6359
diff changeset
777 double d;
16857
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
778
25645
a14111a2a100 Use XCAR, XCDR, XFLOAT_DATA instead of explicit member access.
Ken Raeburn <raeburn@raeburn.org>
parents: 21514
diff changeset
779 IN_FLOAT (d = (*double_round) (XFLOAT_DATA (arg)), name, arg);
16857
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
780 FLOAT_TO_INT (d, arg, name, arg);
6375
212dcd2c06e4 (FLOAT_TO_INT, FLOAT_TO_INT2, range_error2): New macros.
Karl Heuer <kwzh@gnu.org>
parents: 6359
diff changeset
781 }
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
782
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
783 return arg;
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
784 }
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
785
16857
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
786 /* With C's /, the result is implementation-defined if either operand
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
787 is negative, so take care with negative operands in the following
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
788 integer functions. */
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
789
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
790 static EMACS_INT
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
791 ceiling2 (i1, i2)
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
792 EMACS_INT i1, i2;
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
793 {
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
794 return (i2 < 0
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
795 ? (i1 < 0 ? ((-1 - i1) / -i2) + 1 : - (i1 / -i2))
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
796 : (i1 <= 0 ? - (-i1 / i2) : ((i1 - 1) / i2) + 1));
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
797 }
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
798
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
799 static EMACS_INT
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
800 floor2 (i1, i2)
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
801 EMACS_INT i1, i2;
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
802 {
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
803 return (i2 < 0
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
804 ? (i1 <= 0 ? -i1 / -i2 : -1 - ((i1 - 1) / -i2))
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
805 : (i1 < 0 ? -1 - ((-1 - i1) / i2) : i1 / i2));
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
806 }
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
807
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
808 static EMACS_INT
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
809 truncate2 (i1, i2)
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
810 EMACS_INT i1, i2;
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
811 {
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
812 return (i2 < 0
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
813 ? (i1 < 0 ? -i1 / -i2 : - (i1 / -i2))
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
814 : (i1 < 0 ? - (-i1 / i2) : i1 / i2));
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
815 }
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
816
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
817 static EMACS_INT
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
818 round2 (i1, i2)
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
819 EMACS_INT i1, i2;
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
820 {
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
821 /* The C language's division operator gives us one remainder R, but
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
822 we want the remainder R1 on the other side of 0 if R1 is closer
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
823 to 0 than R is; because we want to round to even, we also want R1
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
824 if R and R1 are the same distance from 0 and if C's quotient is
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
825 odd. */
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
826 EMACS_INT q = i1 / i2;
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
827 EMACS_INT r = i1 % i2;
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
828 EMACS_INT abs_r = r < 0 ? -r : r;
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
829 EMACS_INT abs_r1 = (i2 < 0 ? -i2 : i2) - abs_r;
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
830 return q + (abs_r + (q & 1) <= abs_r1 ? 0 : (i2 ^ r) < 0 ? -1 : 1);
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
831 }
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
832
19737
cad2a27d4451 (emacs_rint): Define this,
Richard M. Stallman <rms@gnu.org>
parents: 18627
diff changeset
833 /* The code uses emacs_rint, so that it works to undefine HAVE_RINT
cad2a27d4451 (emacs_rint): Define this,
Richard M. Stallman <rms@gnu.org>
parents: 18627
diff changeset
834 if `rint' exists but does not work right. */
cad2a27d4451 (emacs_rint): Define this,
Richard M. Stallman <rms@gnu.org>
parents: 18627
diff changeset
835 #ifdef HAVE_RINT
cad2a27d4451 (emacs_rint): Define this,
Richard M. Stallman <rms@gnu.org>
parents: 18627
diff changeset
836 #define emacs_rint rint
cad2a27d4451 (emacs_rint): Define this,
Richard M. Stallman <rms@gnu.org>
parents: 18627
diff changeset
837 #else
16929
a6b5ec9a51b4 [!HAVE_RINT] (rint): Convert macro to an actual
Karl Heuer <kwzh@gnu.org>
parents: 16857
diff changeset
838 static double
19737
cad2a27d4451 (emacs_rint): Define this,
Richard M. Stallman <rms@gnu.org>
parents: 18627
diff changeset
839 emacs_rint (d)
16929
a6b5ec9a51b4 [!HAVE_RINT] (rint): Convert macro to an actual
Karl Heuer <kwzh@gnu.org>
parents: 16857
diff changeset
840 double d;
a6b5ec9a51b4 [!HAVE_RINT] (rint): Convert macro to an actual
Karl Heuer <kwzh@gnu.org>
parents: 16857
diff changeset
841 {
16956
bfa27111d35b Whitespace change.
Richard M. Stallman <rms@gnu.org>
parents: 16929
diff changeset
842 return floor (d + 0.5);
16929
a6b5ec9a51b4 [!HAVE_RINT] (rint): Convert macro to an actual
Karl Heuer <kwzh@gnu.org>
parents: 16857
diff changeset
843 }
a6b5ec9a51b4 [!HAVE_RINT] (rint): Convert macro to an actual
Karl Heuer <kwzh@gnu.org>
parents: 16857
diff changeset
844 #endif
a6b5ec9a51b4 [!HAVE_RINT] (rint): Convert macro to an actual
Karl Heuer <kwzh@gnu.org>
parents: 16857
diff changeset
845
16857
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
846 static double
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
847 double_identity (d)
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
848 double d;
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
849 {
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
850 return d;
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
851 }
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
852
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
853 DEFUN ("ceiling", Fceiling, Sceiling, 1, 2, 0,
41028
e4bce1db7f77 (Fceiling, Ffloor): Doc fixes.
Richard M. Stallman <rms@gnu.org>
parents: 41001
diff changeset
854 doc: /* Return the smallest integer no less than ARG.
e4bce1db7f77 (Fceiling, Ffloor): Doc fixes.
Richard M. Stallman <rms@gnu.org>
parents: 41001
diff changeset
855 This rounds the value towards +inf.
41001
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
856 With optional DIVISOR, return the smallest integer no less than ARG/DIVISOR. */)
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
857 (arg, divisor)
16857
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
858 Lisp_Object arg, divisor;
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
859 {
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
860 return rounding_driver (arg, divisor, ceil, ceiling2, "ceiling");
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
861 }
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
862
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
863 DEFUN ("floor", Ffloor, Sfloor, 1, 2, 0,
41028
e4bce1db7f77 (Fceiling, Ffloor): Doc fixes.
Richard M. Stallman <rms@gnu.org>
parents: 41001
diff changeset
864 doc: /* Return the largest integer no greater than ARG.
60651
de22fdff5bad (Ffloor): Doc fix.
Lute Kamstra <lute@gnu.org>
parents: 58986
diff changeset
865 This rounds the value towards -inf.
41001
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
866 With optional DIVISOR, return the largest integer no greater than ARG/DIVISOR. */)
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
867 (arg, divisor)
16857
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
868 Lisp_Object arg, divisor;
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
869 {
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
870 return rounding_driver (arg, divisor, floor, floor2, "floor");
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
871 }
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
872
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
873 DEFUN ("round", Fround, Sround, 1, 2, 0,
41001
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
874 doc: /* Return the nearest integer to ARG.
43179
b8eac678ad15 (Fround): Doc fix.
Eli Zaretskii <eliz@gnu.org>
parents: 41028
diff changeset
875 With optional DIVISOR, return the nearest integer to ARG/DIVISOR.
b8eac678ad15 (Fround): Doc fix.
Eli Zaretskii <eliz@gnu.org>
parents: 41028
diff changeset
876
43181
0f5eacbb65d9 (Fround): A better doc fix.
Eli Zaretskii <eliz@gnu.org>
parents: 43179
diff changeset
877 Rounding a value equidistant between two integers may choose the
0f5eacbb65d9 (Fround): A better doc fix.
Eli Zaretskii <eliz@gnu.org>
parents: 43179
diff changeset
878 integer closer to zero, or it may prefer an even integer, depending on
0f5eacbb65d9 (Fround): A better doc fix.
Eli Zaretskii <eliz@gnu.org>
parents: 43179
diff changeset
879 your machine. For example, \(round 2.5\) can return 3 on some
43189
004354581e54 (Fround): Fix a typo.
Eli Zaretskii <eliz@gnu.org>
parents: 43181
diff changeset
880 systems, but 2 on others. */)
41001
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
881 (arg, divisor)
16857
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
882 Lisp_Object arg, divisor;
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
883 {
19737
cad2a27d4451 (emacs_rint): Define this,
Richard M. Stallman <rms@gnu.org>
parents: 18627
diff changeset
884 return rounding_driver (arg, divisor, emacs_rint, round2, "round");
16857
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
885 }
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
886
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
887 DEFUN ("truncate", Ftruncate, Struncate, 1, 2, 0,
41001
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
888 doc: /* Truncate a floating point number to an int.
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
889 Rounds ARG toward zero.
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
890 With optional DIVISOR, truncate ARG/DIVISOR. */)
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
891 (arg, divisor)
16857
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
892 Lisp_Object arg, divisor;
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
893 {
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
894 return rounding_driver (arg, divisor, double_identity, truncate2,
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
895 "truncate");
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
896 }
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
897
4506
6131dad14f6f (Ffloor): Optional second operand specifies divisor, as in Common Lisp.
Paul Eggert <eggert@twinsun.com>
parents: 4501
diff changeset
898
16786
8907c00c0cc6 <float.h>: Include if STDC_HEADERS.
Paul Eggert <eggert@twinsun.com>
parents: 16220
diff changeset
899 Lisp_Object
8907c00c0cc6 <float.h>: Include if STDC_HEADERS.
Paul Eggert <eggert@twinsun.com>
parents: 16220
diff changeset
900 fmod_float (x, y)
8907c00c0cc6 <float.h>: Include if STDC_HEADERS.
Paul Eggert <eggert@twinsun.com>
parents: 16220
diff changeset
901 register Lisp_Object x, y;
8907c00c0cc6 <float.h>: Include if STDC_HEADERS.
Paul Eggert <eggert@twinsun.com>
parents: 16220
diff changeset
902 {
8907c00c0cc6 <float.h>: Include if STDC_HEADERS.
Paul Eggert <eggert@twinsun.com>
parents: 16220
diff changeset
903 double f1, f2;
8907c00c0cc6 <float.h>: Include if STDC_HEADERS.
Paul Eggert <eggert@twinsun.com>
parents: 16220
diff changeset
904
25645
a14111a2a100 Use XCAR, XCDR, XFLOAT_DATA instead of explicit member access.
Ken Raeburn <raeburn@raeburn.org>
parents: 21514
diff changeset
905 f1 = FLOATP (x) ? XFLOAT_DATA (x) : XINT (x);
a14111a2a100 Use XCAR, XCDR, XFLOAT_DATA instead of explicit member access.
Ken Raeburn <raeburn@raeburn.org>
parents: 21514
diff changeset
906 f2 = FLOATP (y) ? XFLOAT_DATA (y) : XINT (y);
16786
8907c00c0cc6 <float.h>: Include if STDC_HEADERS.
Paul Eggert <eggert@twinsun.com>
parents: 16220
diff changeset
907
8907c00c0cc6 <float.h>: Include if STDC_HEADERS.
Paul Eggert <eggert@twinsun.com>
parents: 16220
diff changeset
908 if (! IEEE_FLOATING_POINT && f2 == 0)
71978
6b2f430f5341 (arith_error, range_error, domain_error): Use xsignal2.
Kim F. Storm <storm@cua.dk>
parents: 68651
diff changeset
909 xsignal0 (Qarith_error);
16786
8907c00c0cc6 <float.h>: Include if STDC_HEADERS.
Paul Eggert <eggert@twinsun.com>
parents: 16220
diff changeset
910
8907c00c0cc6 <float.h>: Include if STDC_HEADERS.
Paul Eggert <eggert@twinsun.com>
parents: 16220
diff changeset
911 /* If the "remainder" comes out with the wrong sign, fix it. */
8907c00c0cc6 <float.h>: Include if STDC_HEADERS.
Paul Eggert <eggert@twinsun.com>
parents: 16220
diff changeset
912 IN_FLOAT2 ((f1 = fmod (f1, f2),
8907c00c0cc6 <float.h>: Include if STDC_HEADERS.
Paul Eggert <eggert@twinsun.com>
parents: 16220
diff changeset
913 f1 = (f2 < 0 ? f1 > 0 : f1 < 0) ? f1 + f2 : f1),
8907c00c0cc6 <float.h>: Include if STDC_HEADERS.
Paul Eggert <eggert@twinsun.com>
parents: 16220
diff changeset
914 "mod", x, y);
8907c00c0cc6 <float.h>: Include if STDC_HEADERS.
Paul Eggert <eggert@twinsun.com>
parents: 16220
diff changeset
915 return make_float (f1);
8907c00c0cc6 <float.h>: Include if STDC_HEADERS.
Paul Eggert <eggert@twinsun.com>
parents: 16220
diff changeset
916 }
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
917
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
918 /* It's not clear these are worth adding. */
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
919
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
920 DEFUN ("fceiling", Ffceiling, Sfceiling, 1, 1, 0,
41001
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
921 doc: /* Return the smallest integer no less than ARG, as a float.
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
922 \(Round toward +inf.\) */)
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
923 (arg)
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
924 register Lisp_Object arg;
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
925 {
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
926 double d = extract_float (arg);
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
927 IN_FLOAT (d = ceil (d), "fceiling", arg);
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
928 return make_float (d);
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
929 }
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
930
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
931 DEFUN ("ffloor", Fffloor, Sffloor, 1, 1, 0,
41001
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
932 doc: /* Return the largest integer no greater than ARG, as a float.
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
933 \(Round towards -inf.\) */)
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
934 (arg)
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
935 register Lisp_Object arg;
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
936 {
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
937 double d = extract_float (arg);
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
938 IN_FLOAT (d = floor (d), "ffloor", arg);
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
939 return make_float (d);
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
940 }
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
941
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
942 DEFUN ("fround", Ffround, Sfround, 1, 1, 0,
41001
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
943 doc: /* Return the nearest integer to ARG, as a float. */)
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
944 (arg)
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
945 register Lisp_Object arg;
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
946 {
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
947 double d = extract_float (arg);
19737
cad2a27d4451 (emacs_rint): Define this,
Richard M. Stallman <rms@gnu.org>
parents: 18627
diff changeset
948 IN_FLOAT (d = emacs_rint (d), "fround", arg);
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
949 return make_float (d);
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
950 }
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
951
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
952 DEFUN ("ftruncate", Fftruncate, Sftruncate, 1, 1, 0,
41001
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
953 doc: /* Truncate a floating point number to an integral float value.
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
954 Rounds the value toward zero. */)
a17c8b15ef1b Change doc-string comments to `new style' [w/`doc:' keyword].
Pavel Janík <Pavel@Janik.cz>
parents: 40656
diff changeset
955 (arg)
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
956 register Lisp_Object arg;
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
957 {
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
958 double d = extract_float (arg);
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
959 if (d >= 0.0)
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
960 IN_FLOAT (d = floor (d), "ftruncate", arg);
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
961 else
5495
87f9165f5b14 [MSDOS]: Don't define HAVE_MATHERR.
Richard M. Stallman <rms@gnu.org>
parents: 4881
diff changeset
962 IN_FLOAT (d = ceil (d), "ftruncate", arg);
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
963 return make_float (d);
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
964 }
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
965
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
966 #ifdef FLOAT_CATCH_SIGILL
621
eca8812e61cd *** empty log message ***
Jim Blandy <jimb@redhat.com>
parents: 485
diff changeset
967 static SIGTYPE
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
968 float_error (signo)
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
969 int signo;
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
970 {
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
971 if (! in_float)
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
972 fatal_error_signal (signo);
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
973
16220
02044b05d8e0 Replaced symbol BSD with BSD_SYSTEM.
Karl Heuer <kwzh@gnu.org>
parents: 16207
diff changeset
974 #ifdef BSD_SYSTEM
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
975 #ifdef BSD4_1
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
976 sigrelse (SIGILL);
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
977 #else /* not BSD4_1 */
638
40b255f55df3 *** empty log message ***
Jim Blandy <jimb@redhat.com>
parents: 624
diff changeset
978 sigsetmask (SIGEMPTYMASK);
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
979 #endif /* not BSD4_1 */
485
8c615e453683 *** empty log message ***
Jim Blandy <jimb@redhat.com>
parents: 102
diff changeset
980 #else
8c615e453683 *** empty log message ***
Jim Blandy <jimb@redhat.com>
parents: 102
diff changeset
981 /* Must reestablish handler each time it is called. */
8c615e453683 *** empty log message ***
Jim Blandy <jimb@redhat.com>
parents: 102
diff changeset
982 signal (SIGILL, float_error);
16220
02044b05d8e0 Replaced symbol BSD with BSD_SYSTEM.
Karl Heuer <kwzh@gnu.org>
parents: 16207
diff changeset
983 #endif /* BSD_SYSTEM */
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
984
58986
59945307b86b * syssignal.h: Declare main_thread.
Jan Djärv <jan.h.d@swipnet.se>
parents: 53002
diff changeset
985 SIGNAL_THREAD_CHECK (signo);
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
986 in_float = 0;
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
987
71978
6b2f430f5341 (arith_error, range_error, domain_error): Use xsignal2.
Kim F. Storm <storm@cua.dk>
parents: 68651
diff changeset
988 xsignal1 (Qarith_error, float_error_arg);
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
989 }
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
990
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
991 /* Another idea was to replace the library function `infnan'
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
992 where SIGILL is signaled. */
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
993
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
994 #endif /* FLOAT_CATCH_SIGILL */
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
995
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
996 #ifdef HAVE_MATHERR
49600
23a1cea22d13 Trailing whitespace deleted.
Juanma Barranquero <lekktu@gmail.com>
parents: 43413
diff changeset
997 int
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
998 matherr (x)
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
999 struct exception *x;
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
1000 {
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
1001 Lisp_Object args;
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
1002 if (! in_float)
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
1003 /* Not called from emacs-lisp float routines; do the default thing. */
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
1004 return 0;
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
1005 if (!strcmp (x->name, "pow"))
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
1006 x->name = "expt";
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
1007
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
1008 args
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
1009 = Fcons (build_string (x->name),
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
1010 Fcons (make_float (x->arg1),
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
1011 ((!strcmp (x->name, "log") || !strcmp (x->name, "pow"))
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
1012 ? Fcons (make_float (x->arg2), Qnil)
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
1013 : Qnil)));
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
1014 switch (x->type)
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
1015 {
71978
6b2f430f5341 (arith_error, range_error, domain_error): Use xsignal2.
Kim F. Storm <storm@cua.dk>
parents: 68651
diff changeset
1016 case DOMAIN: xsignal (Qdomain_error, args); break;
6b2f430f5341 (arith_error, range_error, domain_error): Use xsignal2.
Kim F. Storm <storm@cua.dk>
parents: 68651
diff changeset
1017 case SING: xsignal (Qsingularity_error, args); break;
6b2f430f5341 (arith_error, range_error, domain_error): Use xsignal2.
Kim F. Storm <storm@cua.dk>
parents: 68651
diff changeset
1018 case OVERFLOW: xsignal (Qoverflow_error, args); break;
6b2f430f5341 (arith_error, range_error, domain_error): Use xsignal2.
Kim F. Storm <storm@cua.dk>
parents: 68651
diff changeset
1019 case UNDERFLOW: xsignal (Qunderflow_error, args); break;
6b2f430f5341 (arith_error, range_error, domain_error): Use xsignal2.
Kim F. Storm <storm@cua.dk>
parents: 68651
diff changeset
1020 default: xsignal (Qarith_error, args); break;
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
1021 }
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
1022 return (1); /* don't set errno or print a message */
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
1023 }
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
1024 #endif /* HAVE_MATHERR */
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
1025
21514
fa9ff387d260 Fix -Wimplicit warnings.
Andreas Schwab <schwab@suse.de>
parents: 20122
diff changeset
1026 void
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
1027 init_floatfns ()
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
1028 {
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
1029 #ifdef FLOAT_CATCH_SIGILL
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
1030 signal (SIGILL, float_error);
49600
23a1cea22d13 Trailing whitespace deleted.
Juanma Barranquero <lekktu@gmail.com>
parents: 43413
diff changeset
1031 #endif
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
1032 in_float = 0;
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
1033 }
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
1034
21514
fa9ff387d260 Fix -Wimplicit warnings.
Andreas Schwab <schwab@suse.de>
parents: 20122
diff changeset
1035 void
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
1036 syms_of_floatfns ()
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
1037 {
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
1038 defsubr (&Sacos);
683
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
1039 defsubr (&Sasin);
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
1040 defsubr (&Satan);
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
1041 defsubr (&Scos);
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
1042 defsubr (&Ssin);
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
1043 defsubr (&Stan);
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
1044 #if 0
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
1045 defsubr (&Sacosh);
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
1046 defsubr (&Sasinh);
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
1047 defsubr (&Satanh);
683
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
1048 defsubr (&Scosh);
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
1049 defsubr (&Ssinh);
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
1050 defsubr (&Stanh);
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
1051 defsubr (&Sbessel_y0);
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
1052 defsubr (&Sbessel_y1);
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
1053 defsubr (&Sbessel_yn);
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
1054 defsubr (&Sbessel_j0);
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
1055 defsubr (&Sbessel_j1);
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
1056 defsubr (&Sbessel_jn);
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
1057 defsubr (&Serf);
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
1058 defsubr (&Serfc);
683
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
1059 defsubr (&Slog_gamma);
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
1060 defsubr (&Scube_root);
5595
63bc8a14a073 (Fffloor, Ffceil, Ffround, Fftruncate): New functions.
Richard M. Stallman <rms@gnu.org>
parents: 5495
diff changeset
1061 #endif
2094
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
1062 defsubr (&Sfceiling);
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
1063 defsubr (&Sffloor);
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
1064 defsubr (&Sfround);
c77607f8e32d (FLOAT_CHECK_ERRNO): Define unless NO_FLOAT_CHECK_ERRNO.
Richard M. Stallman <rms@gnu.org>
parents: 1918
diff changeset
1065 defsubr (&Sftruncate);
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
1066 defsubr (&Sexp);
683
7f4d77d29804 *** empty log message ***
Richard M. Stallman <rms@gnu.org>
parents: 638
diff changeset
1067 defsubr (&Sexpt);
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
1068 defsubr (&Slog);
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
1069 defsubr (&Slog10);
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
1070 defsubr (&Ssqrt);
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
1071
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
1072 defsubr (&Sabs);
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
1073 defsubr (&Sfloat);
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
1074 defsubr (&Slogb);
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
1075 defsubr (&Sceiling);
16857
bdafa1f28a64 (rounding_driver): New function for systematic support of
Paul Eggert <eggert@twinsun.com>
parents: 16786
diff changeset
1076 defsubr (&Sfloor);
102
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
1077 defsubr (&Sround);
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
1078 defsubr (&Struncate);
0d3a6b3b64a4 Initial revision
Mike Rowan <mtr@gnu.org>
parents:
diff changeset
1079 }
52401
695cf19ef79e Add arch taglines
Miles Bader <miles@gnu.org>
parents: 49600
diff changeset
1080
695cf19ef79e Add arch taglines
Miles Bader <miles@gnu.org>
parents: 49600
diff changeset
1081 /* arch-tag: be05bf9d-049e-4e31-91b9-e6153d483ae7
695cf19ef79e Add arch taglines
Miles Bader <miles@gnu.org>
parents: 49600
diff changeset
1082 (do not change this comment) */