void vector_substract(double *a, double *b, double *c, unsigned long n) {
  unsigned long int i = 0;
  for (i = 0; i < n; i++) {
    c[i] = a[i] - b[i];
  }
}


x86 assembly (clang, -O0):
<vector_substract>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
movq   $0x0,-0x28(%rbp)
movq   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%rax
cmp    -0x20(%rbp),%rax
jae    116a <vector_substract+0x6a>
mov    -0x8(%rbp),%rax
mov    -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x28(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x28(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x28(%rbp)
jmp    1124 <vector_substract+0x24>
pop    %rbp
ret

x86 assembly (clang, -O2):
<vector_substract>:
test   %rcx,%rcx
je     12c9 <vector_substract+0x1c9>
cmp    $0x4,%rcx
jae    11b3 <vector_substract+0xb3>
xor    %eax,%eax
mov    %rax,%r8
not    %r8
add    %rcx,%r8
mov    %rcx,%r9
and    $0x3,%r9
je     1149 <vector_substract+0x49>
nopw   0x0(%rax,%rax,1)
movsd  (%rdi,%rax,8),%xmm0
subsd  (%rsi,%rax,8),%xmm0
movsd  %xmm0,(%rdx,%rax,8)
add    $0x1,%rax
add    $0xffffffffffffffff,%r9
jne    1130 <vector_substract+0x30>
cmp    $0x3,%r8
jb     12c9 <vector_substract+0x1c9>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movsd  (%rdi,%rax,8),%xmm0
subsd  (%rsi,%rax,8),%xmm0
movsd  %xmm0,(%rdx,%rax,8)
movsd  0x8(%rdi,%rax,8),%xmm0
subsd  0x8(%rsi,%rax,8),%xmm0
movsd  %xmm0,0x8(%rdx,%rax,8)
movsd  0x10(%rdi,%rax,8),%xmm0
subsd  0x10(%rsi,%rax,8),%xmm0
movsd  %xmm0,0x10(%rdx,%rax,8)
movsd  0x18(%rdi,%rax,8),%xmm0
subsd  0x18(%rsi,%rax,8),%xmm0
movsd  %xmm0,0x18(%rdx,%rax,8)
add    $0x4,%rax
cmp    %rax,%rcx
jne    1160 <vector_substract+0x60>
jmp    12c9 <vector_substract+0x1c9>
lea    (%rdx,%rcx,8),%r9
lea    (%rdi,%rcx,8),%rax
cmp    %rdx,%rax
seta   %r10b
lea    (%rsi,%rcx,8),%rax
cmp    %rdi,%r9
seta   %r11b
cmp    %rdx,%rax
seta   %r8b
cmp    %rsi,%r9
seta   %r9b
xor    %eax,%eax
test   %r11b,%r10b
jne    1115 <vector_substract+0x15>
and    %r9b,%r8b
jne    1115 <vector_substract+0x15>
mov    %rcx,%rax
and    $0xfffffffffffffffc,%rax
lea    -0x4(%rax),%r8
mov    %r8,%r9
shr    $0x2,%r9
add    $0x1,%r9
test   %r8,%r8
je     12ca <vector_substract+0x1ca>
mov    %r9,%r10
and    $0xfffffffffffffffe,%r10
xor    %r8d,%r8d
nopl   0x0(%rax,%rax,1)
movupd (%rdi,%r8,8),%xmm0
movupd 0x10(%rdi,%r8,8),%xmm1
movupd (%rsi,%r8,8),%xmm2
subpd  %xmm2,%xmm0
movupd 0x10(%rsi,%r8,8),%xmm2
subpd  %xmm2,%xmm1
movupd %xmm0,(%rdx,%r8,8)
movupd %xmm1,0x10(%rdx,%r8,8)
movupd 0x20(%rdi,%r8,8),%xmm0
movupd 0x30(%rdi,%r8,8),%xmm1
movupd 0x20(%rsi,%r8,8),%xmm2
subpd  %xmm2,%xmm0
movupd 0x30(%rsi,%r8,8),%xmm2
subpd  %xmm2,%xmm1
movupd %xmm0,0x20(%rdx,%r8,8)
movupd %xmm1,0x30(%rdx,%r8,8)
add    $0x8,%r8
add    $0xfffffffffffffffe,%r10
jne    1220 <vector_substract+0x120>
test   $0x1,%r9b
je     12c0 <vector_substract+0x1c0>
movupd (%rdi,%r8,8),%xmm0
movupd 0x10(%rdi,%r8,8),%xmm1
movupd (%rsi,%r8,8),%xmm2
subpd  %xmm2,%xmm0
movupd 0x10(%rsi,%r8,8),%xmm2
subpd  %xmm2,%xmm1
movupd %xmm0,(%rdx,%r8,8)
movupd %xmm1,0x10(%rdx,%r8,8)
cmp    %rcx,%rax
jne    1115 <vector_substract+0x15>
ret
xor    %r8d,%r8d
test   $0x1,%r9b
jne    1291 <vector_substract+0x191>
jmp    12c0 <vector_substract+0x1c0>

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

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