/********************************************************
% 
% Written by:
% -- 
% John L. Weatherwax                2006-08-08
% 
% email: wax@alum.mit.edu
% 
% Please send comments and especially bug reports to the
% above email address.
% 
%-----

*/

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include "mpi.h"
#include "global_grid.h"
#include "sim_consts.h"
#include "createSpaceForForcesToReceive.h"

/* */ 
void createSpaceForForcesToReceive(void){

  /* create NW array */ 
  procInfo.forcesOnExportedPart[NW] = (double*) malloc( 2*procInfo.numToExport[NW]*sizeof(double) );  
  if( !procInfo.forcesOnExportedPart[NW] ){ fprintf(stderr,"No memory!\n"); MPI_Abort(MPI_COMM_WORLD,1); }

  /* create N  array */
  procInfo.forcesOnExportedPart[N]  =  (double*) malloc( 2*procInfo.numToExport[N]*sizeof(double) ); 
  if( !procInfo.forcesOnExportedPart[N] ){ fprintf(stderr,"No memory!\n"); MPI_Abort(MPI_COMM_WORLD,1); } 

  /* create NE array */ 
  procInfo.forcesOnExportedPart[NE] =  (double*) malloc( 2*procInfo.numToExport[NE]*sizeof(double) );  
  if( !procInfo.forcesOnExportedPart[NE] ){ fprintf(stderr,"No memory!\n"); MPI_Abort(MPI_COMM_WORLD,1); }

  /* create W  array */ 
  procInfo.forcesOnExportedPart[W]  =  (double*) malloc( 2*procInfo.numToExport[W]*sizeof(double) );  
  if( !procInfo.forcesOnExportedPart[W] ){ fprintf(stderr,"No memory!\n"); MPI_Abort(MPI_COMM_WORLD,1); }

  /* create I  array ... not needed I don't send to myself ... */ 
  procInfo.forcesOnExportedPart[I] = NULL;

  /* create E  array */ 
  procInfo.forcesOnExportedPart[E]  =  (double*) malloc( 2*procInfo.numToExport[E]*sizeof(double) );  
  if( !procInfo.forcesOnExportedPart[E] ){ fprintf(stderr,"No memory!\n"); MPI_Abort(MPI_COMM_WORLD,1); }

  /* create SW array */ 
  procInfo.forcesOnExportedPart[SW] =  (double*) malloc( 2*procInfo.numToExport[SW]*sizeof(double) );  
  if( !procInfo.forcesOnExportedPart[SW] ){ fprintf(stderr,"No memory!\n"); MPI_Abort(MPI_COMM_WORLD,1); }

  /* create S  array */
  procInfo.forcesOnExportedPart[S]  =  (double*) malloc( 2*procInfo.numToExport[S]*sizeof(double) );  
  if( !procInfo.forcesOnExportedPart[S] ){ fprintf(stderr,"No memory!\n"); MPI_Abort(MPI_COMM_WORLD,1); }

  /* create SE array */ 
  procInfo.forcesOnExportedPart[SE] =  (double*) malloc( 2*procInfo.numToExport[SE]*sizeof(double) );  
  if( !procInfo.forcesOnExportedPart[SE] ){ fprintf(stderr,"No memory!\n"); MPI_Abort(MPI_COMM_WORLD,1); }

}