                        cudaMemcpyHostToDevice));
   gpuErrchk(cudaMemcpy(vars->gpu_z, coords.z, atomNumber * sizeof(double),
                        cudaMemcpyHostToDevice));
-  checkLastErrorCUDA(__FILE__, __LINE__);
 
   // Run the kernel...
   threadsPerBlock = 256;

                          cudaMemcpyDeviceToHost));
  REn = cpu_final_REn;
  LJEn = cpu_final_LJEn;

  cudaDeviceSynchronize();

  cudaFree(gpu_pair1);
  cudaFree(gpu_pair2);

      lambdaCoulomb = DeviceGetLambdaCoulomb(mA, kA, mB, kB, box,
                                             gpu_isFraction, gpu_molIndex,
                                             gpu_kindIndex, gpu_lambdaCoulomb);
      gpu_REn[threadID] = CalcCoulombGPU(distSq, kA, kB,
                                         qi_qj_fact, gpu_rCutLow[0],
                                         gpu_ewald[0], gpu_VDW_Kind[0],
                                         gpu_alpha[box],
