diff options
| author | zimmerma <zimmerma@211d60ee-9f03-0410-a15a-8952a2c7a4e4> | 2013-12-20 06:42:05 +0000 |
|---|---|---|
| committer | zimmerma <zimmerma@211d60ee-9f03-0410-a15a-8952a2c7a4e4> | 2013-12-20 06:42:05 +0000 |
| commit | 5506233919829b5f0ec15f1a4eca8d0f00468d45 (patch) | |
| tree | 14d94aa0273320200216ae669ed8e0c6c7d93469 /src | |
| parent | 5712c925d1ecdcb575ee178a78d1ce8127cb370f (diff) | |
| download | mpc-5506233919829b5f0ec15f1a4eca8d0f00468d45.tar.gz | |
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
Diffstat (limited to 'src')
| -rw-r--r-- | src/sin_cos.c | 5 |
1 files changed, 3 insertions, 2 deletions
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); |