MPFR
|
A C library for Multiple-Precision Floating-point computations with correct Rounding.
GMP를 사용한 실수전용 라이브러리
예제 ¶다음은 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에서 사용된 공간을 지운다.
#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 |
   
|
