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


x86 assembly (clang, -O0):
<update>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
mov    %rcx,-0x18(%rbp)
movsd  %xmm0,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    -0xc(%rbp),%ecx
mov    -0x8(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x8(%rbp)
mov    -0xc(%rbp),%ecx
mov    -0x18(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x18(%rbp)
mov    -0xc(%rbp),%ecx
mov    -0x28(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x28(%rbp)
mov    -0x10(%rbp),%eax
sub    -0xc(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
cmpl   $0x0,-0x2c(%rbp)
je     11c4 <update+0xc4>
mov    -0x18(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x18(%rbp)
movsd  (%rax),%xmm1
movsd  -0x20(%rbp),%xmm0
mov    -0x28(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x28(%rbp)
movsd  (%rax),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x8(%rbp)
movsd  %xmm0,(%rax)
mov    -0x2c(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x2c(%rbp)
jmp    1166 <update+0x66>
pop    %rbp
ret

x86 assembly (clang, -O2):
<update>:
push   %rbp
push   %r15
push   %r14
push   %r12
push   %rbx
sub    %esi,%edx
mov    %edx,%r10d
inc    %r10d
je     1335 <update+0x235>
movslq %esi,%r15
lea    (%r8,%r15,8),%r9
lea    (%rcx,%r15,8),%rsi
lea    (%rdi,%r15,8),%rax
cmp    $0x3,%edx
jae    11e6 <update+0xe6>
lea    -0x1(%r10),%ecx
mov    %r10d,%edx
and    $0x3,%edx
je     1167 <update+0x67>
xor    %edi,%edi
nopl   0x0(%rax)
movsd  (%r9),%xmm1
mulsd  %xmm0,%xmm1
add    $0x8,%r9
addsd  (%rsi),%xmm1
add    $0x8,%rsi
movsd  %xmm1,(%rax)
add    $0x8,%rax
add    $0x1,%edi
cmp    %edi,%edx
jne    1140 <update+0x40>
sub    %edi,%r10d
cmp    $0x3,%ecx
jb     1335 <update+0x235>
mov    %r10d,%ecx
xor    %edx,%edx
cs nopw 0x0(%rax,%rax,1)
nop
movsd  (%r9,%rdx,8),%xmm1
mulsd  %xmm0,%xmm1
addsd  (%rsi,%rdx,8),%xmm1
movsd  %xmm1,(%rax,%rdx,8)
movsd  0x8(%r9,%rdx,8),%xmm1
mulsd  %xmm0,%xmm1
addsd  0x8(%rsi,%rdx,8),%xmm1
movsd  %xmm1,0x8(%rax,%rdx,8)
movsd  0x10(%r9,%rdx,8),%xmm1
mulsd  %xmm0,%xmm1
addsd  0x10(%rsi,%rdx,8),%xmm1
movsd  %xmm1,0x10(%rax,%rdx,8)
movsd  0x18(%r9,%rdx,8),%xmm1
mulsd  %xmm0,%xmm1
addsd  0x18(%rsi,%rdx,8),%xmm1
movsd  %xmm1,0x18(%rax,%rdx,8)
add    $0x4,%rdx
cmp    %edx,%ecx
jne    1180 <update+0x80>
jmp    1335 <update+0x235>
lea    (%r15,%rdx,1),%rbx
lea    (%rdi,%rbx,8),%r11
add    $0x8,%r11
lea    (%rcx,%rbx,8),%r14
add    $0x8,%r14
lea    (%r8,%rbx,8),%rbx
add    $0x8,%rbx
cmp    %r14,%rax
setb   %bpl
cmp    %r11,%rsi
setb   %r14b
cmp    %rbx,%rax
setb   %bl
cmp    %r11,%r9
setb   %r11b
test   %r14b,%bpl
jne    112e <update+0x2e>
and    %r11b,%bl
jne    112e <update+0x2e>
add    $0x1,%rdx
mov    %rdx,%r11
and    $0xfffffffffffffffc,%r11
movapd %xmm0,%xmm1
unpcklpd %xmm0,%xmm1
lea    -0x4(%r11),%rbx
mov    %rbx,%r14
shr    $0x2,%r14
add    $0x1,%r14
test   %rbx,%rbx
je     1352 <update+0x252>
mov    %r14,%rbx
and    $0xfffffffffffffffe,%rbx
lea    (%r8,%r15,8),%r8
add    $0x30,%r8
lea    (%rdi,%r15,8),%r12
add    $0x30,%r12
lea    (%rcx,%r15,8),%rdi
add    $0x30,%rdi
xor    %ecx,%ecx
nopl   0x0(%rax,%rax,1)
movupd -0x30(%rdi,%rcx,8),%xmm2
movupd -0x20(%rdi,%rcx,8),%xmm3
movupd -0x30(%r8,%rcx,8),%xmm4
movupd -0x20(%r8,%rcx,8),%xmm5
mulpd  %xmm1,%xmm4
addpd  %xmm2,%xmm4
mulpd  %xmm1,%xmm5
addpd  %xmm3,%xmm5
movupd %xmm4,-0x30(%r12,%rcx,8)
movupd %xmm5,-0x20(%r12,%rcx,8)
movupd -0x10(%rdi,%rcx,8),%xmm2
movupd (%rdi,%rcx,8),%xmm3
movupd -0x10(%r8,%rcx,8),%xmm4
movupd (%r8,%rcx,8),%xmm5
mulpd  %xmm1,%xmm4
addpd  %xmm2,%xmm4
mulpd  %xmm1,%xmm5
addpd  %xmm3,%xmm5
movupd %xmm4,-0x10(%r12,%rcx,8)
movupd %xmm5,(%r12,%rcx,8)
add    $0x8,%rcx
add    $0xfffffffffffffffe,%rbx
jne    1280 <update+0x180>
test   $0x1,%r14b
je     1330 <update+0x230>
movupd (%rsi,%rcx,8),%xmm2
movupd 0x10(%rsi,%rcx,8),%xmm3
movupd (%r9,%rcx,8),%xmm4
movupd 0x10(%r9,%rcx,8),%xmm5
mulpd  %xmm1,%xmm4
addpd  %xmm2,%xmm4
mulpd  %xmm1,%xmm5
addpd  %xmm3,%xmm5
movupd %xmm4,(%rax,%rcx,8)
movupd %xmm5,0x10(%rax,%rcx,8)
cmp    %r11,%rdx
jne    133e <update+0x23e>
pop    %rbx
pop    %r12
pop    %r14
pop    %r15
pop    %rbp
ret
sub    %r11d,%r10d
lea    (%r9,%r11,8),%r9
lea    (%rax,%r11,8),%rax
lea    (%rsi,%r11,8),%rsi
jmp    112e <update+0x2e>
xor    %ecx,%ecx
test   $0x1,%r14b
jne    12fd <update+0x1fd>
jmp    1330 <update+0x230>

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

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