GMP MPF_ADD의 정밀도 손실. 내 숫자는 어디로 갔습니까?
문제
나는 두 개의 부정적인 수레를 합산합니다.
char * lhs = "-2234.6016114467412141";
char * rhs = "-4939600281397002.2812";
Perl에 따르면 Bignum and 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)
추신 나는 전화를 추가하려고 노력했다 mpf_set_default_prec
64보다 크 (기본값) 값이 있지만 효과가 없습니다.
해결책
64 비트로 만티사를 넘치게하는 것 같습니다. 시도해보십시오 mpf_get_prec(f_res)
확인하는 것은 원하는 정밀도입니다. 전화하지 않으면 mpf_set_default_prec()
MPF VAR (메인 1 행)을 초기화하기 전에.
다른 팁
더 큰 값으로 부팅하십시오. 이 시도:
mpf_set_default_prec(5*1024)
제휴하지 않습니다 StackOverflow