GMP mpf_add的精度下降.我的手指哪里去了? [英] Lost precision on GMP mpf_add. Where have my digits gone?

问题描述

我正在对两个负浮点数求和:

I'm summing two negative floats:

char * lhs = "-2234.6016114467412141";
char * rhs = "-4939600281397002.2812";

根据Perl，使用bignum和Math :: BigFloat，答案是

According to Perl, using bignum and Math::BigFloat, the answer is

-4939600281399236.8828114467412141

但是，根据GMP，使用下面的代码，答案是

However, according to GMP, using the code below, the answer is

-4939600281399236.88281 

我哪里出问题了?剩下的"14467412141"怎么了?

Where have I gone wrong? What happened to the remaining "14467412141"?

#include "stdafx.h"
#include "gmp-static\gmp.h"
#include <stdlib.h>         /* For _MAX_PATH definition */
#include <stdio.h>
#include <malloc.h>
#include <math.h>

#define F(x) mpf_t x; mpf_init( x );

void main(void)
{
    F(f_lhs);
    F(f_rhs);
    F(f_res);

    char * resbuff;

    mp_exp_t exp;

    char * lhs = "-2234.6016114467412141";
    char * rhs = "-4939600281397002.2812";

    int validOp = mpf_set_str( f_lhs, lhs, 10 );
    validOp = mpf_set_str( f_rhs, rhs, 10 );

    mpf_add( f_res, f_lhs, f_rhs );

    resbuff = mpf_get_str( NULL, &exp, 10, 0, f_res );
    printf( "Using mpf_add, %s + %s = %s (exp=%d)\n", lhs, rhs, resbuff, exp );

    free(resbuff);
}

样本输出:

Using mpf_add, -2234.6016114467412141 + -4939600281397002.2812 = -493960028139923688281 (exp=16)

P.S.我尝试将呼叫添加到mpf_set_default_prec，且其值大于64(默认值)，但没有任何效果.

P.S. I have tried adding a call to mpf_set_default_prec with larger than 64 (the default) values, but to no effect.

推荐答案

看起来您的尾数溢出了64位.尝试执行mpf_get_prec(f_res)来检查它是否是您想要的精度.如果未初始化任何mpf var(main的第1行)，请先调用mpf_set_default_prec().

Looks like you are overflowing the mantissa at 64 bits. Try doing mpf_get_prec(f_res) to check it is the precision you want. If not call the mpf_set_default_prec() before you initialize any mpf vars (line 1 of main).

这篇关于GMP mpf_add的精度下降.我的手指哪里去了?的文章就介绍到这了，希望我们推荐的答案对大家有所帮助，也希望大家多多支持IT屋！