src/lmathlib.c

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