A Computational Realization of a Semi-Lagrangian Method for Solving the Advection Equation

Pseudocode 2

__global__  void  vec_foo  (float*  inA,  float  out,  int  size)  {

  for  (int  i=blockIdx.x*blockDim.x+threadIdx.x;

        i<size;  i+=blockDim.x*gridDim.x)

  out[i]  =  do_smb_with(inA[i]);