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

时间:2015-02-26 01:23:38

标签: c gmp

我决定将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代码中看到任何性能改进的可能性吗?

1 个答案:

答案 0 :(得分:2)

打乒乓球。你有两个变量a和b以及一个临时t。您将结果添加到t中,然后将b复制到a和t复制到b。而是在将b添加到a并将b添加到b之间切换。最终结果是a或b,取决于n是奇数还是偶数。