Implement sincos intrinsic to fix exp. (#443)

Author: maleadt

src/device/opencl/math.jl

@@ -90,7 +90,15 @@ for gentype in generic_types
 @device_function rsqrt(x::$gentype) = @builtin_ccall("rsqrt", $gentype, ($gentype,), x)
 
 @device_override Base.sin(x::$gentype) = @builtin_ccall("sin", $gentype, ($gentype,), x)
-# sincos(x::$gentype, $gentype *cosval) = @builtin_ccall("sincos", $gentype, ($gentype, $gentype *), x, cosval)
+@device_override function Base.sincos(x::$gentype)
+    cosval = Ref{$gentype}()
+    sinval = GC.@preserve cosval begin
+        ptr = Base.unsafe_convert(Ptr{$gentype}, cosval)
+        llvm_ptr = reinterpret(LLVMPtr{$gentype, AS.Private}, ptr)
+        @builtin_ccall("sincos", $gentype, ($gentype, LLVMPtr{$gentype, AS.Private}), x, llvm_ptr)
+    end
+    return sinval, cosval[]
+end
 @device_override Base.sinh(x::$gentype) = @builtin_ccall("sinh", $gentype, ($gentype,), x)
 @device_function sinpi(x::$gentype) = @builtin_ccall("sinpi", $gentype, ($gentype,), x)
 

src/device/utils.jl

@@ -63,7 +63,7 @@ macro builtin_ccall(name, ret, argtypes, args...)
             # XXX: where does the V come from?
             "P" * ASstr * "V" * mangle(elt)
         else
-            error("Unknown type $t")
+            error("Unknown type $T")
         end
     end
 

test/device/intrinsics.jl

@@ -37,6 +37,11 @@ end
 
         end
     end
+
+    @testset "exp" begin
+        @test testf(a->exp.(a), Matrix{ComplexF32}([1.0 + 1.0im 1.0 - 1.0im; -1.0 + 1.0im -1.0 - 1.0im]))
+
+    end
 end