                                  int      gridsize, 
                                  int      nbatch,
                                  int*     xpos,
-				 int*     ypos,
-				 int*     zpos,
+                                 int*     ypos,
+                                 int*     zpos,
                                  OutType* kernels,
                                  InType*  d_in,
                                  OutType* d_out,

        int vi_s = batch_no*nbaseline*npol+pol_no;
        int grid_s = batch_no*npol*gridsize*gridsize + pol_no*gridsize*gridsize;
        int vi = 0;
        for(vi = 0; vi < (nbaseline); vi+=npol) {

	    int xl = x[vi];
	    int yl = y[vi];

            // Determine convolution point. This is basically just an
            // optimised way to calculate.
            //int myConvU = myU - u;

						 &args[0], 2*nbaseline*sizeof(int), stream),
				BF_STATUS_INTERNAL_ERROR);
    } else {
    BF_CHECK_CUDA_EXCEPTION(cudaLaunchKernel((void*)romein_kernel<InType,OutType>,
						 grid, block,
						 &args[0], 0, stream),
				BF_STATUS_INTERNAL_ERROR);
