Sheehan Olver

I wonder if we can rewrite besselj? 😉 Exp makes a 4x difference: julia> r=rand(1000000);@time AppleAccelerate.exp(r); 0.003829 seconds (8 allocations: 7.630 MB) julia> r=rand(1000000);@time exp(r); 0.016968 seconds (6 allocations: 7.630...

What's the speed if you call AppleAccellerate.exp!(vec(r))? I think modifying vec(r) also modifies r... Sent from my iPad > On 25 Sep 2015, at 12:59 AM, Richard Mikael Slevinsky [email protected]...

Yes there should. if you search for transform! you’ll see that TensorSpace has one already, but not really anything else. This is only used for ProductFun I think. > On...

Yes . We’ll have to get ApproxFun working in 0.5 first though, which last I checked is held up by a Julia type inference bug > On 25 May 2016,...

`Complex{Float64}` is not (just) a Julia type, it's stored as bits. This is why BLAS is fast for Complex numbers. But so is `Vec{2,Float64}`, with the exact same number of...

(Note that `reinterpret(UInt128,1.0+2.0im)` does not work since `Complex128` is not a `bitstype`, even though `isbits(1.0+2.0im)` is true. )

Yep, I was referring to FixedSizeArrays. It looks from that link that even @SimonDanisch agrees that StaticArrays is the way forward, and after spending a day adding special cases for...

I'd say false, as it's not ≤ all numbers in the interval

After 5 years I've forgotten why I created this issue. Let's just close it.

Just subtract out at infinity