anchor.diff:          ! Construct the symmetry tables.
   
           use basis, only: nbasis, basis_fns, write_basis_fn
  -        use system, only: ndim, system_type, hub_real
  +        use system, only: ndim, system_type, hub_real, hub_k
           use kpoints, only: is_reciprocal_lattice_vector
           use parallel, only: parent
           use utils, only: int_fmt
  
positive.code:                 write (6,'()')
              end if
  
          end if
  
      end subroutine init_symmetry
  
  
negative.code:         integer :: ksum(ndim)
          character(4) :: fmt1
  
          select case(system_type)
          case(hub_real)
  
              ! Symmetry currently not implemented.
              nsym = 1
  
          case(hub_k)
  
              ! Each wavevector corresponds to its own irreducible representation.
              ! The maximum system size considered will be quite small (ie <200)
              ! and the sum of any two wavevectors is another wavevector in the
              ! basis (up to a primitive reciprocal lattice vector).
              ! It is thus feasible to store the nsym^2 product table.
  
              nsym = nbasis/2
              allocate(sym_table(nsym, nsym), stat=ierr)
  