   CUFREE(vars->gpu_rT22);
   CUFREE(vars->gpu_rT23);
   CUFREE(vars->gpu_rT33);
-  CUFREE(vars->gpu_vT11);
-  CUFREE(vars->gpu_vT12);
-  CUFREE(vars->gpu_vT13);
-  CUFREE(vars->gpu_vT22);
-  CUFREE(vars->gpu_vT23);
-  CUFREE(vars->gpu_vT33);
   CUFREE(gpu_particleCharge);
   CUFREE(gpu_final_value);
   CUFREE(d_temp_storage);

  
  if(blockIdx.y == 0) {
    double constant_part = constVal + 1.0 / gpu_hsqrRef[imageID];
    gpu_rT11[imageID] = factor * (1.0 - 2.0 * constant_part *
                                  gpu_kxRef[imageID] * gpu_kxRef[imageID]);
    gpu_rT12[imageID] = factor * (-2.0 * constant_part *
                                  gpu_kxRef[imageID] * gpu_kyRef[imageID]);
    gpu_rT13[imageID] = factor * (-2.0 * constant_part *
                                  gpu_kxRef[imageID] * gpu_kzRef[imageID]);
    gpu_rT22[imageID] = factor * (1.0 - 2.0 * constant_part *
                                  gpu_kyRef[imageID] * gpu_kyRef[imageID]);
    gpu_rT23[imageID] = factor * (-2.0 * constant_part *
                                  gpu_kyRef[imageID] * gpu_kzRef[imageID]);
    gpu_rT33[imageID] = factor * (1.0 - 2.0 * constant_part *
                                  gpu_kzRef[imageID] * gpu_kzRef[imageID]);
  }


                             uint box)
{
  int atomNumber = currentCoords.Count();
  double *gpu_particleCharge;
  double *gpu_final_value;

