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

de-inline mag_mul_2exp_si

Browse source
This commit is contained in:
fredrik 2018-09-17 11:21:11 +02:00
parent ea2c300a1f
commit e9a9293516
3 changed files with 49 additions and 17 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

20
arb/mul_2exp_si.c
View file

@ -11,10 +11,28 @@
#include "arb.h" #include "arb.h"
/* mag_mul_2exp_si is non-inline, but avoid overhead here */
static __inline__ void
_mag_mul_2exp_si(mag_t z, const mag_t x, slong y)
{
if (mag_is_special(x))
{
mag_set(z, x);
}
else
{
if (y >= ADD2_FAST_MIN && y <= ADD2_FAST_MAX)
_fmpz_add_fast(MAG_EXPREF(z), MAG_EXPREF(x), y);
else
fmpz_add_si(MAG_EXPREF(z), MAG_EXPREF(x), y);
MAG_MAN(z) = MAG_MAN(x);
}
}
void void
arb_mul_2exp_si(arb_t y, const arb_t x, slong e) arb_mul_2exp_si(arb_t y, const arb_t x, slong e)
{ {
arf_mul_2exp_si(arb_midref(y), arb_midref(x), e); arf_mul_2exp_si(arb_midref(y), arb_midref(x), e);
mag_mul_2exp_si(arb_radref(y), arb_radref(x), e); _mag_mul_2exp_si(arb_radref(y), arb_radref(x), e);
} }

17
mag.h
View file

@ -319,22 +319,7 @@ mag_inv_lower(mag_t res, const mag_t x)
mag_div_lower(res, t, x); mag_div_lower(res, t, x);
} }
MAG_INLINE void void mag_mul_2exp_si(mag_t z, const mag_t x, slong y);
mag_mul_2exp_si(mag_t z, const mag_t x, slong y)
{
if (mag_is_special(x))
{
mag_set(z, x);
}
else
{
if (y >= ADD2_FAST_MIN && y <= ADD2_FAST_MAX)
_fmpz_add_fast(MAG_EXPREF(z), MAG_EXPREF(x), y);
else
fmpz_add_si(MAG_EXPREF(z), MAG_EXPREF(x), y);
MAG_MAN(z) = MAG_MAN(x);
}
}
MAG_INLINE void MAG_INLINE void
mag_mul_2exp_fmpz(mag_t z, const mag_t x, const fmpz_t y) mag_mul_2exp_fmpz(mag_t z, const mag_t x, const fmpz_t y)

29
mag/mul_2exp_si.c Normal file
View file

@ -0,0 +1,29 @@
/*
Copyright (C) 2014 Fredrik Johansson
This file is part of Arb.
Arb is free software: you can redistribute it and/or modify it under
the terms of the GNU Lesser General Public License (LGPL) as published
by the Free Software Foundation; either version 2.1 of the License, or
(at your option) any later version. See <http://www.gnu.org/licenses/>.
*/
#include "mag.h"
void
mag_mul_2exp_si(mag_t z, const mag_t x, slong y)
{
if (mag_is_special(x))
{
mag_set(z, x);
}
else
{
if (y >= ADD2_FAST_MIN && y <= ADD2_FAST_MAX)
_fmpz_add_fast(MAG_EXPREF(z), MAG_EXPREF(x), y);
else
fmpz_add_si(MAG_EXPREF(z), MAG_EXPREF(x), y);
MAG_MAN(z) = MAG_MAN(x);
}
}