void __attribute__((noipa))
f_sub_double(double *restrict dest, double *restrict cond, double limit,
             double *restrict src1, double *restrict src2, unsigned int n) {
  for (unsigned int i = 0; i < n; ++i) {
    double truev = ((src1[i]) - (src2[i]));
    dest[i] = cond[i] < limit ? truev : src2[i];
  }
}


x86 assembly (clang, -O0):
<f_sub_double>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
movsd  %xmm0,-0x18(%rbp)
mov    %rdx,-0x20(%rbp)
mov    %rcx,-0x28(%rbp)
mov    %r8d,-0x2c(%rbp)
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x2c(%rbp),%eax
jae    11a7 <f_sub_double+0xa7>
mov    -0x20(%rbp),%rax
mov    -0x30(%rbp),%ecx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
mov    -0x30(%rbp),%ecx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x30(%rbp),%ecx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x18(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    1177 <f_sub_double+0x77>
movsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
jmp    1188 <f_sub_double+0x88>
mov    -0x28(%rbp),%rax
mov    -0x30(%rbp),%ecx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x40(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x30(%rbp),%ecx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1124 <f_sub_double+0x24>
pop    %rbp
ret

x86 assembly (clang, -O2):
<f_sub_double>:
test   %r8d,%r8d
je     1134 <f_sub_double+0x34>
mov    %r8d,%r9d
cmp    $0x1,%r8d
jne    1135 <f_sub_double+0x35>
xor    %eax,%eax
test   $0x1,%r9b
je     1134 <f_sub_double+0x34>
movsd  (%rcx,%rax,8),%xmm1
ucomisd (%rsi,%rax,8),%xmm0
jbe    112f <f_sub_double+0x2f>
movsd  (%rdx,%rax,8),%xmm0
subsd  %xmm1,%xmm0
movapd %xmm0,%xmm1
movsd  %xmm1,(%rdi,%rax,8)
ret
mov    %r9d,%r8d
and    $0xfffffffe,%r8d
xor    %eax,%eax
jmp    114f <f_sub_double+0x4f>
movsd  %xmm1,0x8(%rdi,%rax,8)
add    $0x2,%rax
cmp    %rax,%r8
je     1110 <f_sub_double+0x10>
movsd  (%rcx,%rax,8),%xmm1
ucomisd (%rsi,%rax,8),%xmm0
jbe    1168 <f_sub_double+0x68>
movsd  (%rdx,%rax,8),%xmm2
subsd  %xmm1,%xmm2
movapd %xmm2,%xmm1
movsd  %xmm1,(%rdi,%rax,8)
movsd  0x8(%rcx,%rax,8),%xmm1
ucomisd 0x8(%rsi,%rax,8),%xmm0
jbe    1140 <f_sub_double+0x40>
movsd  0x8(%rdx,%rax,8),%xmm2
subsd  %xmm1,%xmm2
movapd %xmm2,%xmm1
jmp    1140 <f_sub_double+0x40>

LLVM IR (-O0):
define void @f_sub_double(double* noalias %0, double* noalias %1, double %2, double* noalias %3, double* noalias %4, i32 %5) {
%7 = alloca double*
%8 = alloca double*
%9 = alloca double
%10 = alloca double*
%11 = alloca double*
%12 = alloca i32
%13 = alloca i32
%14 = alloca double
store double* %0, double** %7
store double* %1, double** %8
store double %2, double* %9
store double* %3, double** %10
store double* %4, double** %11
store i32 %5, i32* %12
store i32 0, i32* %13
br label %15
15: ; preds = %52, %6
%16 = load i32, i32* %13
%17 = load i32, i32* %12
%18 = icmp ult i32 %16, %17
br i1 %18, label %19, label %55
19: ; preds = %15
%20 = load double*, double** %10
%21 = load i32, i32* %13
%22 = zext i32 %21 to i64
%23 = getelementptr inbounds double, double* %20, i64 %22
%24 = load double, double* %23
%25 = load double*, double** %11
%26 = load i32, i32* %13
%27 = zext i32 %26 to i64
%28 = getelementptr inbounds double, double* %25, i64 %27
%29 = load double, double* %28
%30 = fsub double %24, %29
store double %30, double* %14
%31 = load double*, double** %8
%32 = load i32, i32* %13
%33 = zext i32 %32 to i64
%34 = getelementptr inbounds double, double* %31, i64 %33
%35 = load double, double* %34
%36 = load double, double* %9
%37 = fcmp olt double %35, %36
br i1 %37, label %38, label %40
38: ; preds = %19
%39 = load double, double* %14
br label %46
40: ; preds = %19
%41 = load double*, double** %11
%42 = load i32, i32* %13
%43 = zext i32 %42 to i64
%44 = getelementptr inbounds double, double* %41, i64 %43
%45 = load double, double* %44
br label %46
46: ; preds = %40, %38
%47 = phi double [ %39, %38 ], [ %45, %40 ]
%48 = load double*, double** %7
%49 = load i32, i32* %13
%50 = zext i32 %49 to i64
%51 = getelementptr inbounds double, double* %48, i64 %50
store double %47, double* %51
br label %52
52: ; preds = %46
%53 = load i32, i32* %13
%54 = add i32 %53, 1
store i32 %54, i32* %13
br label %15
55: ; preds = %15
ret void
}

LLVM IR (-O2):
define void @f_sub_double(double* noalias nocapture writeonly %0, double* noalias nocapture %1, double %2, double* noalias nocapture %3, double* noalias nocapture %4, i32 %5) {
%7 = icmp eq i32 %5, 0
br i1 %7, label %30, label %8
8: ; preds = %6
%9 = zext i32 %5 to i64
%10 = and i64 %9, 1
%11 = icmp eq i32 %5, 1
br i1 %11, label %14, label %12
12: ; preds = %8
%13 = and i64 %9, 4294967294
br label %31
14: ; preds = %56, %8
%15 = phi i64 [ 0, %8 ], [ %59, %56 ]
%16 = icmp eq i64 %10, 0
br i1 %16, label %30, label %17
17: ; preds = %14
%18 = getelementptr inbounds double, double* %4, i64 %15
%19 = load double, double* %18
%20 = getelementptr inbounds double, double* %1, i64 %15
%21 = load double, double* %20
%22 = fcmp olt double %21, %2
br i1 %22, label %23, label %27
23: ; preds = %17
%24 = getelementptr inbounds double, double* %3, i64 %15
%25 = load double, double* %24
%26 = fsub double %25, %19
br label %27
27: ; preds = %23, %17
%28 = phi double [ %26, %23 ], [ %19, %17 ]
%29 = getelementptr inbounds double, double* %0, i64 %15
store double %28, double* %29
br label %30
30: ; preds = %27, %14, %6
ret void
31: ; preds = %56, %12
%32 = phi i64 [ 0, %12 ], [ %59, %56 ]
%33 = phi i64 [ 0, %12 ], [ %60, %56 ]
%34 = getelementptr inbounds double, double* %4, i64 %32
%35 = load double, double* %34
%36 = getelementptr inbounds double, double* %1, i64 %32
%37 = load double, double* %36
%38 = fcmp olt double %37, %2
br i1 %38, label %39, label %43
39: ; preds = %31
%40 = getelementptr inbounds double, double* %3, i64 %32
%41 = load double, double* %40
%42 = fsub double %41, %35
br label %43
43: ; preds = %31, %39
%44 = phi double [ %42, %39 ], [ %35, %31 ]
%45 = getelementptr inbounds double, double* %0, i64 %32
store double %44, double* %45
%46 = or i64 %32, 1
%47 = getelementptr inbounds double, double* %4, i64 %46
%48 = load double, double* %47
%49 = getelementptr inbounds double, double* %1, i64 %46
%50 = load double, double* %49
%51 = fcmp olt double %50, %2
br i1 %51, label %52, label %56
52: ; preds = %43
%53 = getelementptr inbounds double, double* %3, i64 %46
%54 = load double, double* %53
%55 = fsub double %54, %48
br label %56
56: ; preds = %52, %43
%57 = phi double [ %55, %52 ], [ %48, %43 ]
%58 = getelementptr inbounds double, double* %0, i64 %46
store double %57, double* %58
%59 = add nuw nsw i64 %32, 2
%60 = add i64 %33, 2
%61 = icmp eq i64 %60, %13
br i1 %61, label %14, label %31
}