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

fix pi_chudnovsky on 32-bit

Browse source
This commit is contained in:
Fredrik Johansson 2012-09-03 14:28:01 +02:00
parent 22f69e4758
commit 6dd11555cc
1 changed files with 8 additions and 1 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

9
fmprb/const_pi_chudnovsky.c
View file

@ -42,7 +42,14 @@ chudnovsky_bsplit(fmprb_t G, fmprb_t P, fmprb_t Q, long a, long b, long wp, int
fmprb_mul_si(G, G, 6*b-1, FMPR_PREC_EXACT);
/* p = C^3 * b^3 / 24 */
fmprb_set_si(P, CONST_C * CONST_C * CONST_C / 24); /* xxx: fix for 32-bit */
#if FLINT_BITS == 64
fmprb_set_si(P, CONST_C * CONST_C * CONST_C / 24);
#else
fmprb_set_ui(P, CONST_C / 24);
fmprb_mul_ui(P, P, CONST_C, FMPR_PREC_EXACT);
fmprb_mul_ui(P, P, CONST_C, FMPR_PREC_EXACT);
#endif
fmprb_mul_ui(P, P, b, FMPR_PREC_EXACT);
fmprb_mul_ui(P, P, b, FMPR_PREC_EXACT);
fmprb_mul_ui(P, P, b, FMPR_PREC_EXACT);