 
         tex_drho.normalized = false;
         tex_drho.filterMode = cudaFilterModePoint;
-        error = cudaBindTexture(0, tex_drho, d_drho, sizeof(Scalar4) * eam_data.nrho * eam_data.ntypes * eam_data.ntypes);
+        error = cudaBindTexture(0, tex_drho, d_drho,
+                                sizeof(Scalar4) * eam_data.nrho * eam_data.ntypes * eam_data.ntypes);
         if (error != cudaSuccess) return error;
 
         tex_rphi.normalized = false;
         tex_rphi.filterMode = cudaFilterModePoint;
-        error = cudaBindTexture(0, tex_rphi, d_rphi, sizeof(Scalar4) * (int) (0.5 * eam_data.nr * (eam_data.ntypes + 1) * eam_data.ntypes));
+        error = cudaBindTexture(0, tex_rphi, d_rphi,
+                                sizeof(Scalar4) * (int) (0.5 * eam_data.nr * (eam_data.ntypes + 1) * eam_data.ntypes));
         if (error != cudaSuccess) return error;
 
         tex_drphi.normalized = false;
         tex_drphi.filterMode = cudaFilterModePoint;
-        error = cudaBindTexture(0, tex_drphi, d_drphi, sizeof(Scalar4) * (int) (0.5 * eam_data.nr * (eam_data.ntypes + 1) * eam_data.ntypes));
+        error = cudaBindTexture(0, tex_drphi, d_drphi,
+                                sizeof(Scalar4) * (int) (0.5 * eam_data.nr * (eam_data.ntypes + 1) * eam_data.ntypes));
         if (error != cudaSuccess) return error;
     }
 
     pdata_pos_tex.normalized = false;
     pdata_pos_tex.filterMode = cudaFilterModePoint;
-    error = cudaBindTexture(0, pdata_pos_tex, d_pos, sizeof(Scalar4)*N);
+    error = cudaBindTexture(0, pdata_pos_tex, d_pos, sizeof(Scalar4) * N);
     if (error != cudaSuccess) return error;
 
     tex_dFdP.normalized = false;
     tex_dFdP.filterMode = cudaFilterModePoint;
-    error = cudaBindTexture(0, tex_dFdP, d_dFdP, sizeof(Scalar)*N);
+    error = cudaBindTexture(0, tex_dFdP, d_dFdP, sizeof(Scalar) * N);
     if (error != cudaSuccess) return error;
 
     // run the kernel

            return error;
    }

    if (compute_capability < 350)
    {
        tex_F.normalized = false;
        tex_F.filterMode = cudaFilterModePoint;
        error = cudaBindTexture(0, tex_F, d_F, sizeof(Scalar4) * eam_data.nrho * eam_data.ntypes);

        dim3 grid_2((int) ceil((double) N / (double) run_block_size_2), 1, 1);
        dim3 threads_2(run_block_size_2, 1, 1);

        gpu_kernel_1<0> << < grid_1, threads_1 >> > (d_force,
                d_virial, virial_pitch, N, d_pos, box, d_n_neigh, d_nlist,
                d_head_list, d_F, d_rho, d_rphi, d_dF, d_drho, d_drphi, d_dFdP);
        gpu_kernel_2<0> << < grid_2, threads_2 >> > (d_force,
                d_virial, virial_pitch, N, d_pos, box, d_n_neigh, d_nlist,
                d_head_list, d_F, d_rho, d_rphi, d_dF, d_drho, d_drphi, d_dFdP);
    }

    return cudaSuccess;
}
