The Accelerating Adoption of Julia(lwn.net)

    if(baz()){
      foo();
    } else {
      bar(); // Thread divergence!!
    }

    if(baz()){
      pushIntoFooQueue();
    } else {
      pushIntoBarQueue(); 
      // Thread divergence, but not much of a penalty
    }

    while(fooIsNotEmpty()){ // No thread divergence at all
        task = SIMDPopFoo();
        task.execute();
    }

    while(barIsNotEmpty()){
        task = SIMDPopBar();
        task.execute();
    }

    julia> function my_kernel(xs)
             total = zero(eltype(xs))
             for x in xs
               @fastmath total += x
             end
             total
           end
    julia> @code_native debuginfo=:none my_kernel(rand(10))
    ...
    L96:
     vaddpd 8(%rcx,%rax,8), %ymm0, %ymm0
     vaddpd 40(%rcx,%rax,8), %ymm1, %ymm1
     vaddpd 72(%rcx,%rax,8), %ymm2, %ymm2
     vaddpd 104(%rcx,%rax,8), %ymm3, %ymm3
     addq $16, %rax
     cmpq %rax, %rdi
     jne L96
    ...
    

    def make_loop(f):
        @jit
        def fn(x):
            for i in range(x.shape[0]):
                 x[i] = f(x[i])
        return fn