Question

我已编写C代码来计算https://github.com/gongchengra/hacker/blob/master/c/16_pi.c

上的pi

没有评论的代码如下所示：

#include <stdio.h>
#define NUMBER 100000
int array_divide_number(int *array, int number, int size)
{
    int i,tmp;
    int modulo=0;
    for(i=0;i<size;i++)
    {
    tmp = array[i]+modulo*10;
    array[i] = tmp/number;
    modulo = tmp%number;
    }
    return 0;
}
void print_array(int *array, int size)
{
    int i,last;
    for(last=size-1;last>=0;last--)
    {
    if(array[last] != 0)
    {
        break;
    }
    }
    for(i=0;i<=last;i++)
    {
    printf("%d",array[i]);
    }
    printf("\n");
}
void copy_array(int *source, int *target, int size)
{
    int i;
    for(i=0;i<size;i++)
    {
    target[i] = source[i];
    }
}
void plus_array(int *augend, int *addend, int *sum, int size)
{
    int i;
    for(i=size-1;i>=0;i--)
    {
    sum[i] = augend[i] + addend[i];
    if(sum[i]>9)
    {
        sum[i] = sum[i] % 10;
        sum[i-1]++;
    }
    }
}
void minus_array(int *minuend, int *subtracter, int *answer, int size)
{
    int i;
    for(i=size-1;i>=0;i--)
    {
    if(minuend[i] >= subtracter[i] || i == 0)
    {
        answer[i] = minuend[i] - subtracter[i];
    }
    else
    {
        if(minuend[i-1] == 0)
        {
            minuend[i-2]--;
            minuend[i-1] = 10;
        }
        minuend[i-1]--;
        answer[i] = 10 + minuend[i] - subtracter[i];
    }
    }
}
int main()
{
    int i;
    int flag=1;
    int d5[NUMBER]={0};
    int t5[NUMBER]={0};
    int d239[NUMBER]={0};
    int t239[NUMBER]={0};
    int pi[NUMBER]={0};
    d5[0]=16;
    d239[0]=4;
    array_divide_number(d5,5,NUMBER);
    array_divide_number(d239,239,NUMBER);
    //every iteration will increase three valid digitals
    for(i=1;i<NUMBER*3/2;i+=2)
    {
    copy_array(d5, t5, NUMBER);
    copy_array(d239, t239, NUMBER);
    array_divide_number(t5,i,NUMBER);
    array_divide_number(t239,i,NUMBER);
    if(flag > 0)
    {
        plus_array(pi,t5,pi,NUMBER);
        minus_array(pi,t239,pi,NUMBER);
    }
    else
    {
        minus_array(pi,t5,pi,NUMBER);
        plus_array(pi,t239,pi,NUMBER);
    }
    flag = -1*flag;
    array_divide_number(d5,5*5,NUMBER);
    array_divide_number(d239,239*239,NUMBER);
    }
    print_array(pi, NUMBER);
    return 0;
}

它有效，但需要太长时间：

$time ./16_pi.exe >pi3.log ./16_pi.exe > pi3.log 617.76s user 0.15s system 98% cpu 10:27.94 total

您可以看到，为了计算pi的100000位数，需要10分钟以上。有没有办法优化代码而不改变计算pi的算法（这里使用的算法是pi =（16 / 5-4 / 239）-1 / 3 *（16/5 ^ 3-4 / 239 ^ 3）+ 1/5 *（16/5 ^ 5-4 / 239 ^ 5）+ ...）

Answer 1

优化代码以某种方式改变算法（即使公式相同）。

你正在做一些bignum算术。 Bignum算法既硬又聪明。您应该使用像GMPlib这样的bignum库（它可能也会从特殊的进位传播机器指令中获利;主要收益是聪明的bignum算法算法），它可能比你天真的bignum例程运行得更快。

当然，在基准测试时不要忘记要求编译器进行优化（例如gcc -O3 -Wall -mtune=native）

FWIW，另见Fabrice Bellard's page about Pi。并且GMP有一个page with code about Pi：计算100,000个数字只需要一小部分时间。

优化c代码以计算pi

1 个答案: