This commit is contained in:
Ralf Stephan 2018-03-04 17:46:31 +01:00
parent 6eabee7efa
commit 0673fc935f

View file

@ -20,6 +20,8 @@ acb_acos(acb_t res, const acb_t z, slong prec)
}
else
{
acb_t t;
acb_init(t);
if (arb_is_zero(acb_imagref(z)))
{
arb_t one;
@ -27,15 +29,16 @@ acb_acos(acb_t res, const acb_t z, slong prec)
arb_one(one);
if (arb_gt(acb_realref(z), one))
{
acb_zero(t);
acb_asin(res, z, prec);
acb_sub(res, t, res, prec);
arb_zero(acb_realref(res));
arb_clear(one);
arb_clear(t);
return;
}
arb_clear(one);
}
acb_t t;
acb_init(t);
acb_asin(res, z, prec);
acb_const_pi(t, prec);
acb_mul_2exp_si(t, t, -1);