hiro/arb
1
0
Fork 0
mirror of https://github.com/vale981/arb synced 2025-03-04 17:01:40 -05:00
Code Issues Projects Releases Packages Wiki Activity Actions

small code cleanup

Browse source
This commit is contained in:
fredrik 2018-03-05 13:36:33 +01:00
parent 7afd3bfaf1
commit 97c263a7cc
2 changed files with 22 additions and 31 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

32
acb/acos.c
View file

@ -22,26 +22,22 @@ acb_acos(acb_t res, const acb_t z, slong prec)
{
acb_t t;
acb_init(t);
if (arb_is_zero(acb_imagref(z)))
acb_one(t);
if (arb_is_zero(acb_imagref(z))
&& arb_gt(acb_realref(z), acb_realref(t)))
{
arb_t one;
arb_init(one);
arb_one(one);
if (arb_gt(acb_realref(z), one))
{
acb_asin(res, z, prec);
acb_neg(res, res);
arb_zero(acb_realref(res));
arb_clear(one);
arb_clear(t);
return;
}
arb_clear(one);
/* pure imaginary on (1,inf) */
acb_asin(res, z, prec);
acb_neg(res, res);
arb_zero(acb_realref(res));
}
else
{
acb_asin(res, z, prec);
acb_const_pi(t, prec);
acb_mul_2exp_si(t, t, -1);
acb_sub(res, t, res, prec);
}
acb_asin(res, z, prec);
acb_const_pi(t, prec);
acb_mul_2exp_si(t, t, -1);
acb_sub(res, t, res, prec);
acb_clear(t);
}
}

21
acb/acosh.c
View file

@ -34,24 +34,19 @@ acb_acosh(acb_t res, const acb_t z, slong prec)
if (!arb_is_zero(acb_imagref(z)))
{
acb_log(res, t, prec);
acb_clear(t);
acb_clear(u);
return;
}
arb_t one, abs;
arb_init(one);
arb_init(abs);
arb_one(one);
arb_abs(abs, acb_realref(z));
acb_log(res, t, prec);
if (arb_lt(abs, one))
else
{
arb_zero(acb_realref(res));
/* pure imaginary on (-1,1) */
arb_abs(acb_realref(u), acb_realref(z));
arb_one(acb_imagref(u));
acb_log(res, t, prec);
if (arb_lt(acb_realref(u), acb_imagref(u)))
arb_zero(acb_realref(res));
}
acb_clear(t);
acb_clear(u);
arb_clear(one);
arb_clear(abs);
}
}