Performance of Float32 exponential

[baggepinnen]

See some discussion in this answer and thread

Fast logsumexp

LLVM. GCC will vectorize log/exp/sin/etc with the appropriate optimization flags, but LLVM needs those in addition to -fveclib=SVML or some other vector library. using SIMDPirates, SLEEFPirates, LoopVectorization function logsumexp_simdpirates!(w::Vector{T},we) where T offset = maximum(w) N = length(w) sl = SIMDPirates.Vec{4,T}((0.,0.,0.,0.)) @inbounds @simd for i = 1:4:N wl = SIMDPirates.vload(SIMDPirates.Vec{4,T}, w, i) @pirate wel = SLEEFPirates.exp(…