remove libm.h includes when math.h and float.h are enough
diff --git a/src/math/fdim.c b/src/math/fdim.c
index fb25521..9585460 100644
--- a/src/math/fdim.c
+++ b/src/math/fdim.c
@@ -1,4 +1,4 @@
-#include "libm.h"
+#include <math.h>
double fdim(double x, double y)
{
diff --git a/src/math/fdimf.c b/src/math/fdimf.c
index 5cfeac6..543c364 100644
--- a/src/math/fdimf.c
+++ b/src/math/fdimf.c
@@ -1,4 +1,4 @@
-#include "libm.h"
+#include <math.h>
float fdimf(float x, float y)
{
diff --git a/src/math/fdiml.c b/src/math/fdiml.c
index cda3022..62e29b7 100644
--- a/src/math/fdiml.c
+++ b/src/math/fdiml.c
@@ -1,4 +1,5 @@
-#include "libm.h"
+#include <math.h>
+#include <float.h>
#if LDBL_MANT_DIG == 53 && LDBL_MAX_EXP == 1024
long double fdiml(long double x, long double y)
diff --git a/src/math/fmax.c b/src/math/fmax.c
index 0b6bf6f..94f0caa 100644
--- a/src/math/fmax.c
+++ b/src/math/fmax.c
@@ -1,4 +1,4 @@
-#include "libm.h"
+#include <math.h>
double fmax(double x, double y)
{
diff --git a/src/math/fmaxf.c b/src/math/fmaxf.c
index 7767c30..695d817 100644
--- a/src/math/fmaxf.c
+++ b/src/math/fmaxf.c
@@ -1,4 +1,4 @@
-#include "libm.h"
+#include <math.h>
float fmaxf(float x, float y)
{
diff --git a/src/math/fmaxl.c b/src/math/fmaxl.c
index 8a1e365..4b03158 100644
--- a/src/math/fmaxl.c
+++ b/src/math/fmaxl.c
@@ -1,4 +1,5 @@
-#include "libm.h"
+#include <math.h>
+#include <float.h>
#if LDBL_MANT_DIG == 53 && LDBL_MAX_EXP == 1024
long double fmaxl(long double x, long double y)
diff --git a/src/math/fmin.c b/src/math/fmin.c
index d1f1645..08a8fd1 100644
--- a/src/math/fmin.c
+++ b/src/math/fmin.c
@@ -1,4 +1,4 @@
-#include "libm.h"
+#include <math.h>
double fmin(double x, double y)
{
diff --git a/src/math/fminf.c b/src/math/fminf.c
index 0964cdb..3573c7d 100644
--- a/src/math/fminf.c
+++ b/src/math/fminf.c
@@ -1,4 +1,4 @@
-#include "libm.h"
+#include <math.h>
float fminf(float x, float y)
{
diff --git a/src/math/fminl.c b/src/math/fminl.c
index ae7159a..69bc24a 100644
--- a/src/math/fminl.c
+++ b/src/math/fminl.c
@@ -1,4 +1,5 @@
-#include "libm.h"
+#include <math.h>
+#include <float.h>
#if LDBL_MANT_DIG == 53 && LDBL_MAX_EXP == 1024
long double fminl(long double x, long double y)
diff --git a/src/math/ilogb.c b/src/math/ilogb.c
index c5915a0..0a3a6a4 100644
--- a/src/math/ilogb.c
+++ b/src/math/ilogb.c
@@ -11,7 +11,6 @@
if (u.bits == 0)
return FP_ILOGB0;
/* subnormal x */
- // FIXME: scale up subnormals with a *0x1p53 or find top set bit with a better method
for (e = -0x3ff; u.bits < (uint64_t)1<<63; e--, u.bits<<=1);
return e;
}
diff --git a/src/math/ldexp.c b/src/math/ldexp.c
index 36835db..f4d1cd6 100644
--- a/src/math/ldexp.c
+++ b/src/math/ldexp.c
@@ -1,4 +1,4 @@
-#include "libm.h"
+#include <math.h>
double ldexp(double x, int n)
{
diff --git a/src/math/ldexpf.c b/src/math/ldexpf.c
index f0981ae..3bad5f3 100644
--- a/src/math/ldexpf.c
+++ b/src/math/ldexpf.c
@@ -1,4 +1,4 @@
-#include "libm.h"
+#include <math.h>
float ldexpf(float x, int n)
{
diff --git a/src/math/ldexpl.c b/src/math/ldexpl.c
index 885ff6e..fd145cc 100644
--- a/src/math/ldexpl.c
+++ b/src/math/ldexpl.c
@@ -1,4 +1,4 @@
-#include "libm.h"
+#include <math.h>
long double ldexpl(long double x, int n)
{
diff --git a/src/math/llrintl.c b/src/math/llrintl.c
index 6b0838d..ec2cf86 100644
--- a/src/math/llrintl.c
+++ b/src/math/llrintl.c
@@ -1,4 +1,6 @@
-#include "libm.h"
+#include <math.h>
+#include <float.h>
+
#if LDBL_MANT_DIG == 53 && LDBL_MAX_EXP == 1024
long long llrintl(long double x)
{
diff --git a/src/math/llroundl.c b/src/math/llroundl.c
index 9e2cfdc..107744a 100644
--- a/src/math/llroundl.c
+++ b/src/math/llroundl.c
@@ -1,4 +1,6 @@
-#include "libm.h"
+#include <math.h>
+#include <float.h>
+
#if LDBL_MANT_DIG == 53 && LDBL_MAX_EXP == 1024
long long llroundl(long double x)
{
diff --git a/src/math/lrintl.c b/src/math/lrintl.c
index 7c09653..c076326 100644
--- a/src/math/lrintl.c
+++ b/src/math/lrintl.c
@@ -1,4 +1,6 @@
-#include "libm.h"
+#include <math.h>
+#include <float.h>
+
#if LDBL_MANT_DIG == 53 && LDBL_MAX_EXP == 1024
long lrintl(long double x)
{
diff --git a/src/math/lroundl.c b/src/math/lroundl.c
index 1469127..7b593f7 100644
--- a/src/math/lroundl.c
+++ b/src/math/lroundl.c
@@ -1,4 +1,6 @@
-#include "libm.h"
+#include <math.h>
+#include <float.h>
+
#if LDBL_MANT_DIG == 53 && LDBL_MAX_EXP == 1024
long lroundl(long double x)
{
diff --git a/src/math/nearbyint.c b/src/math/nearbyint.c
index 781769f..714c55c 100644
--- a/src/math/nearbyint.c
+++ b/src/math/nearbyint.c
@@ -1,5 +1,5 @@
#include <fenv.h>
-#include "libm.h"
+#include <math.h>
/*
rint may raise inexact (and it should not alter the fenv otherwise)
diff --git a/src/math/nearbyintf.c b/src/math/nearbyintf.c
index e4bdb26..07df8f5 100644
--- a/src/math/nearbyintf.c
+++ b/src/math/nearbyintf.c
@@ -1,5 +1,5 @@
#include <fenv.h>
-#include "libm.h"
+#include <math.h>
float nearbyintf(float x) {
fenv_t e;
diff --git a/src/math/nearbyintl.c b/src/math/nearbyintl.c
index b58527c..2906f38 100644
--- a/src/math/nearbyintl.c
+++ b/src/math/nearbyintl.c
@@ -1,4 +1,6 @@
-#include "libm.h"
+#include <math.h>
+#include <float.h>
+
#if LDBL_MANT_DIG == 53 && LDBL_MAX_EXP == 1024
long double nearbyintl(long double x)
{
diff --git a/src/math/nexttoward.c b/src/math/nexttoward.c
index 5e12c48..e150eaa 100644
--- a/src/math/nexttoward.c
+++ b/src/math/nexttoward.c
@@ -11,6 +11,7 @@
*/
#include "libm.h"
+
#if LDBL_MANT_DIG == 53 && LDBL_MAX_EXP == 1024
double nexttoward(double x, long double y)
{
diff --git a/src/math/nexttowardl.c b/src/math/nexttowardl.c
index c393ce9..67a6340 100644
--- a/src/math/nexttowardl.c
+++ b/src/math/nexttowardl.c
@@ -1,4 +1,4 @@
-#include "libm.h"
+#include <math.h>
long double nexttowardl(long double x, long double y)
{
diff --git a/src/math/remainderl.c b/src/math/remainderl.c
index b99f938..2a13c1d 100644
--- a/src/math/remainderl.c
+++ b/src/math/remainderl.c
@@ -1,4 +1,5 @@
-#include "libm.h"
+#include <math.h>
+#include <float.h>
#if LDBL_MANT_DIG == 53 && LDBL_MAX_EXP == 1024
long double remainderl(long double x, long double y)
diff --git a/src/math/round.c b/src/math/round.c
index 2137384..bf0b453 100644
--- a/src/math/round.c
+++ b/src/math/round.c
@@ -25,7 +25,7 @@
* THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
*/
-#include "libm.h"
+#include <math.h>
double round(double x)
{
diff --git a/src/math/roundf.c b/src/math/roundf.c
index 3cfd8ae..662a454 100644
--- a/src/math/roundf.c
+++ b/src/math/roundf.c
@@ -25,7 +25,7 @@
* THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
*/
-#include "libm.h"
+#include <math.h>
float roundf(float x)
{
diff --git a/src/math/roundl.c b/src/math/roundl.c
index ce56e8a..99146f0 100644
--- a/src/math/roundl.c
+++ b/src/math/roundl.c
@@ -25,7 +25,9 @@
* THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
*/
-#include "libm.h"
+#include <math.h>
+#include <float.h>
+
#if LDBL_MANT_DIG == 53 && LDBL_MAX_EXP == 1024
long double roundl(long double x)
{
diff --git a/src/math/scalb.c b/src/math/scalb.c
index 7706e9c..a54bdf2 100644
--- a/src/math/scalb.c
+++ b/src/math/scalb.c
@@ -15,7 +15,7 @@
* should use scalbn() instead.
*/
-#include "libm.h"
+#include <math.h>
double scalb(double x, double fn)
{
diff --git a/src/math/scalbf.c b/src/math/scalbf.c
index 0cc091f..d08f2a1 100644
--- a/src/math/scalbf.c
+++ b/src/math/scalbf.c
@@ -13,7 +13,7 @@
* ====================================================
*/
-#include "libm.h"
+#include <math.h>
float scalbf(float x, float fn)
{
diff --git a/src/math/scalbln.c b/src/math/scalbln.c
index 53854fd..e6f3f19 100644
--- a/src/math/scalbln.c
+++ b/src/math/scalbln.c
@@ -1,5 +1,5 @@
#include <limits.h>
-#include "libm.h"
+#include <math.h>
double scalbln(double x, long n)
{
diff --git a/src/math/scalblnf.c b/src/math/scalblnf.c
index 61600f1..d8e8166 100644
--- a/src/math/scalblnf.c
+++ b/src/math/scalblnf.c
@@ -1,5 +1,5 @@
#include <limits.h>
-#include "libm.h"
+#include <math.h>
float scalblnf(float x, long n)
{
diff --git a/src/math/scalblnl.c b/src/math/scalblnl.c
index 82ebbed..854c51c 100644
--- a/src/math/scalblnl.c
+++ b/src/math/scalblnl.c
@@ -1,5 +1,6 @@
#include <limits.h>
-#include "libm.h"
+#include <math.h>
+#include <float.h>
#if LDBL_MANT_DIG == 53 && LDBL_MAX_EXP == 1024
long double scalblnl(long double x, long n)
