al_exchange_dim.c File Reference

Fill the ghost boundaries along selected dimensions. More...

#include "al_hidden.h"

Include dependency graph for al_exchange_dim.c:

Go to the source code of this file.

Functions
int	AL_Exchange_dim (char *buf, int *dims, int sz_ptr)

Variables
SZ *	sz_stack [AL_MAX_ARRAYS]
int	stack_ptr [AL_MAX_ARRAYS]

Detailed Description

Fill the ghost boundaries along selected dimensions.

Fill the ghost boundaries along selected dimensions

Author

A. Malagoli (University of Chicago)

Date

Jul 17, 1999

Definition in file al_exchange_dim.c.

Function Documentation

int AL_Exchange_dim	(	char *	buf,
		int *	dims,
		int	sz_ptr
	)

Fill the ghost boundaries along selected dimensions

Parameters

[in]	buf	pointer to buffer
[in]	dims	if dims[i]=0, do not perform the exchange in this dimension (array if int)
[in]	sz_ptr	integer pointer to the distributed array descriptor

Definition at line 24 of file al_exchange_dim.c.

{
  register int nd;
  int myrank, nproc;
  int ndim, gp, nleft, nright, tag1, tag2;
  int sendb, recvb;
  MPI_Datatype itype;
  MPI_Comm comm;
  MPI_Status status;
  SZ *s;

  /* DIAGNOSTICS
    Check that sz_ptr points to an allocated SZ
  */
  if( stack_ptr[sz_ptr] == AL_STACK_FREE){
    printf("AL_Decompose: wrong SZ pointer\n");
  }

  s = sz_stack[sz_ptr];

  myrank = s->rank;
  nproc = s->size;
  comm = s->comm;
  ndim = s->ndim;

  for(nd=0;nd<ndim;nd++){
    gp = s->bg[nd];
    
    /* If gp=0, do nothing */
    if( gp > 0 && dims[nd] != 0 ){
      nleft = s->left[nd];
      nright = s->right[nd];
      itype = s->type_rl[nd];
      tag1 = s->tag1[nd];

      sendb = s->sendb1[nd];
      recvb = s->recvb1[nd];

      MPI_Sendrecv(&buf[sendb], 1, itype, nleft, tag1,
           &buf[recvb], 1, itype, nright,tag1,
           comm, &status);

      nleft = s->left[nd];
      nright = s->right[nd];
      itype = s->type_lr[nd];
      tag2 = s->tag2[nd];


      sendb = s->sendb2[nd];
      recvb = s->recvb2[nd];

      MPI_Sendrecv(&buf[sendb], 1, itype, nright, tag2,
           &buf[recvb], 1, itype, nleft,tag2,
           comm, &status);
    }
  }

  /* DIAGNOSTICS */
#ifdef DEBUG
  if(myrank==0) printf("AL_Exchange: filled ghost regions\n"); 
#endif

  return (int) AL_SUCCESS;
}

Here is the caller graph for this function:

Variable Documentation

int stack_ptr[AL_MAX_ARRAYS]

Definition at line 26 of file al_szptr_.c.

SZ* sz_stack[AL_MAX_ARRAYS]

Definition at line 18 of file al_szptr_.c.