1993-07-28 09:05:47 +00:00
|
|
|
#ifndef Py_LONGINTREPR_H
|
|
|
|
#define Py_LONGINTREPR_H
|
|
|
|
#ifdef __cplusplus
|
|
|
|
extern "C" {
|
|
|
|
#endif
|
|
|
|
|
1991-05-05 20:09:44 +00:00
|
|
|
|
1992-01-19 16:30:12 +00:00
|
|
|
/* This is published for the benefit of "friend" marshal.c only. */
|
|
|
|
|
1991-05-05 20:09:44 +00:00
|
|
|
/* Parameters of the long integer representation.
|
|
|
|
These shouldn't have to be changed as C should guarantee that a short
|
2000-07-08 04:17:21 +00:00
|
|
|
contains at least 16 bits, but it's made changeable anyway.
|
1991-05-05 20:09:44 +00:00
|
|
|
Note: 'digit' should be able to hold 2*MASK+1, and 'twodigits'
|
|
|
|
should be able to hold the intermediate results in 'mul'
|
2004-08-29 22:16:50 +00:00
|
|
|
(at most (BASE-1)*(2*BASE+1) == MASK*(2*MASK+3)).
|
1991-05-05 20:09:44 +00:00
|
|
|
Also, x_sub assumes that 'digit' is an unsigned type, and overflow
|
|
|
|
is handled by taking the result mod 2**N for some N > SHIFT.
|
SF patch 936813: fast modular exponentiation
This checkin is adapted from part 2 (of 3) of Trevor Perrin's patch set.
BACKWARD INCOMPATIBILITY: SHIFT must now be divisible by 5. AFAIK,
nobody will care. long_pow() could be complicated to worm around that,
if necessary.
long_pow():
- BUGFIX: This leaked the base and power when the power was negative
(and so the computation delegated to float pow).
- Instead of doing right-to-left exponentiation, do left-to-right. This
is more efficient for small bases, which is the common case.
- In addition, if the exponent is large (more than FIVEARY_CUTOFF
digits), precompute [a**i % c for i in range(32)], and go left to
right 5 bits at a time.
l_divmod():
- The signature changed so that callers who don't want the quotient,
or don't want the remainder, can pass NULL in the slot they don't
want. This saves them from having to declare a vrbl for unwanted
stuff, and remembering to decref it.
long_mod(), long_div(), long_classic_div():
- Adjust to new l_divmod() signature, and simplified as a result.
2004-08-30 02:44:38 +00:00
|
|
|
And, at some places it is assumed that MASK fits in an int, as well.
|
|
|
|
long_pow() requires that SHIFT be divisible by 5. */
|
1991-05-05 20:09:44 +00:00
|
|
|
|
|
|
|
typedef unsigned short digit;
|
1991-05-14 12:06:16 +00:00
|
|
|
typedef unsigned int wdigit; /* digit widened to parameter size */
|
2000-07-08 04:17:21 +00:00
|
|
|
#define BASE_TWODIGITS_TYPE long
|
|
|
|
typedef unsigned BASE_TWODIGITS_TYPE twodigits;
|
|
|
|
typedef BASE_TWODIGITS_TYPE stwodigits; /* signed variant of twodigits */
|
1991-05-05 20:09:44 +00:00
|
|
|
|
|
|
|
#define SHIFT 15
|
|
|
|
#define BASE ((digit)1 << SHIFT)
|
|
|
|
#define MASK ((int)(BASE - 1))
|
|
|
|
|
SF patch 936813: fast modular exponentiation
This checkin is adapted from part 2 (of 3) of Trevor Perrin's patch set.
BACKWARD INCOMPATIBILITY: SHIFT must now be divisible by 5. AFAIK,
nobody will care. long_pow() could be complicated to worm around that,
if necessary.
long_pow():
- BUGFIX: This leaked the base and power when the power was negative
(and so the computation delegated to float pow).
- Instead of doing right-to-left exponentiation, do left-to-right. This
is more efficient for small bases, which is the common case.
- In addition, if the exponent is large (more than FIVEARY_CUTOFF
digits), precompute [a**i % c for i in range(32)], and go left to
right 5 bits at a time.
l_divmod():
- The signature changed so that callers who don't want the quotient,
or don't want the remainder, can pass NULL in the slot they don't
want. This saves them from having to declare a vrbl for unwanted
stuff, and remembering to decref it.
long_mod(), long_div(), long_classic_div():
- Adjust to new l_divmod() signature, and simplified as a result.
2004-08-30 02:44:38 +00:00
|
|
|
#if SHIFT % 5 != 0
|
|
|
|
#error "longobject.c requires that SHIFT be divisible by 5"
|
|
|
|
#endif
|
|
|
|
|
1991-05-05 20:09:44 +00:00
|
|
|
/* Long integer representation.
|
|
|
|
The absolute value of a number is equal to
|
|
|
|
SUM(for i=0 through abs(ob_size)-1) ob_digit[i] * 2**(SHIFT*i)
|
|
|
|
Negative numbers are represented with ob_size < 0;
|
|
|
|
zero is represented by ob_size == 0.
|
|
|
|
In a normalized number, ob_digit[abs(ob_size)-1] (the most significant
|
|
|
|
digit) is never zero. Also, in all cases, for all valid i,
|
|
|
|
0 <= ob_digit[i] <= MASK.
|
2000-07-16 12:04:32 +00:00
|
|
|
The allocation function takes care of allocating extra memory
|
2002-03-02 04:33:09 +00:00
|
|
|
so that ob_digit[0] ... ob_digit[abs(ob_size)-1] are actually available.
|
|
|
|
|
|
|
|
CAUTION: Generic code manipulating subtypes of PyVarObject has to
|
|
|
|
aware that longs abuse ob_size's sign bit.
|
|
|
|
*/
|
1991-05-05 20:09:44 +00:00
|
|
|
|
1992-01-19 16:30:12 +00:00
|
|
|
struct _longobject {
|
2002-03-02 04:33:09 +00:00
|
|
|
PyObject_VAR_HEAD
|
1991-05-05 20:09:44 +00:00
|
|
|
digit ob_digit[1];
|
1992-01-19 16:30:12 +00:00
|
|
|
};
|
1991-05-05 20:09:44 +00:00
|
|
|
|
2002-08-12 07:21:58 +00:00
|
|
|
PyAPI_FUNC(PyLongObject *) _PyLong_New(int);
|
1993-07-28 09:05:47 +00:00
|
|
|
|
2001-09-10 20:52:51 +00:00
|
|
|
/* Return a copy of src. */
|
2002-08-12 07:21:58 +00:00
|
|
|
PyAPI_FUNC(PyObject *) _PyLong_Copy(PyLongObject *src);
|
2001-09-10 20:52:51 +00:00
|
|
|
|
1993-07-28 09:05:47 +00:00
|
|
|
#ifdef __cplusplus
|
|
|
|
}
|
|
|
|
#endif
|
|
|
|
#endif /* !Py_LONGINTREPR_H */
|