From 5506233919829b5f0ec15f1a4eca8d0f00468d45 Mon Sep 17 00:00:00 2001 From: zimmerma Date: Fri, 20 Dec 2013 06:42:05 +0000 Subject: speed up mpc_sin for small input (http://lists.gforge.inria.fr/pipermail/mpc-discuss/2013-December/001266.html) git-svn-id: svn://scm.gforge.inria.fr/svn/mpc/trunk@1395 211d60ee-9f03-0410-a15a-8952a2c7a4e4 --- src/sin_cos.c | 5 +++-- 1 file changed, 3 insertions(+), 2 deletions(-) (limited to 'src') diff --git a/src/sin_cos.c b/src/sin_cos.c index 0bb00cf..980e413 100644 --- a/src/sin_cos.c +++ b/src/sin_cos.c @@ -299,7 +299,7 @@ mpc_sin_cos (mpc_ptr rop_sin, mpc_ptr rop_cos, mpc_srcptr op, mpfr_t s, c, sh, ch, sch, csh; mpfr_prec_t prec; int ok; - int inex_re, inex_im, inex_sin, inex_cos; + int inex_re, inex_im, inex_sin, inex_cos, loop = 0; prec = 2; if (rop_sin != NULL) @@ -315,8 +315,9 @@ mpc_sin_cos (mpc_ptr rop_sin, mpc_ptr rop_cos, mpc_srcptr op, mpfr_init2 (csh, 2); do { + loop ++; ok = 1; - prec += mpc_ceil_log2 (prec) + 5; + prec += (loop == 1) ? mpc_ceil_log2 (prec) + 5 : prec / 2; mpfr_set_prec (s, prec); mpfr_set_prec (c, prec); -- cgit v1.2.1