(itom[i].r.x+moldyn->dim.x/2)/lc->x,
(itom[i].r.y+moldyn->dim.y/2)/lc->y,
(itom[i].r.z+moldyn->dim.z/2)/lc->z,
neighbour_i);
(itom[i].r.x+moldyn->dim.x/2)/lc->x,
(itom[i].r.y+moldyn->dim.y/2)/lc->y,
(itom[i].r.z+moldyn->dim.z/2)/lc->z,
neighbour_i);
scale=-0.5*(f_c*(df_r-b*df_a)+df_c*(f_r-b*f_a)); // - in albe formalism
#endif
v3_scale(&force,&(dist_ij),scale);
scale=-0.5*(f_c*(df_r-b*df_a)+df_c*(f_r-b*f_a)); // - in albe formalism
#endif
v3_scale(&force,&(dist_ij),scale);
/* force contribution for atom j */
v3_scale(&force,&force,-1.0); // dri rij = - drj rij
/* force contribution for atom j */
v3_scale(&force,&force,-1.0); // dri rij = - drj rij
albe_v_calc(ai,&force,&(dist_ij));
//virial_calc(ai,&force,&(dist_ij));
albe_v_calc(ai,&force,&(dist_ij));
//virial_calc(ai,&force,&(dist_ij));
- pthread_mutex_unlock(&emutex);
- pthread_mutex_lock(&(amutex[ai->tag]));
+ if(pthread_mutex_unlock(&emutex))
+ perror("[albe fast] mutex unlock (energy)\n");
+ if(pthread_mutex_lock(&(amutex[ai->tag])))
+ perror("[albe fast] mutex lock (4)\n");
albe_v_calc(ai,&force,&dist_ij);
//virial_calc(ai,&force,&dist_ij);
/* force contribution to atom i */
v3_scale(&force,&force,-1.0);
v3_add(&(ai->f),&(ai->f),&force);
albe_v_calc(ai,&force,&dist_ij);
//virial_calc(ai,&force,&dist_ij);
/* force contribution to atom i */
v3_scale(&force,&force,-1.0);
v3_add(&(ai->f),&(ai->f),&force);
albe_v_calc(ai,&force,&dist_ik);
//virial_calc(ai,&force,&dist_ik);
/* force contribution to atom i */
v3_scale(&force,&force,-1.0);
v3_add(&(ai->f),&(ai->f),&force);
albe_v_calc(ai,&force,&dist_ik);
//virial_calc(ai,&force,&dist_ik);
/* force contribution to atom i */
v3_scale(&force,&force,-1.0);
v3_add(&(ai->f),&(ai->f),&force);