Skip to content

Commit

Permalink
Bump Static to 0.4
Browse files Browse the repository at this point in the history
  • Loading branch information
chriselrod committed Oct 25, 2021
1 parent fae0a11 commit 17cfa5d
Show file tree
Hide file tree
Showing 2 changed files with 9 additions and 5 deletions.
4 changes: 2 additions & 2 deletions Project.toml
Original file line number Diff line number Diff line change
@@ -1,7 +1,7 @@
name = "SLEEFPirates"
uuid = "476501e8-09a2-5ece-8869-fb82de89a1fa"
authors = ["chriselrod <[email protected]>"]
version = "0.6.27"
version = "0.6.28"

[deps]
IfElse = "615f187c-cbe4-4ef1-ba3b-2fcf58d6d173"
Expand All @@ -10,7 +10,7 @@ VectorizationBase = "3d5dd08c-fd9d-11e8-17fa-ed2836048c2f"

[compat]
IfElse = "0.1"
Static = "0.2, 0.3"
Static = "0.2, 0.3, 0.4"
VectorizationBase = "0.20.12, 0.21"
julia = "1.5"

Expand Down
10 changes: 7 additions & 3 deletions src/priv.jl
Original file line number Diff line number Diff line change
Expand Up @@ -426,15 +426,19 @@ end
@inline function logk(d::V) where (V <: FloatType)
T = eltype(d)
m, e = splitfloat(d)

# m = significand(d); e = exponent(d)
x = ddiv(dsub2(m, T(1.0)), dadd2(T(1.0), m))

# @show
x2 = dsqu(x)

t = logk_kernel(x2)

s = dmul(MDLN2(T), e)
s = Double(0.6931471803691238 * e, 1.9082149292705877e-10 * e)
# s = dmul(MDLN2(T), e)
# @show Float64(s) s
s = dadd(s, scale(x, T(2.0)))
# @show s
s = dadd(s, dmul(dmul(x2, x), t))
# @show s
return s
end

0 comments on commit 17cfa5d

Please sign in to comment.