VPF::Пересечение массивов с разных процессоров

a11en

Дата 18.1.2011, 11:58 (ссылка)

(нет голосов)

Загрузка ...

Шустрый

Профиль
Группа: Участник
Сообщений: 64
Регистрация: 4.11.2006
Где: Ставрополь

Репутация: нет
Всего: нет

Прошу прощения если тема создана не в том разделе, но более подходящей не нашел.

В общем, есть некоторые массивы на разных процессорах, мне нужно сделать их пересечение, пытался сделать пересечение при помощи пользовательской функции и MPI_Reduce. Вроде само пересечение, если выводить результат именно из функции происходит правильно, но в результирующем векторе результат какой-то не такой.

З.Ы. Я не до конца понимаю как работает MPI_Reduce, видимо от этого все проблемы, но нужно срочно поэтому спрашиваю у вас, а сам пока буду искать описания.

Если есть более правильный способ сделать это буду благодарен. Привожу код своего "творения"

Код


#include "mpi.h"
#include<iostream>
#include<set>
#include <iterator>
#include <algorithm>

using namespace std;

void out(const char *s, const int *begin,const int *end)
{
    cout<<s;
    copy(begin,end,ostream_iterator<int>(cout," "));
    cout<<endl;
}
int *prodend;
void func(int *in, int *inout, int *len, MPI_Datatype *dtype)
{
    int i, prod[4];
    prodend=set_intersection(inout,inout+*len,in,in+*len,prod);
    *inout=*prod;
}
int main(int argc, char* argv[])
{
    MPI_Op intersection;
    MPI_Init(&argc,&argv);
    int i, size, rank,result[8],aS[4];
    MPI_Comm_rank(MPI_COMM_WORLD,&rank);
    if (rank==0)
        for (i=0;i<4;i++)
        aS[i]=1;
    
    else
        for (i=0;i<4;i++)
        aS[i]=i;
    MPI_Comm_size(MPI_COMM_WORLD,&size);
    MPI_Op_create((MPI_User_function *)func,0,&intersection);
    MPI_Reduce(&aS,&result,4,MPI_INT,intersection,0,MPI_COMM_WORLD);
    MPI_Op_free(&intersection);
    
    if (rank==0)
    {
     out("S=",result,result+4);
    }    
    MPI_Finalize();
    return 0;
}

0 Пользователей читают эту тему (0 Гостей и 0 Скрытых Пользователей)
0 Пользователей:
« Предыдущая тема \| C/C++: Общие вопросы \| Следующая тема »