Have interpolate_ll return very negative, but finite ll in case of ou…

…t-of-bounds
brenhinkeller · Feb 1, 2024 · aac28e7 · aac28e7
1 parent 7c7c87c
commit aac28e7
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/src/Utilities.jl b/src/Utilities.jl
@@ -120,15 +120,15 @@
     function interpolate_ll(x::AbstractVector,p::AbstractMatrix{T}) where {T<:Number}
         ll = zero(T)
         @inbounds for i ∈ eachindex(x)
-            ll += linterp_at_index(view(p,:,i), x[i], -Inf)
+            ll += linterp_at_index(view(p,:,i), x[i], -1e9)
         end
         return ll
     end
     function interpolate_ll(x::AbstractVector,ages::AbstractVector{Radiocarbon{T}}) where T
         ll = zero(T)
         @inbounds for i ∈ eachindex(x,ages)
             dist = ages[i].dist
-            ll += linterp_at_index(dist, x[i], -Inf)
+            ll += linterp_at_index(dist, x[i], -1e9)
         end
         return ll
     end