void doubles_subtract(const double *a, const double *b, int len, double *r) {
  for (int i = 0; i < len; i++)
    r[i] = a[i] - b[i];
}


x86 assembly (clang, -O0):
<doubles_subtract>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %rcx,-0x20(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    115b <doubles_subtract+0x5b>
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    111a <doubles_subtract+0x1a>
pop    %rbp
ret

x86 assembly (clang, -O2):
<doubles_subtract>:
test   %edx,%edx
jle    12b7 <doubles_subtract+0x1b7>
mov    %edx,%r8d
cmp    $0x4,%edx
jae    11b3 <doubles_subtract+0xb3>
xor    %edx,%edx
mov    %rdx,%r9
not    %r9
add    %r8,%r9
mov    %r8,%rax
and    $0x3,%rax
je     1149 <doubles_subtract+0x49>
nopl   0x0(%rax,%rax,1)
movsd  (%rdi,%rdx,8),%xmm0
subsd  (%rsi,%rdx,8),%xmm0
movsd  %xmm0,(%rcx,%rdx,8)
add    $0x1,%rdx
add    $0xffffffffffffffff,%rax
jne    1130 <doubles_subtract+0x30>
cmp    $0x3,%r9
jb     12b7 <doubles_subtract+0x1b7>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movsd  (%rdi,%rdx,8),%xmm0
subsd  (%rsi,%rdx,8),%xmm0
movsd  %xmm0,(%rcx,%rdx,8)
movsd  0x8(%rdi,%rdx,8),%xmm0
subsd  0x8(%rsi,%rdx,8),%xmm0
movsd  %xmm0,0x8(%rcx,%rdx,8)
movsd  0x10(%rdi,%rdx,8),%xmm0
subsd  0x10(%rsi,%rdx,8),%xmm0
movsd  %xmm0,0x10(%rcx,%rdx,8)
movsd  0x18(%rdi,%rdx,8),%xmm0
subsd  0x18(%rsi,%rdx,8),%xmm0
movsd  %xmm0,0x18(%rcx,%rdx,8)
add    $0x4,%rdx
cmp    %rdx,%r8
jne    1160 <doubles_subtract+0x60>
jmp    12b7 <doubles_subtract+0x1b7>
lea    (%rcx,%r8,8),%rdx
lea    (%rdi,%r8,8),%rax
cmp    %rcx,%rax
seta   %r10b
lea    (%rsi,%r8,8),%rax
cmp    %rdi,%rdx
seta   %r11b
cmp    %rcx,%rax
seta   %al
cmp    %rsi,%rdx
seta   %r9b
xor    %edx,%edx
test   %r11b,%r10b
jne    1116 <doubles_subtract+0x16>
and    %r9b,%al
jne    1116 <doubles_subtract+0x16>
mov    %r8d,%edx
and    $0xfffffffc,%edx
lea    -0x4(%rdx),%rax
mov    %rax,%r9
shr    $0x2,%r9
add    $0x1,%r9
test   %rax,%rax
je     12b8 <doubles_subtract+0x1b8>
mov    %r9,%r10
and    $0xfffffffffffffffe,%r10
xor    %eax,%eax
cs nopw 0x0(%rax,%rax,1)
nop
movupd (%rdi,%rax,8),%xmm0
movupd 0x10(%rdi,%rax,8),%xmm1
movupd (%rsi,%rax,8),%xmm2
subpd  %xmm2,%xmm0
movupd 0x10(%rsi,%rax,8),%xmm2
subpd  %xmm2,%xmm1
movupd %xmm0,(%rcx,%rax,8)
movupd %xmm1,0x10(%rcx,%rax,8)
movupd 0x20(%rdi,%rax,8),%xmm0
movupd 0x30(%rdi,%rax,8),%xmm1
movupd 0x20(%rsi,%rax,8),%xmm2
subpd  %xmm2,%xmm0
movupd 0x30(%rsi,%rax,8),%xmm2
subpd  %xmm2,%xmm1
movupd %xmm0,0x20(%rcx,%rax,8)
movupd %xmm1,0x30(%rcx,%rax,8)
add    $0x8,%rax
add    $0xfffffffffffffffe,%r10
jne    1220 <doubles_subtract+0x120>
test   $0x1,%r9b
je     12ae <doubles_subtract+0x1ae>
movupd (%rdi,%rax,8),%xmm0
movupd 0x10(%rdi,%rax,8),%xmm1
movupd (%rsi,%rax,8),%xmm2
subpd  %xmm2,%xmm0
movupd 0x10(%rsi,%rax,8),%xmm2
subpd  %xmm2,%xmm1
movupd %xmm0,(%rcx,%rax,8)
movupd %xmm1,0x10(%rcx,%rax,8)
cmp    %r8,%rdx
jne    1116 <doubles_subtract+0x16>
ret
xor    %eax,%eax
test   $0x1,%r9b
jne    1285 <doubles_subtract+0x185>
jmp    12ae <doubles_subtract+0x1ae>

LLVM IR (-O0):
define void @doubles_subtract(double* %0, double* %1, i32 %2, double* %3) {
%5 = alloca double*
%6 = alloca double*
%7 = alloca i32
%8 = alloca double*
%9 = alloca i32
store double* %0, double** %5
store double* %1, double** %6
store i32 %2, i32* %7
store double* %3, double** %8
store i32 0, i32* %9
br label %10
10: ; preds = %30, %4
%11 = load i32, i32* %9
%12 = load i32, i32* %7
%13 = icmp slt i32 %11, %12
br i1 %13, label %14, label %33
14: ; preds = %10
%15 = load double*, double** %5
%16 = load i32, i32* %9
%17 = sext i32 %16 to i64
%18 = getelementptr inbounds double, double* %15, i64 %17
%19 = load double, double* %18
%20 = load double*, double** %6
%21 = load i32, i32* %9
%22 = sext i32 %21 to i64
%23 = getelementptr inbounds double, double* %20, i64 %22
%24 = load double, double* %23
%25 = fsub double %19, %24
%26 = load double*, double** %8
%27 = load i32, i32* %9
%28 = sext i32 %27 to i64
%29 = getelementptr inbounds double, double* %26, i64 %28
store double %25, double* %29
br label %30
30: ; preds = %14
%31 = load i32, i32* %9
%32 = add nsw i32 %31, 1
store i32 %32, i32* %9
br label %10
33: ; preds = %10
ret void
}

LLVM IR (-O2):
define void @doubles_subtract(double* nocapture %0, double* nocapture %1, i32 %2, double* nocapture writeonly %3) {
%5 = icmp sgt i32 %2, 0
br i1 %5, label %6, label %117
6: ; preds = %4
%7 = zext i32 %2 to i64
%8 = icmp ult i32 %2, 4
br i1 %8, label %96, label %9
9: ; preds = %6
%10 = getelementptr double, double* %3, i64 %7
%11 = getelementptr double, double* %0, i64 %7
%12 = getelementptr double, double* %1, i64 %7
%13 = icmp ugt double* %11, %3
%14 = icmp ugt double* %10, %0
%15 = and i1 %13, %14
%16 = icmp ugt double* %12, %3
%17 = icmp ugt double* %10, %1
%18 = and i1 %16, %17
%19 = or i1 %15, %18
br i1 %19, label %96, label %20
20: ; preds = %9
%21 = and i64 %7, 4294967292
%22 = add nsw i64 %21, -4
%23 = lshr exact i64 %22, 2
%24 = add nuw nsw i64 %23, 1
%25 = and i64 %24, 1
%26 = icmp eq i64 %22, 0
br i1 %26, label %72, label %27
27: ; preds = %20
%28 = and i64 %24, 9223372036854775806
br label %29
29: ; preds = %29, %27
%30 = phi i64 [ 0, %27 ], [ %69, %29 ]
%31 = phi i64 [ 0, %27 ], [ %70, %29 ]
%32 = getelementptr inbounds double, double* %0, i64 %30
%33 = bitcast double* %32 to <2 x double>*
%34 = load <2 x double>, <2 x double>* %33
%35 = getelementptr inbounds double, double* %32, i64 2
%36 = bitcast double* %35 to <2 x double>*
%37 = load <2 x double>, <2 x double>* %36
%38 = getelementptr inbounds double, double* %1, i64 %30
%39 = bitcast double* %38 to <2 x double>*
%40 = load <2 x double>, <2 x double>* %39
%41 = getelementptr inbounds double, double* %38, i64 2
%42 = bitcast double* %41 to <2 x double>*
%43 = load <2 x double>, <2 x double>* %42
%44 = fsub <2 x double> %34, %40
%45 = fsub <2 x double> %37, %43
%46 = getelementptr inbounds double, double* %3, i64 %30
%47 = bitcast double* %46 to <2 x double>*
store <2 x double> %44, <2 x double>* %47
%48 = getelementptr inbounds double, double* %46, i64 2
%49 = bitcast double* %48 to <2 x double>*
store <2 x double> %45, <2 x double>* %49
%50 = or i64 %30, 4
%51 = getelementptr inbounds double, double* %0, i64 %50
%52 = bitcast double* %51 to <2 x double>*
%53 = load <2 x double>, <2 x double>* %52
%54 = getelementptr inbounds double, double* %51, i64 2
%55 = bitcast double* %54 to <2 x double>*
%56 = load <2 x double>, <2 x double>* %55
%57 = getelementptr inbounds double, double* %1, i64 %50
%58 = bitcast double* %57 to <2 x double>*
%59 = load <2 x double>, <2 x double>* %58
%60 = getelementptr inbounds double, double* %57, i64 2
%61 = bitcast double* %60 to <2 x double>*
%62 = load <2 x double>, <2 x double>* %61
%63 = fsub <2 x double> %53, %59
%64 = fsub <2 x double> %56, %62
%65 = getelementptr inbounds double, double* %3, i64 %50
%66 = bitcast double* %65 to <2 x double>*
store <2 x double> %63, <2 x double>* %66
%67 = getelementptr inbounds double, double* %65, i64 2
%68 = bitcast double* %67 to <2 x double>*
store <2 x double> %64, <2 x double>* %68
%69 = add nuw i64 %30, 8
%70 = add i64 %31, 2
%71 = icmp eq i64 %70, %28
br i1 %71, label %72, label %29
72: ; preds = %29, %20
%73 = phi i64 [ 0, %20 ], [ %69, %29 ]
%74 = icmp eq i64 %25, 0
br i1 %74, label %94, label %75
75: ; preds = %72
%76 = getelementptr inbounds double, double* %0, i64 %73
%77 = bitcast double* %76 to <2 x double>*
%78 = load <2 x double>, <2 x double>* %77
%79 = getelementptr inbounds double, double* %76, i64 2
%80 = bitcast double* %79 to <2 x double>*
%81 = load <2 x double>, <2 x double>* %80
%82 = getelementptr inbounds double, double* %1, i64 %73
%83 = bitcast double* %82 to <2 x double>*
%84 = load <2 x double>, <2 x double>* %83
%85 = getelementptr inbounds double, double* %82, i64 2
%86 = bitcast double* %85 to <2 x double>*
%87 = load <2 x double>, <2 x double>* %86
%88 = fsub <2 x double> %78, %84
%89 = fsub <2 x double> %81, %87
%90 = getelementptr inbounds double, double* %3, i64 %73
%91 = bitcast double* %90 to <2 x double>*
store <2 x double> %88, <2 x double>* %91
%92 = getelementptr inbounds double, double* %90, i64 2
%93 = bitcast double* %92 to <2 x double>*
store <2 x double> %89, <2 x double>* %93
br label %94
94: ; preds = %72, %75
%95 = icmp eq i64 %21, %7
br i1 %95, label %117, label %96
96: ; preds = %9, %6, %94
%97 = phi i64 [ 0, %9 ], [ 0, %6 ], [ %21, %94 ]
%98 = xor i64 %97, -1
%99 = add nsw i64 %98, %7
%100 = and i64 %7, 3
%101 = icmp eq i64 %100, 0
br i1 %101, label %114, label %102
102: ; preds = %96, %102
%103 = phi i64 [ %111, %102 ], [ %97, %96 ]
%104 = phi i64 [ %112, %102 ], [ 0, %96 ]
%105 = getelementptr inbounds double, double* %0, i64 %103
%106 = load double, double* %105
%107 = getelementptr inbounds double, double* %1, i64 %103
%108 = load double, double* %107
%109 = fsub double %106, %108
%110 = getelementptr inbounds double, double* %3, i64 %103
store double %109, double* %110
%111 = add nuw nsw i64 %103, 1
%112 = add i64 %104, 1
%113 = icmp eq i64 %112, %100
br i1 %113, label %114, label %102
114: ; preds = %102, %96
%115 = phi i64 [ %97, %96 ], [ %111, %102 ]
%116 = icmp ult i64 %99, 3
br i1 %116, label %117, label %118
117: ; preds = %114, %118, %94, %4
ret void
118: ; preds = %114, %118
%119 = phi i64 [ %147, %118 ], [ %115, %114 ]
%120 = getelementptr inbounds double, double* %0, i64 %119
%121 = load double, double* %120
%122 = getelementptr inbounds double, double* %1, i64 %119
%123 = load double, double* %122
%124 = fsub double %121, %123
%125 = getelementptr inbounds double, double* %3, i64 %119
store double %124, double* %125
%126 = add nuw nsw i64 %119, 1
%127 = getelementptr inbounds double, double* %0, i64 %126
%128 = load double, double* %127
%129 = getelementptr inbounds double, double* %1, i64 %126
%130 = load double, double* %129
%131 = fsub double %128, %130
%132 = getelementptr inbounds double, double* %3, i64 %126
store double %131, double* %132
%133 = add nuw nsw i64 %119, 2
%134 = getelementptr inbounds double, double* %0, i64 %133
%135 = load double, double* %134
%136 = getelementptr inbounds double, double* %1, i64 %133
%137 = load double, double* %136
%138 = fsub double %135, %137
%139 = getelementptr inbounds double, double* %3, i64 %133
store double %138, double* %139
%140 = add nuw nsw i64 %119, 3
%141 = getelementptr inbounds double, double* %0, i64 %140
%142 = load double, double* %141
%143 = getelementptr inbounds double, double* %1, i64 %140
%144 = load double, double* %143
%145 = fsub double %142, %144
%146 = getelementptr inbounds double, double* %3, i64 %140
store double %145, double* %146
%147 = add nuw nsw i64 %119, 4
%148 = icmp eq i64 %147, %7
br i1 %148, label %117, label %118
}