调用mpfr的一段示例(极不完整)

来源:互联网 发布:解决网络高峰期问题 编辑:程序博客网 时间:2024/05/20 10:13
  char buffer[2048];  mp_exp_t  e;  unsigned int i;  mpfr_t s, t, u;  mpfr_init2 (t, 200);  mpfr_set_d (t, 1.0, GMP_RNDD);  mpfr_init2 (s, 200);  mpfr_set_d (s, 1.0, GMP_RNDD);  mpfr_init2 (u, 200);  for (i = 1; i <= 100; i++)    {      mpfr_mul_ui (t, t, i, GMP_RNDU);      mpfr_set_d (u, 1.0, GMP_RNDD);      mpfr_div (u, u, t, GMP_RNDD);      mpfr_add (s, s, u, GMP_RNDD);    }  //printf ("Sum is ");  //mpfr_out_str (stdout, 10, 0, s, GMP_RNDD);  //putchar ('\n');  mpfr_sprintf(buffer, "%.9R", s);  Memo1->Text = mpfr_get_str(NULL,&e,10,0, s, GMP_RNDD);  mpfr_clear (s);  mpfr_clear (t);  mpfr_clear (u);

原创粉丝点击