Parallel Numerical Simulations of Three-Dimensional Electromagnetic Radiation with MPI-CUDA Paradigms

Algorithm 3

__global__ void caluH(arguments)

{

  = threadIdx. + blockIdx.  * blockDim.;

  = threadIdx. + blockIdx.  * blockDim.;

  = threadIdx. + blockIdx.  * blockDim.;

 if ( < _size && < _size && < _size)

 {

  int   *  _size *  _size + *  _size + ;

  double tmp = _[];

  _[] =   *  _[] −   *  ((_[ + _size] − _[]) − (_[ + 1] − _[]))/delta;

  _[] =   *  _[] +   *  (  *  _[] −   *  tmp)

 }

}

The host:

dim3 block(BlockDim[0], BlockDim[1], BlockDim[2]);

Grid[0] = /blockDim[0];

Grid[1] = /blockDim[1];

Grid[2] = /blockDim[2];

dim3 grid (Grid[0], Grid[1], Grid[2]);

caluH <<<grid, block>>> (arguments);