anchor.diff:    !!<summary>Lists the HNFs for the given symmetry group.</summary>
     !!<parameter name="A" regular="true">The parent lattice for the
     !!system.</parameter>
  -  !!<parameter name="basis" regular="true">The atomic basis.</parameter>
     !!<parameter name="n_" regular="true">The determinant of the desired
     !!HNF.</parameter>
     !!<parameter name="eps_" regular="true">The floating point
  
positive.code:           if (count_i ==0) then
               count_i = count_i + 1
               max_rmin = this_rmin
               temp_hnfs(:,:,count_i) = unq_hnfs(:,:,j)
            else if (abs(this_rmin-max_rmin)<eps) then
               count_i = count_i + 1
               temp_hnfs(:,:,count_i) = unq_hnfs(:,:,j)             
            else if (this_rmin > max_rmin) then
               count_i = 1
               temp_hnfs = 0
               max_rmin = this_rmin
               temp_hnfs(:,:,count_i) = unq_hnfs(:,:,j)             
            end if
         end do
  
  
negative.code:        eps = eps_
      end if
  
      allocate(d(3,1))
      d = (/0.0_dp,0.0_dp,0.0_dp/)
      hnfs = 0
      n_hnfs = 0
      r_mins = 0
  
      call get_dvector_permutations(A,d,dRPList,LatDim,eps)
      call matrix_inverse(A,inv_lat)
  
      do i=n_min,n_max
         call get_all_HNFs(i,all_hnfs)
  