objfloat.c 4.15 KB
Newer Older
1
2
#include <stdlib.h>
#include <assert.h>
3
#include <math.h>
4
5
6
7

#include "nlr.h"
#include "misc.h"
#include "mpconfig.h"
8
#include "qstr.h"
9
#include "obj.h"
10
#include "parsenum.h"
11
12
13
#include "runtime0.h"

#if MICROPY_ENABLE_FLOAT
14

15
#if MICROPY_FLOAT_IMPL == MICROPY_FLOAT_IMPL_FLOAT
16
#include "formatfloat.h"
17
#endif
18

19
STATIC void float_print(void (*print)(void *env, const char *fmt, ...), void *env, mp_obj_t o_in, mp_print_kind_t kind) {
20
    mp_obj_float_t *o = o_in;
21
22
23
24
25
#if MICROPY_FLOAT_IMPL == MICROPY_FLOAT_IMPL_FLOAT
    char buf[32];
    format_float(o->value, buf, sizeof(buf), 'g', 6, '\0');
    print(env, "%s", buf);
#else
mux's avatar
mux committed
26
    print(env, "%.8g", (double) o->value);
27
#endif
28
29
}

30
STATIC mp_obj_t float_make_new(mp_obj_t type_in, uint n_args, uint n_kw, const mp_obj_t *args) {
31
32
    // TODO check n_kw == 0

33
34
35
36
37
    switch (n_args) {
        case 0:
            return mp_obj_new_float(0);

        case 1:
38
39
40
41
            if (MP_OBJ_IS_STR(args[0])) {
                // a string, parse it
                uint l;
                const char *s = mp_obj_str_get_data(args[0], &l);
42
                return mp_parse_num_decimal(s, l, false, false);
43
            } else if (MP_OBJ_IS_TYPE(args[0], &mp_type_float)) {
44
                // a float, just return it
45
46
                return args[0];
            } else {
47
                // something else, try to cast it to a float
48
49
50
51
                return mp_obj_new_float(mp_obj_get_float(args[0]));
            }

        default:
52
            nlr_jump(mp_obj_new_exception_msg_varg(&mp_type_TypeError, "float takes at most 1 argument, %d given", n_args));
53
54
55
    }
}

56
STATIC mp_obj_t float_unary_op(int op, mp_obj_t o_in) {
57
58
    mp_obj_float_t *o = o_in;
    switch (op) {
Damien George's avatar
Damien George committed
59
60
61
        case MP_UNARY_OP_BOOL: return MP_BOOL(o->value != 0);
        case MP_UNARY_OP_POSITIVE: return o_in;
        case MP_UNARY_OP_NEGATIVE: return mp_obj_new_float(-o->value);
62
63
64
65
        default: return NULL; // op not supported
    }
}

66
STATIC mp_obj_t float_binary_op(int op, mp_obj_t lhs_in, mp_obj_t rhs_in) {
67
    mp_obj_float_t *lhs = lhs_in;
68
    if (MP_OBJ_IS_TYPE(rhs_in, &mp_type_complex)) {
69
70
71
        return mp_obj_complex_binary_op(op, lhs->value, 0, rhs_in);
    } else {
        return mp_obj_float_binary_op(op, lhs->value, rhs_in);
72
73
74
    }
}

75
const mp_obj_type_t mp_type_float = {
76
    { &mp_type_type },
77
    .name = MP_QSTR_float,
78
79
80
81
    .print = float_print,
    .make_new = float_make_new,
    .unary_op = float_unary_op,
    .binary_op = float_binary_op,
82
83
84
85
};

mp_obj_t mp_obj_new_float(mp_float_t value) {
    mp_obj_float_t *o = m_new(mp_obj_float_t, 1);
86
    o->base.type = &mp_type_float;
87
88
89
90
91
    o->value = value;
    return (mp_obj_t)o;
}

mp_float_t mp_obj_float_get(mp_obj_t self_in) {
92
    assert(MP_OBJ_IS_TYPE(self_in, &mp_type_float));
93
94
95
96
    mp_obj_float_t *self = self_in;
    return self->value;
}

97
98
99
mp_obj_t mp_obj_float_binary_op(int op, mp_float_t lhs_val, mp_obj_t rhs_in) {
    mp_float_t rhs_val = mp_obj_get_float(rhs_in); // can be any type, this function will convert to float (if possible)
    switch (op) {
Damien George's avatar
Damien George committed
100
101
102
103
104
105
        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;
106
        /* TODO floor(?) the value
Damien George's avatar
Damien George committed
107
108
        case MP_BINARY_OP_FLOOR_DIVIDE:
        case MP_BINARY_OP_INPLACE_FLOOR_DIVIDE: val = lhs_val / rhs_val; break;
109
        */
Damien George's avatar
Damien George committed
110
111
        case MP_BINARY_OP_TRUE_DIVIDE:
        case MP_BINARY_OP_INPLACE_TRUE_DIVIDE: 
112
113
114
115
116
            lhs_val /= rhs_val; 
            if (isinf(lhs_val)){ // check for division by zero
                nlr_jump(mp_obj_new_exception_msg(&mp_type_ZeroDivisionError, "float division by zero"));
            }
            break;
Damien George's avatar
Damien George committed
117
118
119
120
        case MP_BINARY_OP_LESS: return MP_BOOL(lhs_val < rhs_val);
        case MP_BINARY_OP_MORE: return MP_BOOL(lhs_val > rhs_val);
        case MP_BINARY_OP_LESS_EQUAL: return MP_BOOL(lhs_val <= rhs_val);
        case MP_BINARY_OP_MORE_EQUAL: return MP_BOOL(lhs_val >= rhs_val);
John R. Lenton's avatar
John R. Lenton committed
121

122
123
        default:
            return NULL; // op not supported
124
125
126
127
    }
    return mp_obj_new_float(lhs_val);
}

128
#endif // MICROPY_ENABLE_FLOAT