Perte de précision sur GMP mpf_add. Où sont passés mes chiffres?

Question

Je sommateurs deux flotteurs négatifs:

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

Selon Perl, en utilisant bignum et Math :: BigFloat, la réponse est

-4939600281399236.8828114467412141 

Toutefois, selon les BPF, en utilisant le code ci-dessous, la réponse est

-4939600281399236.88281 

Où est-ce que je me suis trompé? Qu'est-il arrivé au reste "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); 
} 

sortie de l'échantillon:

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

P.S. J'ai essayé d'ajouter un appel à mpf_set_default_prec avec des valeurs supérieures à 64 (par défaut), mais sans aucun effet.

Answer 1

On dirait que vous débordez de la mantisse à 64 bits. Essayez de faire mpf_get_prec(f_res) pour vérifier que c'est la précision que vous voulez. Si ce n'est pas le cas, appelez le mpf_set_default_prec() avant d'initialiser toutes les variables mpf (ligne 1 de la page principale).

Answer 2

Amorce avec une valeur plus grande. Essayez ceci:

mpf_set_default_prec(5*1024)