-/* math.c - these are based on bob smith's csharp routines */\r
-\r
-#define __USE_ISOC99\r
-#include <math.h>\r
-#include <mono/metadata/sysmath.h>\r
-#include <mono/metadata/exception.h>\r
-\r
-#ifndef NAN\r
-# if G_BYTE_ORDER == G_BIG_ENDIAN\r
-# define __nan_bytes { 0x7f, 0xc0, 0, 0 }\r
-# endif\r
-# if G_BYTE_ORDER == G_LITTLE_ENDIAN\r
-# define __nan_bytes { 0, 0, 0xc0, 0x7f }\r
-# endif\r
-\r
-static union { unsigned char __c[4]; float __d; } __nan_union = { __nan_bytes };\r
-# define NAN (__nan_union.__d)\r
-#endif\r
-\r
-#ifndef HUGE_VAL\r
-#define __huge_val_t union { unsigned char __c[8]; double __d; }\r
-# if G_BYTE_ORDER == G_BIG_ENDIAN\r
-# define __HUGE_VAL_bytes { 0x7f, 0xf0, 0, 0, 0, 0, 0, 0 }\r
-# endif\r
-# if G_BYTE_ORDER == G_LITTLE_ENDIAN\r
-# define __HUGE_VAL_bytes { 0, 0, 0, 0, 0, 0, 0xf0, 0x7f }\r
-# endif\r
-static __huge_val_t __huge_val = { __HUGE_VAL_bytes };\r
-# define HUGE_VAL (__huge_val.__d)\r
-#endif\r
-\r
-\r
-gdouble ves_icall_System_Math_Floor (gdouble x) {\r
- MONO_ARCH_SAVE_REGS;\r
- return floor(x);\r
-}\r
-\r
-gdouble ves_icall_System_Math_Round (gdouble x) {\r
- double int_part, dec_part;\r
- MONO_ARCH_SAVE_REGS;\r
- int_part = floor(x);\r
- dec_part = x - int_part;\r
- if (((dec_part == 0.5) &&\r
- ((2.0 * ((int_part / 2.0) - floor(int_part / 2.0))) != 0.0)) ||\r
- (dec_part > 0.5)) {\r
- int_part++;\r
- }\r
- return int_part;\r
-}\r
-\r
-gdouble ves_icall_System_Math_Round2 (gdouble value, gint32 digits) {\r
- double p, int_part, dec_part;\r
- MONO_ARCH_SAVE_REGS;\r
- if (value == HUGE_VAL)\r
- return HUGE_VAL;\r
- if (value == -HUGE_VAL)\r
- return -HUGE_VAL;\r
- if (digits == 0)\r
- return ves_icall_System_Math_Round(value);\r
- p = pow(10, digits);\r
- int_part = floor(value);\r
- dec_part = value - int_part;\r
- dec_part *= 1000000000000000ULL;\r
- dec_part = floor(dec_part);\r
- dec_part /= (1000000000000000ULL / p);\r
- dec_part = ves_icall_System_Math_Round(dec_part);\r
- dec_part /= p;\r
- return int_part + dec_part;\r
-}\r
-\r
-gdouble \r
-ves_icall_System_Math_Sin (gdouble x)\r
-{\r
- MONO_ARCH_SAVE_REGS;\r
-\r
- return sin (x);\r
-}\r
-\r
-gdouble \r
-ves_icall_System_Math_Cos (gdouble x)\r
-{\r
- MONO_ARCH_SAVE_REGS;\r
-\r
- return cos (x);\r
-}\r
-\r
-gdouble \r
-ves_icall_System_Math_Tan (gdouble x)\r
-{\r
- MONO_ARCH_SAVE_REGS;\r
-\r
- return tan (x);\r
-}\r
-\r
-gdouble \r
-ves_icall_System_Math_Sinh (gdouble x)\r
-{\r
- MONO_ARCH_SAVE_REGS;\r
-\r
- return sinh (x);\r
-}\r
-\r
-gdouble \r
-ves_icall_System_Math_Cosh (gdouble x)\r
-{\r
- MONO_ARCH_SAVE_REGS;\r
-\r
- return cosh (x);\r
-}\r
-\r
-gdouble \r
-ves_icall_System_Math_Tanh (gdouble x)\r
-{\r
- MONO_ARCH_SAVE_REGS;\r
-\r
- return tanh (x);\r
-}\r
-\r
-gdouble \r
-ves_icall_System_Math_Acos (gdouble x)\r
-{\r
- MONO_ARCH_SAVE_REGS;\r
-\r
- if (x < -1 || x > 1)\r
- return NAN;\r
-\r
- return acos (x);\r
-}\r
-\r
-gdouble \r
-ves_icall_System_Math_Asin (gdouble x)\r
-{\r
- MONO_ARCH_SAVE_REGS;\r
-\r
- if (x < -1 || x > 1)\r
- return NAN;\r
-\r
- return asin (x);\r
-}\r
-\r
-gdouble \r
-ves_icall_System_Math_Atan (gdouble x)\r
-{\r
- MONO_ARCH_SAVE_REGS;\r
-\r
- return atan (x);\r
-}\r
-\r
-gdouble \r
-ves_icall_System_Math_Atan2 (gdouble y, gdouble x)\r
-{\r
- double result;\r
- MONO_ARCH_SAVE_REGS;\r
-\r
- if ((y == HUGE_VAL && x == HUGE_VAL) ||\r
- (y == HUGE_VAL && x == -HUGE_VAL) ||\r
- (y == -HUGE_VAL && x == HUGE_VAL) ||\r
- (y == -HUGE_VAL && x == -HUGE_VAL)) {\r
- return NAN;\r
- }\r
- result = atan2 (y, x);\r
- return (result == -0)? 0: result;\r
-}\r
-\r
-gdouble \r
-ves_icall_System_Math_Exp (gdouble x)\r
-{\r
- MONO_ARCH_SAVE_REGS;\r
-\r
- return exp (x);\r
-}\r
-\r
-gdouble \r
-ves_icall_System_Math_Log (gdouble x)\r
-{\r
- MONO_ARCH_SAVE_REGS;\r
-\r
- if (x == 0)\r
- return -HUGE_VAL;\r
- else if (x < 0)\r
- return NAN;\r
-\r
- return log (x);\r
-}\r
-\r
-gdouble \r
-ves_icall_System_Math_Log10 (gdouble x)\r
-{\r
- MONO_ARCH_SAVE_REGS;\r
-\r
- if (x == 0)\r
- return -HUGE_VAL;\r
- else if (x < 0)\r
- return NAN;\r
-\r
- return log10 (x);\r
-}\r
-\r
-gdouble \r
-ves_icall_System_Math_Pow (gdouble x, gdouble y)\r
-{\r
- double result;\r
- MONO_ARCH_SAVE_REGS;\r
-\r
- if (isnan(x) || isnan(y)) {\r
- return NAN;\r
- }\r
-\r
- if ((x == 1 || x == -1) && (y == HUGE_VAL || y == -HUGE_VAL)) {\r
- return NAN;\r
- }\r
-\r
- /* This code is for return the same results as MS.NET for certain\r
- * limit values */\r
- if (x < -9007199254740991.0) {\r
- if (y > 9007199254740991.0)\r
- return HUGE_VAL;\r
- if (y < -9007199254740991.0)\r
- return 0;\r
- }\r
-\r
- result = pow (x, y);\r
-\r
- /* This code is for return the same results as MS.NET for certain\r
- * limit values */\r
- if (isnan(result) &&\r
- (x == -1.0) &&\r
- ((y > 9007199254740991.0) || (y < -9007199254740991.0))) {\r
- return 1;\r
- }\r
-\r
- return (result == -0)? 0: result;\r
-}\r
-\r
-gdouble \r
-ves_icall_System_Math_Sqrt (gdouble x)\r
-{\r
- MONO_ARCH_SAVE_REGS;\r
-\r
- if (x < 0)\r
- return NAN;\r
-\r
- return sqrt (x);\r
-}\r
+/* math.c - these are based on bob smith's csharp routines */
+
+#define __USE_ISOC99
+#include <math.h>
+#include <mono/metadata/sysmath.h>
+#include <mono/metadata/exception.h>
+
+#ifndef NAN
+# if G_BYTE_ORDER == G_BIG_ENDIAN
+# define __nan_bytes { 0x7f, 0xc0, 0, 0 }
+# endif
+# if G_BYTE_ORDER == G_LITTLE_ENDIAN
+# define __nan_bytes { 0, 0, 0xc0, 0x7f }
+# endif
+
+static union { unsigned char __c[4]; float __d; } __nan_union = { __nan_bytes };
+# define NAN (__nan_union.__d)
+#endif
+
+#ifndef HUGE_VAL
+#define __huge_val_t union { unsigned char __c[8]; double __d; }
+# if G_BYTE_ORDER == G_BIG_ENDIAN
+# define __HUGE_VAL_bytes { 0x7f, 0xf0, 0, 0, 0, 0, 0, 0 }
+# endif
+# if G_BYTE_ORDER == G_LITTLE_ENDIAN
+# define __HUGE_VAL_bytes { 0, 0, 0, 0, 0, 0, 0xf0, 0x7f }
+# endif
+static __huge_val_t __huge_val = { __HUGE_VAL_bytes };
+# define HUGE_VAL (__huge_val.__d)
+#endif
+
+
+gdouble ves_icall_System_Math_Floor (gdouble x) {
+ MONO_ARCH_SAVE_REGS;
+ return floor(x);
+}
+
+gdouble ves_icall_System_Math_Round (gdouble x) {
+ double int_part, dec_part;
+ MONO_ARCH_SAVE_REGS;
+ int_part = floor(x);
+ dec_part = x - int_part;
+ if (((dec_part == 0.5) &&
+ ((2.0 * ((int_part / 2.0) - floor(int_part / 2.0))) != 0.0)) ||
+ (dec_part > 0.5)) {
+ int_part++;
+ }
+ return int_part;
+}
+
+gdouble ves_icall_System_Math_Round2 (gdouble value, gint32 digits) {
+ double p, int_part, dec_part;
+ MONO_ARCH_SAVE_REGS;
+ if (value == HUGE_VAL)
+ return HUGE_VAL;
+ if (value == -HUGE_VAL)
+ return -HUGE_VAL;
+ if (digits == 0)
+ return ves_icall_System_Math_Round(value);
+ p = pow(10, digits);
+ int_part = floor(value);
+ dec_part = value - int_part;
+ dec_part *= 1000000000000000ULL;
+ dec_part = floor(dec_part);
+ dec_part /= (1000000000000000ULL / p);
+ dec_part = ves_icall_System_Math_Round(dec_part);
+ dec_part /= p;
+ return int_part + dec_part;
+}
+
+gdouble
+ves_icall_System_Math_Sin (gdouble x)
+{
+ MONO_ARCH_SAVE_REGS;
+
+ return sin (x);
+}
+
+gdouble
+ves_icall_System_Math_Cos (gdouble x)
+{
+ MONO_ARCH_SAVE_REGS;
+
+ return cos (x);
+}
+
+gdouble
+ves_icall_System_Math_Tan (gdouble x)
+{
+ MONO_ARCH_SAVE_REGS;
+
+ return tan (x);
+}
+
+gdouble
+ves_icall_System_Math_Sinh (gdouble x)
+{
+ MONO_ARCH_SAVE_REGS;
+
+ return sinh (x);
+}
+
+gdouble
+ves_icall_System_Math_Cosh (gdouble x)
+{
+ MONO_ARCH_SAVE_REGS;
+
+ return cosh (x);
+}
+
+gdouble
+ves_icall_System_Math_Tanh (gdouble x)
+{
+ MONO_ARCH_SAVE_REGS;
+
+ return tanh (x);
+}
+
+gdouble
+ves_icall_System_Math_Acos (gdouble x)
+{
+ MONO_ARCH_SAVE_REGS;
+
+ if (x < -1 || x > 1)
+ return NAN;
+
+ return acos (x);
+}
+
+gdouble
+ves_icall_System_Math_Asin (gdouble x)
+{
+ MONO_ARCH_SAVE_REGS;
+
+ if (x < -1 || x > 1)
+ return NAN;
+
+ return asin (x);
+}
+
+gdouble
+ves_icall_System_Math_Atan (gdouble x)
+{
+ MONO_ARCH_SAVE_REGS;
+
+ return atan (x);
+}
+
+gdouble
+ves_icall_System_Math_Atan2 (gdouble y, gdouble x)
+{
+ double result;
+ MONO_ARCH_SAVE_REGS;
+
+ if ((y == HUGE_VAL && x == HUGE_VAL) ||
+ (y == HUGE_VAL && x == -HUGE_VAL) ||
+ (y == -HUGE_VAL && x == HUGE_VAL) ||
+ (y == -HUGE_VAL && x == -HUGE_VAL)) {
+ return NAN;
+ }
+ result = atan2 (y, x);
+ return (result == -0)? 0: result;
+}
+
+gdouble
+ves_icall_System_Math_Exp (gdouble x)
+{
+ MONO_ARCH_SAVE_REGS;
+
+ return exp (x);
+}
+
+gdouble
+ves_icall_System_Math_Log (gdouble x)
+{
+ MONO_ARCH_SAVE_REGS;
+
+ if (x == 0)
+ return -HUGE_VAL;
+ else if (x < 0)
+ return NAN;
+
+ return log (x);
+}
+
+gdouble
+ves_icall_System_Math_Log10 (gdouble x)
+{
+ MONO_ARCH_SAVE_REGS;
+
+ if (x == 0)
+ return -HUGE_VAL;
+ else if (x < 0)
+ return NAN;
+
+ return log10 (x);
+}
+
+gdouble
+ves_icall_System_Math_Pow (gdouble x, gdouble y)
+{
+ double result;
+ MONO_ARCH_SAVE_REGS;
+
+ if (isnan(x) || isnan(y)) {
+ return NAN;
+ }
+
+ if ((x == 1 || x == -1) && (y == HUGE_VAL || y == -HUGE_VAL)) {
+ return NAN;
+ }
+
+ /* This code is for return the same results as MS.NET for certain
+ * limit values */
+ if (x < -9007199254740991.0) {
+ if (y > 9007199254740991.0)
+ return HUGE_VAL;
+ if (y < -9007199254740991.0)
+ return 0;
+ }
+
+ result = pow (x, y);
+
+ /* This code is for return the same results as MS.NET for certain
+ * limit values */
+ if (isnan(result) &&
+ (x == -1.0) &&
+ ((y > 9007199254740991.0) || (y < -9007199254740991.0))) {
+ return 1;
+ }
+
+ return (result == -0)? 0: result;
+}
+
+gdouble
+ves_icall_System_Math_Sqrt (gdouble x)
+{
+ MONO_ARCH_SAVE_REGS;
+
+ if (x < 0)
+ return NAN;
+
+ return sqrt (x);
+}
projects / mono.git / blobdiff