Fix fastmath trig functions (#497)

PainterQubits · Nov 8, 2021 · 955fc05 · 955fc05 · sostock · Nov 12, 2021
1 parent 1e3b270
commit 955fc05
Show file tree

Hide file tree

Showing 2 changed files with 7 additions and 7 deletions.
diff --git a/Project.toml b/Project.toml
@@ -1,6 +1,6 @@
 name = "Unitful"
 uuid = "1986cc42-f94f-5a68-af5c-568840ba703d"
-version = "1.9.1"
+version = "1.9.2"
 
 [deps]
 ConstructionBase = "187b0558-2788-49d3-abe0-74a17ed4e7c9"

diff --git a/src/fastmath.jl b/src/fastmath.jl
@@ -28,15 +28,17 @@ import Base.FastMath: @fastmath,
     le_fast,
     pow_fast,
     sqrt_fast,
+    cos_fast,
+    sin_fast,
+    tan_fast,
     atan_fast,
     hypot_fast,
     max_fast,
     min_fast,
     minmax_fast,
     cis_fast,
     angle_fast,
-    fast_op,
-    libm
+    fast_op
 
 sub_fast(x::Quantity{T}) where {T <: FloatTypes} = typeof(x)(neg_float_fast(x.val))
 
@@ -166,10 +168,8 @@ sqrt_fast(x::Quantity{T}) where {T <: FloatTypes} =
 for f in (:cos, :sin, :tan)
     f_fast = fast_op[f]
     @eval begin
-        $f_fast(x::DimensionlessQuantity{Float32,U}) where {U} =
-            ccall(($(string(f,"f")),libm), Float32, (Float32,), uconvert(x,NoUnits))
-        $f_fast(x::DimensionlessQuantity{Float64,U}) where {U} =
-            ccall(($(string(f)),libm), Float64, (Float64,), uconvert(x,NoUnits))
+        $(f_fast)(x::DimensionlessQuantity{<:Union{Float32,Float64}}) =
+            $(f_fast)(uconvert(NoUnits, x))
     end
 end