GMP mpf_add的精度損失。我的數字在哪裏？

我總結兩個負花車：GMP mpf_add的精度損失。我的數字在哪裏？

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

據Perl中，使用BIGNUM和數學:: BigFloat，答案是

-4939600281399236.8828114467412141

然而，按GMP，使用下面的代碼，答案是

-4939600281399236.88281

我哪裏出錯了？剩下的「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（默認值），但沒有效果。

來源

2009-01-21 bugmagnet

看起來你已經溢出了64位的尾數。試着做mpf_get_prec(f_res)檢查它是你想要的精度。如果在初始化任何mpf變量（main的第1行）之前未調用mpf_set_default_prec()。

來源

2009-01-21 08:48:38 robottobor

Spot on。我將宏改爲 #define F（x）mpf_t（x）; mpf_init2（（x），128L）; 和一切工作正常。 – bugmagnet 2009-01-21 09:47:08

以較大的值啓動。試試這個：

mpf_set_default_prec(5*1024)

來源

2009-11-28 17:55:25 TANILSON

GMP mpf_add的精度損失。我的數字在哪裏？

回答

相關問題