anchor.diff:          integer :: ierr
   
           if (allocated(sym_table)) deallocate(sym_table, stat=ierr)
  +        if (allocated(inv_sym)) deallocate(inv_sym, stat=ierr)
   
       end subroutine end_symmetry
   
  
positive.code: 
              nsym = nbasis/2
              allocate(sym_table(nsym, nsym), stat=ierr)
  
              fmt1 = int_fmt(nsym)
  
  
negative.code:                         if (is_reciprocal_lattice_vector(ksum - basis_fns(k*2)%l)) then
                              sym_table(i,j) = k
                              sym_table(j,i) = k
                              if (k==1) then
                                  inv_sym(i) = j
                                  inv_sym(j) = i
                              end if
                              exit
                          end if
                      end do
                  end do
              end do
              if (parent) then
                  write (6,'(1X,a20,/,1X,20("-"),2/,1X,a67,2/,1X,a82,/)') &
                      "Symmetry information", &
                      "An index i refers to the wavevector of the i-th alpha spin-orbital.", &
                      "The matrix below gives the result of k_i+k_j to within a reciprocal lattice vector."
                  do i = 1, nsym
                      do j = 1, nsym
                          write (6,'('//fmt1//')', advance='no') sym_table(j,i)
                      end do
                      write (6,'()')
                  end do
                  write (6,'(/,1X,a53,/)') 'The table below gives the inverse of each wavevector.'
                  write (6,'(1X,a5,3X,a7)') 'Index','Inverse'
                  do i = 1, nsym
                      write (6,'(i4,5X,i4)') i,inv_sym(i)
                  end do
                  write (6,'()')
              end if
  
          end if
  
  