diff options
Diffstat (limited to 'src/lj_vmmath.c')
-rw-r--r-- | src/lj_vmmath.c | 75 |
1 files changed, 37 insertions, 38 deletions
diff --git a/src/lj_vmmath.c b/src/lj_vmmath.c index ff41ba28..536199d8 100644 --- a/src/lj_vmmath.c +++ b/src/lj_vmmath.c | |||
@@ -13,16 +13,29 @@ | |||
13 | #include "lj_ir.h" | 13 | #include "lj_ir.h" |
14 | #include "lj_vm.h" | 14 | #include "lj_vm.h" |
15 | 15 | ||
16 | /* -- Helper functions for generated machine code ------------------------- */ | 16 | /* -- Wrapper functions --------------------------------------------------- */ |
17 | 17 | ||
18 | #if LJ_TARGET_X86ORX64 | 18 | #if LJ_TARGET_X86 && __ELF__ && __PIC__ |
19 | /* Wrapper functions to avoid linker issues on OSX. */ | 19 | /* Wrapper functions to deal with the ELF/x86 PIC disaster. */ |
20 | LJ_FUNCA double lj_vm_sinh(double x) { return sinh(x); } | 20 | LJ_FUNCA double lj_wrap_log(double x) { return log(x); } |
21 | LJ_FUNCA double lj_vm_cosh(double x) { return cosh(x); } | 21 | LJ_FUNCA double lj_wrap_log10(double x) { return log10(x); } |
22 | LJ_FUNCA double lj_vm_tanh(double x) { return tanh(x); } | 22 | LJ_FUNCA double lj_wrap_exp(double x) { return exp(x); } |
23 | LJ_FUNCA double lj_wrap_sin(double x) { return sin(x); } | ||
24 | LJ_FUNCA double lj_wrap_cos(double x) { return cos(x); } | ||
25 | LJ_FUNCA double lj_wrap_tan(double x) { return tan(x); } | ||
26 | LJ_FUNCA double lj_wrap_asin(double x) { return asin(x); } | ||
27 | LJ_FUNCA double lj_wrap_acos(double x) { return acos(x); } | ||
28 | LJ_FUNCA double lj_wrap_atan(double x) { return atan(x); } | ||
29 | LJ_FUNCA double lj_wrap_sinh(double x) { return sinh(x); } | ||
30 | LJ_FUNCA double lj_wrap_cosh(double x) { return cosh(x); } | ||
31 | LJ_FUNCA double lj_wrap_tanh(double x) { return tanh(x); } | ||
32 | LJ_FUNCA double lj_wrap_atan2(double x, double y) { return atan2(x, y); } | ||
33 | LJ_FUNCA double lj_wrap_pow(double x, double y) { return pow(x, y); } | ||
34 | LJ_FUNCA double lj_wrap_fmod(double x, double y) { return fmod(x, y); } | ||
23 | #endif | 35 | #endif |
24 | 36 | ||
25 | #if !LJ_TARGET_X86ORX64 | 37 | /* -- Helper functions for generated machine code ------------------------- */ |
38 | |||
26 | double lj_vm_foldarith(double x, double y, int op) | 39 | double lj_vm_foldarith(double x, double y, int op) |
27 | { | 40 | { |
28 | switch (op) { | 41 | switch (op) { |
@@ -35,37 +48,20 @@ double lj_vm_foldarith(double x, double y, int op) | |||
35 | case IR_NEG - IR_ADD: return -x; break; | 48 | case IR_NEG - IR_ADD: return -x; break; |
36 | case IR_ABS - IR_ADD: return fabs(x); break; | 49 | case IR_ABS - IR_ADD: return fabs(x); break; |
37 | #if LJ_HASJIT | 50 | #if LJ_HASJIT |
38 | case IR_ATAN2 - IR_ADD: return atan2(x, y); break; | ||
39 | case IR_LDEXP - IR_ADD: return ldexp(x, (int)y); break; | 51 | case IR_LDEXP - IR_ADD: return ldexp(x, (int)y); break; |
40 | case IR_MIN - IR_ADD: return x > y ? y : x; break; | 52 | case IR_MIN - IR_ADD: return x < y ? x : y; break; |
41 | case IR_MAX - IR_ADD: return x < y ? y : x; break; | 53 | case IR_MAX - IR_ADD: return x > y ? x : y; break; |
42 | #endif | 54 | #endif |
43 | default: return x; | 55 | default: return x; |
44 | } | 56 | } |
45 | } | 57 | } |
46 | #endif | ||
47 | |||
48 | #if LJ_HASJIT | ||
49 | 58 | ||
50 | #ifdef LUAJIT_NO_LOG2 | 59 | #if (LJ_HASJIT && !(LJ_TARGET_ARM || LJ_TARGET_ARM64 || LJ_TARGET_PPC)) || LJ_TARGET_MIPS |
51 | double lj_vm_log2(double a) | ||
52 | { | ||
53 | return log(a) * 1.4426950408889634074; | ||
54 | } | ||
55 | #endif | ||
56 | |||
57 | #ifdef LUAJIT_NO_EXP2 | ||
58 | double lj_vm_exp2(double a) | ||
59 | { | ||
60 | return exp(a * 0.6931471805599453); | ||
61 | } | ||
62 | #endif | ||
63 | |||
64 | #if !(LJ_TARGET_ARM || LJ_TARGET_PPC) | ||
65 | int32_t LJ_FASTCALL lj_vm_modi(int32_t a, int32_t b) | 60 | int32_t LJ_FASTCALL lj_vm_modi(int32_t a, int32_t b) |
66 | { | 61 | { |
67 | uint32_t y, ua, ub; | 62 | uint32_t y, ua, ub; |
68 | lua_assert(b != 0); /* This must be checked before using this function. */ | 63 | /* This must be checked before using this function. */ |
64 | lj_assertX(b != 0, "modulo with zero divisor"); | ||
69 | ua = a < 0 ? (uint32_t)-a : (uint32_t)a; | 65 | ua = a < 0 ? (uint32_t)-a : (uint32_t)a; |
70 | ub = b < 0 ? (uint32_t)-b : (uint32_t)b; | 66 | ub = b < 0 ? (uint32_t)-b : (uint32_t)b; |
71 | y = ua % ub; | 67 | y = ua % ub; |
@@ -75,12 +71,21 @@ int32_t LJ_FASTCALL lj_vm_modi(int32_t a, int32_t b) | |||
75 | } | 71 | } |
76 | #endif | 72 | #endif |
77 | 73 | ||
74 | #if LJ_HASJIT | ||
75 | |||
76 | #ifdef LUAJIT_NO_LOG2 | ||
77 | double lj_vm_log2(double a) | ||
78 | { | ||
79 | return log(a) * 1.4426950408889634074; | ||
80 | } | ||
81 | #endif | ||
82 | |||
78 | #if !LJ_TARGET_X86ORX64 | 83 | #if !LJ_TARGET_X86ORX64 |
79 | /* Unsigned x^k. */ | 84 | /* Unsigned x^k. */ |
80 | static double lj_vm_powui(double x, uint32_t k) | 85 | static double lj_vm_powui(double x, uint32_t k) |
81 | { | 86 | { |
82 | double y; | 87 | double y; |
83 | lua_assert(k != 0); | 88 | lj_assertX(k != 0, "pow with zero exponent"); |
84 | for (; (k & 1) == 0; k >>= 1) x *= x; | 89 | for (; (k & 1) == 0; k >>= 1) x *= x; |
85 | y = x; | 90 | y = x; |
86 | if ((k >>= 1) != 0) { | 91 | if ((k >>= 1) != 0) { |
@@ -107,6 +112,7 @@ double lj_vm_powi(double x, int32_t k) | |||
107 | else | 112 | else |
108 | return 1.0 / lj_vm_powui(x, (uint32_t)-k); | 113 | return 1.0 / lj_vm_powui(x, (uint32_t)-k); |
109 | } | 114 | } |
115 | #endif | ||
110 | 116 | ||
111 | /* Computes fpm(x) for extended math functions. */ | 117 | /* Computes fpm(x) for extended math functions. */ |
112 | double lj_vm_foldfpm(double x, int fpm) | 118 | double lj_vm_foldfpm(double x, int fpm) |
@@ -116,19 +122,12 @@ double lj_vm_foldfpm(double x, int fpm) | |||
116 | case IRFPM_CEIL: return lj_vm_ceil(x); | 122 | case IRFPM_CEIL: return lj_vm_ceil(x); |
117 | case IRFPM_TRUNC: return lj_vm_trunc(x); | 123 | case IRFPM_TRUNC: return lj_vm_trunc(x); |
118 | case IRFPM_SQRT: return sqrt(x); | 124 | case IRFPM_SQRT: return sqrt(x); |
119 | case IRFPM_EXP: return exp(x); | ||
120 | case IRFPM_EXP2: return lj_vm_exp2(x); | ||
121 | case IRFPM_LOG: return log(x); | 125 | case IRFPM_LOG: return log(x); |
122 | case IRFPM_LOG2: return lj_vm_log2(x); | 126 | case IRFPM_LOG2: return lj_vm_log2(x); |
123 | case IRFPM_LOG10: return log10(x); | 127 | default: lj_assertX(0, "bad fpm %d", fpm); |
124 | case IRFPM_SIN: return sin(x); | ||
125 | case IRFPM_COS: return cos(x); | ||
126 | case IRFPM_TAN: return tan(x); | ||
127 | default: lua_assert(0); | ||
128 | } | 128 | } |
129 | return 0; | 129 | return 0; |
130 | } | 130 | } |
131 | #endif | ||
132 | 131 | ||
133 | #if LJ_HASFFI | 132 | #if LJ_HASFFI |
134 | int lj_vm_errno(void) | 133 | int lj_vm_errno(void) |