mpf_add فقدت الدقة على GMP. حيث ذهب بلدي الأرقام؟
سؤال
وأنا تلخيص اثنين يطفو السلبية:
char * lhs = "-2234.6016114467412141";
char * rhs = "-4939600281397002.2812";
ووفقا لبيرل، وذلك باستخدام bignum والرياضيات :: BigFloat، فإن الجواب هو
-4939600281399236.8828114467412141
ولكن، وفقا لبرنامج الرصد العالمي، وذلك باستخدام رمز أدناه، فإن الجواب هو
-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 (الافتراضي) القيم، ولكن دون جدوى.
المحلول
ويبدو أنك تفيض الجزء العشري في 64 بت. حاول القيام mpf_get_prec(f_res)
للتحقق من ذلك هو الدقة التي تريدها. إن لم يكن استدعاء mpf_set_default_prec()
قبل تهيئة أي فار MPF (الخط 1 من الرئيسي).
نصائح أخرى
والتمهيد مع قيمة أكبر. جرب هذا:
mpf_set_default_prec(5*1024)
لا تنتمي إلى StackOverflow