inside_bz_explicit Function

    logical pure function inside_bz_explicit(this, k_vec)
        class(lattice), intent(in) :: this
        integer, intent(in) :: k_vec(sdim)

        integer :: i

        ! oles lattice is defined by four corner points and two lines
        inside_bz_explicit = .false.

        do i = 1, this%get_nsites()
            if (all(k_vec == this%get_k_vec(i))) inside_bz_explicit = .true.
        end do

    end function inside_bz_explicit