GNU Multiple Precision Arithmetic Library
Numerical instability로 고민하고 계시다면 GMP를 써보세요.
예제 ¶
#include <gmp.h>
int main()
{
mpf_t fac;
int i;
mpf_init(fac);
mpf_set_d(fac,(double) 1);
for (i=1;i<40;i++) {
mpf_mul_ui(fac, fac, (unsigned long) i);
}
gmp_printf("%.*Ff\n",40,fac);
}
