2014-05-03 22:27:38 +00:00
|
|
|
/*
|
2017-06-30 07:22:17 +00:00
|
|
|
* This file is part of the MicroPython project, http://micropython.org/
|
2014-05-03 22:27:38 +00:00
|
|
|
*
|
|
|
|
* The MIT License (MIT)
|
|
|
|
*
|
|
|
|
* Copyright (c) 2013, 2014 Damien P. George
|
|
|
|
*
|
|
|
|
* Permission is hereby granted, free of charge, to any person obtaining a copy
|
|
|
|
* of this software and associated documentation files (the "Software"), to deal
|
|
|
|
* in the Software without restriction, including without limitation the rights
|
|
|
|
* to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
|
|
|
|
* copies of the Software, and to permit persons to whom the Software is
|
|
|
|
* furnished to do so, subject to the following conditions:
|
|
|
|
*
|
|
|
|
* The above copyright notice and this permission notice shall be included in
|
|
|
|
* all copies or substantial portions of the Software.
|
|
|
|
*
|
|
|
|
* THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
|
|
|
|
* IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
|
|
|
|
* FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
|
|
|
|
* AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
|
|
|
|
* LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
|
|
|
|
* OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN
|
|
|
|
* THE SOFTWARE.
|
|
|
|
*/
|
|
|
|
|
2013-12-21 18:17:45 +00:00
|
|
|
#include <stdlib.h>
|
2014-03-30 23:12:40 +00:00
|
|
|
#include <stdio.h>
|
|
|
|
#include <string.h>
|
2013-12-21 18:17:45 +00:00
|
|
|
#include <assert.h>
|
|
|
|
|
2015-01-01 20:27:54 +00:00
|
|
|
#include "py/parsenum.h"
|
|
|
|
#include "py/runtime.h"
|
2013-12-21 18:17:45 +00:00
|
|
|
|
2014-06-01 12:32:54 +00:00
|
|
|
#if MICROPY_PY_BUILTINS_FLOAT
|
2014-03-10 13:27:02 +00:00
|
|
|
|
2015-03-14 23:11:25 +00:00
|
|
|
#include <math.h>
|
2015-01-01 20:27:54 +00:00
|
|
|
#include "py/formatfloat.h"
|
2013-12-21 18:17:45 +00:00
|
|
|
|
2015-11-27 17:09:11 +00:00
|
|
|
#if MICROPY_OBJ_REPR != MICROPY_OBJ_REPR_C && MICROPY_OBJ_REPR != MICROPY_OBJ_REPR_D
|
2015-10-17 22:20:57 +00:00
|
|
|
|
2016-05-09 21:11:34 +00:00
|
|
|
// M_E and M_PI are not part of the math.h standard and may not be defined
|
|
|
|
#ifndef M_E
|
|
|
|
#define M_E (2.7182818284590452354)
|
|
|
|
#endif
|
|
|
|
#ifndef M_PI
|
|
|
|
#define M_PI (3.14159265358979323846)
|
|
|
|
#endif
|
|
|
|
|
2015-08-20 22:42:35 +00:00
|
|
|
typedef struct _mp_obj_float_t {
|
|
|
|
mp_obj_base_t base;
|
|
|
|
mp_float_t value;
|
|
|
|
} mp_obj_float_t;
|
|
|
|
|
|
|
|
const mp_obj_float_t mp_const_float_e_obj = {{&mp_type_float}, M_E};
|
|
|
|
const mp_obj_float_t mp_const_float_pi_obj = {{&mp_type_float}, M_PI};
|
|
|
|
|
2015-10-17 22:20:57 +00:00
|
|
|
#endif
|
|
|
|
|
2017-04-06 07:27:33 +00:00
|
|
|
#if MICROPY_FLOAT_HIGH_QUALITY_HASH
|
|
|
|
// must return actual integer value if it fits in mp_int_t
|
|
|
|
mp_int_t mp_float_hash(mp_float_t src) {
|
|
|
|
#if MICROPY_FLOAT_IMPL == MICROPY_FLOAT_IMPL_DOUBLE
|
|
|
|
typedef uint64_t mp_float_uint_t;
|
|
|
|
#elif MICROPY_FLOAT_IMPL == MICROPY_FLOAT_IMPL_FLOAT
|
|
|
|
typedef uint32_t mp_float_uint_t;
|
|
|
|
#endif
|
|
|
|
union {
|
|
|
|
mp_float_t f;
|
|
|
|
#if MP_ENDIANNESS_LITTLE
|
|
|
|
struct { mp_float_uint_t frc:MP_FLOAT_FRAC_BITS, exp:MP_FLOAT_EXP_BITS, sgn:1; } p;
|
|
|
|
#else
|
|
|
|
struct { mp_float_uint_t sgn:1, exp:MP_FLOAT_EXP_BITS, frc:MP_FLOAT_FRAC_BITS; } p;
|
|
|
|
#endif
|
|
|
|
mp_float_uint_t i;
|
|
|
|
} u = {.f = src};
|
|
|
|
|
|
|
|
mp_int_t val;
|
|
|
|
const int adj_exp = (int)u.p.exp - MP_FLOAT_EXP_BIAS;
|
|
|
|
if (adj_exp < 0) {
|
|
|
|
// value < 1; must be sure to handle 0.0 correctly (ie return 0)
|
|
|
|
val = u.i;
|
|
|
|
} else {
|
|
|
|
// if adj_exp is max then: u.p.frc==0 indicates inf, else NaN
|
|
|
|
// else: 1 <= value
|
|
|
|
mp_float_uint_t frc = u.p.frc | ((mp_float_uint_t)1 << MP_FLOAT_FRAC_BITS);
|
|
|
|
|
|
|
|
if (adj_exp <= MP_FLOAT_FRAC_BITS) {
|
|
|
|
// number may have a fraction; xor the integer part with the fractional part
|
|
|
|
val = (frc >> (MP_FLOAT_FRAC_BITS - adj_exp))
|
|
|
|
^ (frc & ((1 << (MP_FLOAT_FRAC_BITS - adj_exp)) - 1));
|
|
|
|
} else if ((unsigned int)adj_exp < BITS_PER_BYTE * sizeof(mp_int_t) - 1) {
|
|
|
|
// the number is a (big) whole integer and will fit in val's signed-width
|
|
|
|
val = (mp_int_t)frc << (adj_exp - MP_FLOAT_FRAC_BITS);
|
|
|
|
} else {
|
|
|
|
// integer part will overflow val's width so just use what bits we can
|
|
|
|
val = frc;
|
|
|
|
}
|
|
|
|
}
|
|
|
|
|
|
|
|
if (u.p.sgn) {
|
|
|
|
val = -val;
|
|
|
|
}
|
|
|
|
|
|
|
|
return val;
|
|
|
|
}
|
|
|
|
#endif
|
|
|
|
|
2015-04-09 22:56:15 +00:00
|
|
|
STATIC void float_print(const mp_print_t *print, mp_obj_t o_in, mp_print_kind_t kind) {
|
2015-01-20 12:47:20 +00:00
|
|
|
(void)kind;
|
2015-10-17 21:57:34 +00:00
|
|
|
mp_float_t o_val = mp_obj_float_get(o_in);
|
2014-03-10 07:10:01 +00:00
|
|
|
#if MICROPY_FLOAT_IMPL == MICROPY_FLOAT_IMPL_FLOAT
|
2014-09-11 21:24:45 +00:00
|
|
|
char buf[16];
|
2017-04-19 08:46:53 +00:00
|
|
|
#if MICROPY_OBJ_REPR == MICROPY_OBJ_REPR_C
|
|
|
|
const int precision = 6;
|
|
|
|
#else
|
2015-05-16 08:54:19 +00:00
|
|
|
const int precision = 7;
|
2017-04-19 08:46:53 +00:00
|
|
|
#endif
|
2014-03-10 07:10:01 +00:00
|
|
|
#else
|
2014-03-30 23:12:40 +00:00
|
|
|
char buf[32];
|
2015-05-16 08:54:19 +00:00
|
|
|
const int precision = 16;
|
|
|
|
#endif
|
2015-10-17 21:57:34 +00:00
|
|
|
mp_format_float(o_val, buf, sizeof(buf), 'g', precision, '\0');
|
2015-04-09 22:56:15 +00:00
|
|
|
mp_print_str(print, buf);
|
2015-04-22 15:51:29 +00:00
|
|
|
if (strchr(buf, '.') == NULL && strchr(buf, 'e') == NULL && strchr(buf, 'n') == NULL) {
|
|
|
|
// Python floats always have decimal point (unless inf or nan)
|
2015-04-09 22:56:15 +00:00
|
|
|
mp_print_str(print, ".0");
|
2014-03-30 23:12:40 +00:00
|
|
|
}
|
2013-12-21 18:17:45 +00:00
|
|
|
}
|
|
|
|
|
2016-01-03 15:55:55 +00:00
|
|
|
STATIC mp_obj_t float_make_new(const mp_obj_type_t *type_in, size_t n_args, size_t n_kw, const mp_obj_t *args) {
|
2015-01-20 12:47:20 +00:00
|
|
|
(void)type_in;
|
2014-05-11 17:37:21 +00:00
|
|
|
mp_arg_check_num(n_args, n_kw, 0, 1, false);
|
2014-01-18 14:10:48 +00:00
|
|
|
|
2014-01-04 20:21:15 +00:00
|
|
|
switch (n_args) {
|
|
|
|
case 0:
|
|
|
|
return mp_obj_new_float(0);
|
|
|
|
|
|
|
|
case 1:
|
2017-11-21 04:01:38 +00:00
|
|
|
default: {
|
|
|
|
mp_buffer_info_t bufinfo;
|
|
|
|
if (mp_get_buffer(args[0], &bufinfo, MP_BUFFER_READ)) {
|
|
|
|
// a textual representation, parse it
|
|
|
|
return mp_parse_num_decimal(bufinfo.buf, bufinfo.len, false, false, NULL);
|
2015-08-20 22:30:12 +00:00
|
|
|
} else if (mp_obj_is_float(args[0])) {
|
2014-03-21 11:45:46 +00:00
|
|
|
// a float, just return it
|
2014-01-04 20:21:15 +00:00
|
|
|
return args[0];
|
|
|
|
} else {
|
2014-03-21 11:45:46 +00:00
|
|
|
// something else, try to cast it to a float
|
2014-01-04 20:21:15 +00:00
|
|
|
return mp_obj_new_float(mp_obj_get_float(args[0]));
|
|
|
|
}
|
2017-11-21 04:01:38 +00:00
|
|
|
}
|
2014-01-04 20:21:15 +00:00
|
|
|
}
|
|
|
|
}
|
|
|
|
|
2017-08-29 03:04:01 +00:00
|
|
|
STATIC mp_obj_t float_unary_op(mp_unary_op_t op, mp_obj_t o_in) {
|
2015-10-17 21:57:34 +00:00
|
|
|
mp_float_t val = mp_obj_float_get(o_in);
|
2013-12-21 18:17:45 +00:00
|
|
|
switch (op) {
|
2015-10-17 21:57:34 +00:00
|
|
|
case MP_UNARY_OP_BOOL: return mp_obj_new_bool(val != 0);
|
2017-04-04 01:57:21 +00:00
|
|
|
case MP_UNARY_OP_HASH: return MP_OBJ_NEW_SMALL_INT(mp_float_hash(val));
|
2014-03-30 12:35:08 +00:00
|
|
|
case MP_UNARY_OP_POSITIVE: return o_in;
|
2015-10-17 21:57:34 +00:00
|
|
|
case MP_UNARY_OP_NEGATIVE: return mp_obj_new_float(-val);
|
2017-09-17 21:06:43 +00:00
|
|
|
case MP_UNARY_OP_ABS: {
|
|
|
|
// TODO check for NaN etc
|
2017-09-18 04:31:03 +00:00
|
|
|
if (val < 0) {
|
|
|
|
return mp_obj_new_float(-val);
|
2017-09-17 21:06:43 +00:00
|
|
|
} else {
|
|
|
|
return o_in;
|
|
|
|
}
|
|
|
|
}
|
2014-05-21 18:42:43 +00:00
|
|
|
default: return MP_OBJ_NULL; // op not supported
|
2013-12-21 18:17:45 +00:00
|
|
|
}
|
|
|
|
}
|
|
|
|
|
2017-08-29 03:04:01 +00:00
|
|
|
STATIC mp_obj_t float_binary_op(mp_binary_op_t op, mp_obj_t lhs_in, mp_obj_t rhs_in) {
|
2015-10-17 21:57:34 +00:00
|
|
|
mp_float_t lhs_val = mp_obj_float_get(lhs_in);
|
2014-06-19 22:48:35 +00:00
|
|
|
#if MICROPY_PY_BUILTINS_COMPLEX
|
2014-03-08 15:24:39 +00:00
|
|
|
if (MP_OBJ_IS_TYPE(rhs_in, &mp_type_complex)) {
|
2015-10-17 21:57:34 +00:00
|
|
|
return mp_obj_complex_binary_op(op, lhs_val, 0, rhs_in);
|
2014-06-19 22:48:35 +00:00
|
|
|
} else
|
|
|
|
#endif
|
|
|
|
{
|
2015-10-17 21:57:34 +00:00
|
|
|
return mp_obj_float_binary_op(op, lhs_val, rhs_in);
|
2013-12-21 18:17:45 +00:00
|
|
|
}
|
|
|
|
}
|
|
|
|
|
2014-03-08 15:24:39 +00:00
|
|
|
const mp_obj_type_t mp_type_float = {
|
2014-02-15 16:10:44 +00:00
|
|
|
{ &mp_type_type },
|
2014-02-15 11:34:50 +00:00
|
|
|
.name = MP_QSTR_float,
|
2014-01-05 20:34:09 +00:00
|
|
|
.print = float_print,
|
|
|
|
.make_new = float_make_new,
|
|
|
|
.unary_op = float_unary_op,
|
|
|
|
.binary_op = float_binary_op,
|
2013-12-21 18:17:45 +00:00
|
|
|
};
|
|
|
|
|
2015-11-27 17:09:11 +00:00
|
|
|
#if MICROPY_OBJ_REPR != MICROPY_OBJ_REPR_C && MICROPY_OBJ_REPR != MICROPY_OBJ_REPR_D
|
2015-10-17 22:20:57 +00:00
|
|
|
|
2013-12-21 18:17:45 +00:00
|
|
|
mp_obj_t mp_obj_new_float(mp_float_t value) {
|
|
|
|
mp_obj_float_t *o = m_new(mp_obj_float_t, 1);
|
2014-03-08 15:24:39 +00:00
|
|
|
o->base.type = &mp_type_float;
|
2013-12-21 18:17:45 +00:00
|
|
|
o->value = value;
|
2015-11-27 17:01:44 +00:00
|
|
|
return MP_OBJ_FROM_PTR(o);
|
2013-12-21 18:17:45 +00:00
|
|
|
}
|
|
|
|
|
|
|
|
mp_float_t mp_obj_float_get(mp_obj_t self_in) {
|
2015-08-20 22:30:12 +00:00
|
|
|
assert(mp_obj_is_float(self_in));
|
2015-11-27 17:01:44 +00:00
|
|
|
mp_obj_float_t *self = MP_OBJ_TO_PTR(self_in);
|
2013-12-21 18:17:45 +00:00
|
|
|
return self->value;
|
|
|
|
}
|
|
|
|
|
2015-10-17 22:20:57 +00:00
|
|
|
#endif
|
|
|
|
|
2015-06-13 21:00:10 +00:00
|
|
|
STATIC void mp_obj_float_divmod(mp_float_t *x, mp_float_t *y) {
|
|
|
|
// logic here follows that of CPython
|
|
|
|
// https://docs.python.org/3/reference/expressions.html#binary-arithmetic-operations
|
|
|
|
// x == (x//y)*y + (x%y)
|
|
|
|
// divmod(x, y) == (x//y, x%y)
|
|
|
|
mp_float_t mod = MICROPY_FLOAT_C_FUN(fmod)(*x, *y);
|
|
|
|
mp_float_t div = (*x - mod) / *y;
|
|
|
|
|
|
|
|
// Python specs require that mod has same sign as second operand
|
|
|
|
if (mod == 0.0) {
|
|
|
|
mod = MICROPY_FLOAT_C_FUN(copysign)(0.0, *y);
|
|
|
|
} else {
|
|
|
|
if ((mod < 0.0) != (*y < 0.0)) {
|
|
|
|
mod += *y;
|
|
|
|
div -= 1.0;
|
|
|
|
}
|
|
|
|
}
|
|
|
|
|
|
|
|
mp_float_t floordiv;
|
|
|
|
if (div == 0.0) {
|
|
|
|
// if division is zero, take the correct sign of zero
|
|
|
|
floordiv = MICROPY_FLOAT_C_FUN(copysign)(0.0, *x / *y);
|
|
|
|
} else {
|
|
|
|
// Python specs require that x == (x//y)*y + (x%y)
|
|
|
|
floordiv = MICROPY_FLOAT_C_FUN(floor)(div);
|
|
|
|
if (div - floordiv > 0.5) {
|
|
|
|
floordiv += 1.0;
|
|
|
|
}
|
|
|
|
}
|
|
|
|
|
|
|
|
// return results
|
|
|
|
*x = floordiv;
|
|
|
|
*y = mod;
|
|
|
|
}
|
|
|
|
|
2017-08-29 03:04:01 +00:00
|
|
|
mp_obj_t mp_obj_float_binary_op(mp_binary_op_t op, mp_float_t lhs_val, mp_obj_t rhs_in) {
|
2017-09-02 18:19:01 +00:00
|
|
|
mp_float_t rhs_val;
|
|
|
|
if (!mp_obj_get_float_maybe(rhs_in, &rhs_val)) {
|
|
|
|
return MP_OBJ_NULL; // op not supported
|
|
|
|
}
|
|
|
|
|
2014-01-06 22:13:00 +00:00
|
|
|
switch (op) {
|
2014-03-30 12:35:08 +00:00
|
|
|
case MP_BINARY_OP_ADD:
|
|
|
|
case MP_BINARY_OP_INPLACE_ADD: lhs_val += rhs_val; break;
|
|
|
|
case MP_BINARY_OP_SUBTRACT:
|
|
|
|
case MP_BINARY_OP_INPLACE_SUBTRACT: lhs_val -= rhs_val; break;
|
|
|
|
case MP_BINARY_OP_MULTIPLY:
|
|
|
|
case MP_BINARY_OP_INPLACE_MULTIPLY: lhs_val *= rhs_val; break;
|
|
|
|
case MP_BINARY_OP_FLOOR_DIVIDE:
|
2014-03-30 23:18:09 +00:00
|
|
|
case MP_BINARY_OP_INPLACE_FLOOR_DIVIDE:
|
2014-04-13 16:45:30 +00:00
|
|
|
if (rhs_val == 0) {
|
|
|
|
zero_division_error:
|
2016-10-17 01:17:37 +00:00
|
|
|
mp_raise_msg(&mp_type_ZeroDivisionError, "division by zero");
|
2014-04-13 16:45:30 +00:00
|
|
|
}
|
2014-09-13 17:43:09 +00:00
|
|
|
// Python specs require that x == (x//y)*y + (x%y) so we must
|
|
|
|
// call divmod to compute the correct floor division, which
|
|
|
|
// returns the floor divide in lhs_val.
|
|
|
|
mp_obj_float_divmod(&lhs_val, &rhs_val);
|
2014-04-13 16:45:30 +00:00
|
|
|
break;
|
2014-03-30 12:35:08 +00:00
|
|
|
case MP_BINARY_OP_TRUE_DIVIDE:
|
2014-04-13 16:45:30 +00:00
|
|
|
case MP_BINARY_OP_INPLACE_TRUE_DIVIDE:
|
|
|
|
if (rhs_val == 0) {
|
|
|
|
goto zero_division_error;
|
2014-03-20 22:40:38 +00:00
|
|
|
}
|
2014-04-13 16:45:30 +00:00
|
|
|
lhs_val /= rhs_val;
|
2014-03-20 22:40:38 +00:00
|
|
|
break;
|
2014-09-13 17:43:09 +00:00
|
|
|
case MP_BINARY_OP_MODULO:
|
|
|
|
case MP_BINARY_OP_INPLACE_MODULO:
|
|
|
|
if (rhs_val == 0) {
|
|
|
|
goto zero_division_error;
|
|
|
|
}
|
|
|
|
lhs_val = MICROPY_FLOAT_C_FUN(fmod)(lhs_val, rhs_val);
|
|
|
|
// Python specs require that mod has same sign as second operand
|
|
|
|
if (lhs_val == 0.0) {
|
|
|
|
lhs_val = MICROPY_FLOAT_C_FUN(copysign)(0.0, rhs_val);
|
|
|
|
} else {
|
|
|
|
if ((lhs_val < 0.0) != (rhs_val < 0.0)) {
|
|
|
|
lhs_val += rhs_val;
|
|
|
|
}
|
|
|
|
}
|
|
|
|
break;
|
2014-04-02 11:26:49 +00:00
|
|
|
case MP_BINARY_OP_POWER:
|
2017-02-02 13:01:37 +00:00
|
|
|
case MP_BINARY_OP_INPLACE_POWER:
|
|
|
|
if (lhs_val == 0 && rhs_val < 0) {
|
|
|
|
goto zero_division_error;
|
|
|
|
}
|
2017-09-26 02:57:51 +00:00
|
|
|
if (lhs_val < 0 && rhs_val != MICROPY_FLOAT_C_FUN(floor)(rhs_val)) {
|
|
|
|
#if MICROPY_PY_BUILTINS_COMPLEX
|
|
|
|
return mp_obj_complex_binary_op(MP_BINARY_OP_POWER, lhs_val, 0, rhs_in);
|
|
|
|
#else
|
|
|
|
mp_raise_ValueError("complex values not supported");
|
|
|
|
#endif
|
|
|
|
}
|
2017-02-02 13:01:37 +00:00
|
|
|
lhs_val = MICROPY_FLOAT_C_FUN(pow)(lhs_val, rhs_val);
|
|
|
|
break;
|
2015-06-13 21:00:10 +00:00
|
|
|
case MP_BINARY_OP_DIVMOD: {
|
|
|
|
if (rhs_val == 0) {
|
|
|
|
goto zero_division_error;
|
|
|
|
}
|
|
|
|
mp_obj_float_divmod(&lhs_val, &rhs_val);
|
|
|
|
mp_obj_t tuple[2] = {
|
|
|
|
mp_obj_new_float(lhs_val),
|
|
|
|
mp_obj_new_float(rhs_val),
|
|
|
|
};
|
|
|
|
return mp_obj_new_tuple(2, tuple);
|
|
|
|
}
|
2015-10-11 09:09:43 +00:00
|
|
|
case MP_BINARY_OP_LESS: return mp_obj_new_bool(lhs_val < rhs_val);
|
|
|
|
case MP_BINARY_OP_MORE: return mp_obj_new_bool(lhs_val > rhs_val);
|
|
|
|
case MP_BINARY_OP_EQUAL: return mp_obj_new_bool(lhs_val == rhs_val);
|
|
|
|
case MP_BINARY_OP_LESS_EQUAL: return mp_obj_new_bool(lhs_val <= rhs_val);
|
|
|
|
case MP_BINARY_OP_MORE_EQUAL: return mp_obj_new_bool(lhs_val >= rhs_val);
|
2014-01-11 00:58:59 +00:00
|
|
|
|
2014-03-30 22:38:25 +00:00
|
|
|
default:
|
2014-05-21 18:42:43 +00:00
|
|
|
return MP_OBJ_NULL; // op not supported
|
2014-01-06 22:13:00 +00:00
|
|
|
}
|
|
|
|
return mp_obj_new_float(lhs_val);
|
|
|
|
}
|
|
|
|
|
2014-06-01 12:32:54 +00:00
|
|
|
#endif // MICROPY_PY_BUILTINS_FLOAT
|