Mercurial > cpython-lru_cache-weakref

annotate Python/pymath.c @ 58664:993f703a61ba

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Merged revisions 77484,77487,77561,77570,77593,77603,77608,77667,77702-77703,77739,77858,77887,77889 via svnmerge from svn+ssh://pythondev@svn.python.org/python/trunk ........ r77484 | skip.montanaro | 2010-01-13 19:12:34 -0600 (Wed, 13 Jan 2010) | 4 lines Update PyEval_EvalFrame to PyEval_EvalFrameEx. This looks to have been done partially before. Also add a comment describing how this might have to work with different versions of the interpreter. ........ r77487 | ezio.melotti | 2010-01-14 05:34:10 -0600 (Thu, 14 Jan 2010) | 1 line Fixed typo ........ r77561 | georg.brandl | 2010-01-17 02:42:30 -0600 (Sun, 17 Jan 2010) | 1 line #7699: improve datetime docs: straightforward linking to strftime/strptime section, mark classmethods as such. ........ r77570 | georg.brandl | 2010-01-17 06:14:42 -0600 (Sun, 17 Jan 2010) | 1 line Add note about usage of STRINGLIB_EMPTY. ........ r77593 | georg.brandl | 2010-01-17 17:33:53 -0600 (Sun, 17 Jan 2010) | 1 line Fix internal reference. ........ r77603 | benjamin.peterson | 2010-01-18 17:07:56 -0600 (Mon, 18 Jan 2010) | 8 lines data descriptors do not override the class dictionary if __get__ is not defined Adjust documentation and add a test to verify this behavior. See http://mail.python.org/pipermail/python-dev/2010-January/095637.html for discussion. ........ r77608 | gregory.p.smith | 2010-01-19 02:19:03 -0600 (Tue, 19 Jan 2010) | 6 lines Do not compile stubs for the sha2 series hashes in the openssl hashlib module when the openssl version is too old to support them. That leads both compiled code bloat and to unittests attempting to test implementations that don't exist for comparison purposes on such platforms. ........ r77667 | mark.dickinson | 2010-01-21 12:32:27 -0600 (Thu, 21 Jan 2010) | 1 line Add two more test_strtod test values. ........ r77702 | georg.brandl | 2010-01-23 02:43:31 -0600 (Sat, 23 Jan 2010) | 1 line #7762: fix refcount annotation of PyUnicode_Tailmatch(). ........ r77703 | georg.brandl | 2010-01-23 02:47:54 -0600 (Sat, 23 Jan 2010) | 1 line #7725: fix referencing issue. ........ r77739 | benjamin.peterson | 2010-01-24 21:52:52 -0600 (Sun, 24 Jan 2010) | 1 line mention from_float() in error message ........ r77858 | georg.brandl | 2010-01-30 11:57:48 -0600 (Sat, 30 Jan 2010) | 1 line #7802: fix invalid example (heh). ........ r77887 | georg.brandl | 2010-01-31 12:51:49 -0600 (Sun, 31 Jan 2010) | 5 lines Fix-up ftplib documentation: move exception descriptions to toplevel, not inside a class remove attribution in "versionadded" spell and grammar check docstring of FTP_TLS ........ r77889 | michael.foord | 2010-01-31 13:59:26 -0600 (Sun, 31 Jan 2010) | 1 line Minor modification to unittest documentation. ........
author Benjamin Peterson <benjamin@python.org>
date Wed, 03 Feb 2010 02:35:45 +0000
parents 9751ddcc1882
children cab14be0ada1
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
46477
b69f443ef480 Merged revisions 62380,62382-62383 via svnmerge from
Christian Heimes <christian@cheimes.de>
parents:
diff changeset
1 #include "Python.h"
b69f443ef480 Merged revisions 62380,62382-62383 via svnmerge from
Christian Heimes <christian@cheimes.de>
parents:
diff changeset
2
51485
f9f2ed79b00f Merged revisions 69459 via svnmerge from
Mark Dickinson <dickinsm@gmail.com>
parents: 49270
diff changeset
3 #ifdef X87_DOUBLE_ROUNDING
f9f2ed79b00f Merged revisions 69459 via svnmerge from
Mark Dickinson <dickinsm@gmail.com>
parents: 49270
diff changeset
4 /* On x86 platforms using an x87 FPU, this function is called from the
f9f2ed79b00f Merged revisions 69459 via svnmerge from
Mark Dickinson <dickinsm@gmail.com>
parents: 49270
diff changeset
5 Py_FORCE_DOUBLE macro (defined in pymath.h) to force a floating-point
f9f2ed79b00f Merged revisions 69459 via svnmerge from
Mark Dickinson <dickinsm@gmail.com>
parents: 49270
diff changeset
6 number out of an 80-bit x87 FPU register and into a 64-bit memory location,
f9f2ed79b00f Merged revisions 69459 via svnmerge from
Mark Dickinson <dickinsm@gmail.com>
parents: 49270
diff changeset
7 thus rounding from extended precision to double precision. */
f9f2ed79b00f Merged revisions 69459 via svnmerge from
Mark Dickinson <dickinsm@gmail.com>
parents: 49270
diff changeset
8 double _Py_force_double(double x)
f9f2ed79b00f Merged revisions 69459 via svnmerge from
Mark Dickinson <dickinsm@gmail.com>
parents: 49270
diff changeset
9 {
f9f2ed79b00f Merged revisions 69459 via svnmerge from
Mark Dickinson <dickinsm@gmail.com>
parents: 49270
diff changeset
10 volatile double y;
f9f2ed79b00f Merged revisions 69459 via svnmerge from
Mark Dickinson <dickinsm@gmail.com>
parents: 49270
diff changeset
11 y = x;
f9f2ed79b00f Merged revisions 69459 via svnmerge from
Mark Dickinson <dickinsm@gmail.com>
parents: 49270
diff changeset
12 return y;
f9f2ed79b00f Merged revisions 69459 via svnmerge from
Mark Dickinson <dickinsm@gmail.com>
parents: 49270
diff changeset
13 }
f9f2ed79b00f Merged revisions 69459 via svnmerge from
Mark Dickinson <dickinsm@gmail.com>
parents: 49270
diff changeset
14 #endif
f9f2ed79b00f Merged revisions 69459 via svnmerge from
Mark Dickinson <dickinsm@gmail.com>
parents: 49270
diff changeset
15
53240
28645a961004 The SSE2 detection and enabling could potentially cause
Mark Dickinson <dickinsm@gmail.com>
parents: 53223
diff changeset
16 #ifdef HAVE_GCC_ASM_FOR_X87
53185
97ae7a709c1d Issue #1580: use short float repr where possible.
Mark Dickinson <dickinsm@gmail.com>
parents: 51485
diff changeset
17
97ae7a709c1d Issue #1580: use short float repr where possible.
Mark Dickinson <dickinsm@gmail.com>
parents: 51485
diff changeset
18 /* inline assembly for getting and setting the 387 FPU control word on
97ae7a709c1d Issue #1580: use short float repr where possible.
Mark Dickinson <dickinsm@gmail.com>
parents: 51485
diff changeset
19 gcc/x86 */
97ae7a709c1d Issue #1580: use short float repr where possible.
Mark Dickinson <dickinsm@gmail.com>
parents: 51485
diff changeset
20
97ae7a709c1d Issue #1580: use short float repr where possible.
Mark Dickinson <dickinsm@gmail.com>
parents: 51485
diff changeset
21 unsigned short _Py_get_387controlword(void) {
97ae7a709c1d Issue #1580: use short float repr where possible.
Mark Dickinson <dickinsm@gmail.com>
parents: 51485
diff changeset
22 unsigned short cw;
97ae7a709c1d Issue #1580: use short float repr where possible.
Mark Dickinson <dickinsm@gmail.com>
parents: 51485
diff changeset
23 __asm__ __volatile__ ("fnstcw %0" : "=m" (cw));
97ae7a709c1d Issue #1580: use short float repr where possible.
Mark Dickinson <dickinsm@gmail.com>
parents: 51485
diff changeset
24 return cw;
97ae7a709c1d Issue #1580: use short float repr where possible.
Mark Dickinson <dickinsm@gmail.com>
parents: 51485
diff changeset
25 }
97ae7a709c1d Issue #1580: use short float repr where possible.
Mark Dickinson <dickinsm@gmail.com>
parents: 51485
diff changeset
26
97ae7a709c1d Issue #1580: use short float repr where possible.
Mark Dickinson <dickinsm@gmail.com>
parents: 51485
diff changeset
27 void _Py_set_387controlword(unsigned short cw) {
97ae7a709c1d Issue #1580: use short float repr where possible.
Mark Dickinson <dickinsm@gmail.com>
parents: 51485
diff changeset
28 __asm__ __volatile__ ("fldcw %0" : : "m" (cw));
97ae7a709c1d Issue #1580: use short float repr where possible.
Mark Dickinson <dickinsm@gmail.com>
parents: 51485
diff changeset
29 }
97ae7a709c1d Issue #1580: use short float repr where possible.
Mark Dickinson <dickinsm@gmail.com>
parents: 51485
diff changeset
30
97ae7a709c1d Issue #1580: use short float repr where possible.
Mark Dickinson <dickinsm@gmail.com>
parents: 51485
diff changeset
31 #endif
97ae7a709c1d Issue #1580: use short float repr where possible.
Mark Dickinson <dickinsm@gmail.com>
parents: 51485
diff changeset
32
97ae7a709c1d Issue #1580: use short float repr where possible.
Mark Dickinson <dickinsm@gmail.com>
parents: 51485
diff changeset
33
46477
b69f443ef480 Merged revisions 62380,62382-62383 via svnmerge from
Christian Heimes <christian@cheimes.de>
parents:
diff changeset
34 #ifndef HAVE_HYPOT
b69f443ef480 Merged revisions 62380,62382-62383 via svnmerge from
Christian Heimes <christian@cheimes.de>
parents:
diff changeset
35 double hypot(double x, double y)
b69f443ef480 Merged revisions 62380,62382-62383 via svnmerge from
Christian Heimes <christian@cheimes.de>
parents:
diff changeset
36 {
b69f443ef480 Merged revisions 62380,62382-62383 via svnmerge from
Christian Heimes <christian@cheimes.de>
parents:
diff changeset
37 double yx;
b69f443ef480 Merged revisions 62380,62382-62383 via svnmerge from
Christian Heimes <christian@cheimes.de>
parents:
diff changeset
38
b69f443ef480 Merged revisions 62380,62382-62383 via svnmerge from
Christian Heimes <christian@cheimes.de>
parents:
diff changeset
39 x = fabs(x);
b69f443ef480 Merged revisions 62380,62382-62383 via svnmerge from
Christian Heimes <christian@cheimes.de>
parents:
diff changeset
40 y = fabs(y);
b69f443ef480 Merged revisions 62380,62382-62383 via svnmerge from
Christian Heimes <christian@cheimes.de>
parents:
diff changeset
41 if (x < y) {
b69f443ef480 Merged revisions 62380,62382-62383 via svnmerge from
Christian Heimes <christian@cheimes.de>
parents:
diff changeset
42 double temp = x;
b69f443ef480 Merged revisions 62380,62382-62383 via svnmerge from
Christian Heimes <christian@cheimes.de>
parents:
diff changeset
43 x = y;
b69f443ef480 Merged revisions 62380,62382-62383 via svnmerge from
Christian Heimes <christian@cheimes.de>
parents:
diff changeset
44 y = temp;
b69f443ef480 Merged revisions 62380,62382-62383 via svnmerge from
Christian Heimes <christian@cheimes.de>
parents:
diff changeset
45 }
b69f443ef480 Merged revisions 62380,62382-62383 via svnmerge from
Christian Heimes <christian@cheimes.de>
parents:
diff changeset
46 if (x == 0.)
b69f443ef480 Merged revisions 62380,62382-62383 via svnmerge from
Christian Heimes <christian@cheimes.de>
parents:
diff changeset
47 return 0.;
b69f443ef480 Merged revisions 62380,62382-62383 via svnmerge from
Christian Heimes <christian@cheimes.de>
parents:
diff changeset
48 else {
b69f443ef480 Merged revisions 62380,62382-62383 via svnmerge from
Christian Heimes <christian@cheimes.de>
parents:
diff changeset
49 yx = y/x;
b69f443ef480 Merged revisions 62380,62382-62383 via svnmerge from
Christian Heimes <christian@cheimes.de>
parents:
diff changeset
50 return x*sqrt(1.+yx*yx);
b69f443ef480 Merged revisions 62380,62382-62383 via svnmerge from
Christian Heimes <christian@cheimes.de>
parents:
diff changeset
51 }
b69f443ef480 Merged revisions 62380,62382-62383 via svnmerge from
Christian Heimes <christian@cheimes.de>
parents:
diff changeset
52 }
b69f443ef480 Merged revisions 62380,62382-62383 via svnmerge from
Christian Heimes <christian@cheimes.de>
parents:
diff changeset
53 #endif /* HAVE_HYPOT */
b69f443ef480 Merged revisions 62380,62382-62383 via svnmerge from
Christian Heimes <christian@cheimes.de>
parents:
diff changeset
54
b69f443ef480 Merged revisions 62380,62382-62383 via svnmerge from
Christian Heimes <christian@cheimes.de>
parents:
diff changeset
55 #ifndef HAVE_COPYSIGN
53223
8209b6c32e61 Merged revisions 71705 via svnmerge from
Mark Dickinson <dickinsm@gmail.com>
parents: 53221
diff changeset
56 double
46477
b69f443ef480 Merged revisions 62380,62382-62383 via svnmerge from
Christian Heimes <christian@cheimes.de>
parents:
diff changeset
57 copysign(double x, double y)
b69f443ef480 Merged revisions 62380,62382-62383 via svnmerge from
Christian Heimes <christian@cheimes.de>
parents:
diff changeset
58 {
b69f443ef480 Merged revisions 62380,62382-62383 via svnmerge from
Christian Heimes <christian@cheimes.de>
parents:
diff changeset
59 /* use atan2 to distinguish -0. from 0. */
b69f443ef480 Merged revisions 62380,62382-62383 via svnmerge from
Christian Heimes <christian@cheimes.de>
parents:
diff changeset
60 if (y > 0. || (y == 0. && atan2(y, -1.) > 0.)) {
b69f443ef480 Merged revisions 62380,62382-62383 via svnmerge from
Christian Heimes <christian@cheimes.de>
parents:
diff changeset
61 return fabs(x);
b69f443ef480 Merged revisions 62380,62382-62383 via svnmerge from
Christian Heimes <christian@cheimes.de>
parents:
diff changeset
62 } else {
b69f443ef480 Merged revisions 62380,62382-62383 via svnmerge from
Christian Heimes <christian@cheimes.de>
parents:
diff changeset
63 return -fabs(x);
b69f443ef480 Merged revisions 62380,62382-62383 via svnmerge from
Christian Heimes <christian@cheimes.de>
parents:
diff changeset
64 }
b69f443ef480 Merged revisions 62380,62382-62383 via svnmerge from
Christian Heimes <christian@cheimes.de>
parents:
diff changeset
65 }
b69f443ef480 Merged revisions 62380,62382-62383 via svnmerge from
Christian Heimes <christian@cheimes.de>
parents:
diff changeset
66 #endif /* HAVE_COPYSIGN */
b69f443ef480 Merged revisions 62380,62382-62383 via svnmerge from
Christian Heimes <christian@cheimes.de>
parents:
diff changeset
67
53221
c5c8680ee0b1 Add check for C99 round function to configure, and define
Mark Dickinson <dickinsm@gmail.com>
parents: 53185
diff changeset
68 #ifndef HAVE_ROUND
c5c8680ee0b1 Add check for C99 round function to configure, and define
Mark Dickinson <dickinsm@gmail.com>
parents: 53185
diff changeset
69 double
c5c8680ee0b1 Add check for C99 round function to configure, and define
Mark Dickinson <dickinsm@gmail.com>
parents: 53185
diff changeset
70 round(double x)
c5c8680ee0b1 Add check for C99 round function to configure, and define
Mark Dickinson <dickinsm@gmail.com>
parents: 53185
diff changeset
71 {
c5c8680ee0b1 Add check for C99 round function to configure, and define
Mark Dickinson <dickinsm@gmail.com>
parents: 53185
diff changeset
72 double absx, y;
c5c8680ee0b1 Add check for C99 round function to configure, and define
Mark Dickinson <dickinsm@gmail.com>
parents: 53185
diff changeset
73 absx = fabs(x);
c5c8680ee0b1 Add check for C99 round function to configure, and define
Mark Dickinson <dickinsm@gmail.com>
parents: 53185
diff changeset
74 y = floor(absx);
c5c8680ee0b1 Add check for C99 round function to configure, and define
Mark Dickinson <dickinsm@gmail.com>
parents: 53185
diff changeset
75 if (absx - y >= 0.5)
c5c8680ee0b1 Add check for C99 round function to configure, and define
Mark Dickinson <dickinsm@gmail.com>
parents: 53185
diff changeset
76 y += 1.0;
c5c8680ee0b1 Add check for C99 round function to configure, and define
Mark Dickinson <dickinsm@gmail.com>
parents: 53185
diff changeset
77 return copysign(y, x);
c5c8680ee0b1 Add check for C99 round function to configure, and define
Mark Dickinson <dickinsm@gmail.com>
parents: 53185
diff changeset
78 }
c5c8680ee0b1 Add check for C99 round function to configure, and define
Mark Dickinson <dickinsm@gmail.com>
parents: 53185
diff changeset
79 #endif /* HAVE_ROUND */