git/typthon
1
0
Fork 0
mirror of https://github.com/johndoe6345789/typthon.git synced 2026-04-25 06:05:05 +00:00
Code Issues Packages Projects Releases Wiki Activity
Files
0c32bc85581fbec1f37844073d4c525feaa344b0
typthon/Objects/floatobject.c
copilot-swe-agent[bot] 0f2e4eb3fd Fix final Py_ patterns - uintptr, arithmetic shift, inc files 
- Fixed Py_uintptr_t → Ty_uintptr_t
- Fixed Py_ARITHMETIC_RIGHT_SHIFT → Ty_ARITHMETIC_RIGHT_SHIFT
- Fixed PyTypeObject, PyNumberMethods, PyAsyncMethods in .inc files

Build is nearing completion.

Co-authored-by: johndoe6345789 <224850594+johndoe6345789@users.noreply.github.com>
2025-12-29 18:38:56 +00:00

2646 lines
73 KiB
C
Raw Blame History

/* Float object implementation */
/* XXX There should be overflow checks here, but it's hard to check
for any kind of float exception without losing portability. */
#include "Python.h"
#include "pycore_abstract.h" // _PyNumber_Index()
#include "pycore_dtoa.h" // _Ty_dg_dtoa()
#include "pycore_floatobject.h" // _TyFloat_FormatAdvancedWriter()
#include "pycore_freelist.h" // _Ty_FREELIST_FREE(), _Ty_FREELIST_POP()
#include "pycore_initconfig.h" // _TyStatus_OK()
#include "pycore_long.h" // _TyLong_GetOne()
#include "pycore_modsupport.h" // _TyArg_NoKwnames()
#include "pycore_object.h" // _TyObject_Init(), _PyDebugAllocatorStats()
#include "pycore_pymath.h" // _PY_SHORT_FLOAT_REPR
#include "pycore_pystate.h" // _TyInterpreterState_GET()
#include "pycore_stackref.h" // PyStackRef_AsPyObjectBorrow()
#include "pycore_structseq.h" // _PyStructSequence_FiniBuiltin()
#include <float.h> // DBL_MAX
#include <stdlib.h> // strtol()
/*[clinic input]
class float "TyObject *" "&TyFloat_Type"
[clinic start generated code]*/
/*[clinic end generated code: output=da39a3ee5e6b4b0d input=dd0003f68f144284]*/
#include "clinic/floatobject.c.h"
double
TyFloat_GetMax(void)
{
return DBL_MAX;
}
double
TyFloat_GetMin(void)
{
return DBL_MIN;
}
static TyTypeObject FloatInfoType;
TyDoc_STRVAR(floatinfo__doc__,
"sys.float_info\n\
\n\
A named tuple holding information about the float type. It contains low level\n\
information about the precision and internal representation. Please study\n\
your system's :file:`float.h` for more information.");
static TyStructSequence_Field floatinfo_fields[] = {
{"max", "DBL_MAX -- maximum representable finite float"},
{"max_exp", "DBL_MAX_EXP -- maximum int e such that radix**(e-1) "
"is representable"},
{"max_10_exp", "DBL_MAX_10_EXP -- maximum int e such that 10**e "
"is representable"},
{"min", "DBL_MIN -- Minimum positive normalized float"},
{"min_exp", "DBL_MIN_EXP -- minimum int e such that radix**(e-1) "
"is a normalized float"},
{"min_10_exp", "DBL_MIN_10_EXP -- minimum int e such that 10**e is "
"a normalized float"},
{"dig", "DBL_DIG -- maximum number of decimal digits that "
"can be faithfully represented in a float"},
{"mant_dig", "DBL_MANT_DIG -- mantissa digits"},
{"epsilon", "DBL_EPSILON -- Difference between 1 and the next "
"representable float"},
{"radix", "FLT_RADIX -- radix of exponent"},
{"rounds", "FLT_ROUNDS -- rounding mode used for arithmetic "
"operations"},
{0}
};
static TyStructSequence_Desc floatinfo_desc = {
"sys.float_info", /* name */
floatinfo__doc__, /* doc */
floatinfo_fields, /* fields */
11
};
TyObject *
TyFloat_GetInfo(void)
{
TyObject* floatinfo;
int pos = 0;
floatinfo = TyStructSequence_New(&FloatInfoType);
if (floatinfo == NULL) {
return NULL;
}
#define SetFlag(CALL) \
do { \
TyObject *flag = (CALL); \
if (flag == NULL) { \
Ty_CLEAR(floatinfo); \
return NULL; \
} \
TyStructSequence_SET_ITEM(floatinfo, pos++, flag); \
} while (0)
#define SetIntFlag(FLAG) SetFlag(TyLong_FromLong((FLAG)))
#define SetDblFlag(FLAG) SetFlag(TyFloat_FromDouble((FLAG)))
SetDblFlag(DBL_MAX);
SetIntFlag(DBL_MAX_EXP);
SetIntFlag(DBL_MAX_10_EXP);
SetDblFlag(DBL_MIN);
SetIntFlag(DBL_MIN_EXP);
SetIntFlag(DBL_MIN_10_EXP);
SetIntFlag(DBL_DIG);
SetIntFlag(DBL_MANT_DIG);
SetDblFlag(DBL_EPSILON);
SetIntFlag(FLT_RADIX);
SetIntFlag(FLT_ROUNDS);
#undef SetIntFlag
#undef SetDblFlag
#undef SetFlag
return floatinfo;
}
TyObject *
TyFloat_FromDouble(double fval)
{
PyFloatObject *op = _Ty_FREELIST_POP(PyFloatObject, floats);
if (op == NULL) {
op = PyObject_Malloc(sizeof(PyFloatObject));
if (!op) {
return TyErr_NoMemory();
}
_TyObject_Init((TyObject*)op, &TyFloat_Type);
}
op->ob_fval = fval;
return (TyObject *) op;
}
_PyStackRef _TyFloat_FromDouble_ConsumeInputs(_PyStackRef left, _PyStackRef right, double value)
{
PyStackRef_CLOSE_SPECIALIZED(left, _TyFloat_ExactDealloc);
PyStackRef_CLOSE_SPECIALIZED(right, _TyFloat_ExactDealloc);
return PyStackRef_FromPyObjectSteal(TyFloat_FromDouble(value));
}
static TyObject *
float_from_string_inner(const char *s, Ty_ssize_t len, void *obj)
{
double x;
const char *end;
const char *last = s + len;
/* strip leading whitespace */
while (s < last && Ty_ISSPACE(*s)) {
s++;
}
if (s == last) {
TyErr_Format(TyExc_ValueError,
"could not convert string to float: "
"%R", obj);
return NULL;
}
/* strip trailing whitespace */
while (s < last - 1 && Ty_ISSPACE(last[-1])) {
last--;
}
/* We don't care about overflow or underflow. If the platform
* supports them, infinities and signed zeroes (on underflow) are
* fine. */
x = TyOS_string_to_double(s, (char **)&end, NULL);
if (end != last) {
TyErr_Format(TyExc_ValueError,
"could not convert string to float: "
"%R", obj);
return NULL;
}
else if (x == -1.0 && TyErr_Occurred()) {
return NULL;
}
else {
return TyFloat_FromDouble(x);
}
}
TyObject *
TyFloat_FromString(TyObject *v)
{
const char *s;
TyObject *s_buffer = NULL;
Ty_ssize_t len;
Ty_buffer view = {NULL, NULL};
TyObject *result = NULL;
if (TyUnicode_Check(v)) {
s_buffer = _TyUnicode_TransformDecimalAndSpaceToASCII(v);
if (s_buffer == NULL)
return NULL;
assert(TyUnicode_IS_ASCII(s_buffer));
/* Simply get a pointer to existing ASCII characters. */
s = TyUnicode_AsUTF8AndSize(s_buffer, &len);
assert(s != NULL);
}
else if (TyBytes_Check(v)) {
s = TyBytes_AS_STRING(v);
len = TyBytes_GET_SIZE(v);
}
else if (TyByteArray_Check(v)) {
s = TyByteArray_AS_STRING(v);
len = TyByteArray_GET_SIZE(v);
}
else if (PyObject_GetBuffer(v, &view, PyBUF_SIMPLE) == 0) {
s = (const char *)view.buf;
len = view.len;
/* Copy to NUL-terminated buffer. */
s_buffer = TyBytes_FromStringAndSize(s, len);
if (s_buffer == NULL) {
PyBuffer_Release(&view);
return NULL;
}
s = TyBytes_AS_STRING(s_buffer);
}
else {
TyErr_Format(TyExc_TypeError,
"float() argument must be a string or a real number, not '%.200s'",
Ty_TYPE(v)->tp_name);
return NULL;
}
result = _Ty_string_to_number_with_underscores(s, len, "float", v, v,
float_from_string_inner);
PyBuffer_Release(&view);
Ty_XDECREF(s_buffer);
return result;
}
void
_TyFloat_ExactDealloc(TyObject *obj)
{
assert(TyFloat_CheckExact(obj));
_Ty_FREELIST_FREE(floats, obj, PyObject_Free);
}
static void
float_dealloc(TyObject *op)
{
assert(TyFloat_Check(op));
if (TyFloat_CheckExact(op))
_TyFloat_ExactDealloc(op);
else
Ty_TYPE(op)->tp_free(op);
}
double
TyFloat_AsDouble(TyObject *op)
{
TyNumberMethods *nb;
TyObject *res;
double val;
if (op == NULL) {
TyErr_BadArgument();
return -1;
}
if (TyFloat_Check(op)) {
return TyFloat_AS_DOUBLE(op);
}
nb = Ty_TYPE(op)->tp_as_number;
if (nb == NULL || nb->nb_float == NULL) {
if (nb && nb->nb_index) {
TyObject *res = _PyNumber_Index(op);
if (!res) {
return -1;
}
double val = TyLong_AsDouble(res);
Ty_DECREF(res);
return val;
}
TyErr_Format(TyExc_TypeError, "must be real number, not %.50s",
Ty_TYPE(op)->tp_name);
return -1;
}
res = (*nb->nb_float) (op);
if (res == NULL) {
return -1;
}
if (!TyFloat_CheckExact(res)) {
if (!TyFloat_Check(res)) {
TyErr_Format(TyExc_TypeError,
"%.50s.__float__ returned non-float (type %.50s)",
Ty_TYPE(op)->tp_name, Ty_TYPE(res)->tp_name);
Ty_DECREF(res);
return -1;
}
if (TyErr_WarnFormat(TyExc_DeprecationWarning, 1,
"%.50s.__float__ returned non-float (type %.50s). "
"The ability to return an instance of a strict subclass of float "
"is deprecated, and may be removed in a future version of Python.",
Ty_TYPE(op)->tp_name, Ty_TYPE(res)->tp_name)) {
Ty_DECREF(res);
return -1;
}
}
val = TyFloat_AS_DOUBLE(res);
Ty_DECREF(res);
return val;
}
/* Macro and helper that convert TyObject obj to a C double and store
the value in dbl. If conversion to double raises an exception, obj is
set to NULL, and the function invoking this macro returns NULL. If
obj is not of float or int type, Ty_NotImplemented is incref'ed,
stored in obj, and returned from the function invoking this macro.
*/
#define CONVERT_TO_DOUBLE(obj, dbl) \
if (TyFloat_Check(obj)) \
dbl = TyFloat_AS_DOUBLE(obj); \
else if (_Ty_convert_int_to_double(&(obj), &(dbl)) < 0) \
return obj;
/* Methods */
int
_Ty_convert_int_to_double(TyObject **v, double *dbl)
{
TyObject *obj = *v;
if (TyLong_Check(obj)) {
*dbl = TyLong_AsDouble(obj);
if (*dbl == -1.0 && TyErr_Occurred()) {
*v = NULL;
return -1;
}
}
else {
*v = Ty_NewRef(Ty_NotImplemented);
return -1;
}
return 0;
}
static TyObject *
float_repr(TyObject *op)
{
PyFloatObject *v = _TyFloat_CAST(op);
TyObject *result;
char *buf;
buf = TyOS_double_to_string(TyFloat_AS_DOUBLE(v),
'r', 0,
Ty_DTSF_ADD_DOT_0,
NULL);
if (!buf)
return TyErr_NoMemory();
result = _TyUnicode_FromASCII(buf, strlen(buf));
TyMem_Free(buf);
return result;
}
/* Comparison is pretty much a nightmare. When comparing float to float,
* we do it as straightforwardly (and long-windedly) as conceivable, so
* that, e.g., Python x == y delivers the same result as the platform
* C x == y when x and/or y is a NaN.
* When mixing float with an integer type, there's no good *uniform* approach.
* Converting the double to an integer obviously doesn't work, since we
* may lose info from fractional bits. Converting the integer to a double
* also has two failure modes: (1) an int may trigger overflow (too
* large to fit in the dynamic range of a C double); (2) even a C long may have
* more bits than fit in a C double (e.g., on a 64-bit box long may have
* 63 bits of precision, but a C double probably has only 53), and then
* we can falsely claim equality when low-order integer bits are lost by
* coercion to double. So this part is painful too.
*/
static TyObject*
float_richcompare(TyObject *v, TyObject *w, int op)
{
double i, j;
int r = 0;
assert(TyFloat_Check(v));
i = TyFloat_AS_DOUBLE(v);
/* Switch on the type of w. Set i and j to doubles to be compared,
* and op to the richcomp to use.
*/
if (TyFloat_Check(w))
j = TyFloat_AS_DOUBLE(w);
else if (!isfinite(i)) {
if (TyLong_Check(w))
/* If i is an infinity, its magnitude exceeds any
* finite integer, so it doesn't matter which int we
* compare i with. If i is a NaN, similarly.
*/
j = 0.0;
else
goto Unimplemented;
}
else if (TyLong_Check(w)) {
int vsign = i == 0.0 ? 0 : i < 0.0 ? -1 : 1;
int wsign;
int exponent;
(void)TyLong_GetSign(w, &wsign);
if (vsign != wsign) {
/* Magnitudes are irrelevant -- the signs alone
* determine the outcome.
*/
i = (double)vsign;
j = (double)wsign;
goto Compare;
}
/* The signs are the same. */
/* Convert w to a double if it fits. In particular, 0 fits. */
int64_t nbits64 = _TyLong_NumBits(w);
assert(nbits64 >= 0);
assert(!TyErr_Occurred());
if (nbits64 > DBL_MAX_EXP) {
/* This Python integer is larger than any finite C double.
* Replace with little doubles
* that give the same outcome -- w is so large that
* its magnitude must exceed the magnitude of any
* finite float.
*/
i = (double)vsign;
assert(wsign != 0);
j = wsign * 2.0;
goto Compare;
}
int nbits = (int)nbits64;
if (nbits <= 48) {
j = TyLong_AsDouble(w);
/* It's impossible that <= 48 bits overflowed. */
assert(j != -1.0 || ! TyErr_Occurred());
goto Compare;
}
assert(wsign != 0); /* else nbits was 0 */
assert(vsign != 0); /* if vsign were 0, then since wsign is
* not 0, we would have taken the
* vsign != wsign branch at the start */
/* We want to work with non-negative numbers. */
if (vsign < 0) {
/* "Multiply both sides" by -1; this also swaps the
* comparator.
*/
i = -i;
op = _Ty_SwappedOp[op];
}
assert(i > 0.0);
(void) frexp(i, &exponent);
/* exponent is the # of bits in v before the radix point;
* we know that nbits (the # of bits in w) > 48 at this point
*/
if (exponent < nbits) {
i = 1.0;
j = 2.0;
goto Compare;
}
if (exponent > nbits) {
i = 2.0;
j = 1.0;
goto Compare;
}
/* v and w have the same number of bits before the radix
* point. Construct two ints that have the same comparison
* outcome.
*/
{
double fracpart;
double intpart;
TyObject *result = NULL;
TyObject *vv = NULL;
TyObject *ww = w;
if (wsign < 0) {
ww = PyNumber_Negative(w);
if (ww == NULL)
goto Error;
}
else
Ty_INCREF(ww);
fracpart = modf(i, &intpart);
vv = TyLong_FromDouble(intpart);
if (vv == NULL)
goto Error;
if (fracpart != 0.0) {
/* Shift left, and or a 1 bit into vv
* to represent the lost fraction.
*/
TyObject *temp;
temp = _TyLong_Lshift(ww, 1);
if (temp == NULL)
goto Error;
Ty_SETREF(ww, temp);
temp = _TyLong_Lshift(vv, 1);
if (temp == NULL)
goto Error;
Ty_SETREF(vv, temp);
temp = PyNumber_Or(vv, _TyLong_GetOne());
if (temp == NULL)
goto Error;
Ty_SETREF(vv, temp);
}
r = PyObject_RichCompareBool(vv, ww, op);
if (r < 0)
goto Error;
result = TyBool_FromLong(r);
Error:
Ty_XDECREF(vv);
Ty_XDECREF(ww);
return result;
}
} /* else if (TyLong_Check(w)) */
else /* w isn't float or int */
goto Unimplemented;
Compare:
switch (op) {
case Py_EQ:
r = i == j;
break;
case Py_NE:
r = i != j;
break;
case Py_LE:
r = i <= j;
break;
case Py_GE:
r = i >= j;
break;
case Py_LT:
r = i < j;
break;
case Py_GT:
r = i > j;
break;
}
return TyBool_FromLong(r);
Unimplemented:
Py_RETURN_NOTIMPLEMENTED;
}
static Ty_hash_t
float_hash(TyObject *op)
{
PyFloatObject *v = _TyFloat_CAST(op);
return _Ty_HashDouble(op, v->ob_fval);
}
static TyObject *
float_add(TyObject *v, TyObject *w)
{
double a,b;
CONVERT_TO_DOUBLE(v, a);
CONVERT_TO_DOUBLE(w, b);
a = a + b;
return TyFloat_FromDouble(a);
}
static TyObject *
float_sub(TyObject *v, TyObject *w)
{
double a,b;
CONVERT_TO_DOUBLE(v, a);
CONVERT_TO_DOUBLE(w, b);
a = a - b;
return TyFloat_FromDouble(a);
}
static TyObject *
float_mul(TyObject *v, TyObject *w)
{
double a,b;
CONVERT_TO_DOUBLE(v, a);
CONVERT_TO_DOUBLE(w, b);
a = a * b;
return TyFloat_FromDouble(a);
}
static TyObject *
float_div(TyObject *v, TyObject *w)
{
double a,b;
CONVERT_TO_DOUBLE(v, a);
CONVERT_TO_DOUBLE(w, b);
if (b == 0.0) {
TyErr_SetString(TyExc_ZeroDivisionError,
"division by zero");
return NULL;
}
a = a / b;
return TyFloat_FromDouble(a);
}
static TyObject *
float_rem(TyObject *v, TyObject *w)
{
double vx, wx;
double mod;
CONVERT_TO_DOUBLE(v, vx);
CONVERT_TO_DOUBLE(w, wx);
if (wx == 0.0) {
TyErr_SetString(TyExc_ZeroDivisionError,
"division by zero");
return NULL;
}
mod = fmod(vx, wx);
if (mod) {
/* ensure the remainder has the same sign as the denominator */
if ((wx < 0) != (mod < 0)) {
mod += wx;
}
}
else {
/* the remainder is zero, and in the presence of signed zeroes
fmod returns different results across platforms; ensure
it has the same sign as the denominator. */
mod = copysign(0.0, wx);
}
return TyFloat_FromDouble(mod);
}
static void
_float_div_mod(double vx, double wx, double *floordiv, double *mod)
{
double div;
*mod = fmod(vx, wx);
/* fmod is typically exact, so vx-mod is *mathematically* an
exact multiple of wx. But this is fp arithmetic, and fp
vx - mod is an approximation; the result is that div may
not be an exact integral value after the division, although
it will always be very close to one.
*/
div = (vx - *mod) / wx;
if (*mod) {
/* ensure the remainder has the same sign as the denominator */
if ((wx < 0) != (*mod < 0)) {
*mod += wx;
div -= 1.0;
}
}
else {
/* the remainder is zero, and in the presence of signed zeroes
fmod returns different results across platforms; ensure
it has the same sign as the denominator. */
*mod = copysign(0.0, wx);
}
/* snap quotient to nearest integral value */
if (div) {
*floordiv = floor(div);
if (div - *floordiv > 0.5) {
*floordiv += 1.0;
}
}
else {
/* div is zero - get the same sign as the true quotient */
*floordiv = copysign(0.0, vx / wx); /* zero w/ sign of vx/wx */
}
}
static TyObject *
float_divmod(TyObject *v, TyObject *w)
{
double vx, wx;
double mod, floordiv;
CONVERT_TO_DOUBLE(v, vx);
CONVERT_TO_DOUBLE(w, wx);
if (wx == 0.0) {
TyErr_SetString(TyExc_ZeroDivisionError, "division by zero");
return NULL;
}
_float_div_mod(vx, wx, &floordiv, &mod);
return Ty_BuildValue("(dd)", floordiv, mod);
}
static TyObject *
float_floor_div(TyObject *v, TyObject *w)
{
double vx, wx;
double mod, floordiv;
CONVERT_TO_DOUBLE(v, vx);
CONVERT_TO_DOUBLE(w, wx);
if (wx == 0.0) {
TyErr_SetString(TyExc_ZeroDivisionError, "division by zero");
return NULL;
}
_float_div_mod(vx, wx, &floordiv, &mod);
return TyFloat_FromDouble(floordiv);
}
/* determine whether x is an odd integer or not; assumes that
x is not an infinity or nan. */
#define DOUBLE_IS_ODD_INTEGER(x) (fmod(fabs(x), 2.0) == 1.0)
static TyObject *
float_pow(TyObject *v, TyObject *w, TyObject *z)
{
double iv, iw, ix;
int negate_result = 0;
if ((TyObject *)z != Ty_None) {
TyErr_SetString(TyExc_TypeError, "pow() 3rd argument not "
"allowed unless all arguments are integers");
return NULL;
}
CONVERT_TO_DOUBLE(v, iv);
CONVERT_TO_DOUBLE(w, iw);
/* Sort out special cases here instead of relying on pow() */
if (iw == 0) { /* v**0 is 1, even 0**0 */
return TyFloat_FromDouble(1.0);
}
if (isnan(iv)) { /* nan**w = nan, unless w == 0 */
return TyFloat_FromDouble(iv);
}
if (isnan(iw)) { /* v**nan = nan, unless v == 1; 1**nan = 1 */
return TyFloat_FromDouble(iv == 1.0 ? 1.0 : iw);
}
if (isinf(iw)) {
/* v**inf is: 0.0 if abs(v) < 1; 1.0 if abs(v) == 1; inf if
* abs(v) > 1 (including case where v infinite)
*
* v**-inf is: inf if abs(v) < 1; 1.0 if abs(v) == 1; 0.0 if
* abs(v) > 1 (including case where v infinite)
*/
iv = fabs(iv);
if (iv == 1.0)
return TyFloat_FromDouble(1.0);
else if ((iw > 0.0) == (iv > 1.0))
return TyFloat_FromDouble(fabs(iw)); /* return inf */
else
return TyFloat_FromDouble(0.0);
}
if (isinf(iv)) {
/* (+-inf)**w is: inf for w positive, 0 for w negative; in
* both cases, we need to add the appropriate sign if w is
* an odd integer.
*/
int iw_is_odd = DOUBLE_IS_ODD_INTEGER(iw);
if (iw > 0.0)
return TyFloat_FromDouble(iw_is_odd ? iv : fabs(iv));
else
return TyFloat_FromDouble(iw_is_odd ?
copysign(0.0, iv) : 0.0);
}
if (iv == 0.0) { /* 0**w is: 0 for w positive, 1 for w zero
(already dealt with above), and an error
if w is negative. */
int iw_is_odd = DOUBLE_IS_ODD_INTEGER(iw);
if (iw < 0.0) {
TyErr_SetString(TyExc_ZeroDivisionError,
"zero to a negative power");
return NULL;
}
/* use correct sign if iw is odd */
return TyFloat_FromDouble(iw_is_odd ? iv : 0.0);
}
if (iv < 0.0) {
/* Whether this is an error is a mess, and bumps into libm
* bugs so we have to figure it out ourselves.
*/
if (iw != floor(iw)) {
/* Negative numbers raised to fractional powers
* become complex.
*/
return TyComplex_Type.tp_as_number->nb_power(v, w, z);
}
/* iw is an exact integer, albeit perhaps a very large
* one. Replace iv by its absolute value and remember
* to negate the pow result if iw is odd.
*/
iv = -iv;
negate_result = DOUBLE_IS_ODD_INTEGER(iw);
}
if (iv == 1.0) { /* 1**w is 1, even 1**inf and 1**nan */
/* (-1) ** large_integer also ends up here. Here's an
* extract from the comments for the previous
* implementation explaining why this special case is
* necessary:
*
* -1 raised to an exact integer should never be exceptional.
* Alas, some libms (chiefly glibc as of early 2003) return
* NaN and set EDOM on pow(-1, large_int) if the int doesn't
* happen to be representable in a *C* integer. That's a
* bug.
*/
return TyFloat_FromDouble(negate_result ? -1.0 : 1.0);
}
/* Now iv and iw are finite, iw is nonzero, and iv is
* positive and not equal to 1.0. We finally allow
* the platform pow to step in and do the rest.
*/
errno = 0;
ix = pow(iv, iw);
_Ty_ADJUST_ERANGE1(ix);
if (negate_result)
ix = -ix;
if (errno != 0) {
/* We don't expect any errno value other than ERANGE, but
* the range of libm bugs appears unbounded.
*/
TyErr_SetFromErrno(errno == ERANGE ? TyExc_OverflowError :
TyExc_ValueError);
return NULL;
}
return TyFloat_FromDouble(ix);
}
#undef DOUBLE_IS_ODD_INTEGER
static TyObject *
float_neg(TyObject *op)
{
PyFloatObject *v = _TyFloat_CAST(op);
return TyFloat_FromDouble(-v->ob_fval);
}
static TyObject *
float_abs(TyObject *op)
{
PyFloatObject *v = _TyFloat_CAST(op);
return TyFloat_FromDouble(fabs(v->ob_fval));
}
static int
float_bool(TyObject *op)
{
PyFloatObject *v = _TyFloat_CAST(op);
return v->ob_fval != 0.0;
}
/*[clinic input]
float.is_integer
Return True if the float is an integer.
[clinic start generated code]*/
static TyObject *
float_is_integer_impl(TyObject *self)
/*[clinic end generated code: output=7112acf95a4d31ea input=311810d3f777e10d]*/
{
double x = TyFloat_AsDouble(self);
TyObject *o;
if (x == -1.0 && TyErr_Occurred())
return NULL;
if (!isfinite(x))
Py_RETURN_FALSE;
errno = 0;
o = (floor(x) == x) ? Ty_True : Ty_False;
if (errno != 0) {
TyErr_SetFromErrno(errno == ERANGE ? TyExc_OverflowError :
TyExc_ValueError);
return NULL;
}
return Ty_NewRef(o);
}
/*[clinic input]
float.__trunc__
Return the Integral closest to x between 0 and x.
[clinic start generated code]*/
static TyObject *
float___trunc___impl(TyObject *self)
/*[clinic end generated code: output=dd3e289dd4c6b538 input=591b9ba0d650fdff]*/
{
return TyLong_FromDouble(TyFloat_AS_DOUBLE(self));
}
/*[clinic input]
float.__floor__
Return the floor as an Integral.
[clinic start generated code]*/
static TyObject *
float___floor___impl(TyObject *self)
/*[clinic end generated code: output=e0551dbaea8c01d1 input=77bb13eb12e268df]*/
{
double x = TyFloat_AS_DOUBLE(self);
return TyLong_FromDouble(floor(x));
}
/*[clinic input]
float.__ceil__
Return the ceiling as an Integral.
[clinic start generated code]*/
static TyObject *
float___ceil___impl(TyObject *self)
/*[clinic end generated code: output=a2fd8858f73736f9 input=79e41ae94aa0a516]*/
{
double x = TyFloat_AS_DOUBLE(self);
return TyLong_FromDouble(ceil(x));
}
/* double_round: rounds a finite double to the closest multiple of
10**-ndigits; here ndigits is within reasonable bounds (typically, -308 <=
ndigits <= 323). Returns a Python float, or sets a Python error and
returns NULL on failure (OverflowError and memory errors are possible). */
#if _PY_SHORT_FLOAT_REPR == 1
/* version of double_round that uses the correctly-rounded string<->double
conversions from Python/dtoa.c */
static TyObject *
double_round(double x, int ndigits) {
double rounded;
Ty_ssize_t buflen, mybuflen=100;
char *buf, *buf_end, shortbuf[100], *mybuf=shortbuf;
int decpt, sign;
TyObject *result = NULL;
_Ty_SET_53BIT_PRECISION_HEADER;
/* round to a decimal string */
_Ty_SET_53BIT_PRECISION_START;
buf = _Ty_dg_dtoa(x, 3, ndigits, &decpt, &sign, &buf_end);
_Ty_SET_53BIT_PRECISION_END;
if (buf == NULL) {
TyErr_NoMemory();
return NULL;
}
/* Get new buffer if shortbuf is too small. Space needed <= buf_end -
buf + 8: (1 extra for '0', 1 for sign, 5 for exp, 1 for '\0'). */
buflen = buf_end - buf;
if (buflen + 8 > mybuflen) {
mybuflen = buflen+8;
mybuf = (char *)TyMem_Malloc(mybuflen);
if (mybuf == NULL) {
TyErr_NoMemory();
goto exit;
}
}
/* copy buf to mybuf, adding exponent, sign and leading 0 */
TyOS_snprintf(mybuf, mybuflen, "%s0%se%d", (sign ? "-" : ""),
buf, decpt - (int)buflen);
/* and convert the resulting string back to a double */
errno = 0;
_Ty_SET_53BIT_PRECISION_START;
rounded = _Ty_dg_strtod(mybuf, NULL);
_Ty_SET_53BIT_PRECISION_END;
if (errno == ERANGE && fabs(rounded) >= 1.)
TyErr_SetString(TyExc_OverflowError,
"rounded value too large to represent");
else
result = TyFloat_FromDouble(rounded);
/* done computing value; now clean up */
if (mybuf != shortbuf)
TyMem_Free(mybuf);
exit:
_Ty_dg_freedtoa(buf);
return result;
}
#else // _PY_SHORT_FLOAT_REPR == 0
/* fallback version, to be used when correctly rounded binary<->decimal
conversions aren't available */
static TyObject *
double_round(double x, int ndigits) {
double pow1, pow2, y, z;
if (ndigits >= 0) {
if (ndigits > 22) {
/* pow1 and pow2 are each safe from overflow, but
pow1*pow2 ~= pow(10.0, ndigits) might overflow */
pow1 = pow(10.0, (double)(ndigits-22));
pow2 = 1e22;
}
else {
pow1 = pow(10.0, (double)ndigits);
pow2 = 1.0;
}
y = (x*pow1)*pow2;
/* if y overflows, then rounded value is exactly x */
if (!isfinite(y))
return TyFloat_FromDouble(x);
}
else {
pow1 = pow(10.0, (double)-ndigits);
pow2 = 1.0; /* unused; silences a gcc compiler warning */
y = x / pow1;
}
z = round(y);
if (fabs(y-z) == 0.5)
/* halfway between two integers; use round-half-even */
z = 2.0*round(y/2.0);
if (ndigits >= 0)
z = (z / pow2) / pow1;
else
z *= pow1;
/* if computation resulted in overflow, raise OverflowError */
if (!isfinite(z)) {
TyErr_SetString(TyExc_OverflowError,
"overflow occurred during round");
return NULL;
}
return TyFloat_FromDouble(z);
}
#endif // _PY_SHORT_FLOAT_REPR == 0
/* round a Python float v to the closest multiple of 10**-ndigits */
/*[clinic input]
float.__round__
ndigits as o_ndigits: object = None
/
Return the Integral closest to x, rounding half toward even.
When an argument is passed, work like built-in round(x, ndigits).
[clinic start generated code]*/
static TyObject *
float___round___impl(TyObject *self, TyObject *o_ndigits)
/*[clinic end generated code: output=374c36aaa0f13980 input=fc0fe25924fbc9ed]*/
{
double x, rounded;
Ty_ssize_t ndigits;
x = TyFloat_AsDouble(self);
if (o_ndigits == Ty_None) {
/* single-argument round or with None ndigits:
* round to nearest integer */
rounded = round(x);
if (fabs(x-rounded) == 0.5)
/* halfway case: round to even */
rounded = 2.0*round(x/2.0);
return TyLong_FromDouble(rounded);
}
/* interpret second argument as a Ty_ssize_t; clips on overflow */
ndigits = PyNumber_AsSsize_t(o_ndigits, NULL);
if (ndigits == -1 && TyErr_Occurred())
return NULL;
/* nans and infinities round to themselves */
if (!isfinite(x))
return TyFloat_FromDouble(x);
/* Deal with extreme values for ndigits. For ndigits > NDIGITS_MAX, x
always rounds to itself. For ndigits < NDIGITS_MIN, x always
rounds to +-0.0. Here 0.30103 is an upper bound for log10(2). */
#define NDIGITS_MAX ((int)((DBL_MANT_DIG-DBL_MIN_EXP) * 0.30103))
#define NDIGITS_MIN (-(int)((DBL_MAX_EXP + 1) * 0.30103))
if (ndigits > NDIGITS_MAX)
/* return x */
return TyFloat_FromDouble(x);
else if (ndigits < NDIGITS_MIN)
/* return 0.0, but with sign of x */
return TyFloat_FromDouble(0.0*x);
else
/* finite x, and ndigits is not unreasonably large */
return double_round(x, (int)ndigits);
#undef NDIGITS_MAX
#undef NDIGITS_MIN
}
static TyObject *
float_float(TyObject *v)
{
if (TyFloat_CheckExact(v)) {
return Ty_NewRef(v);
}
else {
return TyFloat_FromDouble(((PyFloatObject *)v)->ob_fval);
}
}
/*[clinic input]
float.conjugate
Return self, the complex conjugate of any float.
[clinic start generated code]*/
static TyObject *
float_conjugate_impl(TyObject *self)
/*[clinic end generated code: output=8ca292c2479194af input=82ba6f37a9ff91dd]*/
{
return float_float(self);
}
/* turn ASCII hex characters into integer values and vice versa */
static char
char_from_hex(int x)
{
assert(0 <= x && x < 16);
return Ty_hexdigits[x];
}
/* This table maps characters to their hexadecimal values, only
* works with encodings whose lower half is ASCII (like UTF-8).
* '0' maps to 0, ..., '9' maps to 9.
* 'a' and 'A' map to 10, ..., 'f' and 'F' map to 15.
* All other indices map to -1.
*/
static const int
_CHAR_TO_HEX[256] = {
-1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1,
-1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1,
-1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1,
0, 1, 2, 3, 4, 5, 6, 7, 8, 9, -1, -1, -1, -1, -1, -1,
-1, 10, 11, 12, 13, 14, 15, -1, -1, -1, -1, -1, -1, -1, -1, -1,
-1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1,
-1, 10, 11, 12, 13, 14, 15, -1, -1, -1, -1, -1, -1, -1, -1, -1,
-1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1,
-1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1,
-1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1,
-1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1,
-1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1,
-1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1,
-1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1,
-1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1,
-1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1,
};
/* Convert a character to its hexadecimal value, or -1 if it's not a
* valid hexadecimal character, only works with encodings whose lower
* half is ASCII (like UTF-8).
*/
static int
hex_from_char(unsigned char c) {
return _CHAR_TO_HEX[c];
}
/* convert a float to a hexadecimal string */
/* TOHEX_NBITS is DBL_MANT_DIG rounded up to the next integer
of the form 4k+1. */
#define TOHEX_NBITS DBL_MANT_DIG + 3 - (DBL_MANT_DIG+2)%4
/*[clinic input]
float.hex
Return a hexadecimal representation of a floating-point number.
>>> (-0.1).hex()
'-0x1.999999999999ap-4'
>>> 3.14159.hex()
'0x1.921f9f01b866ep+1'
[clinic start generated code]*/
static TyObject *
float_hex_impl(TyObject *self)
/*[clinic end generated code: output=0ebc9836e4d302d4 input=bec1271a33d47e67]*/
{
double x, m;
int e, shift, i, si, esign;
/* Space for 1+(TOHEX_NBITS-1)/4 digits, a decimal point, and the
trailing NUL byte. */
char s[(TOHEX_NBITS-1)/4+3];
CONVERT_TO_DOUBLE(self, x);
if (isnan(x) || isinf(x))
return float_repr(self);
if (x == 0.0) {
if (copysign(1.0, x) == -1.0)
return TyUnicode_FromString("-0x0.0p+0");
else
return TyUnicode_FromString("0x0.0p+0");
}
m = frexp(fabs(x), &e);
shift = 1 - Ty_MAX(DBL_MIN_EXP - e, 0);
m = ldexp(m, shift);
e -= shift;
si = 0;
s[si] = char_from_hex((int)m);
si++;
m -= (int)m;
s[si] = '.';
si++;
for (i=0; i < (TOHEX_NBITS-1)/4; i++) {
m *= 16.0;
s[si] = char_from_hex((int)m);
si++;
m -= (int)m;
}
s[si] = '\0';
if (e < 0) {
esign = (int)'-';
e = -e;
}
else
esign = (int)'+';
if (x < 0.0)
return TyUnicode_FromFormat("-0x%sp%c%d", s, esign, e);
else
return TyUnicode_FromFormat("0x%sp%c%d", s, esign, e);
}
/* Convert a hexadecimal string to a float. */
/*[clinic input]
@classmethod
float.fromhex
string: object
/
Create a floating-point number from a hexadecimal string.
>>> float.fromhex('0x1.ffffp10')
2047.984375
>>> float.fromhex('-0x1p-1074')
-5e-324
[clinic start generated code]*/
static TyObject *
float_fromhex_impl(TyTypeObject *type, TyObject *string)
/*[clinic end generated code: output=c54b4923552e5af5 input=0407bebd354bca89]*/
{
TyObject *result;
double x;
long exp, top_exp, lsb, key_digit;
const char *s, *coeff_start, *s_store, *coeff_end, *exp_start, *s_end;
int half_eps, digit, round_up, negate=0;
Ty_ssize_t length, ndigits, fdigits, i;
/*
* For the sake of simplicity and correctness, we impose an artificial
* limit on ndigits, the total number of hex digits in the coefficient
* The limit is chosen to ensure that, writing exp for the exponent,
*
* (1) if exp > LONG_MAX/2 then the value of the hex string is
* guaranteed to overflow (provided it's nonzero)
*
* (2) if exp < LONG_MIN/2 then the value of the hex string is
* guaranteed to underflow to 0.
*
* (3) if LONG_MIN/2 <= exp <= LONG_MAX/2 then there's no danger of
* overflow in the calculation of exp and top_exp below.
*
* More specifically, ndigits is assumed to satisfy the following
* inequalities:
*
* 4*ndigits <= DBL_MIN_EXP - DBL_MANT_DIG - LONG_MIN/2
* 4*ndigits <= LONG_MAX/2 + 1 - DBL_MAX_EXP
*
* If either of these inequalities is not satisfied, a ValueError is
* raised. Otherwise, write x for the value of the hex string, and
* assume x is nonzero. Then
*
* 2**(exp-4*ndigits) <= |x| < 2**(exp+4*ndigits).
*
* Now if exp > LONG_MAX/2 then:
*
* exp - 4*ndigits >= LONG_MAX/2 + 1 - (LONG_MAX/2 + 1 - DBL_MAX_EXP)
* = DBL_MAX_EXP
*
* so |x| >= 2**DBL_MAX_EXP, which is too large to be stored in C
* double, so overflows. If exp < LONG_MIN/2, then
*
* exp + 4*ndigits <= LONG_MIN/2 - 1 + (
* DBL_MIN_EXP - DBL_MANT_DIG - LONG_MIN/2)
* = DBL_MIN_EXP - DBL_MANT_DIG - 1
*
* and so |x| < 2**(DBL_MIN_EXP-DBL_MANT_DIG-1), hence underflows to 0
* when converted to a C double.
*
* It's easy to show that if LONG_MIN/2 <= exp <= LONG_MAX/2 then both
* exp+4*ndigits and exp-4*ndigits are within the range of a long.
*/
s = TyUnicode_AsUTF8AndSize(string, &length);
if (s == NULL)
return NULL;
s_end = s + length;
/********************
* Parse the string *
********************/
/* leading whitespace */
while (Ty_ISSPACE(*s))
s++;
/* infinities and nans */
x = _Ty_parse_inf_or_nan(s, (char **)&coeff_end);
if (coeff_end != s) {
s = coeff_end;
goto finished;
}
/* optional sign */
if (*s == '-') {
s++;
negate = 1;
}
else if (*s == '+')
s++;
/* [0x] */
s_store = s;
if (*s == '0') {
s++;
if (*s == 'x' || *s == 'X')
s++;
else
s = s_store;
}
/* coefficient: <integer> [. <fraction>] */
coeff_start = s;
while (hex_from_char(*s) >= 0)
s++;
s_store = s;
if (*s == '.') {
s++;
while (hex_from_char(*s) >= 0)
s++;
coeff_end = s-1;
}
else
coeff_end = s;
/* ndigits = total # of hex digits; fdigits = # after point */
ndigits = coeff_end - coeff_start;
fdigits = coeff_end - s_store;
if (ndigits == 0)
goto parse_error;
if (ndigits > Ty_MIN(DBL_MIN_EXP - DBL_MANT_DIG - LONG_MIN/2,
LONG_MAX/2 + 1 - DBL_MAX_EXP)/4)
goto insane_length_error;
/* [p <exponent>] */
if (*s == 'p' || *s == 'P') {
s++;
exp_start = s;
if (*s == '-' || *s == '+')
s++;
if (!('0' <= *s && *s <= '9'))
goto parse_error;
s++;
while ('0' <= *s && *s <= '9')
s++;
exp = strtol(exp_start, NULL, 10);
}
else
exp = 0;
/* for 0 <= j < ndigits, HEX_DIGIT(j) gives the jth most significant digit */
#define HEX_DIGIT(j) hex_from_char(*((j) < fdigits ? \
coeff_end-(j) : \
coeff_end-1-(j)))
/*******************************************
* Compute rounded value of the hex string *
*******************************************/
/* Discard leading zeros, and catch extreme overflow and underflow */
while (ndigits > 0 && HEX_DIGIT(ndigits-1) == 0)
ndigits--;
if (ndigits == 0 || exp < LONG_MIN/2) {
x = 0.0;
goto finished;
}
if (exp > LONG_MAX/2)
goto overflow_error;
/* Adjust exponent for fractional part. */
exp = exp - 4*((long)fdigits);
/* top_exp = 1 more than exponent of most sig. bit of coefficient */
top_exp = exp + 4*((long)ndigits - 1);
for (digit = HEX_DIGIT(ndigits-1); digit != 0; digit /= 2)
top_exp++;
/* catch almost all nonextreme cases of overflow and underflow here */
if (top_exp < DBL_MIN_EXP - DBL_MANT_DIG) {
x = 0.0;
goto finished;
}
if (top_exp > DBL_MAX_EXP)
goto overflow_error;
/* lsb = exponent of least significant bit of the *rounded* value.
This is top_exp - DBL_MANT_DIG unless result is subnormal. */
lsb = Ty_MAX(top_exp, (long)DBL_MIN_EXP) - DBL_MANT_DIG;
x = 0.0;
if (exp >= lsb) {
/* no rounding required */
for (i = ndigits-1; i >= 0; i--)
x = 16.0*x + HEX_DIGIT(i);
x = ldexp(x, (int)(exp));
goto finished;
}
/* rounding required. key_digit is the index of the hex digit
containing the first bit to be rounded away. */
half_eps = 1 << (int)((lsb - exp - 1) % 4);
key_digit = (lsb - exp - 1) / 4;
for (i = ndigits-1; i > key_digit; i--)
x = 16.0*x + HEX_DIGIT(i);
digit = HEX_DIGIT(key_digit);
x = 16.0*x + (double)(digit & (16-2*half_eps));
/* round-half-even: round up if bit lsb-1 is 1 and at least one of
bits lsb, lsb-2, lsb-3, lsb-4, ... is 1. */
if ((digit & half_eps) != 0) {
round_up = 0;
if ((digit & (3*half_eps-1)) != 0 || (half_eps == 8 &&
key_digit+1 < ndigits && (HEX_DIGIT(key_digit+1) & 1) != 0))
round_up = 1;
else
for (i = key_digit-1; i >= 0; i--)
if (HEX_DIGIT(i) != 0) {
round_up = 1;
break;
}
if (round_up) {
x += 2*half_eps;
if (top_exp == DBL_MAX_EXP &&
x == ldexp((double)(2*half_eps), DBL_MANT_DIG))
/* overflow corner case: pre-rounded value <
2**DBL_MAX_EXP; rounded=2**DBL_MAX_EXP. */
goto overflow_error;
}
}
x = ldexp(x, (int)(exp+4*key_digit));
finished:
/* optional trailing whitespace leading to the end of the string */
while (Ty_ISSPACE(*s))
s++;
if (s != s_end)
goto parse_error;
result = TyFloat_FromDouble(negate ? -x : x);
if (type != &TyFloat_Type && result != NULL) {
Ty_SETREF(result, PyObject_CallOneArg((TyObject *)type, result));
}
return result;
overflow_error:
TyErr_SetString(TyExc_OverflowError,
"hexadecimal value too large to represent as a float");
return NULL;
parse_error:
TyErr_SetString(TyExc_ValueError,
"invalid hexadecimal floating-point string");
return NULL;
insane_length_error:
TyErr_SetString(TyExc_ValueError,
"hexadecimal string too long to convert");
return NULL;
}
/*[clinic input]
float.as_integer_ratio
Return a pair of integers, whose ratio is exactly equal to the original float.
The ratio is in lowest terms and has a positive denominator. Raise
OverflowError on infinities and a ValueError on NaNs.
>>> (10.0).as_integer_ratio()
(10, 1)
>>> (0.0).as_integer_ratio()
(0, 1)
>>> (-.25).as_integer_ratio()
(-1, 4)
[clinic start generated code]*/
static TyObject *
float_as_integer_ratio_impl(TyObject *self)
/*[clinic end generated code: output=65f25f0d8d30a712 input=d5ba7765655d75bd]*/
{
double self_double;
double float_part;
int exponent;
int i;
TyObject *py_exponent = NULL;
TyObject *numerator = NULL;
TyObject *denominator = NULL;
TyObject *result_pair = NULL;
TyNumberMethods *long_methods = TyLong_Type.tp_as_number;
CONVERT_TO_DOUBLE(self, self_double);
if (isinf(self_double)) {
TyErr_SetString(TyExc_OverflowError,
"cannot convert Infinity to integer ratio");
return NULL;
}
if (isnan(self_double)) {
TyErr_SetString(TyExc_ValueError,
"cannot convert NaN to integer ratio");
return NULL;
}
float_part = frexp(self_double, &exponent); /* self_double == float_part * 2**exponent exactly */
for (i=0; i<300 && float_part != floor(float_part) ; i++) {
float_part *= 2.0;
exponent--;
}
/* self == float_part * 2**exponent exactly and float_part is integral.
If FLT_RADIX != 2, the 300 steps may leave a tiny fractional part
to be truncated by TyLong_FromDouble(). */
numerator = TyLong_FromDouble(float_part);
if (numerator == NULL)
goto error;
denominator = TyLong_FromLong(1);
if (denominator == NULL)
goto error;
py_exponent = TyLong_FromLong(Ty_ABS(exponent));
if (py_exponent == NULL)
goto error;
/* fold in 2**exponent */
if (exponent > 0) {
Ty_SETREF(numerator,
long_methods->nb_lshift(numerator, py_exponent));
if (numerator == NULL)
goto error;
}
else {
Ty_SETREF(denominator,
long_methods->nb_lshift(denominator, py_exponent));
if (denominator == NULL)
goto error;
}
result_pair = TyTuple_Pack(2, numerator, denominator);
error:
Ty_XDECREF(py_exponent);
Ty_XDECREF(denominator);
Ty_XDECREF(numerator);
return result_pair;
}
static TyObject *
float_subtype_new(TyTypeObject *type, TyObject *x);
/*[clinic input]
@classmethod
float.__new__ as float_new
x: object(c_default="NULL") = 0
/
Convert a string or number to a floating-point number, if possible.
[clinic start generated code]*/
static TyObject *
float_new_impl(TyTypeObject *type, TyObject *x)
/*[clinic end generated code: output=ccf1e8dc460ba6ba input=55909f888aa0c8a6]*/
{
if (type != &TyFloat_Type) {
if (x == NULL) {
x = _TyLong_GetZero();
}
return float_subtype_new(type, x); /* Wimp out */
}
if (x == NULL) {
return TyFloat_FromDouble(0.0);
}
/* If it's a string, but not a string subclass, use
TyFloat_FromString. */
if (TyUnicode_CheckExact(x))
return TyFloat_FromString(x);
return PyNumber_Float(x);
}
/* Wimpy, slow approach to tp_new calls for subtypes of float:
first create a regular float from whatever arguments we got,
then allocate a subtype instance and initialize its ob_fval
from the regular float. The regular float is then thrown away.
*/
static TyObject *
float_subtype_new(TyTypeObject *type, TyObject *x)
{
TyObject *tmp, *newobj;
assert(TyType_IsSubtype(type, &TyFloat_Type));
tmp = float_new_impl(&TyFloat_Type, x);
if (tmp == NULL)
return NULL;
assert(TyFloat_Check(tmp));
newobj = type->tp_alloc(type, 0);
if (newobj == NULL) {
Ty_DECREF(tmp);
return NULL;
}
((PyFloatObject *)newobj)->ob_fval = ((PyFloatObject *)tmp)->ob_fval;
Ty_DECREF(tmp);
return newobj;
}
static TyObject *
float_vectorcall(TyObject *type, TyObject *const *args,
size_t nargsf, TyObject *kwnames)
{
if (!_TyArg_NoKwnames("float", kwnames)) {
return NULL;
}
Ty_ssize_t nargs = PyVectorcall_NARGS(nargsf);
if (!_TyArg_CheckPositional("float", nargs, 0, 1)) {
return NULL;
}
TyObject *x = nargs >= 1 ? args[0] : NULL;
return float_new_impl(_TyType_CAST(type), x);
}
/*[clinic input]
@classmethod
float.from_number
number: object
/
Convert real number to a floating-point number.
[clinic start generated code]*/
static TyObject *
float_from_number_impl(TyTypeObject *type, TyObject *number)
/*[clinic end generated code: output=dda7e4466ab7068d input=1f8424d9bc11866a]*/
{
if (TyFloat_CheckExact(number) && type == &TyFloat_Type) {
Ty_INCREF(number);
return number;
}
double x = TyFloat_AsDouble(number);
if (x == -1.0 && TyErr_Occurred()) {
return NULL;
}
TyObject *result = TyFloat_FromDouble(x);
if (type != &TyFloat_Type && result != NULL) {
Ty_SETREF(result, PyObject_CallOneArg((TyObject *)type, result));
}
return result;
}
/*[clinic input]
float.__getnewargs__
[clinic start generated code]*/
static TyObject *
float___getnewargs___impl(TyObject *self)
/*[clinic end generated code: output=873258c9d206b088 input=002279d1d77891e6]*/
{
return Ty_BuildValue("(d)", ((PyFloatObject *)self)->ob_fval);
}
/* this is for the benefit of the pack/unpack routines below */
typedef enum _py_float_format_type float_format_type;
#define unknown_format _py_float_format_unknown
#define ieee_big_endian_format _py_float_format_ieee_big_endian
#define ieee_little_endian_format _py_float_format_ieee_little_endian
#define float_format (_PyRuntime.float_state.float_format)
#define double_format (_PyRuntime.float_state.double_format)
/*[clinic input]
@classmethod
float.__getformat__
typestr: str
Must be 'double' or 'float'.
/
You probably don't want to use this function.
It exists mainly to be used in Python's test suite.
This function returns whichever of 'unknown', 'IEEE, big-endian' or 'IEEE,
little-endian' best describes the format of floating-point numbers used by the
C type named by typestr.
[clinic start generated code]*/
static TyObject *
float___getformat___impl(TyTypeObject *type, const char *typestr)
/*[clinic end generated code: output=2bfb987228cc9628 input=90d5e246409a246e]*/
{
float_format_type r;
if (strcmp(typestr, "double") == 0) {
r = double_format;
}
else if (strcmp(typestr, "float") == 0) {
r = float_format;
}
else {
TyErr_SetString(TyExc_ValueError,
"__getformat__() argument 1 must be "
"'double' or 'float'");
return NULL;
}
switch (r) {
case unknown_format:
return TyUnicode_FromString("unknown");
case ieee_little_endian_format:
return TyUnicode_FromString("IEEE, little-endian");
case ieee_big_endian_format:
return TyUnicode_FromString("IEEE, big-endian");
default:
TyErr_SetString(TyExc_RuntimeError,
"insane float_format or double_format");
return NULL;
}
}
static TyObject *
float_getreal(TyObject *v, void *Py_UNUSED(closure))
{
return float_float(v);
}
static TyObject *
float_getimag(TyObject *Py_UNUSED(v), void *Py_UNUSED(closure))
{
return TyFloat_FromDouble(0.0);
}
/*[clinic input]
float.__format__
format_spec: unicode
/
Formats the float according to format_spec.
[clinic start generated code]*/
static TyObject *
float___format___impl(TyObject *self, TyObject *format_spec)
/*[clinic end generated code: output=b260e52a47eade56 input=2ece1052211fd0e6]*/
{
_PyUnicodeWriter writer;
int ret;
_PyUnicodeWriter_Init(&writer);
ret = _TyFloat_FormatAdvancedWriter(
&writer,
self,
format_spec, 0, TyUnicode_GET_LENGTH(format_spec));
if (ret == -1) {
_PyUnicodeWriter_Dealloc(&writer);
return NULL;
}
return _PyUnicodeWriter_Finish(&writer);
}
static TyMethodDef float_methods[] = {
FLOAT_FROM_NUMBER_METHODDEF
FLOAT_CONJUGATE_METHODDEF
FLOAT___TRUNC___METHODDEF
FLOAT___FLOOR___METHODDEF
FLOAT___CEIL___METHODDEF
FLOAT___ROUND___METHODDEF
FLOAT_AS_INTEGER_RATIO_METHODDEF
FLOAT_FROMHEX_METHODDEF
FLOAT_HEX_METHODDEF
FLOAT_IS_INTEGER_METHODDEF
FLOAT___GETNEWARGS___METHODDEF
FLOAT___GETFORMAT___METHODDEF
FLOAT___FORMAT___METHODDEF
{NULL, NULL} /* sentinel */
};
static TyGetSetDef float_getset[] = {
{"real",
float_getreal, NULL,
"the real part of a complex number",
NULL},
{"imag",
float_getimag, NULL,
"the imaginary part of a complex number",
NULL},
{NULL} /* Sentinel */
};
static TyNumberMethods float_as_number = {
float_add, /* nb_add */
float_sub, /* nb_subtract */
float_mul, /* nb_multiply */
float_rem, /* nb_remainder */
float_divmod, /* nb_divmod */
float_pow, /* nb_power */
float_neg, /* nb_negative */
float_float, /* nb_positive */
float_abs, /* nb_absolute */
float_bool, /* nb_bool */
0, /* nb_invert */
0, /* nb_lshift */
0, /* nb_rshift */
0, /* nb_and */
0, /* nb_xor */
0, /* nb_or */
float___trunc___impl, /* nb_int */
0, /* nb_reserved */
float_float, /* nb_float */
0, /* nb_inplace_add */
0, /* nb_inplace_subtract */
0, /* nb_inplace_multiply */
0, /* nb_inplace_remainder */
0, /* nb_inplace_power */
0, /* nb_inplace_lshift */
0, /* nb_inplace_rshift */
0, /* nb_inplace_and */
0, /* nb_inplace_xor */
0, /* nb_inplace_or */
float_floor_div, /* nb_floor_divide */
float_div, /* nb_true_divide */
0, /* nb_inplace_floor_divide */
0, /* nb_inplace_true_divide */
};
TyTypeObject TyFloat_Type = {
TyVarObject_HEAD_INIT(&TyType_Type, 0)
"float",
sizeof(PyFloatObject),
0,
float_dealloc, /* tp_dealloc */
0, /* tp_vectorcall_offset */
0, /* tp_getattr */
0, /* tp_setattr */
0, /* tp_as_async */
float_repr, /* tp_repr */
&float_as_number, /* tp_as_number */
0, /* tp_as_sequence */
0, /* tp_as_mapping */
float_hash, /* tp_hash */
0, /* tp_call */
0, /* tp_str */
PyObject_GenericGetAttr, /* tp_getattro */
0, /* tp_setattro */
0, /* tp_as_buffer */
Ty_TPFLAGS_DEFAULT | Ty_TPFLAGS_BASETYPE |
_Ty_TPFLAGS_MATCH_SELF, /* tp_flags */
float_new__doc__, /* tp_doc */
0, /* tp_traverse */
0, /* tp_clear */
float_richcompare, /* tp_richcompare */
0, /* tp_weaklistoffset */
0, /* tp_iter */
0, /* tp_iternext */
float_methods, /* tp_methods */
0, /* tp_members */
float_getset, /* tp_getset */
0, /* tp_base */
0, /* tp_dict */
0, /* tp_descr_get */
0, /* tp_descr_set */
0, /* tp_dictoffset */
0, /* tp_init */
0, /* tp_alloc */
float_new, /* tp_new */
.tp_vectorcall = float_vectorcall,
.tp_version_tag = _Ty_TYPE_VERSION_FLOAT,
};
static void
_init_global_state(void)
{
float_format_type detected_double_format, detected_float_format;
/* We attempt to determine if this machine is using IEEE
floating-point formats by peering at the bits of some
carefully chosen values. If it looks like we are on an
IEEE platform, the float packing/unpacking routines can
just copy bits, if not they resort to arithmetic & shifts
and masks. The shifts & masks approach works on all finite
values, but what happens to infinities, NaNs and signed
zeroes on packing is an accident, and attempting to unpack
a NaN or an infinity will raise an exception.
Note that if we're on some whacked-out platform which uses
IEEE formats but isn't strictly little-endian or big-
endian, we will fall back to the portable shifts & masks
method. */
#if SIZEOF_DOUBLE == 8
{
double x = 9006104071832581.0;
if (memcmp(&x, "\x43\x3f\xff\x01\x02\x03\x04\x05", 8) == 0)
detected_double_format = ieee_big_endian_format;
else if (memcmp(&x, "\x05\x04\x03\x02\x01\xff\x3f\x43", 8) == 0)
detected_double_format = ieee_little_endian_format;
else
detected_double_format = unknown_format;
}
#else
detected_double_format = unknown_format;
#endif
#if SIZEOF_FLOAT == 4
{
float y = 16711938.0;
if (memcmp(&y, "\x4b\x7f\x01\x02", 4) == 0)
detected_float_format = ieee_big_endian_format;
else if (memcmp(&y, "\x02\x01\x7f\x4b", 4) == 0)
detected_float_format = ieee_little_endian_format;
else
detected_float_format = unknown_format;
}
#else
detected_float_format = unknown_format;
#endif
double_format = detected_double_format;
float_format = detected_float_format;
}
void
_TyFloat_InitState(TyInterpreterState *interp)
{
if (!_Ty_IsMainInterpreter(interp)) {
return;
}
_init_global_state();
}
TyStatus
_TyFloat_InitTypes(TyInterpreterState *interp)
{
/* Init float info */
if (_PyStructSequence_InitBuiltin(interp, &FloatInfoType,
&floatinfo_desc) < 0)
{
return _TyStatus_ERR("can't init float info type");
}
return _TyStatus_OK();
}
void
_TyFloat_FiniType(TyInterpreterState *interp)
{
_PyStructSequence_FiniBuiltin(interp, &FloatInfoType);
}
/* Print summary info about the state of the optimized allocator */
void
_TyFloat_DebugMallocStats(FILE *out)
{
_PyDebugAllocatorStats(out,
"free PyFloatObject",
_Ty_FREELIST_SIZE(floats),
sizeof(PyFloatObject));
}
/*----------------------------------------------------------------------------
* TyFloat_{Pack,Unpack}{2,4,8}. See floatobject.h.
* To match the NPY_HALF_ROUND_TIES_TO_EVEN behavior in:
* https://github.com/numpy/numpy/blob/master/numpy/core/src/npymath/halffloat.c
* We use:
* bits = (unsigned short)f; Note the truncation
* if ((f - bits > 0.5) || (f - bits == 0.5 && bits % 2)) {
* bits++;
* }
*/
int
TyFloat_Pack2(double x, char *data, int le)
{
unsigned char *p = (unsigned char *)data;
unsigned char sign;
int e;
double f;
unsigned short bits;
int incr = 1;
if (x == 0.0) {
sign = (copysign(1.0, x) == -1.0);
e = 0;
bits = 0;
}
else if (isinf(x)) {
sign = (x < 0.0);
e = 0x1f;
bits = 0;
}
else if (isnan(x)) {
sign = (copysign(1.0, x) == -1.0);
e = 0x1f;
uint64_t v;
memcpy(&v, &x, sizeof(v));
v &= 0xffc0000000000ULL;
bits = (unsigned short)(v >> 42); /* NaN's type & payload */
}
else {
sign = (x < 0.0);
if (sign) {
x = -x;
}
f = frexp(x, &e);
if (f < 0.5 || f >= 1.0) {
TyErr_SetString(TyExc_SystemError,
"frexp() result out of range");
return -1;
}
/* Normalize f to be in the range [1.0, 2.0) */
f *= 2.0;
e--;
if (e >= 16) {
goto Overflow;
}
else if (e < -25) {
/* |x| < 2**-25. Underflow to zero. */
f = 0.0;
e = 0;
}
else if (e < -14) {
/* |x| < 2**-14. Gradual underflow */
f = ldexp(f, 14 + e);
e = 0;
}
else /* if (!(e == 0 && f == 0.0)) */ {
e += 15;
f -= 1.0; /* Get rid of leading 1 */
}
f *= 1024.0; /* 2**10 */
/* Round to even */
bits = (unsigned short)f; /* Note the truncation */
assert(bits < 1024);
assert(e < 31);
if ((f - bits > 0.5) || ((f - bits == 0.5) && (bits % 2 == 1))) {
++bits;
if (bits == 1024) {
/* The carry propagated out of a string of 10 1 bits. */
bits = 0;
++e;
if (e == 31)
goto Overflow;
}
}
}
bits |= (e << 10) | (sign << 15);
/* Write out result. */
if (le) {
p += 1;
incr = -1;
}
/* First byte */
*p = (unsigned char)((bits >> 8) & 0xFF);
p += incr;
/* Second byte */
*p = (unsigned char)(bits & 0xFF);
return 0;
Overflow:
TyErr_SetString(TyExc_OverflowError,
"float too large to pack with e format");
return -1;
}
int
TyFloat_Pack4(double x, char *data, int le)
{
unsigned char *p = (unsigned char *)data;
if (float_format == unknown_format) {
unsigned char sign;
int e;
double f;
unsigned int fbits;
int incr = 1;
if (le) {
p += 3;
incr = -1;
}
if (x < 0) {
sign = 1;
x = -x;
}
else
sign = 0;
f = frexp(x, &e);
/* Normalize f to be in the range [1.0, 2.0) */
if (0.5 <= f && f < 1.0) {
f *= 2.0;
e--;
}
else if (f == 0.0)
e = 0;
else {
TyErr_SetString(TyExc_SystemError,
"frexp() result out of range");
return -1;
}
if (e >= 128)
goto Overflow;
else if (e < -126) {
/* Gradual underflow */
f = ldexp(f, 126 + e);
e = 0;
}
else if (!(e == 0 && f == 0.0)) {
e += 127;
f -= 1.0; /* Get rid of leading 1 */
}
f *= 8388608.0; /* 2**23 */
fbits = (unsigned int)(f + 0.5); /* Round */
assert(fbits <= 8388608);
if (fbits >> 23) {
/* The carry propagated out of a string of 23 1 bits. */
fbits = 0;
++e;
if (e >= 255)
goto Overflow;
}
/* First byte */
*p = (sign << 7) | (e >> 1);
p += incr;
/* Second byte */
*p = (char) (((e & 1) << 7) | (fbits >> 16));
p += incr;
/* Third byte */
*p = (fbits >> 8) & 0xFF;
p += incr;
/* Fourth byte */
*p = fbits & 0xFF;
/* Done */
return 0;
}
else {
float y = (float)x;
int i, incr = 1;
if (isinf(y) && !isinf(x))
goto Overflow;
/* correct y if x was a sNaN, transformed to qNaN by conversion */
if (isnan(x)) {
uint64_t v;
memcpy(&v, &x, 8);
#ifndef __riscv
if ((v & (1ULL << 51)) == 0) {
uint32_t u32;
memcpy(&u32, &y, 4);
u32 &= ~(1 << 22); /* make sNaN */
memcpy(&y, &u32, 4);
}
#else
uint32_t u32;
memcpy(&u32, &y, 4);
if ((v & (1ULL << 51)) == 0) {
u32 &= ~(1 << 22);
}
/* Workaround RISC-V: "If a NaN value is converted to a
* different floating-point type, the result is the
* canonical NaN of the new type". The canonical NaN here
* is a positive qNaN with zero payload. */
if (v & (1ULL << 63)) {
u32 |= (1 << 31); /* set sign */
}
/* add payload */
u32 -= (u32 & 0x3fffff);
u32 += (uint32_t)((v & 0x7ffffffffffffULL) >> 29);
memcpy(&y, &u32, 4);
#endif
}
unsigned char s[sizeof(float)];
memcpy(s, &y, sizeof(float));
if ((float_format == ieee_little_endian_format && !le)
|| (float_format == ieee_big_endian_format && le)) {
p += 3;
incr = -1;
}
for (i = 0; i < 4; i++) {
*p = s[i];
p += incr;
}
return 0;
}
Overflow:
TyErr_SetString(TyExc_OverflowError,
"float too large to pack with f format");
return -1;
}
int
TyFloat_Pack8(double x, char *data, int le)
{
unsigned char *p = (unsigned char *)data;
if (double_format == unknown_format) {
unsigned char sign;
int e;
double f;
unsigned int fhi, flo;
int incr = 1;
if (le) {
p += 7;
incr = -1;
}
if (x < 0) {
sign = 1;
x = -x;
}
else
sign = 0;
f = frexp(x, &e);
/* Normalize f to be in the range [1.0, 2.0) */
if (0.5 <= f && f < 1.0) {
f *= 2.0;
e--;
}
else if (f == 0.0)
e = 0;
else {
TyErr_SetString(TyExc_SystemError,
"frexp() result out of range");
return -1;
}
if (e >= 1024)
goto Overflow;
else if (e < -1022) {
/* Gradual underflow */
f = ldexp(f, 1022 + e);
e = 0;
}
else if (!(e == 0 && f == 0.0)) {
e += 1023;
f -= 1.0; /* Get rid of leading 1 */
}
/* fhi receives the high 28 bits; flo the low 24 bits (== 52 bits) */
f *= 268435456.0; /* 2**28 */
fhi = (unsigned int)f; /* Truncate */
assert(fhi < 268435456);
f -= (double)fhi;
f *= 16777216.0; /* 2**24 */
flo = (unsigned int)(f + 0.5); /* Round */
assert(flo <= 16777216);
if (flo >> 24) {
/* The carry propagated out of a string of 24 1 bits. */
flo = 0;
++fhi;
if (fhi >> 28) {
/* And it also propagated out of the next 28 bits. */
fhi = 0;
++e;
if (e >= 2047)
goto Overflow;
}
}
/* First byte */
*p = (sign << 7) | (e >> 4);
p += incr;
/* Second byte */
*p = (unsigned char) (((e & 0xF) << 4) | (fhi >> 24));
p += incr;
/* Third byte */
*p = (fhi >> 16) & 0xFF;
p += incr;
/* Fourth byte */
*p = (fhi >> 8) & 0xFF;
p += incr;
/* Fifth byte */
*p = fhi & 0xFF;
p += incr;
/* Sixth byte */
*p = (flo >> 16) & 0xFF;
p += incr;
/* Seventh byte */
*p = (flo >> 8) & 0xFF;
p += incr;
/* Eighth byte */
*p = flo & 0xFF;
/* p += incr; */
/* Done */
return 0;
Overflow:
TyErr_SetString(TyExc_OverflowError,
"float too large to pack with d format");
return -1;
}
else {
unsigned char as_bytes[8];
memcpy(as_bytes, &x, 8);
const unsigned char *s = as_bytes;
int i, incr = 1;
if ((double_format == ieee_little_endian_format && !le)
|| (double_format == ieee_big_endian_format && le)) {
p += 7;
incr = -1;
}
for (i = 0; i < 8; i++) {
*p = *s++;
p += incr;
}
return 0;
}
}
double
TyFloat_Unpack2(const char *data, int le)
{
unsigned char *p = (unsigned char *)data;
unsigned char sign;
int e;
unsigned int f;
double x;
int incr = 1;
if (le) {
p += 1;
incr = -1;
}
/* First byte */
sign = (*p >> 7) & 1;
e = (*p & 0x7C) >> 2;
f = (*p & 0x03) << 8;
p += incr;
/* Second byte */
f |= *p;
if (e == 0x1f) {
if (f == 0) {
/* Infinity */
return sign ? -Ty_INFINITY : Ty_INFINITY;
}
else {
/* NaN */
uint64_t v = sign ? 0xfff0000000000000ULL : 0x7ff0000000000000ULL;
v += (uint64_t)f << 42; /* add NaN's type & payload */
memcpy(&x, &v, sizeof(v));
return x;
}
}
x = (double)f / 1024.0;
if (e == 0) {
e = -14;
}
else {
x += 1.0;
e -= 15;
}
x = ldexp(x, e);
if (sign)
x = -x;
return x;
}
double
TyFloat_Unpack4(const char *data, int le)
{
unsigned char *p = (unsigned char *)data;
if (float_format == unknown_format) {
unsigned char sign;
int e;
unsigned int f;
double x;
int incr = 1;
if (le) {
p += 3;
incr = -1;
}
/* First byte */
sign = (*p >> 7) & 1;
e = (*p & 0x7F) << 1;
p += incr;
/* Second byte */
e |= (*p >> 7) & 1;
f = (*p & 0x7F) << 16;
p += incr;
if (e == 255) {
TyErr_SetString(
TyExc_ValueError,
"can't unpack IEEE 754 special value "
"on non-IEEE platform");
return -1;
}
/* Third byte */
f |= *p << 8;
p += incr;
/* Fourth byte */
f |= *p;
x = (double)f / 8388608.0;
/* XXX This sadly ignores Inf/NaN issues */
if (e == 0)
e = -126;
else {
x += 1.0;
e -= 127;
}
x = ldexp(x, e);
if (sign)
x = -x;
return x;
}
else {
float x;
if ((float_format == ieee_little_endian_format && !le)
|| (float_format == ieee_big_endian_format && le)) {
char buf[4];
char *d = &buf[3];
int i;
for (i = 0; i < 4; i++) {
*d-- = *p++;
}
memcpy(&x, buf, 4);
}
else {
memcpy(&x, p, 4);
}
/* return sNaN double if x was sNaN float */
if (isnan(x)) {
uint32_t v;
memcpy(&v, &x, 4);
#ifndef __riscv
if ((v & (1 << 22)) == 0) {
double y = x; /* will make qNaN double */
uint64_t u64;
memcpy(&u64, &y, 8);
u64 &= ~(1ULL << 51); /* make sNaN */
memcpy(&y, &u64, 8);
return y;
}
#else
double y = x;
uint64_t u64;
memcpy(&u64, &y, 8);
if ((v & (1 << 22)) == 0) {
u64 &= ~(1ULL << 51);
}
/* Workaround RISC-V, see TyFloat_Pack4() */
if (v & (1 << 31)) {
u64 |= (1ULL << 63); /* set sign */
}
/* add payload */
u64 -= (u64 & 0x7ffffffffffffULL);
u64 += ((v & 0x3fffffULL) << 29);
memcpy(&y, &u64, 8);
return y;
#endif
}
return x;
}
}
double
TyFloat_Unpack8(const char *data, int le)
{
unsigned char *p = (unsigned char *)data;
if (double_format == unknown_format) {
unsigned char sign;
int e;
unsigned int fhi, flo;
double x;
int incr = 1;
if (le) {
p += 7;
incr = -1;
}
/* First byte */
sign = (*p >> 7) & 1;
e = (*p & 0x7F) << 4;
p += incr;
/* Second byte */
e |= (*p >> 4) & 0xF;
fhi = (*p & 0xF) << 24;
p += incr;
if (e == 2047) {
TyErr_SetString(
TyExc_ValueError,
"can't unpack IEEE 754 special value "
"on non-IEEE platform");
return -1.0;
}
/* Third byte */
fhi |= *p << 16;
p += incr;
/* Fourth byte */
fhi |= *p << 8;
p += incr;
/* Fifth byte */
fhi |= *p;
p += incr;
/* Sixth byte */
flo = *p << 16;
p += incr;
/* Seventh byte */
flo |= *p << 8;
p += incr;
/* Eighth byte */
flo |= *p;
x = (double)fhi + (double)flo / 16777216.0; /* 2**24 */
x /= 268435456.0; /* 2**28 */
if (e == 0)
e = -1022;
else {
x += 1.0;
e -= 1023;
}
x = ldexp(x, e);
if (sign)
x = -x;
return x;
}
else {
double x;
if ((double_format == ieee_little_endian_format && !le)
|| (double_format == ieee_big_endian_format && le)) {
char buf[8];
char *d = &buf[7];
int i;
for (i = 0; i < 8; i++) {
*d-- = *p++;
}
memcpy(&x, buf, 8);
}
else {
memcpy(&x, p, 8);
}
return x;
}
}
Reference in New Issue View Git Blame Copy Permalink