Matrix-Multiplikation unter Verwendung mehrerer threads?

Ich soll 2 multiplizieren von Matrizen unter Verwendung von Fäden. Zwei Dinge: ich bekomme immer 0 ist, wenn ich das Programm starte. Ich bekomme auch Fehler Meldung(für jeden, der es sagt "warning: passing argument 1 of 'printMatrix' from incompatible pointer type" auf den Fett gedruckten Zeilen(wo ich versuche zu drucken, die Ausgabe). Auch der Hinweis, der erste block, das ist Fett, ich, das war mein Versuch das problem zu lösen. Ich denke, ich bin in der Nähe, aber ich kann nicht. Kann jemand helfen? Vielen Dank 🙂
Die Ausgabe sieht wie folgt aus:
A=
1 4
2 5
3 6
B=
8 7 6
5 4 3
A*B=
0 0 0
0 0 0
0 0 0

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

#define M 3
#define K 2
#define N 3

struct v
{
    int i; //row
    int j; //column
};

int A[M][K] = {{1,4},{2,5},{3,6}};
int B[K][N] = {{8,7,6},{5,4,3}};
int C[M][N];

void *workerThread(void *data)
{
    int i=((struct v*)data)->i;
    int j=((struct v*)data)->j;
    int accumulator = 0;

/*this is where you should calculate the assigned Cell. You will need to use the row(i) of
            A and column[j] of B. Accumulate the result in accumulator */
    **int k;
    for(k=0; k<k; k++)
    {
            accumulator = A[i][k]*B[k][j];
    }
    C[i][j]=accumulator;
    pthread_exit(NULL);**
 }

 void printMatrix(int *matrixIn, int rows, int columns)
{
    int *matrix = matrixIn;
    int i,j;
    for (i=0;i<rows;i++)
    {
}

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

    pthread_t threads[M*N];
    int i,j;
    int counter = 0;
    int numThreadsCreated = 0;

    /*the following 5 lines demonstrates how to create 1 thread to calculate C[0][0], you
            will need to create a loop for all of C's cells*/
    struct v *data = (struct v *)malloc(sizeof(struct v));
    data->i = 0; //assign the row of C for thread to calculate
    data->j = 0; //assign the column of C for thread to calculate
    pthread_create(&threads[0], NULL, workerThread, data);
    numThreadsCreated++;

    /*wait for all the threads to finish before printing out the matrices*/
    for(j=0; j < numThreadsCreated; j++)
    {
            pthread_join( threads[j], NULL);
    }

    printf("A=\n");
    **printMatrix(A,3,2);**
    printf("B=\n");
    **printMatrix(B,2,3);**
    printf("A*B=\n");
    **printMatrix(C,M,N);**
    pthread_exit(NULL);
}
Diese Zeile ist fehlerhaft--->for(k=0; k<k; k++)!
Wie so? Dies war die Linie, die ich wurde empfohlen von jemand anderem. @shekharsuman
Überprüfen Sie meine Antwort für mehr Klarheit und das gleiche versuchen und dann weitere Kommentar über dein problem!!!
Es sei denn, Ihr M und N sind riesig, das ist nicht zu viel Leistung gewinnen. Sie sollten die Suche "parallel-matrix-Multiplikation" um zu sehen, verschiedene alternative Ansätze, um das problem beim Umgang mit großen Matrizen.

InformationsquelleAutor Jay | 2014-10-21

  1. 1

    Ihrem Programm zu haben scheint, implementiert eine falsche Codierung Algorithmus zur matrix-Multiplikation.

    Folgende Stück code, das scheint absurd :-

    for(k=0; k<k; k++)
{
        accumulator = A[i][k]*B[k][j];
}
C[i][j]=accumulator;    //  Not a code for matrix multiplication...

    Sollten Sie etwas umsetzen, wie :-

    for(i=0;i<M;i++){
 for(j=0;j<N;j++){
  accumulator=0;
   for(int something=0;something<K;something++){
    accumulator=accumulator+A[i][something]*B[something][j];
   }
  C[i][j]=accumulator;
  accumulator=0;
 }
}
    Ich habe gerade versucht dieses. Ich ersetzte meinen vorherigen code mit der von Ihnen empfohlenen. Es gibt mir diese Fehlermeldung ("Fehler: 'k" schwarzarbeit (erste Benutzung in dieser Funktion)" & ("Fehler: (jeder nicht deklarierte Bezeichner wird nur einmal berichtet, für jede Funktion erscheint in") @shekharsuman
    Sie müssen die Kapital-ODER Blockschrift K anstelle von Kleinbuchstaben k!!! Sie haben #define K=3. Wie kann es Fehler geben? Bitte überprüfen Sie Ihre Fälle!
    Die Fehlermeldungen beziehen sich auf diese Zeile: for (int irgendwas = 0; etwas<k; irgendwas++)
    Gut, es ist immer noch, dass Sie mir die selben 3 Fehlermeldungen für den Druck der Matrizen, und auch eine andere Fehlermeldung für die Letzte Zeile(end bracket), die sagt: Fehler: expected declaration oder Erklärung am Ende der Eingabe
    Es hat geklappt! ... Ich denke, gibt es mir immer noch, die Fehlermeldung, sondern es multipliziert sich.

    InformationsquelleAutor Am_I_Helpful

  2. 0

    FWIW... für eine kleine matrix, die Kosten der ab-threads, und was-nicht-ist wahrscheinlich, die meinen, das ist alles eine Verschwendung von Zeit... und für eine große matrix, die Sie zu sein scheinen, beginnend M*N threads !!

    Da können wir davon ausgehen, dass dies compute-bound, dort ist nicht viel Punkt ab mehr threads als CPUs zu verwenden. Dann wird jeder thread muss wählen Sie den nächsten result-Element zu berechnen, und zu stoppen, wenn das Ergebnis voll ist. Eine Möglichkeit, dies zu tun wäre, um zu übergeben "control" - Struktur, um alle pthreads, entlang der Linien von:

      struct control
  {
     pthread_mutex_t   mutex ;
     int  q ;
  } ;

    und dann die pthread-Funktion loop:

      while (1)
    {
       int q, i, j ;

       pthread_mutex_lock(&ctrl->mutex) ;
       q = ctrl->q ;
       ctrl->q += 1 ;
       pthread_mutex_unlock(&ctrl->mutex) ;

       if (q >= (N * M))
         break ;

       i = q / K ;
       j = q % K ;

       ....
    } ;

    Eine Atomare fetch und hinzufügen 1 q würde den job tun, wie gut. (NB: es wird vorausgesetzt, dass N*M + thread-count <= INT_MAX !)

    Könnte man, natürlich, teilen Sie die N * M Ergebnisse statisch über die Anzahl der pthreads Sie starten, so dass jedem Durchgang mit einer anderen anfänglichen q und die Anzahl der Ergebnis posten, für die es zuständig ist (denken Sie an das Geschäft mit der Rest nach der Division N * M durch die Anzahl der pthreads !). Dies vermeidet alle brauchen Sie für die Pthread-Interaktion -- also kein mutex und keine atomics erforderlich.

    InformationsquelleAutor

  3. 0
    #include<stdio.h>
#include<stdlib.h>
#include<pthread.h>
#include<unistd.h>

pthread_mutex_t lock;

typedef struct
{

    int mat1;

    int mat2;
}address;

int matrix1[3][3],matrix2[3][3],result[3][3];

void *matrixOne()

{

    int i,j,data=0;

    for(i=0;i<3;i++)

    {

        for(j=0;j<3;j++)

        {

            matrix1[i][j]=data+1;

            //         sleep(2);

            data=data+1;

        }

    }

    //printf matrix

    printf("\n****matrix 1 created By thread-1***\n\n");

    for(i=0;i<3;i++)
    {
        printf("\n");
        for(j=0;j<3;j++)

        {

            printf("Element matrix[%d][%d] of matrixOne = %d\n",i,j,matrix1[i][j]);

            sleep(1);

        }
    }

    printf("\n**********1st matrix*********\n");

    for(i=0;i<3;i++)
    {
        printf("\n");
        for(j=0;j<3;j++)

        {

            printf(" %d ",matrix1[i][j]);

        }
    }



}
void *matrixTwo()

{

    int i,j,data=9;

    for(i=0;i<3;i++)

    {

        for(j=0;j<3;j++)

        {

            matrix2[i][j]=data;

            //         sleep(2);

            data=data-1;

        }

    }

    //printf matrix

    printf("\n\t\t\t\t\t****matrix 2 created By thread-2***\n");

    for(i=0;i<3;i++)
    {
        printf("\n");
        for(j=0;j<3;j++)

        {

            printf("\t\t\t\t\tElement matrix[%d][%d] of matrixTwo = %d\n",i,j,matrix2[i][j]);

            sleep(1);

        }
    }

    printf("\n**********2nd matrix*********\n");

    for(i=0;i<3;i++)
    {
        printf("\n");
        for(j=0;j<3;j++)

        {

            printf(" %d ",matrix2[i][j]);

        }
    }

    printf("\n");

}

void *multiply(void *data)

{
    pthread_mutex_lock(&lock);
    int i,j,k,sum=0;

    address *data1=(address*)data;

    for(i=0;i<1;i++)

    {

        for(j=0;j<3;j++)

        {

            sum=0;
            data1->mat2=0;
            for(k=0;k<3;k++)

            {

                //printf("%d\n",data1->mat1[i][k]);

                //             sum=sum+matrix1[i][k]*matrix2[k][j];

                sum=sum+(matrix1[data1->mat1][k])*(matrix2[data1->mat2][j]);

                ++(data1->mat2);

            }

            result[data1->mat1][j]=sum;
            //data1->mat2=0;

        }



    }

    printf("\n**********thread %d multiplication of matrix*********\n",data1->mat1);

    for(i=0;i<1;i++)
    {
        printf("\n");
        for(j=0;j<3;j++)

        {

            printf(" %d ",result[data1->mat1][j]);

        }
    }

    printf("\n");


    pthread_mutex_unlock(&lock);
}



void main()

{      

    int i,j;

    //     address add[3];


    address *add = (address*) malloc(3 * sizeof(address));


    pthread_t thread1,thread2,thread3[3];


    pthread_create(&thread1,NULL,matrixOne,NULL);

    //pthread_join(thread1,NULL);

    pthread_create(&thread2,NULL,matrixTwo,NULL);

    pthread_join(thread1,NULL);

    pthread_join(thread2,NULL);

    for(i=0;i<3;i++)

    {
        // add[i].mat1=*matrix1+i*3*4;

        // add[i].mat2=*matrix2+i*4;

        add[i].mat1=i;

        add[i].mat2=i;



        pthread_create(&thread3[i],NULL,multiply,&add[i]);

    }

    for(i=0;i<3;i++)

    {

        pthread_join(thread3[i],NULL);

    }

    printf("\n\n********** multiplication of matrix*********\n");

    for(i=0;i<3;i++)
    {
        printf("\n");
        for(j=0;j<3;j++)

        {

            printf(" %d ",result[i][j]);

        }
    }


    printf("\n\n");


}
    gcc Programmname.c -lpthread

    InformationsquelleAutor Sandeep Kumar

Schreibe einen Kommentar