diff options
author | zimmerma <zimmerma@280ebfd0-de03-0410-8827-d642c229c3f4> | 1999-07-01 09:39:56 +0000 |
---|---|---|
committer | zimmerma <zimmerma@280ebfd0-de03-0410-8827-d642c229c3f4> | 1999-07-01 09:39:56 +0000 |
commit | b4b31ef7efc218ab482af439bbc72e862d771937 (patch) | |
tree | 93eba7027dff0c3a1a2bbbd51aa1609bf815604c /sqrt.c | |
parent | 2e82db13421ca2b670bf9f4e00a0ca244d892f9e (diff) | |
download | mpfr-b4b31ef7efc218ab482af439bbc72e862d771937.tar.gz |
initialize temporary variables with current precision
git-svn-id: svn://scm.gforge.inria.fr/svn/mpfr/trunk@231 280ebfd0-de03-0410-8827-d642c229c3f4
Diffstat (limited to 'sqrt.c')
-rw-r--r-- | sqrt.c | 2 |
1 files changed, 1 insertions, 1 deletions
@@ -31,7 +31,7 @@ mpfr_sqrt(X, a, rnd_mode) q = p + (int) ceil(log(4.0*ceil(log((double)q)/log(2.0))+2.0)/log(2.0)); err = q-p; /* the error is at most 2^err ulp */ q = (q/mp_bits_per_limb)*mp_bits_per_limb; /* adjust to entire limb */ - mpfr_init(t); mpfr_init(u); + mpfr_init(t, q+mp_bits_per_limb); mpfr_init(u, q+mp_bits_per_limb); do { q += mp_bits_per_limb; if (q>3*p+mp_bits_per_limb) { |