simplify y

Author: alyst

src/dual.jl

@@ -803,13 +803,12 @@ end
 # creates the copy of x and applies fvalue!(values) to its values, and fpartial!(partial, ix) to its partials
 function _map_dual_components(fvalue!, fpartial!, x::AbstractArray{DT}) where DT <: Dual{<:Any, T} where T
     N = npartials(DT)
-    @assert Base.elsize(x) == sizeof(T) * (N + 1) # check that accessing x as Vector{T} is possible
     res = similar(x) # result
     t = similar(x, T) # temporary Array{T} for fvalue!/fpartial! application
-    ystride = N + 1
-    # y allows res to be accessed as Vector{T}
-    y = Base.unsafe_wrap(Array, Base.unsafe_convert(Ptr{T}, res),
-                         length(res) * ystride, own=false)
+    # y allows res to be accessed as Array{T}
+    y = reinterpret(reshape, T, res)
+    @assert size(y) == (N + 1, length(res))
+    ystride = size(y, 1)
 
     # calculate res values
     @inbounds for (j, v) in enumerate(x)