diff options
-rw-r--r-- | src/atan.c | 2 | ||||
-rw-r--r-- | src/const_euler.c | 2 | ||||
-rw-r--r-- | src/eint.c | 2 | ||||
-rw-r--r-- | src/erandom.c | 2 | ||||
-rw-r--r-- | src/fpif.c | 13 | ||||
-rw-r--r-- | src/li2.c | 2 | ||||
-rw-r--r-- | src/lngamma.c | 2 | ||||
-rw-r--r-- | src/mpfr-impl.h | 3 | ||||
-rw-r--r-- | src/mpfr.h | 14 | ||||
-rw-r--r-- | src/nrandom.c | 2 | ||||
-rw-r--r-- | src/random_deviate.c | 2 | ||||
-rw-r--r-- | src/random_deviate.h | 2 | ||||
-rw-r--r-- | src/rndna.c | 4 | ||||
-rw-r--r-- | src/set_sj.c | 4 | ||||
-rw-r--r-- | src/set_str.c | 2 | ||||
-rw-r--r-- | src/set_uj.c | 4 | ||||
-rw-r--r-- | src/sin_cos.c | 4 | ||||
-rw-r--r-- | src/strtofr.c | 6 | ||||
-rw-r--r-- | src/vasprintf.c | 2 | ||||
-rw-r--r-- | src/zeta.c | 12 | ||||
-rw-r--r-- | tests/tcmp2.c | 2 | ||||
-rw-r--r-- | tests/tdiv.c | 2 | ||||
-rw-r--r-- | tests/teq.c | 2 | ||||
-rw-r--r-- | tests/terandom_chisq.c | 6 | ||||
-rw-r--r-- | tests/tfmma.c | 2 | ||||
-rw-r--r-- | tests/tfmod.c | 6 | ||||
-rw-r--r-- | tests/tfprintf.c | 2 | ||||
-rw-r--r-- | tests/tgamma.c | 2 | ||||
-rw-r--r-- | tests/tnrandom_chisq.c | 6 | ||||
-rw-r--r-- | tests/tprintf.c | 2 | ||||
-rw-r--r-- | tests/trint.c | 2 | ||||
-rw-r--r-- | tests/tsinh_cosh.c | 4 | ||||
-rw-r--r-- | tests/tsqr.c | 4 | ||||
-rw-r--r-- | tests/tsum.c | 4 |
34 files changed, 68 insertions, 64 deletions
diff --git a/src/atan.c b/src/atan.c index 0f63416c9..9a1459081 100644 --- a/src/atan.c +++ b/src/atan.c @@ -56,7 +56,7 @@ static const mp_limb_t atan_table[][3] = { }; static void -set_table (mpfr_t y, const mp_limb_t x[3]) +set_table (mpfr_ptr y, const mp_limb_t x[3]) { mpfr_prec_t p = MPFR_PREC(y); mp_size_t n = MPFR_PREC2LIMBS(p); diff --git a/src/const_euler.c b/src/const_euler.c index 9e3b88ac5..4672e2d6e 100644 --- a/src/const_euler.c +++ b/src/const_euler.c @@ -181,7 +181,7 @@ mpfr_const_euler_bs_2 (mpz_t P, mpz_t Q, mpz_t T, } int -mpfr_const_euler_internal (mpfr_t x, mpfr_rnd_t rnd) +mpfr_const_euler_internal (mpfr_ptr x, mpfr_rnd_t rnd) { mpfr_const_euler_bs_t sum; mpz_t t, u, v; diff --git a/src/eint.c b/src/eint.c index 4df65821d..edd852fce 100644 --- a/src/eint.c +++ b/src/eint.c @@ -36,7 +36,7 @@ https://www.gnu.org/licenses/ or write to the Free Software Foundation, Inc., Return PREC(y) when the truncated series does not converge. */ static mpfr_exp_t -mpfr_eint_aux (mpfr_t y, mpfr_srcptr x) +mpfr_eint_aux (mpfr_ptr y, mpfr_srcptr x) { mpfr_t eps; /* dynamic (absolute) error bound on t */ mpfr_t erru, errs; diff --git a/src/erandom.c b/src/erandom.c index 86c9c4ddc..d9249d29f 100644 --- a/src/erandom.c +++ b/src/erandom.c @@ -80,7 +80,7 @@ E (mpfr_random_deviate_t x, gmp_randstate_t r, /* return an exponential random deviate with mean 1 as a MPFR */ int -mpfr_erandom (mpfr_t z, gmp_randstate_t r, mpfr_rnd_t rnd) +mpfr_erandom (mpfr_ptr z, gmp_randstate_t r, mpfr_rnd_t rnd) { mpfr_random_deviate_t x, p, q; int inex; diff --git a/src/fpif.c b/src/fpif.c index df775dbd7..67f98b5aa 100644 --- a/src/fpif.c +++ b/src/fpif.c @@ -291,7 +291,8 @@ mpfr_fpif_read_precision_from_file (FILE *fh) * until one has integer types larger than 128 bits). */ static unsigned char* -mpfr_fpif_store_exponent (unsigned char *buffer, size_t *buffer_size, mpfr_t x) +mpfr_fpif_store_exponent (unsigned char *buffer, size_t *buffer_size, + mpfr_ptr x) { unsigned char *result; mpfr_uexp_t uexp; @@ -372,7 +373,7 @@ mpfr_fpif_store_exponent (unsigned char *buffer, size_t *buffer_size, mpfr_t x) * than 128 bits). */ static int -mpfr_fpif_read_exponent_from_file (mpfr_t x, FILE * fh) +mpfr_fpif_read_exponent_from_file (mpfr_ptr x, FILE * fh) { mpfr_exp_t exponent; mpfr_uexp_t uexp; @@ -456,7 +457,7 @@ mpfr_fpif_read_exponent_from_file (mpfr_t x, FILE * fh) * format */ static unsigned char* -mpfr_fpif_store_limbs (unsigned char *buffer, size_t *buffer_size, mpfr_t x) +mpfr_fpif_store_limbs (unsigned char *buffer, size_t *buffer_size, mpfr_ptr x) { unsigned char *result; mpfr_prec_t precision; @@ -492,7 +493,7 @@ mpfr_fpif_store_limbs (unsigned char *buffer, size_t *buffer_size, mpfr_t x) * Assume buffer is not NULL. */ static void -mpfr_fpif_read_limbs (mpfr_t x, unsigned char *buffer, size_t nb_byte) +mpfr_fpif_read_limbs (mpfr_ptr x, unsigned char *buffer, size_t nb_byte) { size_t mp_bytes_per_limb; size_t nb_partial_byte; @@ -522,7 +523,7 @@ mpfr_fpif_read_limbs (mpfr_t x, unsigned char *buffer, size_t nb_byte) * return 0 if successful */ int -mpfr_fpif_export (FILE *fh, mpfr_t x) +mpfr_fpif_export (FILE *fh, mpfr_ptr x) { int status; unsigned char *buf; @@ -582,7 +583,7 @@ mpfr_fpif_export (FILE *fh, mpfr_t x) * Return 0 if the import was successful. */ int -mpfr_fpif_import (mpfr_t x, FILE *fh) +mpfr_fpif_import (mpfr_ptr x, FILE *fh) { int status; mpfr_prec_t precision; @@ -31,7 +31,7 @@ https://www.gnu.org/licenses/ or write to the Free Software Foundation, Inc., for determinating the relative error. */ static int -li2_series (mpfr_t sum, mpfr_srcptr z, mpfr_rnd_t rnd_mode) +li2_series (mpfr_ptr sum, mpfr_srcptr z, mpfr_rnd_t rnd_mode) { int i; mpfr_t s, u, v, w; diff --git a/src/lngamma.c b/src/lngamma.c index 8edc2ddc1..029af42c7 100644 --- a/src/lngamma.c +++ b/src/lngamma.c @@ -31,7 +31,7 @@ https://www.gnu.org/licenses/ or write to the Free Software Foundation, Inc., precision should be >= 4. */ static void -mpfr_gamma_alpha (mpfr_t s, mpfr_prec_t p) +mpfr_gamma_alpha (mpfr_ptr s, mpfr_prec_t p) { MPFR_LOG_FUNC (("p=%Pu", p), diff --git a/src/mpfr-impl.h b/src/mpfr-impl.h index c9c792a21..90bc2d238 100644 --- a/src/mpfr-impl.h +++ b/src/mpfr-impl.h @@ -2474,7 +2474,8 @@ __MPFR_DECLSPEC void mpfr_rand_raw (mpfr_limb_ptr, gmp_randstate_t, __MPFR_DECLSPEC mpz_srcptr mpfr_bernoulli_cache (unsigned long); __MPFR_DECLSPEC void mpfr_bernoulli_freecache (void); -__MPFR_DECLSPEC int mpfr_sincos_fast (mpfr_t, mpfr_t, mpfr_srcptr, mpfr_rnd_t); +__MPFR_DECLSPEC int mpfr_sincos_fast (mpfr_ptr, mpfr_ptr, mpfr_srcptr, + mpfr_rnd_t); __MPFR_DECLSPEC double mpfr_scale2 (double, int); diff --git a/src/mpfr.h b/src/mpfr.h index ef5e8ee42..0d6643e33 100644 --- a/src/mpfr.h +++ b/src/mpfr.h @@ -781,8 +781,8 @@ __MPFR_DECLSPEC int mpfr_subnormalize (mpfr_ptr, int, mpfr_rnd_t); __MPFR_DECLSPEC int mpfr_strtofr (mpfr_ptr, const char *, char **, int, mpfr_rnd_t); -__MPFR_DECLSPEC void mpfr_round_nearest_away_begin (mpfr_t); -__MPFR_DECLSPEC int mpfr_round_nearest_away_end (mpfr_t, int); +__MPFR_DECLSPEC void mpfr_round_nearest_away_begin (mpfr_ptr); +__MPFR_DECLSPEC int mpfr_round_nearest_away_end (mpfr_ptr, int); __MPFR_DECLSPEC size_t mpfr_custom_get_size (mpfr_prec_t); __MPFR_DECLSPEC void mpfr_custom_init (void *, mpfr_prec_t); @@ -1080,10 +1080,12 @@ extern "C" { #define mpfr_set_uj_2exp __gmpfr_set_uj_2exp #define mpfr_get_sj __gmpfr_mpfr_get_sj #define mpfr_get_uj __gmpfr_mpfr_get_uj -__MPFR_DECLSPEC int mpfr_set_sj (mpfr_t, intmax_t, mpfr_rnd_t); -__MPFR_DECLSPEC int mpfr_set_sj_2exp (mpfr_t, intmax_t, intmax_t, mpfr_rnd_t); -__MPFR_DECLSPEC int mpfr_set_uj (mpfr_t, uintmax_t, mpfr_rnd_t); -__MPFR_DECLSPEC int mpfr_set_uj_2exp (mpfr_t, uintmax_t, intmax_t, mpfr_rnd_t); +__MPFR_DECLSPEC int mpfr_set_sj (mpfr_ptr, intmax_t, mpfr_rnd_t); +__MPFR_DECLSPEC int mpfr_set_sj_2exp (mpfr_ptr, intmax_t, intmax_t, + mpfr_rnd_t); +__MPFR_DECLSPEC int mpfr_set_uj (mpfr_ptr, uintmax_t, mpfr_rnd_t); +__MPFR_DECLSPEC int mpfr_set_uj_2exp (mpfr_ptr, uintmax_t, intmax_t, + mpfr_rnd_t); __MPFR_DECLSPEC intmax_t mpfr_get_sj (mpfr_srcptr, mpfr_rnd_t); __MPFR_DECLSPEC uintmax_t mpfr_get_uj (mpfr_srcptr, mpfr_rnd_t); diff --git a/src/nrandom.c b/src/nrandom.c index b081c87ad..990987ee6 100644 --- a/src/nrandom.c +++ b/src/nrandom.c @@ -155,7 +155,7 @@ B (unsigned long k, mpfr_random_deviate_t x, gmp_randstate_t r, /* return a normal random deviate with mean 0 and variance 1 as a MPFR */ int -mpfr_nrandom (mpfr_t z, gmp_randstate_t r, mpfr_rnd_t rnd) +mpfr_nrandom (mpfr_ptr z, gmp_randstate_t r, mpfr_rnd_t rnd) { mpfr_random_deviate_t x, p, q; int inex; diff --git a/src/random_deviate.c b/src/random_deviate.c index 72531e774..3d875c6ea 100644 --- a/src/random_deviate.c +++ b/src/random_deviate.c @@ -280,7 +280,7 @@ mpfr_random_deviate_less (mpfr_random_deviate_t x, mpfr_random_deviate_t y, /* set mpfr_t z = (neg ? -1 : 1) * (n + x) */ int mpfr_random_deviate_value (int neg, unsigned long n, - mpfr_random_deviate_t x, mpfr_t z, + mpfr_random_deviate_t x, mpfr_ptr z, gmp_randstate_t r, mpfr_rnd_t rnd) { /* r is used to add as many bits as necessary to match the precision of z */ diff --git a/src/random_deviate.h b/src/random_deviate.h index fe9246236..faa196eee 100644 --- a/src/random_deviate.h +++ b/src/random_deviate.h @@ -76,7 +76,7 @@ __MPFR_DECLSPEC int /* set mpfr_t z = (neg ? -1 : 1) * (n + x) */ __MPFR_DECLSPEC int mpfr_random_deviate_value (int, unsigned long, - mpfr_random_deviate_ptr, mpfr_t, + mpfr_random_deviate_ptr, mpfr_ptr, gmp_randstate_t, mpfr_rnd_t); #if defined(__cplusplus) diff --git a/src/rndna.c b/src/rndna.c index 23997bf0e..b76517739 100644 --- a/src/rndna.c +++ b/src/rndna.c @@ -61,7 +61,7 @@ typedef enum { and prepares rop to give it one more bit of precision and to save its old value within it. */ void -mpfr_round_nearest_away_begin (mpfr_t rop) +mpfr_round_nearest_away_begin (mpfr_ptr rop) { mpfr_t tmp; mp_size_t xsize; @@ -129,7 +129,7 @@ mpfr_round_nearest_away_begin (mpfr_t rop) copying it back the result of the applied function and performing additional roundings. */ int -mpfr_round_nearest_away_end (mpfr_t rop, int inex) +mpfr_round_nearest_away_end (mpfr_ptr rop, int inex) { mpfr_t tmp; mp_size_t xsize; diff --git a/src/set_sj.c b/src/set_sj.c index 488125076..a4454928d 100644 --- a/src/set_sj.c +++ b/src/set_sj.c @@ -26,13 +26,13 @@ https://www.gnu.org/licenses/ or write to the Free Software Foundation, Inc., #ifdef _MPFR_H_HAVE_INTMAX_T int -mpfr_set_sj (mpfr_t x, intmax_t j, mpfr_rnd_t rnd) +mpfr_set_sj (mpfr_ptr x, intmax_t j, mpfr_rnd_t rnd) { return mpfr_set_sj_2exp (x, j, 0, rnd); } int -mpfr_set_sj_2exp (mpfr_t x, intmax_t j, intmax_t e, mpfr_rnd_t rnd) +mpfr_set_sj_2exp (mpfr_ptr x, intmax_t j, intmax_t e, mpfr_rnd_t rnd) { if (j >= 0) return mpfr_set_uj_2exp (x, j, e, rnd); diff --git a/src/set_str.c b/src/set_str.c index 408051dbb..7876fdaf0 100644 --- a/src/set_str.c +++ b/src/set_str.c @@ -23,7 +23,7 @@ https://www.gnu.org/licenses/ or write to the Free Software Foundation, Inc., #include "mpfr-impl.h" int -mpfr_set_str (mpfr_t x, const char *str, int base, mpfr_rnd_t rnd) +mpfr_set_str (mpfr_ptr x, const char *str, int base, mpfr_rnd_t rnd) { char *p; diff --git a/src/set_uj.c b/src/set_uj.c index 4dab72e19..d69d3d391 100644 --- a/src/set_uj.c +++ b/src/set_uj.c @@ -29,13 +29,13 @@ https://www.gnu.org/licenses/ or write to the Free Software Foundation, Inc., #define uintmaxpml (sizeof(uintmax_t) / sizeof(mp_limb_t)) int -mpfr_set_uj (mpfr_t x, uintmax_t j, mpfr_rnd_t rnd) +mpfr_set_uj (mpfr_ptr x, uintmax_t j, mpfr_rnd_t rnd) { return mpfr_set_uj_2exp (x, j, 0, rnd); } int -mpfr_set_uj_2exp (mpfr_t x, uintmax_t j, intmax_t e, mpfr_rnd_t rnd) +mpfr_set_uj_2exp (mpfr_ptr x, uintmax_t j, intmax_t e, mpfr_rnd_t rnd) { int cnt, inex; mp_size_t i, k; diff --git a/src/sin_cos.c b/src/sin_cos.c index daf3040ee..00bc18858 100644 --- a/src/sin_cos.c +++ b/src/sin_cos.c @@ -463,7 +463,7 @@ sin_bs_aux (mpz_t Q0, mpz_t S0, mpz_t C0, mpz_srcptr p, mpfr_prec_t r, Return err such that the relative error is bounded by 2^err ulps. */ static int -sincos_aux (mpfr_t s, mpfr_t c, mpfr_srcptr x, mpfr_rnd_t rnd_mode) +sincos_aux (mpfr_ptr s, mpfr_ptr c, mpfr_srcptr x, mpfr_rnd_t rnd_mode) { mpfr_prec_t prec_s, sh; mpz_t Q, S, C, Q2, S2, C2, y; @@ -577,7 +577,7 @@ sincos_aux (mpfr_t s, mpfr_t c, mpfr_srcptr x, mpfr_rnd_t rnd_mode) Assumes s differs from c. */ int -mpfr_sincos_fast (mpfr_t s, mpfr_t c, mpfr_srcptr x, mpfr_rnd_t rnd) +mpfr_sincos_fast (mpfr_ptr s, mpfr_ptr c, mpfr_srcptr x, mpfr_rnd_t rnd) { int inexs, inexc; mpfr_t x_red, ts, tc; diff --git a/src/strtofr.c b/src/strtofr.c index c6aeb7391..f45fd5516 100644 --- a/src/strtofr.c +++ b/src/strtofr.c @@ -226,7 +226,7 @@ fast_casecmp (const char *s1, const char *s2) BUT if it returns 0 (NAN or INF), the ternary value is also '0' (ie NAN and INF are exact) */ static int -parse_string (mpfr_t x, struct parsed_string *pstr, +parse_string (mpfr_ptr x, struct parsed_string *pstr, const char **string, int base) { const char *str = *string; @@ -451,7 +451,7 @@ parse_string (mpfr_t x, struct parsed_string *pstr, and the precision of x. Returns the ternary value. */ static int -parsed_string_to_mpfr (mpfr_t x, struct parsed_string *pstr, mpfr_rnd_t rnd) +parsed_string_to_mpfr (mpfr_ptr x, struct parsed_string *pstr, mpfr_rnd_t rnd) { mpfr_prec_t precx, prec, ysize_bits, pstr_size; mpfr_exp_t exp; @@ -934,7 +934,7 @@ free_parsed_string (struct parsed_string *pstr) } int -mpfr_strtofr (mpfr_t x, const char *string, char **end, int base, +mpfr_strtofr (mpfr_ptr x, const char *string, char **end, int base, mpfr_rnd_t rnd) { int res; diff --git a/src/vasprintf.c b/src/vasprintf.c index 6173f149c..9c8e533ea 100644 --- a/src/vasprintf.c +++ b/src/vasprintf.c @@ -963,7 +963,7 @@ floor_log10 (mpfr_srcptr x) #define NDIGITS 8 MPFR_RETURNS_NONNULL static char * -mpfr_get_str_wrapper (mpfr_exp_t *exp, int base, size_t n, const mpfr_t op, +mpfr_get_str_wrapper (mpfr_exp_t *exp, int base, size_t n, mpfr_srcptr op, const struct printf_spec spec) { size_t ndigits; diff --git a/src/zeta.c b/src/zeta.c index ec2d88de7..99382a270 100644 --- a/src/zeta.c +++ b/src/zeta.c @@ -35,7 +35,7 @@ https://www.gnu.org/licenses/ or write to the Free Software Foundation, Inc., sum(tc[i]*product((s+2j)*(s+2j-1)/n^2,j=1..i-1), i=1..p)*s*n^(-s-1) */ static void -mpfr_zeta_part_b (mpfr_t b, mpfr_srcptr s, int n, int p, mpfr_t *tc) +mpfr_zeta_part_b (mpfr_ptr b, mpfr_srcptr s, int n, int p, mpfr_t *tc) { mpfr_t s1, d, u; unsigned long n2; @@ -130,7 +130,7 @@ mpfr_zeta_c (int p, mpfr_t *tc) n - an integer Output: sum - a floating-point number approximating sum(1/i^s, i=1..n-1) */ static void -mpfr_zeta_part_a (mpfr_t sum, mpfr_srcptr s, int n) +mpfr_zeta_part_a (mpfr_ptr sum, mpfr_srcptr s, int n) { mpfr_t u, s1; int i; @@ -158,7 +158,7 @@ mpfr_zeta_part_a (mpfr_t sum, mpfr_srcptr s, int n) Output: z - Zeta(s) rounded to the precision of z with direction rnd_mode */ static int -mpfr_zeta_pos (mpfr_t z, mpfr_srcptr s, mpfr_rnd_t rnd_mode) +mpfr_zeta_pos (mpfr_ptr z, mpfr_srcptr s, mpfr_rnd_t rnd_mode) { mpfr_t b, c, z_pre, f, s1; double beta, sd, dnep; @@ -356,8 +356,8 @@ compute_add (mpfr_srcptr s, mpfr_prec_t precz) At input, p is Pi rounded down. The comments in the code are for rnd = RNDD. */ static void -mpfr_reflection_overflow (mpfr_t z, mpfr_t s1, const mpfr_t s, mpfr_t y, - mpfr_t p, mpfr_rnd_t rnd) +mpfr_reflection_overflow (mpfr_ptr z, mpfr_ptr s1, mpfr_srcptr s, mpfr_ptr y, + mpfr_ptr p, mpfr_rnd_t rnd) { mpz_t sint; @@ -432,7 +432,7 @@ mpfr_reflection_overflow (mpfr_t z, mpfr_t s1, const mpfr_t s, mpfr_t y, } int -mpfr_zeta (mpfr_t z, mpfr_srcptr s, mpfr_rnd_t rnd_mode) +mpfr_zeta (mpfr_ptr z, mpfr_srcptr s, mpfr_rnd_t rnd_mode) { mpfr_t z_pre, s1, y, p; long add; diff --git a/tests/tcmp2.c b/tests/tcmp2.c index 71e4bde56..858324dd4 100644 --- a/tests/tcmp2.c +++ b/tests/tcmp2.c @@ -24,7 +24,7 @@ https://www.gnu.org/licenses/ or write to the Free Software Foundation, Inc., /* set bit n of x to b, where bit 0 is the most significant one */ static void -set_bit (mpfr_t x, unsigned int n, int b) +set_bit (mpfr_ptr x, unsigned int n, int b) { unsigned l; mp_size_t xn; diff --git a/tests/tdiv.c b/tests/tdiv.c index 0beff3fda..736bc0130 100644 --- a/tests/tdiv.c +++ b/tests/tdiv.c @@ -369,7 +369,7 @@ check_convergence (void) /* given y = o(x/u), x, u, find the inexact flag by multiplying y by u */ static int -get_inexact (mpfr_t y, mpfr_t x, mpfr_t u) +get_inexact (mpfr_ptr y, mpfr_ptr x, mpfr_ptr u) { mpfr_t xx; int inex; diff --git a/tests/teq.c b/tests/teq.c index 81d2f4e5d..fbd59dcdb 100644 --- a/tests/teq.c +++ b/tests/teq.c @@ -23,7 +23,7 @@ https://www.gnu.org/licenses/ or write to the Free Software Foundation, Inc., #include "mpfr-test.h" static void -teq (mpfr_t x) +teq (mpfr_ptr x) { mpfr_t y; unsigned long k, px, mx; diff --git a/tests/terandom_chisq.c b/tests/terandom_chisq.c index c2f9af4ab..70a5f536a 100644 --- a/tests/terandom_chisq.c +++ b/tests/terandom_chisq.c @@ -26,7 +26,7 @@ https://www.gnu.org/licenses/ or write to the Free Software Foundation, Inc., * exponential distribution. We only take differences of this function so the * offset doesn't matter; here Phi(0) = 0. */ static void -exponential_cumulative (mpfr_t z, mpfr_t x, mpfr_rnd_t rnd) +exponential_cumulative (mpfr_ptr z, mpfr_ptr x, mpfr_rnd_t rnd) { mpfr_neg (z, x, rnd); mpfr_expm1 (z, z, rnd); @@ -43,7 +43,7 @@ exponential_cumulative (mpfr_t z, mpfr_t x, mpfr_rnd_t rnd) * TAOCP, Vol 2, 3.3.1, Table 1. It more accurate than the similar formula, * DLMF 8.11.10. */ static void -chisq_prob (mpfr_t q, long nu, mpfr_t chisqp) +chisq_prob (mpfr_ptr q, long nu, mpfr_ptr chisqp) { mpfr_t t; mpfr_rnd_t rnd; @@ -170,7 +170,7 @@ test_erandom_chisq_cont (long num, mpfr_prec_t prec, int nu, * this function. low precision means prec = 2, 3, or 4. High values of * precision will result in integer overflow. */ static long -sequential (mpfr_t x) +sequential (mpfr_ptr x) { long expt, prec; diff --git a/tests/tfmma.c b/tests/tfmma.c index 34567bd5a..97efc540e 100644 --- a/tests/tfmma.c +++ b/tests/tfmma.c @@ -24,7 +24,7 @@ https://www.gnu.org/licenses/ or write to the Free Software Foundation, Inc., /* check both mpfr_fmma and mpfr_fmms */ static void -random_test (mpfr_t a, mpfr_t b, mpfr_t c, mpfr_t d, mpfr_rnd_t rnd) +random_test (mpfr_ptr a, mpfr_ptr b, mpfr_ptr c, mpfr_ptr d, mpfr_rnd_t rnd) { mpfr_t ref, res, ab, cd; int inex_ref, inex_res; diff --git a/tests/tfmod.c b/tests/tfmod.c index 2fd24ed42..3148047ea 100644 --- a/tests/tfmod.c +++ b/tests/tfmod.c @@ -62,8 +62,8 @@ slow_fmod (mpfr_ptr r, mpfr_srcptr x, mpfr_srcptr y, mpfr_rnd_t rnd) } static void -test_failed (mpfr_t erem, mpfr_t grem, int eret, int gret, mpfr_t x, mpfr_t y, - mpfr_rnd_t rnd) +test_failed (mpfr_ptr erem, mpfr_ptr grem, int eret, int gret, + mpfr_ptr x, mpfr_ptr y, mpfr_rnd_t rnd) { printf ("error: mpfr_fmod (r, x, y, rnd)\n x = "); mpfr_out_str (stdout, 10, 0, x, MPFR_RNDD); @@ -83,7 +83,7 @@ test_failed (mpfr_t erem, mpfr_t grem, int eret, int gret, mpfr_t x, mpfr_t y, } static void -check (mpfr_t r0, mpfr_t x, mpfr_t y, mpfr_rnd_t rnd) +check (mpfr_ptr r0, mpfr_ptr x, mpfr_ptr y, mpfr_rnd_t rnd) { int inex0, inex1; mpfr_t r1; diff --git a/tests/tfprintf.c b/tests/tfprintf.c index a52f19f68..4b3bc113d 100644 --- a/tests/tfprintf.c +++ b/tests/tfprintf.c @@ -65,7 +65,7 @@ https://www.gnu.org/licenses/ or write to the Free Software Foundation, Inc., const int prec_max_printf = 5000; static void -check (FILE *fout, const char *fmt, mpfr_t x) +check (FILE *fout, const char *fmt, mpfr_ptr x) { if (mpfr_fprintf (fout, fmt, x) == -1) { diff --git a/tests/tgamma.c b/tests/tgamma.c index bb43e351d..04af9d367 100644 --- a/tests/tgamma.c +++ b/tests/tgamma.c @@ -890,7 +890,7 @@ tiny (int stop) computing with a working precision p2. Assume that x is not an integer <= 2. */ static void -exp_lgamma (mpfr_t x, mpfr_prec_t p1, mpfr_prec_t p2) +exp_lgamma (mpfr_ptr x, mpfr_prec_t p1, mpfr_prec_t p2) { mpfr_t yd, yu, zd, zu; int inexd, inexu, sign; diff --git a/tests/tnrandom_chisq.c b/tests/tnrandom_chisq.c index 0dcea31ad..122e413de 100644 --- a/tests/tnrandom_chisq.c +++ b/tests/tnrandom_chisq.c @@ -26,7 +26,7 @@ https://www.gnu.org/licenses/ or write to the Free Software Foundation, Inc., * for the normal distribution. We only take differences of this function so * the offset doesn't matter; here Phi(0) = 0. */ static void -normal_cumulative (mpfr_t z, mpfr_t x, mpfr_rnd_t rnd) +normal_cumulative (mpfr_ptr z, mpfr_ptr x, mpfr_rnd_t rnd) { mpfr_sqrt_ui (z, 2, rnd); mpfr_div (z, x, z, rnd); @@ -44,7 +44,7 @@ normal_cumulative (mpfr_t z, mpfr_t x, mpfr_rnd_t rnd) * TAOCP, Vol 2, 3.3.1, Table 1. It more accurate than the similar formula, * DLMF 8.11.10. */ static void -chisq_prob (mpfr_t q, long nu, mpfr_t chisqp) +chisq_prob (mpfr_ptr q, long nu, mpfr_ptr chisqp) { mpfr_t t; mpfr_rnd_t rnd; @@ -166,7 +166,7 @@ test_nrandom_chisq_cont (long num, mpfr_prec_t prec, int nu, * this function. low precision means prec = 2, 3, or 4. High values of * precision will result in integer overflow. */ static long -sequential (mpfr_t x) +sequential (mpfr_ptr x) { long expt, prec; diff --git a/tests/tprintf.c b/tests/tprintf.c index 4ea0ef620..2a724dd75 100644 --- a/tests/tprintf.c +++ b/tests/tprintf.c @@ -74,7 +74,7 @@ const int prec_max_printf = 5000; int stdout_redirect; static void -check (const char *fmt, mpfr_t x) +check (const char *fmt, mpfr_ptr x) { if (mpfr_printf (fmt, x) == -1) { diff --git a/tests/trint.c b/tests/trint.c index 2c136e670..1ab5246e2 100644 --- a/tests/trint.c +++ b/tests/trint.c @@ -367,7 +367,7 @@ test_against_libc (void) #endif static void -err (const char *str, mp_size_t s, mpfr_t x, mpfr_t y, mpfr_prec_t p, +err (const char *str, mp_size_t s, mpfr_ptr x, mpfr_ptr y, mpfr_prec_t p, mpfr_rnd_t r, int trint, int inexact) { printf ("Error: %s\ns = %u, p = %u, r = %s, trint = %d, inexact = %d\nx = ", diff --git a/tests/tsinh_cosh.c b/tests/tsinh_cosh.c index 3cce163df..712c24099 100644 --- a/tests/tsinh_cosh.c +++ b/tests/tsinh_cosh.c @@ -23,7 +23,7 @@ https://www.gnu.org/licenses/ or write to the Free Software Foundation, Inc., #include "mpfr-test.h" static void -failed (mpfr_t x, mpfr_t esh, mpfr_t gsh, mpfr_t ech, mpfr_t gch) +failed (mpfr_ptr x, mpfr_ptr esh, mpfr_ptr gsh, mpfr_ptr ech, mpfr_ptr gch) { printf ("error : mpfr_sinh_cosh (x) x = "); mpfr_out_str (stdout, 10, 0, x, MPFR_RNDD); @@ -43,7 +43,7 @@ failed (mpfr_t x, mpfr_t esh, mpfr_t gsh, mpfr_t ech, mpfr_t gch) /* check against sinh, cosh */ static void -check (mpfr_t x, mpfr_rnd_t rnd) +check (mpfr_ptr x, mpfr_rnd_t rnd) { mpfr_t s, c, sx, cx; int isc, is, ic; diff --git a/tests/tsqr.c b/tests/tsqr.c index 4907d6433..4cafcc617 100644 --- a/tests/tsqr.c +++ b/tests/tsqr.c @@ -33,7 +33,7 @@ inexact_sign (int x) static void error1 (mpfr_rnd_t rnd, mpfr_prec_t prec, - mpfr_t in, mpfr_t outmul, mpfr_t outsqr) + mpfr_t in, mpfr_ptr outmul, mpfr_ptr outsqr) { printf("ERROR: for %s and prec=%lu\nINPUT=", mpfr_print_rnd_mode(rnd), (unsigned long) prec); @@ -44,7 +44,7 @@ error1 (mpfr_rnd_t rnd, mpfr_prec_t prec, } static void -error2 (mpfr_rnd_t rnd, mpfr_prec_t prec, mpfr_t in, mpfr_t out, +error2 (mpfr_rnd_t rnd, mpfr_prec_t prec, mpfr_ptr in, mpfr_ptr out, int inexactmul, int inexactsqr) { printf("ERROR: for %s and prec=%lu\nINPUT=", mpfr_print_rnd_mode(rnd), diff --git a/tests/tsum.c b/tests/tsum.c index 866465d2c..29561f433 100644 --- a/tests/tsum.c +++ b/tests/tsum.c @@ -59,7 +59,7 @@ get_prec_max (mpfr_t *t, int n) } static void -get_exact_sum (mpfr_t sum, mpfr_t *tab, int n) +get_exact_sum (mpfr_ptr sum, mpfr_t *tab, int n) { int i; @@ -1198,7 +1198,7 @@ check_coverage (void) } static int -mpfr_sum_naive (mpfr_t s, mpfr_t *x, int n, mpfr_rnd_t rnd) +mpfr_sum_naive (mpfr_ptr s, mpfr_t *x, int n, mpfr_rnd_t rnd) { int ret, i; switch (n) |