int process_2b_bonds(t_moldyn *moldyn,void *data,
int (*process)(t_moldyn *moldyn,t_atom *itom,t_atom *jtom,
void *data,u8 bc));
+int process_neighbours(t_moldyn *moldyn,void *data,t_atom *atom,
+ int (*process)(t_moldyn *moldyn,t_atom *atom,t_atom *natom,
+ void *data,u8 bc));
+
int get_line(int fd,char *line,int max);
int pair_correlation_init(t_moldyn *moldyn,double dr);