OpenXM_contrib/gmp/mpf/get_d.c - diff

Return to get_d.c CVS log

Up to [local] / OpenXM_contrib / gmp / mpf

version 1.1.1.2, 2000/09/09 14:13:14

version 1.1.1.3, 2003/08/25 16:06:35

Line 1

/* double mpf_get_d (mpf_t src) -- Return the double approximation to SRC.

This file is part of the GNU MP Library.

Line 23 MA 02111-1307, USA. */

#include "gmp-impl.h"

double

#if __STDC__

mpf_get_d (mpf_srcptr src)

#else

mpf_get_d (src)

mpf_srcptr src;

#endif

{

double res;

mp_size_t size, i, n_limbs_to_use;

Line 48 mpf_get_d (src)

Line 43 mpf_get_d (src)

for (i = 2; i <= n_limbs_to_use; i++)

res = res * MP_BASE_AS_DOUBLE + qp[size - i];

res = __gmp_scale2 (res, (EXP(src) - n_limbs_to_use) * BITS_PER_MP_LIMB);

res = __gmp_scale2 (res, (EXP(src) - n_limbs_to_use) * GMP_NUMB_BITS);

return negative ? -res : res;

}