为什么我的C程序比等效的Haskell慢？

Question

我决定将Haskell的Fibonacci程序的标准严格尾递归版本与使用C语言编写的程序的性能进行比较，使用GMP进行比较，以使结果大到适合一个词（在Haskell中，我使用多精度Integer类型）。 我将省略Haskell程序，因为这是有关C和GMP的问题。 C实现是这样的：

#include <stdio.h>
#include <stdlib.h>
#include <gmp.h>

void fib(unsigned int n){
    mpz_t a, b, t;

    mpz_init_set_ui(a, 0);
    mpz_init_set_ui(b, 1);
    mpz_init(t);
    for(; n > 1; n --){
        mpz_add(t, a, b);
        mpz_set(a, b);
        mpz_set(b, t);
    }
    //mpz_out_str(stdout, 10, b);
}


int main(int argc, char **argv){
    unsigned long n, f;
    if(argc != 2){
        printf("Usage: fibc <number>\n");
        return 1;
    }

    fib(atol(argv[1]));

    return 0;
}

请注意，我注释掉了输出值的行，这大约需要一秒钟（我在Haskell版本中保留了此行为）。

结果是：

time ./fibhs 1000000
./fibhs 1000000  5.77s user 0.05s system 99% cpu 5.831 total

time ./fibc 1000000
./fibc 1000000  11.19s user 0.00s system 100% cpu 11.194 total

我认为我一定使用GMP错误。 谁能在C代码中看到任何性能改进的可能性吗？

Answer 1

打乒乓球。 您有两个变量a和b和一个临时t。 您将结果相加并将其放入t，然后将b复制到a并将t复制到b。 而是在将b添加到a和将a添加到b之间交替。 根据n是奇数还是偶数，最终结果是a或b。