+ if(!(atom[i].attr&ATOM_ATTR_3BP)||
+ !(btom->attr&ATOM_ATTR_3BP))
+ continue;
+
+ link_cell_neighbour_index(moldyn,
+ (btom->r.x+moldyn->dim.x/2)/lc->x,
+ (btom->r.y+moldyn->dim.y/2)/lc->y,
+ (btom->r.z+moldyn->dim.z/2)/lc->z,
+ neighbourk);
+
+ for(k=0;k<lc->countn;k++) {
+
+ thisk=&(neighbourk[k]);
+ list_reset(thisk);
+
+ if(thisk->start==NULL)
+ continue;
+
+ bck=(k<lc->dnlc)?0:1;
+
+ do {
+
+ ktom=thisk->current->data;
+
+ if(!(ktom->attr&ATOM_ATTR_3BP))
+ continue;
+
+ if(ktom==btom)
+ continue;
+
+ if(ktom==&(atom[i]))
+ continue;
+
+ moldyn->pf_func3b(moldyn,&(atom[i]),btom,ktom,bck);
+
+ } while(list_next(thisk)!=\
+ L_NO_NEXT_ELEMENT);
+
+ } while(list_next(this)!=L_NO_NEXT_ELEMENT);
+ }