scale=-0.5*(f_c*(df_r+b*df_a)+df_c*(f_r+b*df_a));
v3_scale(&force,&(exchange->dist_ij),scale);
v3_add(&(ai->f),&(ai->f),&force);
scale=-0.5*(f_c*(df_r+b*df_a)+df_c*(f_r+b*df_a));
v3_scale(&force,&(exchange->dist_ij),scale);
v3_add(&(ai->f),&(ai->f),&force);
- v3_sub(&(aj->f),&(aj->f),&force); // dri rij = - drj rij
+ v3_scale(&force,&force,-1.0); // dri rij = - drj rij
+ v3_add(&(aj->f),&(aj->f),&force);