A C library for '''M'''ultiple-'''P'''recision '''F'''loating-point computations with correct '''R'''ounding.

http://www.mpfr.org/mpfr500.png

http://www.mpfr.org/

[GMP]를 사용한 실수전용 라이브러리
== 예제 ==
from http://www.mpfr.org/sample.html

다음은 200-bit 정밀도(precision)로 1+1/1!+1/2!+...+1/100!를 구하는 프로그램입니다.

{{{mpfr_t s, t, u;}}}
 세개의 s, t, u 실수를 정의
{{{mpfr_init2 (t, 200);}}}
 200-bit 정밀도로 초기화
{{{mpfr_set_d (t, 1.0, GMP_RNDD);}}}
 변수 t에 double형 변수 1.0을 집어넣되, -Infinity쪽(downward rounding) 반올림을 하도록 한다. (이 경우 반올림이 일어나지 않는데, 1은 200비트에서 정확하게 1로 표현되기 때문.
{{{mpfr_mul_ui (t, t, i, GMP_RNDU);}}}
 t에 unsigned int인 i를 곱해서 t에 넣는다. 결과값은 +Infinity쪽 (upward rounding)
{{{mpfr_div (u, u, t, GMP_RNDD);}}}
 u를 t로 나눈다. 결과는 -Infinity쪽으로 반올림하고 u에 저장한다.
{{{mpfr_out_str (stdout, 10, 0, s, GMP_RNDD);}}}
 10을 base로 하는 출력을 s로 -Infinity쪽으로 반올림. 세번째 인자 0은 프린트될 숫자를 자동으로 설정하도록 한다.

{{{mpfr_clear}}}
 MPFR에서 사용된 공간을 지운다.

{{{#!vim c
#include <stdio.h>

#include <gmp.h>
#include <mpfr.h>

int main (void)
{
  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_clear (s);
  mpfr_clear (t);
  mpfr_clear (u);
  return 0;
}

}}}

결과
{{{
$ ./sample
Sum is 2.718281828459045235360287471352662497757247093699959574966913
}}}