Форум программистов, компьютерный форум CyberForum.ru

Умножение матрицы на вектор с использованием MPI - C++

Восстановить пароль Регистрация
 
pozdyaevans
Сообщений: n/a
23.12.2013, 22:12     Умножение матрицы на вектор с использованием MPI #1
Добрый!! Народ, помогите, пожалуйста!! имеется код программы, но она работает только для нулевого процессора. Для первого процессора она вылетает. Ошибку никак не могу найти.....
C++
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
#include <cstdlib>
#include <cstdio>
#include <cmath>
#include <cstring>
#include "mpi.h"
 
const int N=100;
const int a=1;
const int b=9;
 
double *x, **z, *y, *buf;
 
int rank, size;
int m;
MPI_Status status;
 
int random(int min, int max){
return rand()%(max-min+1)+min;}
 
int main(int argc, char** argv){
    MPI_Init(&argc, &argv);
    MPI_Comm_size(MPI_COMM_WORLD, &size);
    MPI_Comm_rank(MPI_COMM_WORLD, &rank);
    m=N/size;
    x = new double [N];
    z = new double* [m];
    y = new double [m];
    buf= new double [m];
    for(int i=0; i<m; i++){
        z[i] = new double [N];}
    if(rank==0){
        for(int j=0; j<N; j++){
                x[j]=random(a,b);
            }
        for(int p=1; p<size; p++){
            MPI_Send(&x, N, MPI_DOUBLE, p, 2, MPI_COMM_WORLD);
        }
        for(int p=1; p<size; p++){
            for(int i=0; i<m; i++){
                //z[i] = new double [N];
                for(int j=0; j<N; j++){
                    z[i][j]=random(a,b);
                    buf[i*N+j]=z[i][j];
                }
            }
            MPI_Send(&buf, m*N, MPI_DOUBLE, p, 1, MPI_COMM_WORLD);
        }
    }
 
    printf("1");
    if(rank!=0){
        printf("2.1");
        MPI_Recv(&x, N, MPI_DOUBLE, 0,2, MPI_COMM_WORLD, &status);
        MPI_Recv(&buf, m*N, MPI_DOUBLE, 0, 1, MPI_COMM_WORLD, &status);
        for(int i=0; i<m; i++){
            //z[i] = new double [N];
            for(int j=0; j<N; j++){
                z[i][j]=random(a,b);
                z[i][j]=buf[i*N+j];
            }
        }
        printf("2.2");
    }
    
    printf("3");
    for(int i=0; i<m; i++){
        y[i]=0.0;
        for(int j=0; j<N; j++){
        y[i]+=z[i][j]*x[j];}
        //printf("y[i]= %d\n", y[i]);
    }
 
    if(rank==0){
        FILE* fp=fopen("result.txt", "w");
        for(int p=1; p<size; p++){
            MPI_Recv(&y[p*m], m, MPI_DOUBLE, p, 1, MPI_COMM_WORLD, &status);
        }
        for(int i=0; i<N; i++){
            fprintf(fp, "y[%d]=%f\n", i, y[i]);
                //fclose(fp);
            }
        fclose(fp);
    }
    else {MPI_Send(&buf, m, MPI_DOUBLE, 0, 1, MPI_COMM_WORLD);
    }
 
    MPI_Finalize();
     return 0;
}
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
23.12.2013, 22:12     Умножение матрицы на вектор с использованием MPI
Посмотрите здесь:

C++ Умножение матрицы на вектор
Умножение матрицы на вектор C++
C++ Ленточное умножение матриц с использованием MPI
Умножение матрицы на вектор-столбец C++
Умножение матрицы на вектор C++
C++ Умножение матрицы на вектор с использованием функций
Умножение матрицы на вектор C++
C++ Умножение матрицы на вектор

Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:
После регистрации реклама в сообщениях будет скрыта и будут доступны все возможности форума.
Ответ Создать тему
Опции темы

Текущее время: 15:41. Часовой пояс GMT +3.
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2016, vBulletin Solutions, Inc.
Рейтинг@Mail.ru