cheng/mpir
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Change MPIR v3+ tuning for mpn_get_str and mpn_set_str to match GMP5

Browse Source
This commit is contained in:
(no author) 2010-02-26 10:59:11 +00:00
parent 88f3536359
commit 5e5f7e57a5
3 changed files with 90 additions and 14 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
build.vc9/lib_speed/lib_speed.vcproj
View File

@ -496,6 +496,10 @@
RelativePath="..\..\tune\set_strb.c"
>
</File>
<File
RelativePath="..\..\tune\set_strp.c"
>
</File>
<File
RelativePath="..\..\tune\set_strs.c"
>

10
gmp-impl.h
View File

@ -4301,9 +4301,13 @@ extern mp_size_t get_str_dc_threshold;
#define GET_STR_PRECOMPUTE_THRESHOLD get_str_precompute_threshold
extern mp_size_t get_str_precompute_threshold;
#undef SET_STR_THRESHOLD
#define SET_STR_THRESHOLD set_str_threshold
extern mp_size_t SET_STR_THRESHOLD;
#undef SET_STR_DC_THRESHOLD
#define SET_STR_DC_THRESHOLD set_str_dc_threshold
extern mp_size_t set_str_dc_threshold;
#undef SET_STR_PRECOMPUTE_THRESHOLD
#define SET_STR_PRECOMPUTE_THRESHOLD set_str_precompute_threshold
extern mp_size_t set_str_precompute_threshold;
#undef FFT_TABLE_ATTRS
#define FFT_TABLE_ATTRS

90
tune/tuneup.c
View File

@ -224,7 +224,8 @@ mp_size_t mod_1_unnorm_threshold = MP_SIZE_T_MAX;
mp_size_t divrem_2_threshold = MP_SIZE_T_MAX;
mp_size_t get_str_dc_threshold = MP_SIZE_T_MAX;
mp_size_t get_str_precompute_threshold = MP_SIZE_T_MAX;
mp_size_t set_str_threshold = MP_SIZE_T_MAX;
mp_size_t set_str_dc_threshold = MP_SIZE_T_MAX;
mp_size_t set_str_precompute_threshold = MP_SIZE_T_MAX;
mp_size_t rootrem_threshold = MP_SIZE_T_MAX;
mp_size_t divrem_hensel_qr_1_threshold = MP_SIZE_T_MAX;
mp_size_t rsh_divrem_hensel_qr_1_threshold = MP_SIZE_T_MAX;
@ -1809,7 +1810,7 @@ tune_get_str (gmp_randstate_t rands)
GET_STR_PRECOMPUTE_THRESHOLD = 0;
param.name = "GET_STR_DC_THRESHOLD";
param.function = speed_mpn_get_str;
param.min_size = 2;
param.min_size = 4;
param.max_size = GET_STR_THRESHOLD_LIMIT;
one (&get_str_dc_threshold, rands, &param);
}
@ -1823,20 +1824,87 @@ tune_get_str (gmp_randstate_t rands)
}
}
double
speed_mpn_pre_set_str (struct speed_params *s)
{
unsigned char *str;
mp_ptr wp;
mp_size_t wn;
unsigned i;
int base;
double t;
mp_ptr powtab_mem, tp;
powers_t powtab[GMP_LIMB_BITS];
mp_size_t un;
int chars_per_limb;
TMP_DECL;
SPEED_RESTRICT_COND (s->size >= 1);
base = s->r == 0 ? 10 : s->r;
SPEED_RESTRICT_COND (base >= 2 && base <= 256);
TMP_MARK;
str = TMP_ALLOC (s->size);
for (i = 0; i < s->size; i++)
str[i] = s->xp[i] % base;
wn = ((mp_size_t) (s->size / mp_bases[base].chars_per_bit_exactly))
/ GMP_LIMB_BITS + 2;
SPEED_TMP_ALLOC_LIMBS (wp, wn, s->align_wp);
/* use this during development to check wn is big enough */
/*
ASSERT_ALWAYS (mpn_set_str (wp, str, s->size, base) <= wn);
*/
speed_operand_src (s, (mp_ptr) str, s->size/BYTES_PER_MP_LIMB);
speed_operand_dst (s, wp, wn);
speed_cache_fill (s);
chars_per_limb = mp_bases[base].chars_per_limb;
un = s->size / chars_per_limb + 1;
powtab_mem = TMP_BALLOC_LIMBS (mpn_dc_set_str_powtab_alloc (un));
mpn_set_str_compute_powtab (powtab, powtab_mem, un, base);
tp = TMP_BALLOC_LIMBS (mpn_dc_set_str_itch (un));
speed_starttime ();
i = s->reps;
do
{
mpn_pre_set_str (wp, str, s->size, powtab, tp);
}
while (--i != 0);
t = speed_endtime ();
TMP_FREE;
return t;
}
void
tune_set_str (gmp_randstate_t rands)
{
static struct param_t param;
s.r = 10; /* decimal */
param.step_factor = 0.1;
param.name = "SET_STR_THRESHOLD";
param.function = speed_mpn_set_str_basecase;
param.function2 = speed_mpn_set_str_subquad;
param.min_size = 100;
param.max_size = 150000;
one (&set_str_threshold, rands, &param);
{
static struct param_t param;
SET_STR_PRECOMPUTE_THRESHOLD = 0;
param.step_factor = 0.01;
param.name = "SET_STR_DC_THRESHOLD";
param.function = speed_mpn_pre_set_str;
param.min_size = 100;
param.max_size = 50000;
one (&set_str_dc_threshold, rands, &param);
}
{
static struct param_t param;
param.step_factor = 0.02;
param.name = "SET_STR_PRECOMPUTE_THRESHOLD";
param.function = speed_mpn_set_str;
param.min_size = SET_STR_DC_THRESHOLD;
param.max_size = 100000;
one (&set_str_precompute_threshold, rands, &param);
}
}