Загубена точност на GMP mpf_add. Къде изчезнаха цифрите ми?

Сумирам две отрицателни плаващи стойности:

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

Според Perl, използвайки bignum и Math::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 (по подразбиране) стойности, но без ефект.


c precision gmp
person bugmagnet    schedule 21.01.2009    source източник

Отговори (2)


arrow_upward
3
arrow_downward

Изглежда, че препълвате мантисата при 64 бита. Опитайте да направите mpf_get_prec(f_res), за да проверите дали е желаната точност. Ако не, извикайте mpf_set_default_prec(), преди да инициализирате всякакви mpf vars (ред 1 от main).

person robottobor    schedule 21.01.2009
comment
На място. Промених моя макрос на #define F(x) mpf_t (x); mpf_init2((x), 128L); и всичко работеше добре. - person bugmagnet; 21.01.2009

arrow_upward
0
arrow_downward

Зареждане с по-голяма стойност. Опитайте тази:

mpf_set_default_prec(5*1024)
person TANILSON    schedule 28.11.2009