Remove cache structs and use Tuples, allocating only when necessary (#21)

Author: DanielVandH

Project.toml

@@ -1,7 +1,7 @@
 name = "AdaptivePredicates"
 uuid = "35492f91-a3bd-45ad-95db-fcad7dcfedb7"
 authors = ["Valentin Churavy <v.churavy@gmail.com> and Daniel VandenHeuvel <danj.vandenheuvel@gmail.com"]
-version = "1.1.1"
+version = "1.2.0"
 
 [compat]
 julia = "1.6"

src/AdaptivePredicates.jl

@@ -1,7 +1,12 @@
 module AdaptivePredicates
 
+@static if isdefined(Base, :Memory)
+    const Vec{T} = Memory{T}
+else
+    const Vec{T} = Vector{T}
+end
+
 include("init.jl")
-include("caches.jl")
 include("macros.jl")
 include("arithmetic.jl")
 include("predicates.jl")
@@ -10,4 +15,8 @@ export orient2, orient3, incircle, insphere
 export orient2p, orient3p, incirclep, inspherep
 export orient2fast, orient3fast, incirclefast, inspherefast
 
+if VERSION ≥ v"1.11.0-DEV.469"
+    eval(Meta.parse("public orient3adapt_cache, incircleadapt_cache, insphereexact_cache"))
+end
+
 end # module

src/arithmetic.jl

@@ -1,6 +1,6 @@
-@inline setindex!!(tup::Tuple, value, index) = @inbounds Base.setindex(tup, value, index)
+@inline setindex!!(tup::Tuple, value, index) = @inbounds index > length(tup) ? tup : Base.setindex(tup, value, index)
 @inline setindex!!(vec::AbstractVector, value, index) = @inbounds Base.setindex!(vec, value, index)
-@inline safe_getindex(e, eindex, elen) = eindex ≤ elen ? @inbounds(e[eindex]) : zero(eltype(e)) # Shewchuk's code is relying on undefined behaviour from out-of-bounds access where we call this. We need to be careful.
+@inline safe_getindex(e, eindex, elen) = (eindex ≤ elen && eindex ≤ length(e)) ? @inbounds(e[eindex]) : zero(eltype(e)) # Shewchuk's code is relying on undefined behaviour from out-of-bounds access where we call this. We need to be careful.
 
 function grow_expansion(elen, e, b, h)
     @inbounds begin
@@ -197,8 +197,8 @@ function fast_expansion_sum_zeroelim(elen, e, flen, f, h)
         fnow = f[1]
         eindex = findex = 1
         if (fnow > enow) == (fnow > -enow)
-            Q = enow
-            eindex += 1
+            Q = enow 
+            eindex += 1 
             enow = safe_getindex(e, eindex, elen)
         else
             Q = fnow