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

      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],

      vars->gpu_lambdaCoulomb,
      vars->gpu_isFraction,
      box);
  cudaDeviceSynchronize();
  checkLastErrorCUDA(__FILE__, __LINE__);

  // ReduceSum
