MPI_Reduce for array in C

Question

I have recently started trying to learn MPI to code in C. I am trying to write very small test codes to make sure I know what I'm doing as I go along. Unfortunately, I seem to be having an issue with one of them, which multiplies a matrix with a vector and outputs the resulting vector. Specifically, when I call:

MPI_Reduce(c, myc, 3, MPI_DOUBLE, MPI_SUM, 0, MPI_COMM_WORLD);

where myc is the part of the vector that is calculated by each processor, my final result is that c[i] = 0 for all i. The code that calculates myc is correct (checked using one processor and outputting myc instead of c). I assume I'm doing something very stupid here, but can't figure out what.

tpg2114 · Accepted Answer · 2011-12-21 20:04:05Z

5

So you want your answer in c and not myc?

The syntax is:

int MPI_Reduce(void *sendbuf, void *recvbuf, int count, MPI_Datatype datatype, 
           MPI_Op op, int root, MPI_Comm comm)

so you are sending c into myc with your example.

Also note that recvbuf is only valid on the root node, which in your case is node 0. All the others will have garbage in recvbuf.

answered Dec 21, 2011 at 20:04

tpg2114

15.3k6 gold badges45 silver badges58 bronze badges

Sign up to request clarification or add additional context in comments.

2 Comments

wolfPack88 Over a year ago

Sigh... I really am an idiot. Thanks for pointing out my idiocy. It works now.

tpg2114 Over a year ago

We've all been there, no worries!

Collectives™ on Stack Overflow

MPI_Reduce for array in C

1 Answer 1

2 Comments

Your Answer

Hot Network Questions

Collectives™ on Stack Overflow

1 Answer 1

2 Comments

Your Answer

Sign up or log in

Post as a guest

Related