t_albe_mult_params *p;
+ // set cutoff before parameters (actually only necessary for some pots)
+ if(moldyn->cutoff==0.0) {
+ printf("[albe] WARNING: no cutoff!\n");
+ return -1;
+ }
+
/* alloc mem for potential parameters */
moldyn->pot_params=malloc(sizeof(t_albe_mult_params));
if(moldyn->pot_params==NULL) {
v3_add(&(aj->f),&(aj->f),&force);
/* virial */
- virial_calc(aj,&force,&(exchange->dist_ij));
+ virial_calc(ai,&force,&(exchange->dist_ij));
#ifdef DEBUG
if(moldyn->time>DSTART&&moldyn->time<DEND) {
}
#endif
+ /* virial */
+ virial_calc(ai,&force,&dist_ij);
+
/* force contribution to atom i */
v3_scale(&force,&force,-1.0);
v3_add(&(ai->f),&(ai->f),&force);
- /* virial */
- virial_calc(ai,&force,&dist_ij);
-
/* derivative wrt k */
v3_scale(&force,&dist_ik,-1.0*dfcg); // dri rik = - drk rik
v3_scale(&tmp,&dcosdrk,fcdg);
}
#endif
+ /* virial */
+ virial_calc(ai,&force,&dist_ik);
+
/* force contribution to atom i */
v3_scale(&force,&force,-1.0);
v3_add(&(ai->f),&(ai->f),&force);
- /* virial */
- virial_calc(ai,&force,&dist_ik);
-
/* increase k counter */
exchange->kcount++;