nth2

#include<stdio.h>
#include<pthread.h>
#include<time.h>

float *result;//result
int n;//the number which every thread should count

void p(void *i)
{
pthread_mutex_t mutex;//lock
pthread_mutex_init(&mutex, NULL);
pthread_mutex_lock(&mutex);

int a=*((int *)i);
int begin=0;
int end=0;
float x=0,y=0;
begin=a*n+1;
end=(a+1)*n;
int j=0;
float total=0;
for(j=begin;j<=end;j++)
{
if(j%3==0) x=y=1.0;
if(j%3==1) x=y=-1.0;
if(j%3==2) x=y=0.0;
total+=x*y;
}

result[a]=total;

pthread_mutex_unlock(&mutex);//realease the lock
}


int main(int args,char*argv[])
{



int threadnumber;//the number of threadnumber
int length;//

scanf("%d",&threadnumber);
scanf("%d",&length);



pthread_t pt[threadnumber];
int ret;
int i=0;

n=length/threadnumber;//每个线程计算的个数


result=(float*)malloc(sizeof(float)*threadnumber);



int a[threadnumber];
for(i=0;i<threadnumber;i++)
a[i]=i;

for(i=0;i<=threadnumber-2;i++)//开始线程
{
ret=pthread_create(&pt[i],NULL,(void*)p,&a[i]);
if(ret!=0)printf("error\n");
}


int h=0;
float total=0,x=0,y=0;
for(h=(threadnumber-1)*n+1;h<=length;h++)
{
if(h%3==0) x=y=1.0;
if(h%3==1) x=y=-1.0;
if(h%3==2) x=y=0.0;
total+=x*y;

}
pthread_join(pt[0],NULL);
pthread_join(pt[threadnumber-2],NULL);

pthread_mutex_t mutex;
pthread_mutex_init(&mutex, NULL);
pthread_mutex_lock(&mutex);
result[threadnumber-1]=total+1;



float last=0;//count the total
for(i=0;i<threadnumber;i++)
{
last+=result[i];
}
pthread_mutex_unlock(&mutex);

printf("the last result is %f \n ",last);




return (0);
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值