Question

我已经编写了这两个版本的代码来计算两个数组上的dotproduct操作。每个长度为256.这是非常简单的顺序代码：

#include <stdlib.h>
#include <stdio.h>

int main(int argc, char* argv[]){
    double sum;
    double a[256], b[256];
    int n = 256, i;
    for (i=0; i<n; i++){
        a[i] = i * 0.5;
        b[i] = i * 2.0; 
    }
    sum = 0;
    for (i=1; i<=n; i++){
        sum = sum + a[i]*b[i];
    }
    printf ("sum = %f\n", sum);
}//main

答案是5559680

但并行代码：

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

#define NUMTHRDS 4

double sum;
double a[256], b[256];
int status;
int n=256;
pthread_t thds[NUMTHRDS];
pthread_mutex_t mutexsum;

void* dotprod(void *arg){
    int myid, i, my_first, my_last;
    double sum_local;

    myid = (int)arg;
    my_first = myid * n/NUMTHRDS;
    my_last = (myid + 1) * n/NUMTHRDS;

    sum_local = 0;
    for (i=my_first; i<=my_last; i++){
        sum_local = sum_local + a[i]*b[i];
    }

    pthread_mutex_lock(&mutexsum);
    sum = sum + sum_local;
    pthread_mutex_unlock(&mutexsum);

    pthread_exit((void*)0);
}//dotprod

int main(int argc, char* argv[]){
    int i;
    pthread_attr_t attr;
    for (i=0; i<n; i++){
        a[i] = i * 0.5;
        b[i] = i * 2.0;
    }

    pthread_mutex_init(&mutexsum, NULL);
    pthread_attr_init(&attr);
    pthread_attr_setdetachstate(&attr, PTHREAD_CREATE_JOINABLE);

    for (i=0; i<NUMTHRDS; i++){
        pthread_create(&thds[i], &attr, dotprod, (void*)i);
    }

    pthread_attr_destroy(&attr);

    for(i=0; i<NUMTHRDS; i++){
        pthread_join(thds[i], (void **)&status);
    }

    printf("sum = %f \n", sum);
    pthread_mutex_destroy(&mutexsum);
    pthread_exit(NULL);

    return 0;   
}//main

答案是5617024

我完全混淆了这个差异是什么？

Answer 1

关闭一个错误。

for (i=1; i<=n; i++){

for (i=0; i<n; i++) {

和

for (i=my_first; i<=my_last; i++){

for (i=my_first; i<my_last; i++){

在第一个程序中，您要添加[256]和b [256]，它不在数组的末尾。很可能这些值是0，所以你得到了正确的答案。

在第二个程序中，您计算两次数组的某些部分：64,128,192和仍然在索引256中添加。

始终检查循环的边界条件，尤其是数组访问时。

并行和顺序dotproduct程序的结果不同

1 个答案: