Версия для печати темы
Нажмите сюда для просмотра этой темы в оригинальном формате
Форум программистов > C/C++: Общие вопросы > Пересечение массивов с разных процессоров


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

В общем, есть некоторые массивы на разных процессорах, мне нужно сделать их пересечение, пытался сделать пересечение  при помощи пользовательской функции и 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;
}

Powered by Invision Power Board (http://www.invisionboard.com)
© Invision Power Services (http://www.invisionpower.com)