MPI multidimensional arrays

[Carlis]

Hi
Im passing some multidimensional al arrays over MPI. However, Im struggling with some hangups.

Essentially I do

static cont int NK=16
static cont int NT=400
static cont int NK=S=4
static cont int NK=IR=2

double SIGMA [NK][NK][NT][IR];

...
cout<< MPI_RANK<< " begin "<<endl;
MPI_Bcast(&SIGMA, NK*NK*NT*S*IR, MPI_DOUBLE, 0, MPI_COMM_WORLD);
cout<< MPI_RANK<< " end "<<endl;
...

Now, I have done this before and never had any trouble. But now, some of my processes get stuck. They do print 'begin', but not 'end'.
Also, some of the processes even finish and quit entirely, while others just wait forever.
Is it something wrong with this MPI call?

Best
//
Johan

 

[TheRyuu]

static cont int NK=16
static cont int NT=400
static cont int NK=S=4
static cont int NK=IR=2

Please use code tags. Is this a typo or are you actually declaring NK several times over?

 

[Carlis]

Oh thats a typo. IR=2 and nothing else. Sorry

 

[Schmide]

Looking it up, it seems the default stack size for MPI is 128k, your arrays are 6.6MB. You are probably overflowing the stack. If you are working with such large arrays, you should probably look at passing by reference and or better memory management techniques.

Edit: Never mind. I see you pass by reference. I blame the lack of code tags.