anchor.diff:  
               nsym = nbasis/2
               allocate(sym_table(nsym, nsym), stat=ierr)
  +            allocate(inv_sym(nsym), stat=ierr)
   
               fmt1 = int_fmt(nsym)
   
  
positive.code:         integer :: ierr
  
          if (allocated(sym_table)) deallocate(sym_table, stat=ierr)
  
      end subroutine end_symmetry
  
  
negative.code: ! sym_table(i,j) = k means that k_i + k_j = k_k to within a primitive reciprocal lattice vector.
  integer, allocatable :: sym_table(:,:) ! (nsym, nsym)
  
  ! inv_sym(i) = j means that k_i + k_j = 0 (ie k_j is the inverse of k_i).
  integer, allocatable :: inv_sym(:) ! nsym
  
  contains
  
      subroutine init_symmetry()
  
          ! Construct the symmetry tables.
  
          use basis, only: nbasis, basis_fns
          use system, only: ndim, system_type, hub_real
  