Move libphoenix math library implementation

Modify libphoenix Makefile, move libphoenix math library implementation, add definitions of float functions, decouple strtod.c from math library

JIRA: RTOS-1132

Author: Mikolaj Matalowski

Makefile

@@ -46,7 +46,6 @@ include ctype/Makefile
 include err/Makefile
 include errno/Makefile
 include locale/Makefile
-include math/Makefile
 include misc/Makefile
 include net/Makefile
 include netinet/Makefile
@@ -65,6 +64,7 @@ include time/Makefile
 include unistd/Makefile
 include wchar/Makefile
 include ubsan/Makefile
+include libm/Makefile
 
 #include test/Makefile
 
@@ -82,7 +82,8 @@ install: install-headers install-libs
 install-headers: $(SRCHEADERS)
 	@echo INSTALL "$(HEADERS_INSTALL_DIR)/*"; \
 	mkdir -p "$(HEADERS_INSTALL_DIR)"; \
-	cp -a include/* "$(HEADERS_INSTALL_DIR)";
+	cp -a include/* "$(HEADERS_INSTALL_DIR)"; \
+	cp -a libm/libmcs/libm/include/* "$(HEADERS_INSTALL_DIR)"
 
 # TODO: remove `rm crt0.o` when we will be sure it's not a symlink to libphoenix.a anymore
 install-libs: $(LIB_TARGETS)

math/common.c renamed to libm/math/common.c

@@ -19,7 +19,6 @@
 #include <float.h>
 #include "common.h"
 
-
 void normalizeSub(double *x, int *exp)
 {
 	conv_t *conv = (conv_t *)x;

math/common.h renamed to libm/math/common.h

@@ -0,0 +1,15 @@
+#include <math.h>
+
+int __fpclassifyf(float x)
+{
+	return __builtin_fpclassify(FP_NAN, FP_INFINITE, FP_NORMAL, FP_SUBNORMAL, FP_ZERO, x);
+}
+
+
+int __fpclassifyd(double x)
+{
+	return __builtin_fpclassify(FP_NAN, FP_INFINITE, FP_NORMAL, FP_SUBNORMAL, FP_ZERO, x);
+}
+
+const float __inff = __builtin_inff();
+const double __infd = __builtin_inf();

libm/math/compatibility.c

@@ -47,6 +47,12 @@ double frexp(double x, int *exp)
 }
 
 
+float frexpf(float x, int *exp)
+{
+	return (float)frexp((double)x, exp);
+}
+
+
 double ldexp(double x, int exp)
 {
 	if (isnan(x) != 0) {
@@ -84,6 +90,12 @@ double ldexp(double x, int exp)
 }
 
 
+float ldexpf(float x, int exp)
+{
+	return (float)ldexp((double)x, exp);
+}
+
+
 double log(double x)
 {
 	double tmp, pow, res;
@@ -129,6 +141,12 @@ double log(double x)
 }
 
 
+float logf(float x)
+{
+	return (float)log((double)x);
+}
+
+
 double log2(double x)
 {
 	return (log(x) / M_LN2);
@@ -142,6 +160,12 @@ double log10(double x)
 }
 
 
+float log10f(float x)
+{
+	return (float)log10((double)x);
+}
+
+
 double modf(double x, double *intpart)
 {
 	conv_t *conv = (conv_t *)&x;
@@ -232,6 +256,12 @@ double exp(double x)
 }
 
 
+float expf(float x)
+{
+	return (float)exp((double)x);
+}
+
+
 double ceil(double x)
 {
 #ifdef __IEEE754_CEIL
@@ -321,6 +351,12 @@ double fmod(double number, double denom)
 }
 
 
+float fmodf(float x, float y)
+{
+	return (float)fmod((double)x, (double)y);
+}
+
+
 double round(double x)
 {
 #ifdef __IEEE754_ROUND

math/complex.c renamed to libm/math/complex.c

@@ -41,6 +41,12 @@ double cosh(double x)
 }
 
 
+float coshf(float x)
+{
+	return (float)cosh((double)x);
+}
+
+
 double sinh(double x)
 {
 	double y, f;
@@ -72,6 +78,12 @@ double sinh(double x)
 }
 
 
+float sinhf(float x)
+{
+	return (double)sin((double)x);
+}
+
+
 double tanh(double x)
 {
 	if (isnan(x) != 0) {
@@ -88,3 +100,9 @@ double tanh(double x)
 	/* cosh is never equal to zero */
 	return (sinh(x) / cosh(x));
 }
+
+
+float tanhf(float x)
+{
+	return (float)tanh((double)x);
+}

math/exp.c renamed to libm/math/exp.c

@@ -18,7 +18,7 @@
 
 
 #include <stdarg.h>
-#include <math/consts.h>
+#include <consts.h>
 
 
 #ifdef __cplusplus
@@ -60,126 +60,126 @@ extern "C" {
 
 
 /* Returns the cosine of an angle of x radians. */
-extern double cos(double x);
+double cos(double x);
 
 
 /* Returns the sine of an angle of x radians. */
-extern double sin(double x);
+double sin(double x);
 
 
 /* Returns the tangent of an angle of x radians. */
-extern double tan(double x);
+double tan(double x);
 
 
 /* Returns the principal value of the arc cosine of x, expressed in radians. */
-extern double acos(double x);
+double acos(double x);
 
 
 /* Returns the principal value of the arc sine of x, expressed in radians. */
-extern double asin(double x);
+double asin(double x);
 
 
 /* Returns the principal value of the arc tan of x, expressed in radians. */
-extern double atan(double x);
+double atan(double x);
 
 
 /* Returns the principal value of the arc tangent of y/x, expressed in radians. */
-extern double atan2(double y, double x);
+double atan2(double y, double x);
 
 
 /* Hyperbolic functions */
 
 
 /* Returns the hyperbolic cosine of x. */
-extern double cosh(double x);
+double cosh(double x);
 
 
 /* Returns the hyperbolic sine of x. */
-extern double sinh(double x);
+double sinh(double x);
 
 
 /* Returns the hyperbolic tangent of x. */
-extern double tanh(double x);
+double tanh(double x);
 
 
 /* Exponential and logarithmic functions */
 
 
 /* Returns the base-e exponential function of x, which is e raised to the power x: e^x. */
-extern double exp(double x);
+double exp(double x);
 
 
 /* Breaks the floating point number x into its binary significand
  * (a floating point with an absolute value between 0.5(included) and 1.0(excluded))
  * and an integral exponent for 2. */
-extern double frexp(double x, int *exp);
+double frexp(double x, int *exp);
 
 
 /* Returns the result of multiplying x (the significand) by 2 raised to the power of exp (the exponent). */
-extern double ldexp(double x, int exp);
+double ldexp(double x, int exp);
 
 
 /* Returns the natural logarithm of x. */
-extern double log(double x);
+double log(double x);
 
 
 /* Returns the common (base-2) logarithm of x. */
-extern double log2(double x);
+double log2(double x);
 
 
 /* Returns the common (base-10) logarithm of x. */
-extern double log10(double x);
+double log10(double x);
 
 
 /* Breaks x into an integral and a fractional part. */
-extern double modf(double x, double *intpart);
-extern float modff(float x, float *intpart);
+double modf(double x, double *intpart);
+float modff(float x, float *intpart);
 
 
 /* Power functions */
 
 
 /* Returns base raised to the power exponent. */
-extern double pow(double base, double exponent);
+double pow(double base, double exponent);
 
 
 /* Returns the square root of x. */
-extern double sqrt(double x);
+double sqrt(double x);
 
 
 /* Rounding and remainder functions */
 
 
 /* Rounds x upward, returning the smallest integral value that is not less than x. */
-extern double ceil(double x);
-extern float ceilf(float x);
+double ceil(double x);
+float ceilf(float x);
 
 
 /* Rounds x downward, returning the largest integral value that is not greater than x. */
-extern double floor(double x);
-extern float floorf(float x);
+double floor(double x);
+float floorf(float x);
 
 
 /* Returns the floating-point remainder of numer/denom (rounded towards zero). */
-extern double fmod(double numer, double denom);
+double fmod(double numer, double denom);
 
 
 /* Return the integral value nearest to x */
-extern double round(double x);
-extern float roundf(float x);
+double round(double x);
+float roundf(float x);
 
 
 /* Rounds x toward zero, returning the nearest integral value that is not larger in magnitude than x. */
-extern double trunc(double x);
-extern float truncf(float x);
+double trunc(double x);
+float truncf(float x);
 
 
 /* Miscellaneous */
 
 
 /* Returns the absolute value of x: |x|. */
-extern double fabs(double x);
-extern float fabsf(float x);
+double fabs(double x);
+float fabsf(float x);
 
 
 /* C99 extensions */
@@ -203,25 +203,6 @@ float sqrtf(float x);
 float fmodf(float num, float denom);
 
 
-#define cosf(x)              ((float)cos(x))
-#define sinf(x)              ((float)sin(x))
-#define tanf(x)              ((float)tan(x))
-#define acosf(x)             ((float)acos(x))
-#define asinf(x)             ((float)asin(x))
-#define atanf(x)             ((float)atan(x))
-#define atan2f(y, x)         ((float)atan2(y, x))
-#define coshf(x)             ((float)cosh(x))
-#define sinhf(x)             ((float)sinh(x))
-#define tanhf(x)             ((float)tanh(x))
-#define expf(x)              ((float)exp(x))
-#define frexpf(x, exp)       ((float)frexp(x, exp))
-#define ldexpf(x, exp)       ((float)ldexp(x, exp))
-#define logf(x)              ((float)log(x))
-#define log10f(x)            ((float)log10(x))
-#define powf(base, exponent) ((float)pow(base, exponent))
-#define fmodf(num, denom)    ((float)fmod(num, denom))
-
-
 #ifdef __cplusplus
 }
 #endif