ba8b5c17 |
1 | /* |
2 | ** $Id: lmathlib.c,v 1.83 2013/03/07 18:21:32 roberto Exp $ |
3 | ** Standard mathematical library |
4 | ** See Copyright Notice in lua.h |
5 | */ |
6 | |
7 | |
8 | #include <stdlib.h> |
9 | #include <math.h> |
10 | |
11 | #define lmathlib_c |
12 | #define LUA_LIB |
13 | |
14 | #include "lua.h" |
15 | |
16 | #include "lauxlib.h" |
17 | #include "lualib.h" |
18 | |
19 | |
20 | #undef PI |
21 | #define PI ((lua_Number)(3.1415926535897932384626433832795)) |
22 | #define RADIANS_PER_DEGREE ((lua_Number)(PI/180.0)) |
23 | |
24 | |
25 | |
26 | static int math_abs (lua_State *L) { |
27 | lua_pushnumber(L, l_mathop(fabs)(luaL_checknumber(L, 1))); |
28 | return 1; |
29 | } |
30 | |
31 | static int math_sin (lua_State *L) { |
32 | lua_pushnumber(L, l_mathop(sin)(luaL_checknumber(L, 1))); |
33 | return 1; |
34 | } |
35 | |
36 | static int math_sinh (lua_State *L) { |
37 | lua_pushnumber(L, l_mathop(sinh)(luaL_checknumber(L, 1))); |
38 | return 1; |
39 | } |
40 | |
41 | static int math_cos (lua_State *L) { |
42 | lua_pushnumber(L, l_mathop(cos)(luaL_checknumber(L, 1))); |
43 | return 1; |
44 | } |
45 | |
46 | static int math_cosh (lua_State *L) { |
47 | lua_pushnumber(L, l_mathop(cosh)(luaL_checknumber(L, 1))); |
48 | return 1; |
49 | } |
50 | |
51 | static int math_tan (lua_State *L) { |
52 | lua_pushnumber(L, l_mathop(tan)(luaL_checknumber(L, 1))); |
53 | return 1; |
54 | } |
55 | |
56 | static int math_tanh (lua_State *L) { |
57 | lua_pushnumber(L, l_mathop(tanh)(luaL_checknumber(L, 1))); |
58 | return 1; |
59 | } |
60 | |
61 | static int math_asin (lua_State *L) { |
62 | lua_pushnumber(L, l_mathop(asin)(luaL_checknumber(L, 1))); |
63 | return 1; |
64 | } |
65 | |
66 | static int math_acos (lua_State *L) { |
67 | lua_pushnumber(L, l_mathop(acos)(luaL_checknumber(L, 1))); |
68 | return 1; |
69 | } |
70 | |
71 | static int math_atan (lua_State *L) { |
72 | lua_pushnumber(L, l_mathop(atan)(luaL_checknumber(L, 1))); |
73 | return 1; |
74 | } |
75 | |
76 | static int math_atan2 (lua_State *L) { |
77 | lua_pushnumber(L, l_mathop(atan2)(luaL_checknumber(L, 1), |
78 | luaL_checknumber(L, 2))); |
79 | return 1; |
80 | } |
81 | |
82 | static int math_ceil (lua_State *L) { |
83 | lua_pushnumber(L, l_mathop(ceil)(luaL_checknumber(L, 1))); |
84 | return 1; |
85 | } |
86 | |
87 | static int math_floor (lua_State *L) { |
88 | lua_pushnumber(L, l_mathop(floor)(luaL_checknumber(L, 1))); |
89 | return 1; |
90 | } |
91 | |
92 | static int math_fmod (lua_State *L) { |
93 | lua_pushnumber(L, l_mathop(fmod)(luaL_checknumber(L, 1), |
94 | luaL_checknumber(L, 2))); |
95 | return 1; |
96 | } |
97 | |
98 | static int math_modf (lua_State *L) { |
99 | lua_Number ip; |
100 | lua_Number fp = l_mathop(modf)(luaL_checknumber(L, 1), &ip); |
101 | lua_pushnumber(L, ip); |
102 | lua_pushnumber(L, fp); |
103 | return 2; |
104 | } |
105 | |
106 | static int math_sqrt (lua_State *L) { |
107 | lua_pushnumber(L, l_mathop(sqrt)(luaL_checknumber(L, 1))); |
108 | return 1; |
109 | } |
110 | |
111 | static int math_pow (lua_State *L) { |
112 | lua_Number x = luaL_checknumber(L, 1); |
113 | lua_Number y = luaL_checknumber(L, 2); |
114 | lua_pushnumber(L, l_mathop(pow)(x, y)); |
115 | return 1; |
116 | } |
117 | |
118 | static int math_log (lua_State *L) { |
119 | lua_Number x = luaL_checknumber(L, 1); |
120 | lua_Number res; |
121 | if (lua_isnoneornil(L, 2)) |
122 | res = l_mathop(log)(x); |
123 | else { |
124 | lua_Number base = luaL_checknumber(L, 2); |
125 | if (base == (lua_Number)10.0) res = l_mathop(log10)(x); |
126 | else res = l_mathop(log)(x)/l_mathop(log)(base); |
127 | } |
128 | lua_pushnumber(L, res); |
129 | return 1; |
130 | } |
131 | |
132 | #if defined(LUA_COMPAT_LOG10) |
133 | static int math_log10 (lua_State *L) { |
134 | lua_pushnumber(L, l_mathop(log10)(luaL_checknumber(L, 1))); |
135 | return 1; |
136 | } |
137 | #endif |
138 | |
139 | static int math_exp (lua_State *L) { |
140 | lua_pushnumber(L, l_mathop(exp)(luaL_checknumber(L, 1))); |
141 | return 1; |
142 | } |
143 | |
144 | static int math_deg (lua_State *L) { |
145 | lua_pushnumber(L, luaL_checknumber(L, 1)/RADIANS_PER_DEGREE); |
146 | return 1; |
147 | } |
148 | |
149 | static int math_rad (lua_State *L) { |
150 | lua_pushnumber(L, luaL_checknumber(L, 1)*RADIANS_PER_DEGREE); |
151 | return 1; |
152 | } |
153 | |
154 | static int math_frexp (lua_State *L) { |
155 | int e; |
156 | lua_pushnumber(L, l_mathop(frexp)(luaL_checknumber(L, 1), &e)); |
157 | lua_pushinteger(L, e); |
158 | return 2; |
159 | } |
160 | |
161 | static int math_ldexp (lua_State *L) { |
162 | lua_Number x = luaL_checknumber(L, 1); |
163 | int ep = luaL_checkint(L, 2); |
164 | lua_pushnumber(L, l_mathop(ldexp)(x, ep)); |
165 | return 1; |
166 | } |
167 | |
168 | |
169 | |
170 | static int math_min (lua_State *L) { |
171 | int n = lua_gettop(L); /* number of arguments */ |
172 | lua_Number dmin = luaL_checknumber(L, 1); |
173 | int i; |
174 | for (i=2; i<=n; i++) { |
175 | lua_Number d = luaL_checknumber(L, i); |
176 | if (d < dmin) |
177 | dmin = d; |
178 | } |
179 | lua_pushnumber(L, dmin); |
180 | return 1; |
181 | } |
182 | |
183 | |
184 | static int math_max (lua_State *L) { |
185 | int n = lua_gettop(L); /* number of arguments */ |
186 | lua_Number dmax = luaL_checknumber(L, 1); |
187 | int i; |
188 | for (i=2; i<=n; i++) { |
189 | lua_Number d = luaL_checknumber(L, i); |
190 | if (d > dmax) |
191 | dmax = d; |
192 | } |
193 | lua_pushnumber(L, dmax); |
194 | return 1; |
195 | } |
196 | |
197 | |
198 | static int math_random (lua_State *L) { |
199 | /* the `%' avoids the (rare) case of r==1, and is needed also because on |
200 | some systems (SunOS!) `rand()' may return a value larger than RAND_MAX */ |
201 | lua_Number r = (lua_Number)(rand()%RAND_MAX) / (lua_Number)RAND_MAX; |
202 | switch (lua_gettop(L)) { /* check number of arguments */ |
203 | case 0: { /* no arguments */ |
204 | lua_pushnumber(L, r); /* Number between 0 and 1 */ |
205 | break; |
206 | } |
207 | case 1: { /* only upper limit */ |
208 | lua_Number u = luaL_checknumber(L, 1); |
209 | luaL_argcheck(L, (lua_Number)1.0 <= u, 1, "interval is empty"); |
210 | lua_pushnumber(L, l_mathop(floor)(r*u) + (lua_Number)(1.0)); /* [1, u] */ |
211 | break; |
212 | } |
213 | case 2: { /* lower and upper limits */ |
214 | lua_Number l = luaL_checknumber(L, 1); |
215 | lua_Number u = luaL_checknumber(L, 2); |
216 | luaL_argcheck(L, l <= u, 2, "interval is empty"); |
217 | lua_pushnumber(L, l_mathop(floor)(r*(u-l+1)) + l); /* [l, u] */ |
218 | break; |
219 | } |
220 | default: return luaL_error(L, "wrong number of arguments"); |
221 | } |
222 | return 1; |
223 | } |
224 | |
225 | |
226 | static int math_randomseed (lua_State *L) { |
227 | srand(luaL_checkunsigned(L, 1)); |
228 | (void)rand(); /* discard first value to avoid undesirable correlations */ |
229 | return 0; |
230 | } |
231 | |
232 | |
233 | static const luaL_Reg mathlib[] = { |
234 | {"abs", math_abs}, |
235 | {"acos", math_acos}, |
236 | {"asin", math_asin}, |
237 | {"atan2", math_atan2}, |
238 | {"atan", math_atan}, |
239 | {"ceil", math_ceil}, |
240 | {"cosh", math_cosh}, |
241 | {"cos", math_cos}, |
242 | {"deg", math_deg}, |
243 | {"exp", math_exp}, |
244 | {"floor", math_floor}, |
245 | {"fmod", math_fmod}, |
246 | {"frexp", math_frexp}, |
247 | {"ldexp", math_ldexp}, |
248 | #if defined(LUA_COMPAT_LOG10) |
249 | {"log10", math_log10}, |
250 | #endif |
251 | {"log", math_log}, |
252 | {"max", math_max}, |
253 | {"min", math_min}, |
254 | {"modf", math_modf}, |
255 | {"pow", math_pow}, |
256 | {"rad", math_rad}, |
257 | {"random", math_random}, |
258 | {"randomseed", math_randomseed}, |
259 | {"sinh", math_sinh}, |
260 | {"sin", math_sin}, |
261 | {"sqrt", math_sqrt}, |
262 | {"tanh", math_tanh}, |
263 | {"tan", math_tan}, |
264 | {NULL, NULL} |
265 | }; |
266 | |
267 | |
268 | /* |
269 | ** Open math library |
270 | */ |
271 | LUAMOD_API int luaopen_math (lua_State *L) { |
272 | luaL_newlib(L, mathlib); |
273 | lua_pushnumber(L, PI); |
274 | lua_setfield(L, -2, "pi"); |
275 | lua_pushnumber(L, HUGE_VAL); |
276 | lua_setfield(L, -2, "huge"); |
277 | return 1; |
278 | } |
279 | |