Include standard headers instead of using builtins

src/acb/lambertw.c

@@ -9,19 +9,9 @@
     (at your option) any later version.  See <https://www.gnu.org/licenses/>.
 */
 
+#include <math.h>
 #include "acb.h"
 
-#ifdef __GNUC__
-# define atan2 __builtin_atan2
-# define cos __builtin_cos
-# define fabs __builtin_fabs
-# define log __builtin_log
-# define pow __builtin_pow
-# define sin __builtin_sin
-#else
-# include <math.h>
-#endif
-
 /* Check if z crosses a branch cut. */
 int
 acb_lambertw_branch_crossing(const acb_t z, const acb_t ez1, const fmpz_t k)

src/acb/test/t-agm1.c

@@ -9,15 +9,10 @@
     (at your option) any later version.  See <https://www.gnu.org/licenses/>.
 */
 
+#include <math.h>
 #include "test_helpers.h"
 #include "acb.h"
 
-#ifdef __GNUC__
-# define fabs __builtin_fabs
-#else
-# include <math.h>
-#endif
-
 #define EPS 1e-13
 #define NUM_DERIVS 4
 #define NUM_TESTS 11

src/acb_calc/integrate_taylor.c

@@ -9,16 +9,11 @@
     (at your option) any later version.  See <https://www.gnu.org/licenses/>.
 */
 
+#include <math.h>
 #include "acb_poly.h"
 #include "arb_calc.h"
 #include "acb_calc.h"
 
-#ifdef __GNUC__
-# define log __builtin_log
-#else
-# include <math.h>
-#endif
-
 int
 acb_calc_integrate_taylor(acb_t res,
     acb_calc_func_t func, void * param,

src/acb_dirichlet/euler_product_real_ui.c

@@ -9,18 +9,12 @@
     (at your option) any later version.  See <https://www.gnu.org/licenses/>.
 */
 
+#include <math.h>
 #include "thread_pool.h"
 #include "thread_support.h"
 #include "arb.h"
 #include "acb_dirichlet.h"
 
-#ifdef __GNUC__
-# define log __builtin_log
-# define sqrt __builtin_sqrt
-#else
-# include <math.h>
-#endif
-
 #define ONE_OVER_LOG2 1.4426950408889634
 
 typedef struct

src/acb_dirichlet/hurwitz_precomp_choose_param.c

@@ -9,17 +9,10 @@
     (at your option) any later version.  See <https://www.gnu.org/licenses/>.
 */
 
+#include <math.h>
 #include "acb.h"
 #include "acb_dirichlet.h"
 
-#ifdef __GNUC__
-# define fabs __builtin_fabs
-# define log __builtin_log
-# define sqrt __builtin_sqrt
-#else
-# include <math.h>
-#endif
-
 void
 acb_dirichlet_hurwitz_precomp_choose_param(ulong * _A, ulong *_K, ulong *_N,
         const acb_t s, double num_eval, slong prec)

src/acb_dirichlet/l.c

@@ -10,15 +10,10 @@
     (at your option) any later version.  See <https://www.gnu.org/licenses/>.
 */
 
+#include <math.h>
 #include "acb.h"
 #include "acb_dirichlet.h"
 
-#ifdef __GNUC__
-# define log __builtin_log
-#else
-# include <math.h>
-#endif
-
 void
 acb_dirichlet_l_general(acb_t res, const acb_t s,
     const dirichlet_group_t G, const dirichlet_char_t chi, slong prec)

src/acb_dirichlet/l_euler_product.c

@@ -9,17 +9,10 @@
     (at your option) any later version.  See <https://www.gnu.org/licenses/>.
 */
 
+#include <math.h>
 #include "acb.h"
 #include "acb_dirichlet.h"
 
-#ifdef __GNUC__
-# define log __builtin_log
-# define pow __builtin_pow
-# define sqrt __builtin_sqrt
-#else
-# include <math.h>
-#endif
-
 #define ONE_OVER_LOG2 1.4426950408889634
 
 void

src/acb_dirichlet/l_fmpq.c

@@ -9,18 +9,13 @@
     (at your option) any later version.  See <https://www.gnu.org/licenses/>.
 */
 
+#include <math.h>
 #include "fmpq.h"
 #include "arb_mat.h"
 #include "arb_hypgeom.h"
 #include "acb.h"
 #include "acb_dirichlet.h"
 
-#ifdef __GNUC__
-# define fabs __builtin_fabs
-#else
-# include <math.h>
-#endif
-
 int
 acb_dirichlet_l_fmpq_use_afe(ulong q, const fmpq_t s, slong prec)
 {

src/acb_dirichlet/l_fmpq_afe.c

@@ -9,18 +9,13 @@
     (at your option) any later version.  See <https://www.gnu.org/licenses/>.
 */
 
+#include <math.h>
 #include "fmpq.h"
 #include "arb_mat.h"
 #include "arb_hypgeom.h"
 #include "acb.h"
 #include "acb_dirichlet.h"
 
-#ifdef __GNUC__
-# define log __builtin_log
-#else
-# include <math.h>
-#endif
-
 #define VERBOSE 0
 
 static double

src/acb_dirichlet/platt_local_hardy_z_zeros.c

@@ -11,16 +11,11 @@
     (at your option) any later version.  See <https://www.gnu.org/licenses/>.
 */
 
+#include <math.h>
 #include "acb.h"
 #include "acb_dirichlet.h"
 #include "arb_calc.h"
 
-#ifdef __GNUC__
-# define exp __builtin_exp
-#else
-# include <math.h>
-#endif
-
 static const slong LOOPCOUNT = 4;
 
 /*

src/acb_dirichlet/stieltjes.c

@@ -9,22 +9,11 @@
     (at your option) any later version.  See <https://www.gnu.org/licenses/>.
 */
 
+#include <math.h>
 #include "acb_poly.h"
 #include "acb_dirichlet.h"
 #include "acb_calc.h"
 
-#ifdef __GNUC__
-# define atan2 __builtin_atan2
-# define cos __builtin_cos
-# define exp __builtin_exp
-# define log __builtin_log
-# define pow __builtin_pow
-# define sqrt __builtin_sqrt
-# define tan __builtin_tan
-#else
-# include <math.h>
-#endif
-
 /* Bound the quadratic Taylor error term. */
 static void
 stieltjes_bound_quadratic_term(arb_t B, const acb_t z,

src/acb_dirichlet/test/t-euler_product_real_ui.c

@@ -9,15 +9,10 @@
     (at your option) any later version.  See <https://www.gnu.org/licenses/>.
 */
 
+#include <math.h>
 #include "test_helpers.h"
 #include "acb_dirichlet.h"
 
-#ifdef __GNUC__
-# define log __builtin_log
-#else
-# include <math.h>
-#endif
-
 const signed char chi[8][6] = {
   {1, 1},
   {2, 0, 1},

src/acb_dirichlet/test/t-l_jet.c

@@ -9,15 +9,10 @@
     (at your option) any later version.  See <https://www.gnu.org/licenses/>.
 */
 
+#include <math.h>
 #include "test_helpers.h"
 #include "acb_dirichlet.h"
 
-#ifdef __GNUC__
-# define fabs __builtin_fabs
-#else
-# include <math.h>
-#endif
-
 /* Laurent expansions at s = 1 of first 10 principal L-functions */
 
 /* with mpmath:

src/acb_dirichlet/zeta.c

@@ -9,15 +9,10 @@
     (at your option) any later version.  See <https://www.gnu.org/licenses/>.
 */
 
+#include <math.h>
 #include "acb_poly.h"
 #include "acb_dirichlet.h"
 
-#ifdef __GNUC__
-# define sqrt __builtin_sqrt
-#else
-# include <math.h>
-#endif
-
 void acb_zeta_si(acb_t z, slong s, slong prec);
 
 static void

src/acb_dirichlet/zeta_jet.c

@@ -9,15 +9,10 @@
     (at your option) any later version.  See <https://www.gnu.org/licenses/>.
 */
 
+#include <math.h>
 #include "acb_poly.h"
 #include "acb_dirichlet.h"
 
-#ifdef __GNUC__
-# define sqrt __builtin_sqrt
-#else
-# include <math.h>
-#endif
-
 void
 _acb_dirichlet_zeta_jet(acb_t t, const acb_t h, int deflate, slong len, slong prec)
 {

src/acb_dirichlet/zeta_rs_r.c

@@ -9,16 +9,10 @@
     (at your option) any later version.  See <https://www.gnu.org/licenses/>.
 */
 
+#include <math.h>
 #include "acb.h"
 #include "acb_dirichlet.h"
 
-#ifdef __GNUC__
-# define log __builtin_log
-# define pow __builtin_pow
-#else
-# include <math.h>
-#endif
-
 void
 acb_dirichlet_zeta_rs_r(acb_t res, const acb_t s, slong K, slong prec)
 {

src/acb_elliptic/rf.c

@@ -9,15 +9,10 @@
     (at your option) any later version.  See <https://www.gnu.org/licenses/>.
 */
 
+#include <math.h>
 #include "acb.h"
 #include "acb_elliptic.h"
 
-#ifdef __GNUC__
-# define pow __builtin_pow
-#else
-# include <math.h>
-#endif
-
 static const unsigned short den_ratio_tab[512] = {
     1,1,10,7,12,11,26,1,136,19,2,23,20,1,58,31,
     16,1,74,1,164,43,2,47,56,1,106,1,4,59,122,1,

src/acb_elliptic/rj.c

@@ -9,16 +9,11 @@
     (at your option) any later version.  See <https://www.gnu.org/licenses/>.
 */
 
+#include <math.h>
 #include "acb.h"
 #include "acb_elliptic.h"
 #include "acb_calc.h"
 
-#ifdef __GNUC__
-# define pow __builtin_pow
-#else
-# include <math.h>
-#endif
-
 static const unsigned short den_ratio_tab[512] = {
     1,1,14,3,44,13,10,17,152,1,46,1,12,29,62,1,
     16,37,2,41,172,1,94,7,8,53,2,1,236,61,2,1,

src/acb_elliptic/test/t-pi.c

@@ -9,17 +9,12 @@
     (at your option) any later version.  See <https://www.gnu.org/licenses/>.
 */
 
+#include <math.h>
 #include "test_helpers.h"
 #include "acb.h"
 #include "acb_elliptic.h"
 #include "acb_modular.h"
 
-#ifdef __GNUC__
-# define fabs __builtin_fabs
-#else
-# include <math.h>
-#endif
-
 static const double testdata_pi[17][6] = {
   {-2.0, 0.0, 0.0, 0.0, 0.9068996821171089253, 0.0},
   {-2.0, 0.0, 2.0, 0.0, 0.90023699434287315449, -0.5446596026314193538},

src/acb_elliptic/test/t-rg.c

@@ -9,16 +9,11 @@
     (at your option) any later version.  See <https://www.gnu.org/licenses/>.
 */
 
+#include <math.h>
 #include "test_helpers.h"
 #include "acb.h"
 #include "acb_elliptic.h"
 
-#ifdef __GNUC__
-# define fabs __builtin_fabs
-#else
-# include <math.h>
-#endif
-
 /* Test input from Carlson's paper and checked with mpmath. */
 
 static const double testdata_rg[7][8] = {

src/acb_elliptic/test/t-rj.c

@@ -9,16 +9,11 @@
     (at your option) any later version.  See <https://www.gnu.org/licenses/>.
 */
 
+#include <math.h>
 #include "test_helpers.h"
 #include "acb.h"
 #include "acb_elliptic.h"
 
-#ifdef __GNUC__
-# define fabs __builtin_fabs
-#else
-# include <math.h>
-#endif
-
 /* Test input from Carlson's paper and checked with mpmath. */
 
 static const double testdata_rj[16][10] = {