#include "global.h"

void UpdateNeighbourlist(x,N,size,periodic,neighbourlist,rskincut)

struct vector   *x;
int N,*neighbourlist;
struct vector   size;
struct ivector  periodic;
double   rskincut;
  

// Subroutine constructs a Verlet neighbour list in the
// following format:
//
// nngbr1          Number of neighbours for atom 1
// ingbr1,1        Index of neighbour 1 of atom 1
// ingbr2,1        Index of neighbour 2 of atom 1
// ....
// ingbrnngbr1,1   Index of neighbour nngbr1 of atom 1
// nngbr2          Number of neighbours for atom 2
// ...
//
// Note that the first index is 1, not 0 !!
{

  int nneighbourstot;


  // Stuff removed

  printf("%s %13.5f %s\n","Neighbour list update found",
       1.0*nneighbourstot/N,"neighbours per atom");
  
}