real *f,
real *fshift);
-enum { coultRF, coultTAB, coultTAB_TWIN, coultNR };
+enum {
+ coultRF, coultTAB, coultTAB_TWIN, coultNR
+};
p_nbk_func_ener p_nbk_c_ener[coultNR] =
-{ nbnxn_kernel_ref_rf_ener,
- nbnxn_kernel_ref_tab_ener,
- nbnxn_kernel_ref_tab_twin_ener };
+{
+ nbnxn_kernel_ref_rf_ener,
+ nbnxn_kernel_ref_tab_ener,
+ nbnxn_kernel_ref_tab_twin_ener
+};
p_nbk_func_ener p_nbk_c_energrp[coultNR] =
-{ nbnxn_kernel_ref_rf_energrp,
- nbnxn_kernel_ref_tab_energrp,
- nbnxn_kernel_ref_tab_twin_energrp};
+{
+ nbnxn_kernel_ref_rf_energrp,
+ nbnxn_kernel_ref_tab_energrp,
+ nbnxn_kernel_ref_tab_twin_energrp
+};
p_nbk_func_noener p_nbk_c_noener[coultNR] =
-{ nbnxn_kernel_ref_rf_noener,
- nbnxn_kernel_ref_tab_noener,
- nbnxn_kernel_ref_tab_twin_noener };
+{
+ nbnxn_kernel_ref_rf_noener,
+ nbnxn_kernel_ref_tab_noener,
+ nbnxn_kernel_ref_tab_twin_noener
+};
void
nbnxn_kernel_ref(const nbnxn_pairlist_set_t *nbl_list,
const nbnxn_atomdata_t *nbat,
const interaction_const_t *ic,
rvec *shift_vec,
- int force_flags,
- int clearF,
+ int force_flags,
+ int clearF,
real *fshift,
real *Vc,
real *Vvdw)
{
- int nnbl;
+ int nnbl;
nbnxn_pairlist_t **nbl;
- int coult;
- int nb;
+ int coult;
+ int nb;
nnbl = nbl_list->nnbl;
nbl = nbl_list->nbl;
}
#pragma omp parallel for schedule(static) num_threads(gmx_omp_nthreads_get(emntNonbonded))
- for(nb=0; nb<nnbl; nb++)
+ for (nb = 0; nb < nnbl; nb++)
{
nbnxn_atomdata_output_t *out;
- real *fshift_p;
+ real *fshift_p;
out = &nbat->out[nb];
if (clearF == enbvClearFYes)
{
- clear_f(nbat,nb,out->f);
+ clear_f(nbat, nb, out->f);
}
if ((force_flags & GMX_FORCE_VIRIAL) && nnbl == 1)
if (!(force_flags & GMX_FORCE_ENERGY))
{
/* Don't calculate energies */
- p_nbk_c_noener[coult](nbl[nb],nbat,
+ p_nbk_c_noener[coult](nbl[nb], nbat,
ic,
shift_vec,
out->f,
out->Vvdw[0] = 0;
out->Vc[0] = 0;
- p_nbk_c_ener[coult](nbl[nb],nbat,
+ p_nbk_c_ener[coult](nbl[nb], nbat,
ic,
shift_vec,
out->f,
/* Calculate energy group contributions */
int i;
- for(i=0; i<out->nV; i++)
+ for (i = 0; i < out->nV; i++)
{
out->Vvdw[i] = 0;
}
- for(i=0; i<out->nV; i++)
+ for (i = 0; i < out->nV; i++)
{
out->Vc[i] = 0;
}
- p_nbk_c_energrp[coult](nbl[nb],nbat,
+ p_nbk_c_energrp[coult](nbl[nb], nbat,
ic,
shift_vec,
out->f,
if (force_flags & GMX_FORCE_ENERGY)
{
- reduce_energies_over_lists(nbat,nnbl,Vvdw,Vc);
+ reduce_energies_over_lists(nbat, nnbl, Vvdw, Vc);
}
}