1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
|
/*
* Copyright © 2017 Keith Packard <keithp@keithp.com>
*
* This program is free software; you can redistribute it and/or modify
* it under the terms of the GNU General Public License as published by
* the Free Software Foundation, either version 2 of the License, or
* (at your option) any later version.
*
* This program is distributed in the hope that it will be useful, but
* WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
* General Public License for more details.
*/
#include "ao_scheme.h"
#include <math.h>
#ifdef AO_SCHEME_FEATURE_FLOAT
static void float_mark(void *addr)
{
(void) addr;
}
static int float_size(void *addr)
{
if (!addr)
return 0;
return sizeof (struct ao_scheme_float);
}
static void float_move(void *addr)
{
(void) addr;
}
const struct ao_scheme_type ao_scheme_float_type = {
.mark = float_mark,
.size = float_size,
.move = float_move,
.name = "float",
};
#ifndef FLOAT_FORMAT
#define FLOAT_FORMAT "%g"
#endif
void
ao_scheme_float_write(ao_poly p)
{
struct ao_scheme_float *f = ao_scheme_poly_float(p);
float v = f->value;
if (isnanf(v))
printf("+nan.0");
else if (isinff(v)) {
if (v < 0)
printf("-");
else
printf("+");
printf("inf.0");
} else
printf (FLOAT_FORMAT, v);
}
float
ao_scheme_poly_number(ao_poly p)
{
switch (ao_scheme_poly_base_type(p)) {
case AO_SCHEME_INT:
return ao_scheme_poly_int(p);
case AO_SCHEME_OTHER:
switch (ao_scheme_other_type(ao_scheme_poly_other(p))) {
case AO_SCHEME_BIGINT:
return ao_scheme_bigint_int(ao_scheme_poly_bigint(p)->value);
case AO_SCHEME_FLOAT:
return ao_scheme_poly_float(p)->value;
}
}
return NAN;
}
ao_poly
ao_scheme_float_get(float value)
{
struct ao_scheme_float *f;
f = ao_scheme_alloc(sizeof (struct ao_scheme_float));
f->type = AO_SCHEME_FLOAT;
f->value = value;
return ao_scheme_float_poly(f);
}
ao_poly
ao_scheme_do_inexactp(struct ao_scheme_cons *cons)
{
if (!ao_scheme_check_argc(_ao_scheme_atom_led, cons, 1, 1))
return AO_SCHEME_NIL;
if (ao_scheme_poly_type(ao_scheme_arg(cons, 0)) == AO_SCHEME_FLOAT)
return _ao_scheme_bool_true;
return _ao_scheme_bool_false;
}
ao_poly
ao_scheme_do_finitep(struct ao_scheme_cons *cons)
{
ao_poly value;
float f;
if (!ao_scheme_check_argc(_ao_scheme_atom_led, cons, 1, 1))
return AO_SCHEME_NIL;
value = ao_scheme_arg(cons, 0);
switch (ao_scheme_poly_type(value)) {
case AO_SCHEME_INT:
case AO_SCHEME_BIGINT:
return _ao_scheme_bool_true;
case AO_SCHEME_FLOAT:
f = ao_scheme_poly_float(value)->value;
if (!isnan(f) && !isinf(f))
return _ao_scheme_bool_true;
}
return _ao_scheme_bool_false;
}
ao_poly
ao_scheme_do_infinitep(struct ao_scheme_cons *cons)
{
ao_poly value;
float f;
if (!ao_scheme_check_argc(_ao_scheme_atom_led, cons, 1, 1))
return AO_SCHEME_NIL;
value = ao_scheme_arg(cons, 0);
switch (ao_scheme_poly_type(value)) {
case AO_SCHEME_FLOAT:
f = ao_scheme_poly_float(value)->value;
if (isinf(f))
return _ao_scheme_bool_true;
}
return _ao_scheme_bool_false;
}
ao_poly
ao_scheme_do_sqrt(struct ao_scheme_cons *cons)
{
ao_poly value;
if (!ao_scheme_check_argc(_ao_scheme_atom_sqrt, cons, 1, 1))
return AO_SCHEME_NIL;
value = ao_scheme_arg(cons, 0);
if (!ao_scheme_number_typep(ao_scheme_poly_type(value)))
return ao_scheme_error(AO_SCHEME_INVALID, "%s: non-numeric", ao_scheme_poly_atom(_ao_scheme_atom_sqrt)->name);
return ao_scheme_float_get(sqrtf(ao_scheme_poly_number(value)));
}
#endif
|
