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

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

Войти
Регистрация
Восстановить пароль
Другие темы раздела
C++ Преобразовать строку http://www.cyberforum.ru/cpp-beginners/thread1053998.html
Преобразовать строку, удалив из нее все запятые, предшествующие первой точке (точек может быть одна, несколько, или ни одной)
C++ Создание статической библиотеки из нескольких файлов с несколькими заголовочными файлами В общем, имеется несколько классов, как водится - каждый имеет .h/.cpp. Решил собрать все это дело в библиотеку, но вот беда - все равно приходится таскать помимо lib все хидеры. Как можно этого... http://www.cyberforum.ru/cpp-beginners/thread1053995.html
Ошибка в коде C++
1 Составить программу, которая описывает структуру в соответствии с вариантом и записывает введенные пользователем данные в бинарный файл (использовать потоки). Процесс введения прерывается...
C++ Ошибки в программе поиска символа в файле
в файле должен найти открывающиеся скобки ( и закрывающиеся ) и сравнить количество #include <iostream> #include <stdlib.h> #include <string> using namespace std; int main() { int k...
C++ файлы, чтение запись идет не так как нужно http://www.cyberforum.ru/cpp-beginners/thread1053941.html
и так, мне очень нужна ваша помощь! в кратце эта программа считывает файл в буфер, потом сравнивает 1 значение из буфера, и если оно подходит изменяет на нужное. и вобщем всё бы хорошо, даже...
C++ Ошибка: левый операнд должен быть левосторонним значением Ребят, проблема такая: редактирую код быстрой сортировки, выдает ошибку (см. выше), подскажите пожалуйста в чем может быть проблема int Quick() { for(int stolbec=0;stolbec!=sizeY ;stolbec++) {... подробнее

Показать сообщение отдельно
pozdyaevans
Сообщений: n/a

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

23.12.2013, 22:12. Просмотров 934. Ответов 0
Метки (Все метки)

Добрый!! Народ, помогите, пожалуйста!! имеется код программы, но она работает только для нулевого процессора. Для первого процессора она вылетает. Ошибку никак не могу найти.....
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;
}
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
 
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2017, vBulletin Solutions, Inc.
Рейтинг@Mail.ru