anchor.diff:                           iio(:,3,i)*kpr(3)
                       ! map to reciprocal unit cell
                       call mapto01(epslat,kprt)
  -                    !***call r3frac(epslat,kprt,iv)
                       ! from lattice coordinates to internal coordinates
                       ! for unshifted mesh
                       kprt(:)=(kprt(:)*div(:)*divsh-shift(:))/dble(divsh)
  
positive.code:                     nkp(:)=nint(kprt(:))
                      ! determine fractional part of rotated k-point
                      call mapto01(epslat,kprt)
                      !***call r3frac(epslat,kprt,iv)
                      ! if fractional part is present discard k-point
                      if (any(kprt.gt.epslat)) nkp(:)=-1
                      nkpid=idkp(nkp)
  
negative.code:     end do
    end subroutine mapto01
  !</sag>
  end subroutine reduk
  !EOC
  