用Pthreads进程共享内存编程

更多的可以参考我的博客，也在陆续更新ing
http://www.hspweb.cn/

1、用VIM编写创建一个进程的代码，每隔一秒递增输出1-5。

#include <pthread.h>
#include <stdio.h>
void* fun(void* rank){
    int i;
    int my_rank=(int)rank;
    for(i=1;i<=5;i++){
        printf("Thread %d:%d\n",my_rank,i);
        sleep(1);
    }
    pthread_exit(0);
}
int main(){
  pthread_t tid1;   //the 2 thread identifier
  int num1=0;

  //create the thread
  pthread_create(&tid1,NULL,fun,(void *)num1);

  //now wait for the thread to exit
  pthread_join(tid1,NULL);
}
image

2、用VIM编写创建两个进程的代码，每隔一秒递增输出1-5。

#include <pthread.h>
#include <stdio.h>
void* fun(void* rank){
    int i;
    int my_rank=(int)rank;
    for(i=1;i<=5;i++){
        printf("Thread %d:%d\n",my_rank,i);
        sleep(1);
    }
    pthread_exit(0);
}
int main(){
  pthread_t tid1,tid2;  //the 2 thread identifier
  int num1=0,num2=1;

  //create the thread
  pthread_create(&tid1,NULL,fun,(void *)num1);
  pthread_create(&tid2,NULL,fun,(void *)num2);
  //now wait for the thread to exit
  pthread_join(tid1,NULL);
  pthread_join(tid2,NULL);
}
image

3、用VIM编写创建两个进程计算π的代码，观察竞争条件。

#include <stdio.h>
#include <stdlib.h>
#include <pthread.h>
#include <math.h>
int n,thread_count;
double sum=0.0;

void* Thread_sum(void* rank){

    long my_rank=(long)rank;
    double factor;
    long long i;
    long long my_n=n/thread_count;
    long long my_first_i=my_n*my_rank;
    long long my_last_i=my_first_i+my_n;

    if(my_first_i % 2==0 ){
        factor=1.0;
    }else{
        factor=-1.0;
    }
    for(i=my_first_i;i<my_last_i;i++,factor=-factor){
        sum+=factor/(2*i+1);
    }
    return NULL;
}

int main(int argc,char* argv[]){
    long long i;
    //scanf("%d%d",&thread_count,&n);
    thread_count=atoi(argv[1]);
    n=atoi(argv[2]);
    //printf("%d\n%d\n",thread_count,n);
    //printf("%d\n%d\n%d\n",argc,argv[0],argv[1],argv[2]);
    printf("With n=%d terms\n",n);
    pthread_t thread_handles[thread_count];
    for(i=0;i<thread_count;i++){
        pthread_create(&thread_handles[i],NULL,Thread_sum,(void *)i);
    }
    for(i=0;i<thread_count;i++){
        pthread_join(thread_handles[i],NULL);
    }
    printf("our estimate of  pi = %.15lf\n",4*sum);
    sum=0;
    pthread_create(&thread_handles[0],NULL,Thread_sum,(void *)0);
    pthread_join(thread_handles[0],NULL);
    printf("Single thread of pi = %.15lf\n",4*sum);
    printf("                 pi = %.15lf\n",4*atan(1.0));
    return 0; 
}
image

4、用忙等待互斥（严格轮换法）解决计算π的竞争条件。

#include <stdio.h>
#include <stdlib.h>
#include <pthread.h>
#include <math.h>
#include "timer.h"
int n,thread_count;
double sum=0.0;
int flag=0;
void* Thread_sum(void* rank){
    long my_rank=(long)rank;
    double factor,my_sum=0.0;
    long long i;
    long long my_n=n/thread_count;
    long long my_first_i=my_n*my_rank;
    long long my_last_i=my_first_i+my_n;

    if(my_first_i % 2==0 ){
        factor=1.0;
    }else{
        factor=-1.0;
    }
    for(i=my_first_i;i<my_last_i;i++,factor=-factor){
        my_sum+=factor/(2*i+1);
    }
    while(flag!=my_rank);
    sum+=my_sum;
    flag=(flag+1)%thread_count;
    // pthread_exit(NULL);
    return NULL;
}
int main(int argc,char* argv[]){
    double start, finish, elapsed;
    long long i;
    //scanf("%d%d",&thread_count,&n);
    thread_count=atoi(argv[1]);
    n=atoi(argv[2]);
    //printf("%d\n%d\n",thread_count,n);
    //printf("%d\n%d\n%d\n",argc,argv[0],argv[1],argv[2]);
    printf("With n=%d terms\n",n);
    pthread_t thread_handles[thread_count];
    GET_TIME(start);
    for(i=0;i<thread_count;i++){
        pthread_create(&thread_handles[i],NULL,Thread_sum,(void *)i);
    }
    for(i=0;i<thread_count;i++){
        pthread_join(thread_handles[i],NULL);
    }
    GET_TIME(finish);
    printf("    Multi-threaded estimate of  pi = %.15lf\n",4*sum);
    elapsed = finish - start;
    printf("    Elapsed time = %e seconds\n", elapsed);

    sum=0;
    GET_TIME(start);
    pthread_create(&thread_handles[0],NULL,Thread_sum,(void *)0);
    pthread_join(thread_handles[0],NULL);
    GET_TIME(finish);
    printf("    Single-threaded estimate of pi = %.15lf\n",4*sum);
    elapsed = finish - start;
    printf("    Elapsed time= %e seconds\n",elapsed);
    printf("    Math library estimate of    pi = %.15lf\n",4*atan(1.0));
    return 0;
}
image

5、用Peterson解法解决计算π的竞争条件。

#include <stdio.h>
#include <stdlib.h>
#include <pthread.h>
#include <math.h>
#include "timer.h"
int n,thread_count;
double sum=0.0;

int turn;
int interested[2];

void* Thread_sum(void* rank){
    long my_rank=(long)rank;
    double factor,my_sum=0.0;
    long long i;
    long long my_n=n/thread_count;
    long long my_first_i=my_n*my_rank;
    long long my_last_i=my_first_i+my_n;

    int other;
    other=(my_rank+1)%2;
    interested[my_rank]=1;

    if(my_first_i % 2==0 ){
        factor=1.0;
    }else{
        factor=-1.0;
    }
    for(i=my_first_i;i<my_last_i;i++,factor=-factor){
        my_sum+=factor/(2*i+1);
    }
    turn=my_rank;
    while(turn==my_rank && interested[other]==1);
    sum+=my_sum;
    interested[my_rank]=0;
    // pthread_exit(NULL);
    return NULL;
}
int main(int argc,char* argv[]){
    double start, finish, elapsed;
    long long i;
    //scanf("%d%d",&thread_count,&n);
    thread_count=atoi(argv[1]);
    n=atoi(argv[2]);
    //printf("%d\n%d\n",thread_count,n);
    //printf("%d\n%d\n%d\n",argc,argv[0],argv[1],argv[2]);
    printf("With n=%d terms\n",n);
    pthread_t thread_handles[thread_count];
    GET_TIME(start);
    for(i=0;i<thread_count;i++){
        pthread_create(&thread_handles[i],NULL,Thread_sum,(void *)i);
    }
    for(i=0;i<thread_count;i++){
        pthread_join(thread_handles[i],NULL);
    }
    GET_TIME(finish);
    printf("     Peterson estimate of       pi = %.15lf\n",4*sum);
    elapsed = finish - start;
    printf("    Elapsed time = %e seconds\n", elapsed);

    sum=0;
    GET_TIME(start);
    pthread_create(&thread_handles[0],NULL,Thread_sum,(void *)0);
    pthread_join(thread_handles[0],NULL);
    GET_TIME(finish);
    printf("    Single-threaded estimate of pi = %.15lf\n",4*sum);
    elapsed = finish - start;
    printf("    Elapsed time= %e seconds\n",elapsed);
    printf("    Math library estimate of    pi = %.15lf\n",4*atan(1.0));
    return 0;
}
image

6、用互斥量解决计算π的竞争条件。

#include <stdio.h>
#include <stdlib.h>
#include <pthread.h>
#include <math.h>
#include "timer.h"
int n,thread_count;
double sum=0.0;
pthread_mutex_t mutex;

void* Thread_sum(void* rank){
    long my_rank=(long)rank;
    double factor,my_sum=0.0;
    long long i;
    long long my_n=n/thread_count;
    long long my_first_i=my_n*my_rank;
    long long my_last_i=my_first_i+my_n;

    if(my_first_i % 2==0 ){
        factor=1.0;
    }else{
        factor=-1.0;
    }
    for(i=my_first_i;i<my_last_i;i++,factor=-factor){
        my_sum+=factor/(2*i+1);
    }
    pthread_mutex_lock(&mutex);
    sum+=my_sum;
    pthread_mutex_unlock(&mutex);

    // pthread_exit(NULL);
    return NULL;
}
int main(int argc,char* argv[]){
    double start, finish, elapsed;
    long long i;
    //scanf("%d%d",&thread_count,&n);
    thread_count=atoi(argv[1]);
    n=atoi(argv[2]);
    //printf("%d\n%d\n",thread_count,n);
    //printf("%d\n%d\n%d\n",argc,argv[0],argv[1],argv[2]);
    printf("With n=%d terms\n",n);
    pthread_t thread_handles[thread_count];
    GET_TIME(start);
     for(i=0;i<thread_count;i++){
        pthread_create(&thread_handles[i],NULL,Thread_sum,(void *)i);
    }
    for(i=0;i<thread_count;i++){
        pthread_join(thread_handles[i],NULL);
    }
    GET_TIME(finish);
    printf("                   estimate of  pi = %.15lf\n",4*sum);
    elapsed = finish - start;
    printf("    Elapsed time = %e seconds\n", elapsed);

    sum=0;
    GET_TIME(start);
    pthread_create(&thread_handles[0],NULL,Thread_sum,(void *)0);
    pthread_join(thread_handles[0],NULL);
    GET_TIME(finish);
    printf("    Single-threaded estimate of pi = %.15lf\n",4*sum);
    elapsed = finish - start;
    printf("    Elapsed time= %e seconds\n",elapsed);
    printf("    Math library estimate of    pi = %.15lf\n",4*atan(1.0));
    return 0;
}
image

7、用信号量实现互斥解决计算π的竞争条件。

#include <stdio.h>
#include <stdlib.h>
#include <pthread.h>
#include <semaphore.h>
#include <math.h>
#include "timer.h"
int n,thread_count;
double sum=0.0;

sem_t sem1,sem2;

void* Thread_sum(void* rank){
    long my_rank=(long)rank;
    double factor,my_sum=0.0;
    long long i;
    long long my_n=n/thread_count;
    long long my_first_i=my_n*my_rank;
    long long my_last_i=my_first_i+my_n;

    if(my_first_i % 2==0 ){
        factor=1.0;
    }else{
        factor=-1.0;
    }
    for(i=my_first_i;i<my_last_i;i++,factor=-factor){
        my_sum+=factor/(2*i+1);
    }
    if(my_rank==0)
        sem_wait(&sem1);
    else
        sem_wait(&sem2);
    sum+=my_sum;
    if(my_rank == 0)
        sem_post(&sem2);
    else
        sem_post(&sem1);
    // pthread_exit(NULL);
    return NULL;
}
int main(int argc,char* argv[]){
    double start, finish, elapsed;
    long long i;
    //scanf("%d%d",&thread_count,&n);
    thread_count=atoi(argv[1]);
    n=atoi(argv[2]);
    //printf("%d\n%d\n",thread_count,n);
    //printf("%d\n%d\n%d\n",argc,argv[0],argv[1],argv[2]);
    printf("With n=%d terms\n",n);
    pthread_t thread_handles[thread_count];
    GET_TIME(start);
    sem_init(&sem1,0,1);
    sem_init(&sem2,0,0);
    for(i=0;i<thread_count;i++){
        pthread_create(&thread_handles[i],NULL,Thread_sum,(void *)i);
    }
    for(i=0;i<thread_count;i++){
        pthread_join(thread_handles[i],NULL);
    }
    sem_destroy(&sem1);
    sem_destroy(&sem2);
    GET_TIME(finish);
    printf("                   estimate of  pi = %.15lf\n",4*sum);
    elapsed = finish - start;
    printf("    Elapsed time = %e seconds\n", elapsed);

    sum=0;
    GET_TIME(start);
    pthread_create(&thread_handles[0],NULL,Thread_sum,(void *)0);
    pthread_join(thread_handles[0],NULL);
    GET_TIME(finish);
    printf("    Single-threaded estimate of pi = %.15lf\n",4*sum);
    elapsed = finish - start;
    printf("    Elapsed time= %e seconds\n",elapsed);
    printf("    Math library estimate of    pi = %.15lf\n",4*atan(1.0));
    return 0;
}
image

8、修改信号量的初值，调换上述2中的线程顺序。

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

sem_t sem1,sem2;
void* fun(void* rank){
    int i;
    int my_rank=(int)rank;
    for(i=1;i<=5;i++){
      if(my_rank==0)
        sem_wait(&sem1);
      else
        sem_wait(&sem2);
        printf("Thread %d:%d\n",my_rank,i);
        sleep(1);
      if(my_rank==0)
        sem_post(&sem2);
      else
        sem_post(&sem1);
    }
    pthread_exit(0);
}
int main(){
  pthread_t tid1,tid2;  //the 2 thread identifier
  int num1=0,num2=1;

  sem_init(&sem1,0,1);
  sem_init(&sem2,0,0);
  //create the thread
  pthread_create(&tid1,NULL,fun,(void *)num1);
  pthread_create(&tid2,NULL,fun,(void *)num2);
  //now wait for the thread to exit
  pthread_join(tid1,NULL);
  pthread_join(tid2,NULL);

  sem_destroy(&sem1);
  sem_destroy(&sem2);
}
image
最后编辑于 ：
©著作权归作者所有,转载或内容合作请联系作者
  • 人面猴
    序言：七十年代末，一起剥皮案震惊了整个滨河市，随后出现的几起案子，更是在滨河造成了极大的恐慌，老刑警刘岩，带你破解...
    沈念sama阅读 212,185评论 6赞 493
  • 死咒
    序言：滨河连续发生了三起死亡事件，死亡现场离奇诡异，居然都是意外死亡，警方通过查阅死者的电脑和手机，发现死者居然都...
    沈念sama阅读 90,445评论 3赞 385
  • 救了他两次的神仙让他今天三更去死
    文/潘晓璐 我一进店门，熙熙楼的掌柜王于贵愁眉苦脸地迎上来，“玉大人，你说我怎么就摊上这事。” “怎么了？”我有些...
    开封第一讲书人阅读 157,684评论 0赞 348
  • 道士缉凶录：失踪的卖姜人
    文/不坏的土叔 我叫张陵，是天一观的道长。 经常有香客问我，道长，这世上最难降的妖魔是什么？ 我笑而不...
    开封第一讲书人阅读 56,564评论 1赞 284
  • 港岛之恋（遗憾婚礼）
    正文 为了忘掉前任，我火速办了婚礼，结果婚礼上，老公的妹妹穿的比我还像新娘。我一直安慰自己，他们只是感情好，可当我...
    茶点故事阅读 65,681评论 6赞 386
  • 恶毒庶女顶嫁案：这布局不是一般人想出来的
    文/花漫 我一把揭开白布。 她就那样静静地躺着，像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上，一...
    开封第一讲书人阅读 49,874评论 1赞 290
  • 城市分裂传说
    那天，我揣着相机与录音，去河边找鬼。 笑死，一个胖子当着我的面吹牛，可吹牛的内容都是我干的。 我是一名探鬼主播，决...
    沈念sama阅读 39,025评论 3赞 408
  • 双鸳鸯连环套：你想象不到人心有多黑
    文/苍兰香墨 我猛地睁开眼，长吁一口气：“原来是场噩梦啊……” “哼！你这毒妇竟也来了？” 一声冷哼从身侧响起，我...
    开封第一讲书人阅读 37,761评论 0赞 268
  • 万荣杀人案实录
    序言：老挝万荣一对情侣失踪，失踪者是张志新（化名）和其女友刘颖，没想到半个月后，有当地人在树林里发现了一具尸体，经...
    沈念sama阅读 44,217评论 1赞 303
  • 护林员之死
    正文 独居荒郊野岭守林人离奇死亡，尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 36,545评论 2赞 327
  • 白月光启示录
    正文 我和宋清朗相恋三年，在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 38,694评论 1赞 341
  • 活死人
    序言：一个原本活蹦乱跳的男人离奇死亡，死状恐怖，灵堂内的尸体忽然破棺而出，到底是诈尸还是另有隐情，我是刑警宁泽，带...
    沈念sama阅读 34,351评论 4赞 332
  • 日本核电站爆炸内幕
    正文 年R本政府宣布，位于F岛的核电站，受9级特大地震影响，放射性物质发生泄漏。R本人自食恶果不足惜，却给世界环境...
    茶点故事阅读 39,988评论 3赞 315
  • 男人毒药：我在死后第九天来索命
    文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹，春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 30,778评论 0赞 21
  • 一桩弑父案，背后竟有这般阴谋
    文/苍兰香墨 我抬头看了看天上的太阳。三九已至，却和暖如春，着一层夹袄步出监牢的瞬间，已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 32,007评论 1赞 266
  • 情欲美人皮
    我被黑心中介骗来泰国打工， 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留，地道东北人。 一个月前我还...
    沈念sama阅读 46,427评论 2赞 360
  • 代替公主和亲
    正文 我出身青楼，却偏偏与公主长得像，于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子，可洞房花烛夜当晚...
    茶点故事阅读 43,580评论 2赞 349