****************************************
Short helpfile for action TOPOLOGY_MATRIX
****************************************
The input trajectory can be in any of the following formats: 

                  NODES - The list of atoms for which you would like to calculate the 
                          contact matrix. The atoms involved must be specified as a list of 
                          labels of \ref mcolv or labels of a \ref multicolvarfunction 
                          actions. If you would just like to use the atomic positions you can 
                          use a \ref DENSITY command to specify a group of atoms. 
                          Specifying your atomic positions using labels of other \ref mcolv or 
                          \ref multicolvarfunction commands is useful, however, as you can 
                          then exploit a much wider variety of functions of the contact 
                          matrix as described in \ref contactmatrix. For more information on 
                          how to specify lists of atoms see \ref Group 
                  ATOMS - . For more information on how to specify lists of atoms see 
                          \ref Group 

The following arguments are compulsory: 

      DENSITY_THRESHOLD - 
                  SIGMA - the width of the function to be used for kernel density 
                          estimation 
                 KERNEL - ( default=gaussian ) the type of kernel function to be used 

In addition you may use the following options: 

  NUMERICAL_DERIVATIVES - ( default=off ) calculate the derivatives for these 
                          quantities numerically 
                  NOPBC - ( default=off ) ignore the periodic boundary conditions 
                          when calculating distances 
                 SERIAL - ( default=off ) do the calculation in serial. Do not use 
                          MPI 
                TIMINGS - ( default=off ) output information on the timings of the 
                          various parts of the calculation 
                HIGHMEM - ( default=off ) use a more memory intensive version of this 
                          collective variable 
                 SWITCH - This keyword is used if you want to employ an alternative 
                          to the continuous switching function defined above. The 
                          following provides information on the \ref switchingfunction that are 
                          available.. You can use multiple instances of this keyword i.e. 
                          SWITCH1, SWITCH2, SWITCH3... 
                 RADIUS - . You can use multiple instances of this keyword i.e. 
                          RADIUS1, RADIUS2, RADIUS3... 
        CYLINDER_SWITCH - a switching function on the distance from the center of the 
                          cylinder. You can use multiple instances of this keyword i.e. 
                          CYLINDER_SWITCH1, CYLINDER_SWITCH2, CYLINDER_SWITCH3... 
               BIN_SIZE - the size to use for the bins. You can use multiple 
                          instances of this keyword i.e. BIN_SIZE1, BIN_SIZE2, BIN_SIZE3... 


