Compute factorial(1000) with libtommath

factorial_tommath.c

#include <stdio.h>
#include "tommath.h"

int main()
{
    mp_int a, b, c;

    mp_init(&a);
    mp_init(&b);
    mp_init(&c);

    mp_set_int(&a, 1);

    int n = 0;
    while (n <= 1000) {
        //char buf[4096];
        //mp_toradix_n(&a, buf, 10, sizeof(buf)-1);
        //printf("fac(%d) = %s\n", n, buf);
        mp_set_int(&b, ++n);
        mp_mul(&a, &b, &c);
        mp_copy(&c, &a);
    }

    mp_clear(&a);
    mp_clear(&b);
    mp_clear(&c);

    return 0;
}