1 /*
2 Copyright (C) 2014 Fredrik Johansson
3
4 This file is part of Arb.
5
6 Arb is free software: you can redistribute it and/or modify it under
7 the terms of the GNU Lesser General Public License (LGPL) as published
8 by the Free Software Foundation; either version 2.1 of the License, or
9 (at your option) any later version. See <http://www.gnu.org/licenses/>.
10 */
11
12 #include "acb.h"
13
14 void
acb_get_rad_ubound_arf(arf_t u,const acb_t z,slong prec)15 acb_get_rad_ubound_arf(arf_t u, const acb_t z, slong prec)
16 {
17 /* fixme: this bound is very sloppy */
18
19 if (mag_cmp(arb_radref(acb_realref(z)), arb_radref(acb_imagref(z))) >= 0)
20 arf_set_mag(u, arb_radref(acb_realref(z)));
21 else
22 arf_set_mag(u, arb_radref(acb_imagref(z)));
23
24 arf_mul_2exp_si(u, u, 1);
25 }
26