long -> slong acb.h.

This commit is contained in:
William Hart 2015-11-06 11:11:35 +00:00
parent 0aba1fa6f4
commit 24217696d7

348
acb.h
View file

@ -69,9 +69,9 @@ acb_clear(acb_t x)
}
ACB_INLINE acb_ptr
_acb_vec_init(long n)
_acb_vec_init(slong n)
{
long i;
slong i;
acb_ptr v = (acb_ptr) flint_malloc(sizeof(acb_struct) * n);
for (i = 0; i < n; i++)
@ -81,9 +81,9 @@ _acb_vec_init(long n)
}
ACB_INLINE void
_acb_vec_clear(acb_ptr v, long n)
_acb_vec_clear(acb_ptr v, slong n)
{
long i;
slong i;
for (i = 0; i < n; i++)
acb_clear(v + i);
flint_free(v);
@ -157,14 +157,14 @@ acb_set(acb_t z, const acb_t x)
}
ACB_INLINE void
acb_set_round(acb_t z, const acb_t x, long prec)
acb_set_round(acb_t z, const acb_t x, slong prec)
{
arb_set_round(acb_realref(z), acb_realref(x), prec);
arb_set_round(acb_imagref(z), acb_imagref(x), prec);
}
ACB_INLINE void
acb_neg_round(acb_t z, const acb_t x, long prec)
acb_neg_round(acb_t z, const acb_t x, slong prec)
{
arb_neg_round(acb_realref(z), acb_realref(x), prec);
arb_neg_round(acb_imagref(z), acb_imagref(x), prec);
@ -248,14 +248,14 @@ acb_set_d(acb_t z, double c)
}
ACB_INLINE void
acb_set_si(acb_t z, long c)
acb_set_si(acb_t z, slong c)
{
arb_set_si(acb_realref(z), c);
arb_zero(acb_imagref(z));
}
ACB_INLINE void
acb_set_si_si(acb_t z, long x, long y)
acb_set_si_si(acb_t z, slong x, slong y)
{
arb_set_si(acb_realref(z), x);
arb_set_si(acb_imagref(z), y);
@ -283,7 +283,7 @@ acb_set_fmpz_fmpz(acb_t z, const fmpz_t x, const fmpz_t y)
}
ACB_INLINE void
acb_set_round_fmpz(acb_t z, const fmpz_t y, long prec)
acb_set_round_fmpz(acb_t z, const fmpz_t y, slong prec)
{
arb_set_round_fmpz(acb_realref(z), y, prec);
arb_zero(acb_imagref(z));
@ -294,7 +294,7 @@ int acb_contains_int(const acb_t x);
int acb_get_unique_fmpz(fmpz_t z, const acb_t x);
ACB_INLINE void
acb_set_fmpq(acb_t z, const fmpq_t c, long prec)
acb_set_fmpq(acb_t z, const fmpq_t c, slong prec)
{
arb_set_fmpq(acb_realref(z), c, prec);
arb_zero(acb_imagref(z));
@ -315,7 +315,7 @@ acb_set_arb_arb(acb_t z, const arb_t x, const arb_t y)
}
ACB_INLINE void
acb_set_round_arb(acb_t z, const arb_t x, long prec)
acb_set_round_arb(acb_t z, const arb_t x, slong prec)
{
arb_set_round(acb_realref(z), x, prec);
arb_zero(acb_imagref(z));
@ -347,7 +347,7 @@ void acb_get_mag(mag_t z, const acb_t x);
void acb_get_mag_lower(mag_t z, const acb_t x);
ACB_INLINE void
acb_get_abs_ubound_arf(arf_t u, const acb_t z, long prec)
acb_get_abs_ubound_arf(arf_t u, const acb_t z, slong prec)
{
if (arb_is_zero(acb_imagref(z)))
{
@ -375,7 +375,7 @@ acb_get_abs_ubound_arf(arf_t u, const acb_t z, long prec)
}
ACB_INLINE void
acb_get_abs_lbound_arf(arf_t u, const acb_t z, long prec)
acb_get_abs_lbound_arf(arf_t u, const acb_t z, slong prec)
{
if (arb_is_zero(acb_imagref(z)))
{
@ -403,7 +403,7 @@ acb_get_abs_lbound_arf(arf_t u, const acb_t z, long prec)
}
ACB_INLINE void
acb_get_rad_ubound_arf(arf_t u, const acb_t z, long prec)
acb_get_rad_ubound_arf(arf_t u, const acb_t z, slong prec)
{
/* fixme: this bound is very sloppy */
@ -415,73 +415,73 @@ acb_get_rad_ubound_arf(arf_t u, const acb_t z, long prec)
arf_mul_2exp_si(u, u, 1);
}
void acb_arg(arb_t r, const acb_t z, long prec);
void acb_arg(arb_t r, const acb_t z, slong prec);
ACB_INLINE void
acb_add(acb_t z, const acb_t x, const acb_t y, long prec)
acb_add(acb_t z, const acb_t x, const acb_t y, slong prec)
{
arb_add(acb_realref(z), acb_realref(x), acb_realref(y), prec);
arb_add(acb_imagref(z), acb_imagref(x), acb_imagref(y), prec);
}
ACB_INLINE void
acb_sub(acb_t z, const acb_t x, const acb_t y, long prec)
acb_sub(acb_t z, const acb_t x, const acb_t y, slong prec)
{
arb_sub(acb_realref(z), acb_realref(x), acb_realref(y), prec);
arb_sub(acb_imagref(z), acb_imagref(x), acb_imagref(y), prec);
}
ACB_INLINE void
acb_add_si(acb_t z, const acb_t x, ulong c, long prec)
acb_add_si(acb_t z, const acb_t x, ulong c, slong prec)
{
arb_add_si(acb_realref(z), acb_realref(x), c, prec);
arb_set_round(acb_imagref(z), acb_imagref(x), prec);
}
ACB_INLINE void
acb_add_ui(acb_t z, const acb_t x, ulong c, long prec)
acb_add_ui(acb_t z, const acb_t x, ulong c, slong prec)
{
arb_add_ui(acb_realref(z), acb_realref(x), c, prec);
arb_set_round(acb_imagref(z), acb_imagref(x), prec);
}
ACB_INLINE void
acb_sub_si(acb_t z, const acb_t x, ulong c, long prec)
acb_sub_si(acb_t z, const acb_t x, ulong c, slong prec)
{
arb_sub_si(acb_realref(z), acb_realref(x), c, prec);
arb_set_round(acb_imagref(z), acb_imagref(x), prec);
}
ACB_INLINE void
acb_sub_ui(acb_t z, const acb_t x, ulong c, long prec)
acb_sub_ui(acb_t z, const acb_t x, ulong c, slong prec)
{
arb_sub_ui(acb_realref(z), acb_realref(x), c, prec);
arb_set_round(acb_imagref(z), acb_imagref(x), prec);
}
ACB_INLINE void
acb_add_fmpz(acb_t z, const acb_t x, const fmpz_t y, long prec)
acb_add_fmpz(acb_t z, const acb_t x, const fmpz_t y, slong prec)
{
arb_add_fmpz(acb_realref(z), acb_realref(x), y, prec);
arb_set_round(acb_imagref(z), acb_imagref(x), prec);
}
ACB_INLINE void
acb_add_arb(acb_t z, const acb_t x, const arb_t y, long prec)
acb_add_arb(acb_t z, const acb_t x, const arb_t y, slong prec)
{
arb_add(acb_realref(z), acb_realref(x), y, prec);
arb_set_round(acb_imagref(z), acb_imagref(x), prec);
}
ACB_INLINE void
acb_sub_fmpz(acb_t z, const acb_t x, const fmpz_t y, long prec)
acb_sub_fmpz(acb_t z, const acb_t x, const fmpz_t y, slong prec)
{
arb_sub_fmpz(acb_realref(z), acb_realref(x), y, prec);
arb_set_round(acb_imagref(z), acb_imagref(x), prec);
}
ACB_INLINE void
acb_sub_arb(acb_t z, const acb_t x, const arb_t y, long prec)
acb_sub_arb(acb_t z, const acb_t x, const arb_t y, slong prec)
{
arb_sub(acb_realref(z), acb_realref(x), y, prec);
arb_set_round(acb_imagref(z), acb_imagref(x), prec);
@ -502,34 +502,34 @@ acb_conj(acb_t z, const acb_t x)
}
ACB_INLINE void
acb_abs(arb_t u, const acb_t z, long prec)
acb_abs(arb_t u, const acb_t z, slong prec)
{
arb_hypot(u, acb_realref(z), acb_imagref(z), prec);
}
ACB_INLINE void
acb_mul_ui(acb_t z, const acb_t x, ulong y, long prec)
acb_mul_ui(acb_t z, const acb_t x, ulong y, slong prec)
{
arb_mul_ui(acb_realref(z), acb_realref(x), y, prec);
arb_mul_ui(acb_imagref(z), acb_imagref(x), y, prec);
}
ACB_INLINE void
acb_mul_si(acb_t z, const acb_t x, long y, long prec)
acb_mul_si(acb_t z, const acb_t x, slong y, slong prec)
{
arb_mul_si(acb_realref(z), acb_realref(x), y, prec);
arb_mul_si(acb_imagref(z), acb_imagref(x), y, prec);
}
ACB_INLINE void
acb_mul_fmpz(acb_t z, const acb_t x, const fmpz_t y, long prec)
acb_mul_fmpz(acb_t z, const acb_t x, const fmpz_t y, slong prec)
{
arb_mul_fmpz(acb_realref(z), acb_realref(x), y, prec);
arb_mul_fmpz(acb_imagref(z), acb_imagref(x), y, prec);
}
ACB_INLINE void
acb_mul_arb(acb_t z, const acb_t x, const arb_t y, long prec)
acb_mul_arb(acb_t z, const acb_t x, const arb_t y, slong prec)
{
arb_mul(acb_realref(z), acb_realref(x), y, prec);
arb_mul(acb_imagref(z), acb_imagref(x), y, prec);
@ -565,12 +565,12 @@ acb_div_onei(acb_t z, const acb_t x)
}
}
void acb_mul(acb_t z, const acb_t x, const acb_t y, long prec);
void acb_mul(acb_t z, const acb_t x, const acb_t y, slong prec);
void acb_mul_naive(acb_t z, const acb_t x, const acb_t y, long prec);
void acb_mul_naive(acb_t z, const acb_t x, const acb_t y, slong prec);
ACB_INLINE void
acb_mul_2exp_si(acb_t z, const acb_t x, long e)
acb_mul_2exp_si(acb_t z, const acb_t x, slong e)
{
arb_mul_2exp_si(acb_realref(z), acb_realref(x), e);
arb_mul_2exp_si(acb_imagref(z), acb_imagref(x), e);
@ -583,132 +583,132 @@ acb_mul_2exp_fmpz(acb_t z, const acb_t x, const fmpz_t c)
arb_mul_2exp_fmpz(acb_imagref(z), acb_imagref(x), c);
}
void acb_addmul(acb_t z, const acb_t x, const acb_t y, long prec);
void acb_addmul(acb_t z, const acb_t x, const acb_t y, slong prec);
void acb_submul(acb_t z, const acb_t x, const acb_t y, long prec);
void acb_submul(acb_t z, const acb_t x, const acb_t y, slong prec);
ACB_INLINE void
acb_addmul_ui(acb_t z, const acb_t x, ulong y, long prec)
acb_addmul_ui(acb_t z, const acb_t x, ulong y, slong prec)
{
arb_addmul_ui(acb_realref(z), acb_realref(x), y, prec);
arb_addmul_ui(acb_imagref(z), acb_imagref(x), y, prec);
}
ACB_INLINE void
acb_addmul_si(acb_t z, const acb_t x, long y, long prec)
acb_addmul_si(acb_t z, const acb_t x, slong y, slong prec)
{
arb_addmul_si(acb_realref(z), acb_realref(x), y, prec);
arb_addmul_si(acb_imagref(z), acb_imagref(x), y, prec);
}
ACB_INLINE void
acb_submul_ui(acb_t z, const acb_t x, ulong y, long prec)
acb_submul_ui(acb_t z, const acb_t x, ulong y, slong prec)
{
arb_submul_ui(acb_realref(z), acb_realref(x), y, prec);
arb_submul_ui(acb_imagref(z), acb_imagref(x), y, prec);
}
ACB_INLINE void
acb_submul_si(acb_t z, const acb_t x, long y, long prec)
acb_submul_si(acb_t z, const acb_t x, slong y, slong prec)
{
arb_submul_si(acb_realref(z), acb_realref(x), y, prec);
arb_submul_si(acb_imagref(z), acb_imagref(x), y, prec);
}
ACB_INLINE void
acb_addmul_fmpz(acb_t z, const acb_t x, const fmpz_t y, long prec)
acb_addmul_fmpz(acb_t z, const acb_t x, const fmpz_t y, slong prec)
{
arb_addmul_fmpz(acb_realref(z), acb_realref(x), y, prec);
arb_addmul_fmpz(acb_imagref(z), acb_imagref(x), y, prec);
}
ACB_INLINE void
acb_submul_fmpz(acb_t z, const acb_t x, const fmpz_t y, long prec)
acb_submul_fmpz(acb_t z, const acb_t x, const fmpz_t y, slong prec)
{
arb_submul_fmpz(acb_realref(z), acb_realref(x), y, prec);
arb_submul_fmpz(acb_imagref(z), acb_imagref(x), y, prec);
}
ACB_INLINE void
acb_addmul_arb(acb_t z, const acb_t x, const arb_t y, long prec)
acb_addmul_arb(acb_t z, const acb_t x, const arb_t y, slong prec)
{
arb_addmul(acb_realref(z), acb_realref(x), y, prec);
arb_addmul(acb_imagref(z), acb_imagref(x), y, prec);
}
ACB_INLINE void
acb_submul_arb(acb_t z, const acb_t x, const arb_t y, long prec)
acb_submul_arb(acb_t z, const acb_t x, const arb_t y, slong prec)
{
arb_submul(acb_realref(z), acb_realref(x), y, prec);
arb_submul(acb_imagref(z), acb_imagref(x), y, prec);
}
void acb_inv(acb_t z, const acb_t x, long prec);
void acb_inv(acb_t z, const acb_t x, slong prec);
void acb_div(acb_t z, const acb_t x, const acb_t y, long prec);
void acb_div(acb_t z, const acb_t x, const acb_t y, slong prec);
ACB_INLINE void
acb_div_ui(acb_t z, const acb_t x, ulong c, long prec)
acb_div_ui(acb_t z, const acb_t x, ulong c, slong prec)
{
arb_div_ui(acb_realref(z), acb_realref(x), c, prec);
arb_div_ui(acb_imagref(z), acb_imagref(x), c, prec);
}
ACB_INLINE void
acb_div_si(acb_t z, const acb_t x, long c, long prec)
acb_div_si(acb_t z, const acb_t x, slong c, slong prec)
{
arb_div_si(acb_realref(z), acb_realref(x), c, prec);
arb_div_si(acb_imagref(z), acb_imagref(x), c, prec);
}
ACB_INLINE void
acb_div_arb(acb_t z, const acb_t x, const arb_t c, long prec)
acb_div_arb(acb_t z, const acb_t x, const arb_t c, slong prec)
{
arb_div(acb_realref(z), acb_realref(x), c, prec);
arb_div(acb_imagref(z), acb_imagref(x), c, prec);
}
ACB_INLINE void
acb_div_fmpz(acb_t z, const acb_t x, const fmpz_t c, long prec)
acb_div_fmpz(acb_t z, const acb_t x, const fmpz_t c, slong prec)
{
arb_div_fmpz(acb_realref(z), acb_realref(x), c, prec);
arb_div_fmpz(acb_imagref(z), acb_imagref(x), c, prec);
}
void acb_cube(acb_t y, const acb_t x, long prec);
void acb_pow_fmpz(acb_t y, const acb_t b, const fmpz_t e, long prec);
void acb_pow_ui(acb_t y, const acb_t b, ulong e, long prec);
void acb_pow_si(acb_t y, const acb_t b, long e, long prec);
void acb_cube(acb_t y, const acb_t x, slong prec);
void acb_pow_fmpz(acb_t y, const acb_t b, const fmpz_t e, slong prec);
void acb_pow_ui(acb_t y, const acb_t b, ulong e, slong prec);
void acb_pow_si(acb_t y, const acb_t b, slong e, slong prec);
ACB_INLINE void
acb_const_pi(acb_t x, long prec)
acb_const_pi(acb_t x, slong prec)
{
arb_const_pi(acb_realref(x), prec);
arb_zero(acb_imagref(x));
}
void acb_log(acb_t r, const acb_t z, long prec);
void acb_log1p(acb_t r, const acb_t z, long prec);
void acb_log(acb_t r, const acb_t z, slong prec);
void acb_log1p(acb_t r, const acb_t z, slong prec);
void acb_exp(acb_t r, const acb_t z, long prec);
void acb_exp_pi_i(acb_t r, const acb_t z, long prec);
void acb_exp_invexp(acb_t r, acb_t s, const acb_t z, long prec);
void acb_exp(acb_t r, const acb_t z, slong prec);
void acb_exp_pi_i(acb_t r, const acb_t z, slong prec);
void acb_exp_invexp(acb_t r, acb_t s, const acb_t z, slong prec);
void acb_sin(acb_t r, const acb_t z, long prec);
void acb_cos(acb_t r, const acb_t z, long prec);
void acb_sin_cos(acb_t s, acb_t c, const acb_t z, long prec);
void acb_tan(acb_t r, const acb_t z, long prec);
void acb_cot(acb_t r, const acb_t z, long prec);
void acb_sin(acb_t r, const acb_t z, slong prec);
void acb_cos(acb_t r, const acb_t z, slong prec);
void acb_sin_cos(acb_t s, acb_t c, const acb_t z, slong prec);
void acb_tan(acb_t r, const acb_t z, slong prec);
void acb_cot(acb_t r, const acb_t z, slong prec);
void acb_asin(acb_t r, const acb_t z, long prec);
void acb_acos(acb_t r, const acb_t z, long prec);
void acb_atan(acb_t r, const acb_t z, long prec);
void acb_asinh(acb_t r, const acb_t z, long prec);
void acb_acosh(acb_t r, const acb_t z, long prec);
void acb_atanh(acb_t r, const acb_t z, long prec);
void acb_asin(acb_t r, const acb_t z, slong prec);
void acb_acos(acb_t r, const acb_t z, slong prec);
void acb_atan(acb_t r, const acb_t z, slong prec);
void acb_asinh(acb_t r, const acb_t z, slong prec);
void acb_acosh(acb_t r, const acb_t z, slong prec);
void acb_atanh(acb_t r, const acb_t z, slong prec);
ACB_INLINE void
acb_sinh(acb_t y, const acb_t x, long prec)
acb_sinh(acb_t y, const acb_t x, slong prec)
{
acb_mul_onei(y, x);
acb_sin(y, y, prec);
@ -716,14 +716,14 @@ acb_sinh(acb_t y, const acb_t x, long prec)
}
ACB_INLINE void
acb_cosh(acb_t y, const acb_t x, long prec)
acb_cosh(acb_t y, const acb_t x, slong prec)
{
acb_mul_onei(y, x);
acb_cos(y, y, prec);
}
ACB_INLINE void
acb_sinh_cosh(acb_t y, acb_t z, const acb_t x, long prec)
acb_sinh_cosh(acb_t y, acb_t z, const acb_t x, slong prec)
{
acb_mul_onei(y, x);
acb_sin_cos(y, z, y, prec);
@ -731,7 +731,7 @@ acb_sinh_cosh(acb_t y, acb_t z, const acb_t x, long prec)
}
ACB_INLINE void
acb_tanh(acb_t y, const acb_t x, long prec)
acb_tanh(acb_t y, const acb_t x, slong prec)
{
acb_mul_onei(y, x);
acb_tan(y, y, prec);
@ -739,67 +739,67 @@ acb_tanh(acb_t y, const acb_t x, long prec)
}
ACB_INLINE void
acb_coth(acb_t y, const acb_t x, long prec)
acb_coth(acb_t y, const acb_t x, slong prec)
{
acb_mul_onei(y, x);
acb_cot(y, y, prec);
acb_mul_onei(y, y);
}
void acb_sin_pi(acb_t r, const acb_t z, long prec);
void acb_cos_pi(acb_t r, const acb_t z, long prec);
void acb_sin_cos_pi(acb_t s, acb_t c, const acb_t z, long prec);
void acb_tan_pi(acb_t r, const acb_t z, long prec);
void acb_cot_pi(acb_t r, const acb_t z, long prec);
void acb_sin_pi(acb_t r, const acb_t z, slong prec);
void acb_cos_pi(acb_t r, const acb_t z, slong prec);
void acb_sin_cos_pi(acb_t s, acb_t c, const acb_t z, slong prec);
void acb_tan_pi(acb_t r, const acb_t z, slong prec);
void acb_cot_pi(acb_t r, const acb_t z, slong prec);
void acb_pow_arb(acb_t z, const acb_t x, const arb_t y, long prec);
void acb_pow(acb_t r, const acb_t x, const acb_t y, long prec);
void acb_pow_arb(acb_t z, const acb_t x, const arb_t y, slong prec);
void acb_pow(acb_t r, const acb_t x, const acb_t y, slong prec);
void acb_sqrt(acb_t y, const acb_t x, long prec);
void acb_rsqrt(acb_t y, const acb_t x, long prec);
void acb_sqrt(acb_t y, const acb_t x, slong prec);
void acb_rsqrt(acb_t y, const acb_t x, slong prec);
void acb_chebyshev_t_ui(acb_t a, ulong n, const acb_t x, long prec);
void acb_chebyshev_t2_ui(acb_t a, acb_t b, ulong n, const acb_t x, long prec);
void acb_chebyshev_u_ui(acb_t a, ulong n, const acb_t x, long prec);
void acb_chebyshev_u2_ui(acb_t a, acb_t b, ulong n, const acb_t x, long prec);
void acb_chebyshev_t_ui(acb_t a, ulong n, const acb_t x, slong prec);
void acb_chebyshev_t2_ui(acb_t a, acb_t b, ulong n, const acb_t x, slong prec);
void acb_chebyshev_u_ui(acb_t a, ulong n, const acb_t x, slong prec);
void acb_chebyshev_u2_ui(acb_t a, acb_t b, ulong n, const acb_t x, slong prec);
void acb_rising_ui_bs(acb_t y, const acb_t x, ulong n, long prec);
void acb_rising_ui_rs(acb_t y, const acb_t x, ulong n, ulong m, long prec);
void acb_rising_ui_rec(acb_t y, const acb_t x, ulong n, long prec);
void acb_rising_ui(acb_t z, const acb_t x, ulong n, long prec);
void acb_rising(acb_t z, const acb_t x, const acb_t n, long prec);
void acb_rising_ui_bs(acb_t y, const acb_t x, ulong n, slong prec);
void acb_rising_ui_rs(acb_t y, const acb_t x, ulong n, ulong m, slong prec);
void acb_rising_ui_rec(acb_t y, const acb_t x, ulong n, slong prec);
void acb_rising_ui(acb_t z, const acb_t x, ulong n, slong prec);
void acb_rising(acb_t z, const acb_t x, const acb_t n, slong prec);
void acb_rising2_ui_bs(acb_t u, acb_t v, const acb_t x, ulong n, long prec);
void acb_rising2_ui_rs(acb_t u, acb_t v, const acb_t x, ulong n, ulong m, long prec);
void acb_rising2_ui(acb_t u, acb_t v, const acb_t x, ulong n, long prec);
void acb_rising2_ui_bs(acb_t u, acb_t v, const acb_t x, ulong n, slong prec);
void acb_rising2_ui_rs(acb_t u, acb_t v, const acb_t x, ulong n, ulong m, slong prec);
void acb_rising2_ui(acb_t u, acb_t v, const acb_t x, ulong n, slong prec);
void acb_rising_ui_get_mag(mag_t bound, const acb_t s, ulong n);
void acb_gamma(acb_t y, const acb_t x, long prec);
void acb_rgamma(acb_t y, const acb_t x, long prec);
void acb_lgamma(acb_t y, const acb_t x, long prec);
void acb_log_sin_pi(acb_t res, const acb_t z, long prec);
void acb_digamma(acb_t y, const acb_t x, long prec);
void acb_zeta(acb_t z, const acb_t s, long prec);
void acb_hurwitz_zeta(acb_t z, const acb_t s, const acb_t a, long prec);
void acb_polygamma(acb_t res, const acb_t s, const acb_t z, long prec);
void acb_gamma(acb_t y, const acb_t x, slong prec);
void acb_rgamma(acb_t y, const acb_t x, slong prec);
void acb_lgamma(acb_t y, const acb_t x, slong prec);
void acb_log_sin_pi(acb_t res, const acb_t z, slong prec);
void acb_digamma(acb_t y, const acb_t x, slong prec);
void acb_zeta(acb_t z, const acb_t s, slong prec);
void acb_hurwitz_zeta(acb_t z, const acb_t s, const acb_t a, slong prec);
void acb_polygamma(acb_t res, const acb_t s, const acb_t z, slong prec);
void acb_log_barnes_g(acb_t res, const acb_t z, long prec);
void acb_barnes_g(acb_t res, const acb_t z, long prec);
void acb_log_barnes_g(acb_t res, const acb_t z, slong prec);
void acb_barnes_g(acb_t res, const acb_t z, slong prec);
void acb_polylog(acb_t w, const acb_t s, const acb_t z, long prec);
void acb_polylog_si(acb_t w, long s, const acb_t z, long prec);
void acb_polylog(acb_t w, const acb_t s, const acb_t z, slong prec);
void acb_polylog_si(acb_t w, slong s, const acb_t z, slong prec);
void acb_agm1(acb_t m, const acb_t z, long prec);
void acb_agm1_cpx(acb_ptr m, const acb_t z, long len, long prec);
void acb_agm1(acb_t m, const acb_t z, slong prec);
void acb_agm1_cpx(acb_ptr m, const acb_t z, slong len, slong prec);
/*
TBD
void acb_invroot_newton(acb_t r, const acb_t a, ulong m, const acb_t r0, long startprec, long prec);
void acb_root_exp(acb_t r, const acb_t a, long m, long index, long prec);
void acb_root_newton(acb_t r, const acb_t a, long m, long index, long prec);
void acb_root(acb_t r, const acb_t a, long m, long index, long prec);
void acb_invroot_newton(acb_t r, const acb_t a, ulong m, const acb_t r0, slong startprec, slong prec);
void acb_root_exp(acb_t r, const acb_t a, slong m, slong index, slong prec);
void acb_root_newton(acb_t r, const acb_t a, slong m, slong index, slong prec);
void acb_root(acb_t r, const acb_t a, slong m, slong index, slong prec);
*/
ACB_INLINE int
@ -816,23 +816,23 @@ acb_indeterminate(acb_t x)
}
ACB_INLINE acb_ptr
_acb_vec_entry_ptr(acb_ptr vec, long i)
_acb_vec_entry_ptr(acb_ptr vec, slong i)
{
return vec + i;
}
ACB_INLINE void
_acb_vec_zero(acb_ptr A, long n)
_acb_vec_zero(acb_ptr A, slong n)
{
long i;
slong i;
for (i = 0; i < n; i++)
acb_zero(A + i);
}
ACB_INLINE int
_acb_vec_is_zero(acb_srcptr vec, long len)
_acb_vec_is_zero(acb_srcptr vec, slong len)
{
long i;
slong i;
for (i = 0; i < len; i++)
if (!acb_is_zero(vec + i))
return 0;
@ -840,121 +840,121 @@ _acb_vec_is_zero(acb_srcptr vec, long len)
}
ACB_INLINE void
_acb_vec_set(acb_ptr res, acb_srcptr vec, long len)
_acb_vec_set(acb_ptr res, acb_srcptr vec, slong len)
{
long i;
slong i;
for (i = 0; i < len; i++)
acb_set(res + i, vec + i);
}
ACB_INLINE void
_acb_vec_set_round(acb_ptr res, acb_srcptr vec, long len, long prec)
_acb_vec_set_round(acb_ptr res, acb_srcptr vec, slong len, slong prec)
{
long i;
slong i;
for (i = 0; i < len; i++)
acb_set_round(res + i, vec + i, prec);
}
ACB_INLINE void
_acb_vec_neg(acb_ptr res, acb_srcptr vec, long len)
_acb_vec_neg(acb_ptr res, acb_srcptr vec, slong len)
{
long i;
slong i;
for (i = 0; i < len; i++)
acb_neg(res + i, vec + i);
}
ACB_INLINE void
_acb_vec_add(acb_ptr res, acb_srcptr vec1, acb_srcptr vec2, long len, long prec)
_acb_vec_add(acb_ptr res, acb_srcptr vec1, acb_srcptr vec2, slong len, slong prec)
{
long i;
slong i;
for (i = 0; i < len; i++)
acb_add(res + i, vec1 + i, vec2 + i, prec);
}
ACB_INLINE void
_acb_vec_sub(acb_ptr res, acb_srcptr vec1, acb_srcptr vec2, long len, long prec)
_acb_vec_sub(acb_ptr res, acb_srcptr vec1, acb_srcptr vec2, slong len, slong prec)
{
long i;
slong i;
for (i = 0; i < len; i++)
acb_sub(res + i, vec1 + i, vec2 + i, prec);
}
ACB_INLINE void
_acb_vec_scalar_submul(acb_ptr res, acb_srcptr vec, long len, const acb_t c, long prec)
_acb_vec_scalar_submul(acb_ptr res, acb_srcptr vec, slong len, const acb_t c, slong prec)
{
long i;
slong i;
for (i = 0; i < len; i++)
acb_submul(res + i, vec + i, c, prec);
}
ACB_INLINE void
_acb_vec_scalar_addmul(acb_ptr res, acb_srcptr vec, long len, const acb_t c, long prec)
_acb_vec_scalar_addmul(acb_ptr res, acb_srcptr vec, slong len, const acb_t c, slong prec)
{
long i;
slong i;
for (i = 0; i < len; i++)
acb_addmul(res + i, vec + i, c, prec);
}
ACB_INLINE void
_acb_vec_scalar_mul(acb_ptr res, acb_srcptr vec, long len, const acb_t c, long prec)
_acb_vec_scalar_mul(acb_ptr res, acb_srcptr vec, slong len, const acb_t c, slong prec)
{
long i;
slong i;
for (i = 0; i < len; i++)
acb_mul(res + i, vec + i, c, prec);
}
ACB_INLINE void
_acb_vec_scalar_mul_ui(acb_ptr res, acb_srcptr vec, long len, ulong c, long prec)
_acb_vec_scalar_mul_ui(acb_ptr res, acb_srcptr vec, slong len, ulong c, slong prec)
{
long i;
slong i;
for (i = 0; i < len; i++)
acb_mul_ui(res + i, vec + i, c, prec);
}
ACB_INLINE void
_acb_vec_scalar_mul_2exp_si(acb_ptr res, acb_srcptr vec, long len, long c)
_acb_vec_scalar_mul_2exp_si(acb_ptr res, acb_srcptr vec, slong len, slong c)
{
long i;
slong i;
for (i = 0; i < len; i++)
acb_mul_2exp_si(res + i, vec + i, c);
}
ACB_INLINE void
_acb_vec_scalar_mul_onei(acb_ptr res, acb_srcptr vec, long len)
_acb_vec_scalar_mul_onei(acb_ptr res, acb_srcptr vec, slong len)
{
long i;
slong i;
for (i = 0; i < len; i++)
acb_mul_onei(res + i, vec + i);
}
ACB_INLINE void
_acb_vec_scalar_div_ui(acb_ptr res, acb_srcptr vec, long len, ulong c, long prec)
_acb_vec_scalar_div_ui(acb_ptr res, acb_srcptr vec, slong len, ulong c, slong prec)
{
long i;
slong i;
for (i = 0; i < len; i++)
acb_div_ui(res + i, vec + i, c, prec);
}
ACB_INLINE void
_acb_vec_scalar_div(acb_ptr res, acb_srcptr vec, long len, const acb_t c, long prec)
_acb_vec_scalar_div(acb_ptr res, acb_srcptr vec, slong len, const acb_t c, slong prec)
{
long i;
slong i;
for (i = 0; i < len; i++)
acb_div(res + i, vec + i, c, prec);
}
ACB_INLINE void
_acb_vec_scalar_mul_arb(acb_ptr res, acb_srcptr vec, long len, const arb_t c, long prec)
_acb_vec_scalar_mul_arb(acb_ptr res, acb_srcptr vec, slong len, const arb_t c, slong prec)
{
long i;
slong i;
for (i = 0; i < len; i++)
acb_mul_arb(res + i, vec + i, c, prec);
}
ACB_INLINE void
_acb_vec_scalar_div_arb(acb_ptr res, acb_srcptr vec, long len, const arb_t c, long prec)
_acb_vec_scalar_div_arb(acb_ptr res, acb_srcptr vec, slong len, const arb_t c, slong prec)
{
long i;
slong i;
for (i = 0; i < len; i++)
{
arb_div(acb_realref(res + i), acb_realref(vec + i), c, prec);
@ -963,17 +963,17 @@ _acb_vec_scalar_div_arb(acb_ptr res, acb_srcptr vec, long len, const arb_t c, lo
}
ACB_INLINE void
_acb_vec_scalar_mul_fmpz(acb_ptr res, acb_srcptr vec, long len, const fmpz_t c, long prec)
_acb_vec_scalar_mul_fmpz(acb_ptr res, acb_srcptr vec, slong len, const fmpz_t c, slong prec)
{
long i;
slong i;
for (i = 0; i < len; i++)
acb_mul_fmpz(res + i, vec + i, c, prec);
}
ACB_INLINE void
_acb_vec_scalar_div_fmpz(acb_ptr res, acb_srcptr vec, long len, const fmpz_t c, long prec)
_acb_vec_scalar_div_fmpz(acb_ptr res, acb_srcptr vec, slong len, const fmpz_t c, slong prec)
{
long i;
slong i;
for (i = 0; i < len; i++)
acb_div_fmpz(res + i, vec + i, c, prec);
}
@ -988,15 +988,15 @@ acb_print(const acb_t x)
printf(")");
}
void acb_printd(const acb_t z, long digits);
void acb_printd(const acb_t z, slong digits);
void acb_randtest(acb_t z, flint_rand_t state, long prec, long mag_bits);
void acb_randtest(acb_t z, flint_rand_t state, slong prec, slong mag_bits);
void acb_randtest_special(acb_t z, flint_rand_t state, long prec, long mag_bits);
void acb_randtest_special(acb_t z, flint_rand_t state, slong prec, slong mag_bits);
void acb_randtest_precise(acb_t z, flint_rand_t state, long prec, long mag_bits);
void acb_randtest_precise(acb_t z, flint_rand_t state, slong prec, slong mag_bits);
void acb_randtest_param(acb_t z, flint_rand_t state, long prec, long mag_bits);
void acb_randtest_param(acb_t z, flint_rand_t state, slong prec, slong mag_bits);
long acb_rel_error_bits(const acb_t x);
@ -1010,7 +1010,7 @@ acb_rel_accuracy_bits(const acb_t x)
ACB_INLINE long
acb_bits(const acb_t x)
{
long b1, b2;
slong b1, b2;
b1 = arb_bits(acb_realref(x));
b2 = arb_bits(acb_imagref(x));
return FLINT_MAX(b1, b2);
@ -1023,9 +1023,9 @@ acb_is_real(const acb_t x)
}
ACB_INLINE int
_acb_vec_is_real(acb_srcptr v, long len)
_acb_vec_is_real(acb_srcptr v, slong len)
{
long i;
slong i;
for (i = 0; i < len; i++)
{
@ -1037,15 +1037,15 @@ _acb_vec_is_real(acb_srcptr v, long len)
}
ACB_INLINE long
_acb_vec_bits(acb_srcptr vec, long len)
_acb_vec_bits(acb_srcptr vec, slong len)
{
return _arb_vec_bits((arb_srcptr) vec, 2 * len);
}
ACB_INLINE void
_acb_vec_set_powers(acb_ptr xs, const acb_t x, long len, long prec)
_acb_vec_set_powers(acb_ptr xs, const acb_t x, slong len, slong prec)
{
long i;
slong i;
for (i = 0; i < len; i++)
{
@ -1061,17 +1061,17 @@ _acb_vec_set_powers(acb_ptr xs, const acb_t x, long len, long prec)
}
ACB_INLINE void
_acb_vec_add_error_arf_vec(acb_ptr res, arf_srcptr err, long len)
_acb_vec_add_error_arf_vec(acb_ptr res, arf_srcptr err, slong len)
{
long i;
slong i;
for (i = 0; i < len; i++)
acb_add_error_arf(res + i, err + i);
}
ACB_INLINE void
_acb_vec_add_error_mag_vec(acb_ptr res, mag_srcptr err, long len)
_acb_vec_add_error_mag_vec(acb_ptr res, mag_srcptr err, slong len)
{
long i;
slong i;
for (i = 0; i < len; i++)
{
mag_add(arb_radref(acb_realref(res + i)),
@ -1082,23 +1082,23 @@ _acb_vec_add_error_mag_vec(acb_ptr res, mag_srcptr err, long len)
}
ACB_INLINE void
_acb_vec_indeterminate(acb_ptr vec, long len)
_acb_vec_indeterminate(acb_ptr vec, slong len)
{
_arb_vec_indeterminate((arb_ptr) vec, 2 * len);
}
ACB_INLINE void
_acb_vec_trim(acb_ptr res, acb_srcptr vec, long len)
_acb_vec_trim(acb_ptr res, acb_srcptr vec, slong len)
{
long i;
slong i;
for (i = 0; i < len; i++)
acb_trim(res + i, vec + i);
}
ACB_INLINE int
_acb_vec_get_unique_fmpz_vec(fmpz * res, acb_srcptr vec, long len)
_acb_vec_get_unique_fmpz_vec(fmpz * res, acb_srcptr vec, slong len)
{
long i;
slong i;
for (i = 0; i < len; i++)
if (!acb_get_unique_fmpz(res + i, vec + i))
return 0;
@ -1106,7 +1106,7 @@ _acb_vec_get_unique_fmpz_vec(fmpz * res, acb_srcptr vec, long len)
}
/* sort complex numbers in a nice-to-display order */
void _acb_vec_sort_pretty(acb_ptr vec, long len);
void _acb_vec_sort_pretty(acb_ptr vec, slong len);
#ifdef __cplusplus
}