Linux下计算pai

最新推荐文章于 2023-05-26 22:21:24 发布

甲壳

最新推荐文章于 2023-05-26 22:21:24 发布

阅读量756

点赞数

分类专栏： Linux 文章标签： linux c pai

Linux 专栏收录该内容

1 篇文章 0 订阅

订阅专栏

双线程计算pai

#include <stdio.h>
#include <stdlib.h>
#include <pthread.h>
#include <sys/time.h>
#define N 1000000000
double pai = 0.0, sum = 0.0;
double step = 1.0/(double)N;
void *count_thread(void)
{
    int i;
    double x, temp;
    for(i = 0; i < N / 2; i ++)
    {
        x = (i + 0.5) * step;
        temp = (4.0 / (1.0 + x *x));
        pai += temp;
    }
}
int main()
{
    pthread_t child_thread;
    int i;
    struct timeval start_time, end_time;
        double x, temp, tempsum;
    double spend_time;
    gettimeofday(&start_time, NULL);
    pthread_create(&child_thread, NULL, count_thread, NULL);
        for(i = N / 2; i < N; i++)
        {
        x = (i + 0.5) * step;
        temp = (4.0 / (1.0 + x * x));
            tempsum += temp;
        }
    sum += tempsum;
    pthread_join(child_thread,NULL);
        pai = pai + sum;
    printf("PAI = %12.9f\n", pai * step);
    gettimeofday(&end_time, 0);
    spend_time = 1000000*(end_time.tv_sec - start_time.tv_sec)+end_time.tv_usec - start_time.tv_usec;
    spend_time /= 1000000;
    printf("It takes %f seconds to calculate PI.\n",spend_time);
    exit(0);
}

单线程计算pai

#include <stdio.h>
#include <stdlib.h>
#include <pthread.h>
#include <sys/time.h>
#define N 1000000000
int main()
{
    double pai = 0.0;
    double step, x, temp, spend_time;
    int i;
    struct timeval start_time, end_time;
    gettimeofday(&start_time, NULL);
    step = 1 / (double)N;
    for(i = 0; i < N; i++)
    {
    x = (i + 0.5) * step;
    temp = (4.0 / (1.0 + x * x));
    pai += temp;
    }
    pai *= step;
    printf("PAI = %.9f\n",pai);
    gettimeofday(&end_time, 0);
    spend_time = 1000000*(end_time.tv_sec - start_time.tv_sec)+end_time.tv_usec - start_time.tv_usec;
    spend_time /=1000000;
    printf("It takes %f seconds to calculate PI.\n",spend_time);
    exit(0);
}