void scadd(double *vec1, int beg, int end, double fac, double *vec2) {
  int i;
  vec1 = vec1 + beg;
  vec2 = vec2 + beg;
  for (i = end - beg + 1; i; i--) {
    (*vec1++) += fac * (*vec2++);
  }
}


x86 assembly (clang, -O0):
<scadd>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
movsd  %xmm0,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    -0x8(%rbp),%rax
movslq -0xc(%rbp),%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x8(%rbp)
mov    -0x20(%rbp),%rax
movslq -0xc(%rbp),%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x20(%rbp)
mov    -0x10(%rbp),%eax
sub    -0xc(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
cmpl   $0x0,-0x24(%rbp)
je     1198 <scadd+0x98>
movsd  -0x18(%rbp),%xmm0
mov    -0x20(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x20(%rbp)
movsd  (%rax),%xmm2
mov    -0x8(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x8(%rbp)
movsd  (%rax),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax)
mov    -0x24(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x24(%rbp)
jmp    1149 <scadd+0x49>
pop    %rbp
ret

x86 assembly (clang, -O2):
<scadd>:
push   %r14
push   %rbx
sub    %esi,%edx
mov    %edx,%r8d
inc    %r8d
je     12ed <scadd+0x1ed>
movslq %esi,%r14
lea    (%rcx,%r14,8),%rax
lea    (%rdi,%r14,8),%rsi
cmp    $0x3,%edx
jb     124a <scadd+0x14a>
lea    (%r14,%rdx,1),%r9
lea    (%rcx,%r9,8),%rbx
add    $0x8,%rbx
cmp    %rbx,%rsi
jae    1147 <scadd+0x47>
lea    (%rdi,%r9,8),%rbx
add    $0x8,%rbx
cmp    %rbx,%rax
jb     124a <scadd+0x14a>
add    $0x1,%rdx
mov    %rdx,%r9
and    $0xfffffffffffffffc,%r9
movapd %xmm0,%xmm1
unpcklpd %xmm0,%xmm1
lea    -0x4(%r9),%rbx
mov    %rbx,%r10
shr    $0x2,%r10
add    $0x1,%r10
test   %rbx,%rbx
je     12f1 <scadd+0x1f1>
mov    %r10,%r11
and    $0xfffffffffffffffe,%r11
lea    (%rcx,%r14,8),%rbx
add    $0x30,%rbx
lea    (%rdi,%r14,8),%rdi
add    $0x30,%rdi
xor    %ecx,%ecx
nopl   0x0(%rax,%rax,1)
movupd -0x30(%rbx,%rcx,8),%xmm2
movupd -0x20(%rbx,%rcx,8),%xmm3
mulpd  %xmm1,%xmm2
movupd -0x30(%rdi,%rcx,8),%xmm4
addpd  %xmm2,%xmm4
movupd -0x20(%rdi,%rcx,8),%xmm2
movupd -0x10(%rdi,%rcx,8),%xmm5
movupd (%rdi,%rcx,8),%xmm6
mulpd  %xmm1,%xmm3
addpd  %xmm2,%xmm3
movupd %xmm4,-0x30(%rdi,%rcx,8)
movupd %xmm3,-0x20(%rdi,%rcx,8)
movupd -0x10(%rbx,%rcx,8),%xmm2
movupd (%rbx,%rcx,8),%xmm3
mulpd  %xmm1,%xmm2
addpd  %xmm5,%xmm2
mulpd  %xmm1,%xmm3
addpd  %xmm6,%xmm3
movupd %xmm2,-0x10(%rdi,%rcx,8)
movupd %xmm3,(%rdi,%rcx,8)
add    $0x8,%rcx
add    $0xfffffffffffffffe,%r11
jne    1190 <scadd+0x90>
test   $0x1,%r10b
je     1236 <scadd+0x136>
movupd (%rax,%rcx,8),%xmm2
movupd 0x10(%rax,%rcx,8),%xmm3
movupd (%rsi,%rcx,8),%xmm4
movupd 0x10(%rsi,%rcx,8),%xmm5
mulpd  %xmm1,%xmm2
addpd  %xmm4,%xmm2
mulpd  %xmm1,%xmm3
addpd  %xmm5,%xmm3
movupd %xmm2,(%rsi,%rcx,8)
movupd %xmm3,0x10(%rsi,%rcx,8)
cmp    %r9,%rdx
je     12ed <scadd+0x1ed>
sub    %r9d,%r8d
lea    (%rax,%r9,8),%rax
lea    (%rsi,%r9,8),%rsi
lea    -0x1(%r8),%ecx
mov    %r8d,%edx
and    $0x3,%edx
je     1282 <scadd+0x182>
xor    %edi,%edi
nopl   0x0(%rax,%rax,1)
movsd  (%rax),%xmm1
mulsd  %xmm0,%xmm1
addsd  (%rsi),%xmm1
movsd  %xmm1,(%rsi)
add    $0x8,%rsi
add    $0x8,%rax
add    $0x1,%edi
cmp    %edi,%edx
jne    1260 <scadd+0x160>
sub    %edi,%r8d
cmp    $0x3,%ecx
jb     12ed <scadd+0x1ed>
mov    %r8d,%ecx
xor    %edx,%edx
nopl   0x0(%rax)
movsd  (%rax,%rdx,8),%xmm1
mulsd  %xmm0,%xmm1
addsd  (%rsi,%rdx,8),%xmm1
movsd  %xmm1,(%rsi,%rdx,8)
movsd  0x8(%rax,%rdx,8),%xmm1
mulsd  %xmm0,%xmm1
addsd  0x8(%rsi,%rdx,8),%xmm1
movsd  %xmm1,0x8(%rsi,%rdx,8)
movsd  0x10(%rax,%rdx,8),%xmm1
mulsd  %xmm0,%xmm1
addsd  0x10(%rsi,%rdx,8),%xmm1
movsd  %xmm1,0x10(%rsi,%rdx,8)
movsd  0x18(%rax,%rdx,8),%xmm1
mulsd  %xmm0,%xmm1
addsd  0x18(%rsi,%rdx,8),%xmm1
movsd  %xmm1,0x18(%rsi,%rdx,8)
add    $0x4,%rdx
cmp    %edx,%ecx
jne    1290 <scadd+0x190>
pop    %rbx
pop    %r14
ret
xor    %ecx,%ecx
test   $0x1,%r10b
jne    1205 <scadd+0x105>
jmp    1236 <scadd+0x136>

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

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