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 (file)
--- 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)
 {
 
 double fdim(double x, double y)
 {

diff --git a/src/math/fdimf.c b/src/math/fdimf.c
index 5cfeac6..543c364 100644 (file)
--- 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)
 {
 
 float fdimf(float x, float y)
 {

diff --git a/src/math/fdiml.c b/src/math/fdiml.c
index cda3022..62e29b7 100644 (file)
--- 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)
 
 #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 (file)
--- 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)
 {
 
 double fmax(double x, double y)
 {

diff --git a/src/math/fmaxf.c b/src/math/fmaxf.c
index 7767c30..695d817 100644 (file)
--- 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)
 {
 
 float fmaxf(float x, float y)
 {

diff --git a/src/math/fmaxl.c b/src/math/fmaxl.c
index 8a1e365..4b03158 100644 (file)
--- 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)
 
 #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 (file)
--- 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)
 {
 
 double fmin(double x, double y)
 {

diff --git a/src/math/fminf.c b/src/math/fminf.c
index 0964cdb..3573c7d 100644 (file)
--- 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)
 {
 
 float fminf(float x, float y)
 {

diff --git a/src/math/fminl.c b/src/math/fminl.c
index ae7159a..69bc24a 100644 (file)
--- 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)
 
 #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 (file)
--- a/src/math/ilogb.c
+++ b/src/math/ilogb.c
@@ -11,7 +11,6 @@ int ilogb(double x)
                if (u.bits == 0)
                        return FP_ILOGB0;
                /* subnormal x */
                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;
        }
                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 (file)
--- 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)
 {
 
 double ldexp(double x, int n)
 {

diff --git a/src/math/ldexpf.c b/src/math/ldexpf.c
index f0981ae..3bad5f3 100644 (file)
--- 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)
 {
 
 float ldexpf(float x, int n)
 {

diff --git a/src/math/ldexpl.c b/src/math/ldexpl.c
index 885ff6e..fd145cc 100644 (file)
--- 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)
 {
 
 long double ldexpl(long double x, int n)
 {

diff --git a/src/math/llrintl.c b/src/math/llrintl.c
index 6b0838d..ec2cf86 100644 (file)
--- 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)
 {
 #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 (file)
--- 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)
 {
 #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 (file)
--- 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)
 {
 #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 (file)
--- 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)
 {
 #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 (file)
--- a/src/math/nearbyint.c
+++ b/src/math/nearbyint.c
@@ -1,5 +1,5 @@
 #include <fenv.h>
 #include <fenv.h>

-#include "libm.h"
+#include <math.h>

 
 /*
 rint may raise inexact (and it should not alter the fenv otherwise)
 
 /*
 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 (file)
--- a/src/math/nearbyintf.c
+++ b/src/math/nearbyintf.c
@@ -1,5 +1,5 @@
 #include <fenv.h>
 #include <fenv.h>

-#include "libm.h"
+#include <math.h>

 
 float nearbyintf(float x) {
        fenv_t e;
 
 float nearbyintf(float x) {
        fenv_t e;

diff --git a/src/math/nearbyintl.c b/src/math/nearbyintl.c
index b58527c..2906f38 100644 (file)
--- 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)
 {
 #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 (file)
--- a/src/math/nexttoward.c
+++ b/src/math/nexttoward.c
@@ -11,6 +11,7 @@
  */
 
 #include "libm.h"
  */
 
 #include "libm.h"

+

 #if LDBL_MANT_DIG == 53 && LDBL_MAX_EXP == 1024
 double nexttoward(double x, long double y)
 {
 #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 (file)
--- 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)
 {
 
 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 (file)
--- 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)
 
 #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 (file)
--- 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.
  */
 
  * THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
  */
 

-#include "libm.h"
+#include <math.h>

 
 double round(double x)
 {
 
 double round(double x)
 {

diff --git a/src/math/roundf.c b/src/math/roundf.c
index 3cfd8ae..662a454 100644 (file)
--- 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.
  */
 
  * THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
  */
 

-#include "libm.h"
+#include <math.h>

 
 float roundf(float x)
 {
 
 float roundf(float x)
 {

diff --git a/src/math/roundl.c b/src/math/roundl.c
index ce56e8a..99146f0 100644 (file)
--- 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.
  */
 
  * 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)
 {
 #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 (file)
--- a/src/math/scalb.c
+++ b/src/math/scalb.c
@@ -15,7 +15,7 @@
  * should use scalbn() instead.
  */
 
  * should use scalbn() instead.
  */
 

-#include "libm.h"
+#include <math.h>

 
 double scalb(double x, double fn)
 {
 
 double scalb(double x, double fn)
 {

diff --git a/src/math/scalbf.c b/src/math/scalbf.c
index 0cc091f..d08f2a1 100644 (file)
--- 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)
 {
 
 float scalbf(float x, float fn)
 {

diff --git a/src/math/scalbln.c b/src/math/scalbln.c
index 53854fd..e6f3f19 100644 (file)
--- a/src/math/scalbln.c
+++ b/src/math/scalbln.c
@@ -1,5 +1,5 @@
 #include <limits.h>
 #include <limits.h>

-#include "libm.h"
+#include <math.h>

 
 double scalbln(double x, long n)
 {
 
 double scalbln(double x, long n)
 {

diff --git a/src/math/scalblnf.c b/src/math/scalblnf.c
index 61600f1..d8e8166 100644 (file)
--- a/src/math/scalblnf.c
+++ b/src/math/scalblnf.c
@@ -1,5 +1,5 @@
 #include <limits.h>
 #include <limits.h>

-#include "libm.h"
+#include <math.h>

 
 float scalblnf(float x, long n)
 {
 
 float scalblnf(float x, long n)
 {

diff --git a/src/math/scalblnl.c b/src/math/scalblnl.c
index 82ebbed..854c51c 100644 (file)
--- a/src/math/scalblnl.c
+++ b/src/math/scalblnl.c
@@ -1,5 +1,6 @@
 #include <limits.h>
 #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)
 
 #if LDBL_MANT_DIG == 53 && LDBL_MAX_EXP == 1024
 long double scalblnl(long double x, long n)