         int vi = 0;
         for(vi = 0; vi < (nbaseline); vi+=npol) {
 
-	    int xl = xdata[vi];
-	    int yl = ydata[vi];
+	    int xl = x[vi];
+	    int yl = y[vi];
 
             // Determine convolution point. This is basically just an
             // optimised way to calculate.

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

class BFromein_impl {

    return Complex<RealType>(real_res, imag_res);
}


template<typename InType, typename OutType>
__global__ void romein_kernel(int                         nbaseline,
			      int                         npol,
