void gather_array_float_max(float *restrict out, const float *restrict in,
                            unsigned int n) {
  for (; n; --n) {
    float q = *in++, *p = out++;
    if (q > *p)
      *p = q;
  }
}


x86 assembly (clang, -O0):
<gather_array_float_max>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
cmpl   $0x0,-0x14(%rbp)
je     1176 <gather_array_float_max+0x76>
mov    -0x10(%rbp),%rax
mov    %rax,%rcx
add    $0x4,%rcx
mov    %rcx,-0x10(%rbp)
movss  (%rax),%xmm0
movss  %xmm0,-0x18(%rbp)
mov    -0x8(%rbp),%rax
mov    %rax,%rcx
add    $0x4,%rcx
mov    %rcx,-0x8(%rbp)
mov    %rax,-0x20(%rbp)
movss  -0x18(%rbp),%xmm0
mov    -0x20(%rbp),%rax
ucomiss (%rax),%xmm0
jbe    1163 <gather_array_float_max+0x63>
movss  -0x18(%rbp),%xmm0
mov    -0x20(%rbp),%rax
movss  %xmm0,(%rax)
jmp    1168 <gather_array_float_max+0x68>
mov    -0x14(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x14(%rbp)
jmp    110f <gather_array_float_max+0xf>
pop    %rbp
ret

x86 assembly (clang, -O2):
<gather_array_float_max>:
test   %edx,%edx
je     1208 <gather_array_float_max+0x108>
lea    -0x1(%rdx),%r8d
cmp    $0x7,%r8d
jae    111d <gather_array_float_max+0x1d>
mov    %rdi,%r11
mov    %rsi,%r10
jmp    1209 <gather_array_float_max+0x109>
add    $0x1,%r8
mov    %r8,%r9
and    $0xfffffffffffffff8,%r9
lea    (%rdi,%r9,4),%r11
lea    (%rsi,%r9,4),%r10
sub    %r9d,%edx
xor    %ecx,%ecx
jmp    114d <gather_array_float_max+0x4d>
nopw   0x0(%rax,%rax,1)
add    $0x8,%rcx
cmp    %rcx,%r9
je     1203 <gather_array_float_max+0x103>
movups (%rsi,%rcx,4),%xmm0
movups (%rdi,%rcx,4),%xmm1
ucomiss %xmm1,%xmm0
jbe    115f <gather_array_float_max+0x5f>
movss  %xmm0,(%rdi,%rcx,4)
cmpltps %xmm0,%xmm1
pextrw $0x2,%xmm1,%eax
test   $0x1,%al
je     1179 <gather_array_float_max+0x79>
movaps %xmm0,%xmm2
shufps $0x55,%xmm0,%xmm2
movss  %xmm2,0x4(%rdi,%rcx,4)
pextrw $0x4,%xmm1,%eax
test   $0x1,%al
je     118f <gather_array_float_max+0x8f>
movaps %xmm0,%xmm2
unpckhpd %xmm0,%xmm2
movss  %xmm2,0x8(%rdi,%rcx,4)
pextrw $0x6,%xmm1,%eax
test   $0x1,%al
je     11a2 <gather_array_float_max+0xa2>
shufps $0xff,%xmm0,%xmm0
movss  %xmm0,0xc(%rdi,%rcx,4)
movups 0x10(%rsi,%rcx,4),%xmm0
movups 0x10(%rdi,%rcx,4),%xmm1
ucomiss %xmm1,%xmm0
jbe    11b7 <gather_array_float_max+0xb7>
movss  %xmm0,0x10(%rdi,%rcx,4)
cmpltps %xmm0,%xmm1
pextrw $0x2,%xmm1,%eax
test   $0x1,%al
je     11d1 <gather_array_float_max+0xd1>
movaps %xmm0,%xmm2
shufps $0x55,%xmm0,%xmm2
movss  %xmm2,0x14(%rdi,%rcx,4)
pextrw $0x4,%xmm1,%eax
test   $0x1,%al
je     11e7 <gather_array_float_max+0xe7>
movaps %xmm0,%xmm2
unpckhpd %xmm0,%xmm2
movss  %xmm2,0x18(%rdi,%rcx,4)
pextrw $0x6,%xmm1,%eax
test   $0x1,%al
je     1140 <gather_array_float_max+0x40>
shufps $0xff,%xmm0,%xmm0
movss  %xmm0,0x1c(%rdi,%rcx,4)
jmp    1140 <gather_array_float_max+0x40>
cmp    %r9,%r8
jne    1209 <gather_array_float_max+0x109>
ret
mov    %edx,%ecx
xor    %edx,%edx
jmp    1218 <gather_array_float_max+0x118>
nop
add    $0x1,%rdx
cmp    %edx,%ecx
je     1208 <gather_array_float_max+0x108>
movss  (%r10,%rdx,4),%xmm0
ucomiss (%r11,%rdx,4),%xmm0
jbe    1210 <gather_array_float_max+0x110>
movss  %xmm0,(%r11,%rdx,4)
jmp    1210 <gather_array_float_max+0x110>

LLVM IR (-O0):
define void @gather_array_float_max(float* noalias %0, float* noalias %1, i32 %2) {
%4 = alloca float*
%5 = alloca float*
%6 = alloca i32
%7 = alloca float
%8 = alloca float*
store float* %0, float** %4
store float* %1, float** %5
store i32 %2, i32* %6
br label %9
9: ; preds = %26, %3
%10 = load i32, i32* %6
%11 = icmp ne i32 %10, 0
br i1 %11, label %12, label %29
12: ; preds = %9
%13 = load float*, float** %5
%14 = getelementptr inbounds float, float* %13, i32 1
store float* %14, float** %5
%15 = load float, float* %13
store float %15, float* %7
%16 = load float*, float** %4
%17 = getelementptr inbounds float, float* %16, i32 1
store float* %17, float** %4
store float* %16, float** %8
%18 = load float, float* %7
%19 = load float*, float** %8
%20 = load float, float* %19
%21 = fcmp ogt float %18, %20
br i1 %21, label %22, label %25
22: ; preds = %12
%23 = load float, float* %7
%24 = load float*, float** %8
store float %23, float* %24
br label %25
25: ; preds = %22, %12
br label %26
26: ; preds = %25
%27 = load i32, i32* %6
%28 = add i32 %27, -1
store i32 %28, i32* %6
br label %9
29: ; preds = %9
ret void
}

LLVM IR (-O2):
define void @gather_array_float_max(float* noalias nocapture %0, float* noalias nocapture %1, i32 %2) {
%4 = icmp eq i32 %2, 0
br i1 %4, label %99, label %5
5: ; preds = %3
%6 = add i32 %2, -1
%7 = zext i32 %6 to i64
%8 = add nuw nsw i64 %7, 1
%9 = icmp ult i32 %6, 7
br i1 %9, label %82, label %10
10: ; preds = %5
%11 = and i64 %8, 8589934584
%12 = getelementptr float, float* %0, i64 %11
%13 = getelementptr float, float* %1, i64 %11
%14 = trunc i64 %11 to i32
%15 = sub i32 %2, %14
br label %16
16: ; preds = %77, %10
%17 = phi i64 [ 0, %10 ], [ %78, %77 ]
%18 = getelementptr float, float* %0, i64 %17
%19 = getelementptr float, float* %1, i64 %17
%20 = bitcast float* %19 to <4 x float>*
%21 = load <4 x float>, <4 x float>* %20
%22 = getelementptr float, float* %19, i64 4
%23 = bitcast float* %22 to <4 x float>*
%24 = load <4 x float>, <4 x float>* %23
%25 = bitcast float* %18 to <4 x float>*
%26 = load <4 x float>, <4 x float>* %25
%27 = getelementptr float, float* %18, i64 4
%28 = bitcast float* %27 to <4 x float>*
%29 = load <4 x float>, <4 x float>* %28
%30 = fcmp ogt <4 x float> %21, %26
%31 = fcmp ogt <4 x float> %24, %29
%32 = extractelement <4 x i1> %30, i64 0
br i1 %32, label %33, label %35
33: ; preds = %16
%34 = extractelement <4 x float> %21, i64 0
store float %34, float* %18
br label %35
35: ; preds = %33, %16
%36 = extractelement <4 x i1> %30, i64 1
br i1 %36, label %37, label %41
37: ; preds = %35
%38 = or i64 %17, 1
%39 = getelementptr float, float* %0, i64 %38
%40 = extractelement <4 x float> %21, i64 1
store float %40, float* %39
br label %41
41: ; preds = %37, %35
%42 = extractelement <4 x i1> %30, i64 2
br i1 %42, label %43, label %47
43: ; preds = %41
%44 = or i64 %17, 2
%45 = getelementptr float, float* %0, i64 %44
%46 = extractelement <4 x float> %21, i64 2
store float %46, float* %45
br label %47
47: ; preds = %43, %41
%48 = extractelement <4 x i1> %30, i64 3
br i1 %48, label %49, label %53
49: ; preds = %47
%50 = or i64 %17, 3
%51 = getelementptr float, float* %0, i64 %50
%52 = extractelement <4 x float> %21, i64 3
store float %52, float* %51
br label %53
53: ; preds = %49, %47
%54 = extractelement <4 x i1> %31, i64 0
br i1 %54, label %55, label %59
55: ; preds = %53
%56 = or i64 %17, 4
%57 = getelementptr float, float* %0, i64 %56
%58 = extractelement <4 x float> %24, i64 0
store float %58, float* %57
br label %59
59: ; preds = %55, %53
%60 = extractelement <4 x i1> %31, i64 1
br i1 %60, label %61, label %65
61: ; preds = %59
%62 = or i64 %17, 5
%63 = getelementptr float, float* %0, i64 %62
%64 = extractelement <4 x float> %24, i64 1
store float %64, float* %63
br label %65
65: ; preds = %61, %59
%66 = extractelement <4 x i1> %31, i64 2
br i1 %66, label %67, label %71
67: ; preds = %65
%68 = or i64 %17, 6
%69 = getelementptr float, float* %0, i64 %68
%70 = extractelement <4 x float> %24, i64 2
store float %70, float* %69
br label %71
71: ; preds = %67, %65
%72 = extractelement <4 x i1> %31, i64 3
br i1 %72, label %73, label %77
73: ; preds = %71
%74 = or i64 %17, 7
%75 = getelementptr float, float* %0, i64 %74
%76 = extractelement <4 x float> %24, i64 3
store float %76, float* %75
br label %77
77: ; preds = %73, %71
%78 = add nuw i64 %17, 8
%79 = icmp eq i64 %78, %11
br i1 %79, label %80, label %16
80: ; preds = %77
%81 = icmp eq i64 %8, %11
br i1 %81, label %99, label %82
82: ; preds = %5, %80
%83 = phi float* [ %0, %5 ], [ %12, %80 ]
%84 = phi float* [ %1, %5 ], [ %13, %80 ]
%85 = phi i32 [ %2, %5 ], [ %15, %80 ]
br label %86
86: ; preds = %82, %96
%87 = phi float* [ %92, %96 ], [ %83, %82 ]
%88 = phi float* [ %90, %96 ], [ %84, %82 ]
%89 = phi i32 [ %97, %96 ], [ %85, %82 ]
%90 = getelementptr inbounds float, float* %88, i64 1
%91 = load float, float* %88
%92 = getelementptr inbounds float, float* %87, i64 1
%93 = load float, float* %87
%94 = fcmp ogt float %91, %93
br i1 %94, label %95, label %96
95: ; preds = %86
store float %91, float* %87
br label %96
96: ; preds = %95, %86
%97 = add i32 %89, -1
%98 = icmp eq i32 %97, 0
br i1 %98, label %99, label %86
99: ; preds = %96, %80, %3
ret void
}