summaryrefslogtreecommitdiff
path: root/atan.c
diff options
context:
space:
mode:
authorvlefevre <vlefevre@280ebfd0-de03-0410-8827-d642c229c3f4>2007-06-20 13:22:49 +0000
committervlefevre <vlefevre@280ebfd0-de03-0410-8827-d642c229c3f4>2007-06-20 13:22:49 +0000
commit6ace8c2b5075f9aa849469ec73435e6e78713f26 (patch)
tree8d6d8aca017a947b461ac5bb5937a6f85edb6e82 /atan.c
parentd5913d021b41388498583ad4cd87eb00a586b37c (diff)
downloadmpfr-6ace8c2b5075f9aa849469ec73435e6e78713f26.tar.gz
Untabified and removed trailing spaces.
git-svn-id: svn://scm.gforge.inria.fr/svn/mpfr/trunk@4553 280ebfd0-de03-0410-8827-d642c229c3f4
Diffstat (limited to 'atan.c')
-rw-r--r--atan.c66
1 files changed, 33 insertions, 33 deletions
diff --git a/atan.c b/atan.c
index 755c89332..59c8ed9af 100644
--- a/atan.c
+++ b/atan.c
@@ -98,31 +98,31 @@ mpfr_atan_aux (mpfr_ptr y, mpz_ptr p, long r, int m, mpz_t *tab)
/* p <> 1: precompute ptoj table */
mpz_set (ptoj[0], p);
for (im = 1 ; im <= m ; im ++)
- mpz_mul (ptoj[im], ptoj[im - 1], ptoj[im - 1]);
+ mpz_mul (ptoj[im], ptoj[im - 1], ptoj[im - 1]);
/* main loop */
n = 1UL << m;
/* the ith term being X^i/(2i+1) with X=p/2^r, we can stop when
p^i/2^(r*i) < 2^(-precy), i.e. r*i > precy + log2(p^i) */
for (i = k = done = 0; (i < n) && (done == 0); i += 2, k ++)
- {
+ {
/* initialize both S[k],Q[k] and S[k+1],Q[k+1] */
- mpz_set_ui (Q[k+1], 2 * i + 3); /* Q(i+1,i+2) */
- mpz_mul_ui (S[k+1], p, 2 * i + 1); /* S(i+1,i+2) */
- mpz_mul_2exp (S[k], Q[k+1], r);
- mpz_sub (S[k], S[k], S[k+1]); /* S(i,i+2) */
- mpz_mul_ui (Q[k], Q[k+1], 2 * i + 1); /* Q(i,i+2) */
+ mpz_set_ui (Q[k+1], 2 * i + 3); /* Q(i+1,i+2) */
+ mpz_mul_ui (S[k+1], p, 2 * i + 1); /* S(i+1,i+2) */
+ mpz_mul_2exp (S[k], Q[k+1], r);
+ mpz_sub (S[k], S[k], S[k+1]); /* S(i,i+2) */
+ mpz_mul_ui (Q[k], Q[k+1], 2 * i + 1); /* Q(i,i+2) */
log2_nb_terms[k] = 1; /* S[k]/Q[k] corresponds to 2 terms */
- for (j = (i + 2) >> 1, l = 1; (j & 1) == 0; l ++, j >>= 1, k --)
- {
+ for (j = (i + 2) >> 1, l = 1; (j & 1) == 0; l ++, j >>= 1, k --)
+ {
/* invariant: S[k-1]/Q[k-1] and S[k]/Q[k] correspond
to 2^l terms each. We combine them into S[k-1]/Q[k-1] */
- MPFR_ASSERTD (k > 0);
- mpz_mul (S[k], S[k], Q[k-1]);
- mpz_mul (S[k], S[k], ptoj[l]);
- mpz_mul (S[k-1], S[k-1], Q[k]);
- mpz_mul_2exp (S[k-1], S[k-1], r << l);
- mpz_add (S[k-1], S[k-1], S[k]);
- mpz_mul (Q[k-1], Q[k-1], Q[k]);
+ MPFR_ASSERTD (k > 0);
+ mpz_mul (S[k], S[k], Q[k-1]);
+ mpz_mul (S[k], S[k], ptoj[l]);
+ mpz_mul (S[k-1], S[k-1], Q[k]);
+ mpz_mul_2exp (S[k-1], S[k-1], r << l);
+ mpz_add (S[k-1], S[k-1], S[k]);
+ mpz_mul (Q[k-1], Q[k-1], Q[k]);
log2_nb_terms[k-1] = l + 1;
/* now S[k-1]/Q[k-1] corresponds to 2^(l+1) terms */
MPFR_MPZ_SIZEINBASE2(mult, ptoj[l+1]);
@@ -131,31 +131,31 @@ mpfr_atan_aux (mpfr_ptr y, mpz_ptr p, long r, int m, mpz_t *tab)
accu[k-1] = (k == 1) ? mult : accu[k-2] + mult;
if (accu[k-1] > precy)
done = 1;
- }
- }
+ }
+ }
}
else /* special case p=1: the ith term being X^i/(2i+1) with X=1/2^r,
we can stop when r*i > precy i.e. i > precy/r */
{
n = 1UL << m;
for (i = k = 0; (i < n) && (i <= precy / r); i += 2, k ++)
- {
- mpz_set_ui (Q[k + 1], 2 * i + 3);
- mpz_mul_2exp (S[k], Q[k+1], r);
- mpz_sub_ui (S[k], S[k], 1 + 2 * i);
- mpz_mul_ui (Q[k], Q[k + 1], 1 + 2 * i);
+ {
+ mpz_set_ui (Q[k + 1], 2 * i + 3);
+ mpz_mul_2exp (S[k], Q[k+1], r);
+ mpz_sub_ui (S[k], S[k], 1 + 2 * i);
+ mpz_mul_ui (Q[k], Q[k + 1], 1 + 2 * i);
log2_nb_terms[k] = 1; /* S[k]/Q[k] corresponds to 2 terms */
- for (j = (i + 2) >> 1, l = 1; (j & 1) == 0; l++, j >>= 1, k --)
- {
- MPFR_ASSERTD (k > 0);
- mpz_mul (S[k], S[k], Q[k-1]);
- mpz_mul (S[k-1], S[k-1], Q[k]);
- mpz_mul_2exp (S[k-1], S[k-1], r << l);
- mpz_add (S[k-1], S[k-1], S[k]);
- mpz_mul (Q[k-1], Q[k-1], Q[k]);
+ for (j = (i + 2) >> 1, l = 1; (j & 1) == 0; l++, j >>= 1, k --)
+ {
+ MPFR_ASSERTD (k > 0);
+ mpz_mul (S[k], S[k], Q[k-1]);
+ mpz_mul (S[k-1], S[k-1], Q[k]);
+ mpz_mul_2exp (S[k-1], S[k-1], r << l);
+ mpz_add (S[k-1], S[k-1], S[k]);
+ mpz_mul (Q[k-1], Q[k-1], Q[k]);
log2_nb_terms[k-1] = l + 1;
- }
- }
+ }
+ }
}
/* we need to combine S[0]/Q[0]...S[k-1]/Q[k-1] */