From 17cfa5dd726661fbbe5dedb77a48d0f83093b74e Mon Sep 17 00:00:00 2001
From: Chris Elrod <elrodc@gmail.com>
Date: Sun, 24 Oct 2021 20:00:04 -0400
Subject: [PATCH] Bump Static to 0.4

---
 Project.toml |  4 ++--
 src/priv.jl  | 10 +++++++---
 2 files changed, 9 insertions(+), 5 deletions(-)

diff --git a/Project.toml b/Project.toml
index 80072ec..abf8744 100644
--- a/Project.toml
+++ b/Project.toml
@@ -1,7 +1,7 @@
 name = "SLEEFPirates"
 uuid = "476501e8-09a2-5ece-8869-fb82de89a1fa"
 authors = ["chriselrod <elrodc@gmail.com>"]
-version = "0.6.27"
+version = "0.6.28"
 
 [deps]
 IfElse = "615f187c-cbe4-4ef1-ba3b-2fcf58d6d173"
@@ -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"
 
diff --git a/src/priv.jl b/src/priv.jl
index 0a3929e..0382d1b 100644
--- a/src/priv.jl
+++ b/src/priv.jl
@@ -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