MPFR库计算精度的工作原理
MPFR库计算精度的工作原理
提问于 2020-04-04 08:56:13
我承认我在理解MPFR库时遇到了困难，我正在尝试计算大数的平方根，但我不确定如何定义mpft_t root变量的精度，也不确定最好的舍入方式是什么。
我的代码如下：
#include <stdio.h>
#include <stdlib.h>
#include <mpfr.h>
int main(){
mpfr_t bigNumber, bigNumber2, root2;
unsigned long int size = 1000000;
mpfr_init2(bigNumber, size);
mpfr_init2(bigNumber2, size);
mpfr_init2(root2, size);
mpfr_ui_pow_ui(bigNumber, 8, 20000, MPFR_RNDZ);
mpfr_sqrt(root2, bigNumber, MPFR_RNDA);
mpfr_pow_ui(bigNumber2, root2, 2, MPFR_RNDA);
return 0;
}但是，无论为mpz_t root设置的精度如何，结果通常都不令人满意。
mpfr_pow_ui (bigNumber2, root, 2, MPFR_RNDA)的结果通常不同于mpfr_t bigNumber，我指的不是浮点值，而是整个部分。
如何找出执行某个计算所需的精度？
此计算的最佳舍入形式是什么？
为什么会发生这种不准确的情况？
回答 1
Stack Overflow用户
回答已采纳
发布于 2021-02-19 02:40:02
代码工作正常(在GCC 9.3.0和MPFR 4.1.0上)，并返回原始数字。
#include <stdio.h>
#include <stdlib.h>
#include <mpfr.h>
int main()
{
mpfr_t bigNumber, bigNumber2, root2;
unsigned long int size = 100000;
mpfr_init2(bigNumber, size);
mpfr_init2(bigNumber2, size);
mpfr_init2(root2, size);
mpfr_ui_pow_ui(bigNumber, 8, 20000, MPFR_RNDD);
mpfr_out_str (stdout, 10, 0, bigNumber, MPFR_RNDD);
putchar ('\n');
putchar ('\n');
putchar ('\n');
mpfr_sqrt(root2, bigNumber, MPFR_RNDD);
mpfr_out_str (stdout, 10, 0, root2, MPFR_RNDD);
putchar ('\n');
putchar ('\n');
putchar ('\n');
mpfr_pow_ui(bigNumber2, root2, 2, MPFR_RNDD);
mpfr_out_str (stdout, 10, 0, bigNumber2, MPFR_RNDD);
putchar ('\n');
putchar ('\n');
putchar ('\n');
putchar ('\n');
return 0;
}页面原文内容由Stack Overflow提供。腾讯云小微IT领域专用引擎提供翻译支持
原文链接：
https://stackoverflow.com/questions/61022525
复制相关文章
相似问题
腾讯云开发者
