Performance improvements for PointVector, RayVector, SubObjectIterator (

#3369)

Add type assertions to their getindex methods to improve type stability.
Similarly ensure that the return value of `size` is recognized as Tuple{Int}.

Silly microbenchmark, before:

    julia> v = PointVector{ZZRingElem}(matrix(ZZ, [1 2 3 ; ]));

    julia> @Btime prod(v)
      1.217 μs (21 allocations: 336 bytes)

After:

    julia> @Btime prod(v);
      288.449 ns (5 allocations: 80 bytes)

src/PolyhedralGeometry/iterators.jl

@@ -13,7 +13,7 @@ for (T, _t) in ((:PointVector, :point_vector), (:RayVector, :ray_vector))
     end
 
     Base.IndexStyle(::Type{<:$T}) = IndexLinear()
-    Base.getindex(po::$T, i::Base.Integer) = po.p[1, i]
+    Base.getindex(po::$T{U}, i::Base.Integer) where U = po.p[1, i]::U
 
     function Base.setindex!(po::$T, val, i::Base.Integer)
       @boundscheck 1 <= length(po) <= i
@@ -23,7 +23,7 @@ for (T, _t) in ((:PointVector, :point_vector), (:RayVector, :ray_vector))
 
     Base.firstindex(::$T) = 1
     Base.lastindex(iter::$T) = length(iter)
-    Base.size(po::$T) = (size(po.p, 2),)
+    Base.size(po::$T) = (size(po.p, 2)::Int,)
 
     coefficient_field(po::$T) = base_ring(po.p)
 
@@ -237,7 +237,7 @@ Base.IndexStyle(::Type{<:SubObjectIterator}) = IndexLinear()
 
 function Base.getindex(iter::SubObjectIterator{T}, i::Base.Integer) where T
     @boundscheck 1 <= i && i <= iter.n
-    return iter.Acc(T, iter.Obj, i; iter.options...)
+    return iter.Acc(T, iter.Obj, i; iter.options...)::T
 end
 
 Base.firstindex(::SubObjectIterator) = 1