void update_float(float *vec1, int beg, int end, float *vec2, float fac,
                  float *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_float>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
mov    %rcx,-0x18(%rbp)
movss  %xmm0,-0x1c(%rbp)
mov    %r8,-0x28(%rbp)
mov    -0xc(%rbp),%ecx
mov    -0x8(%rbp),%rax
movslq %ecx,%rcx
shl    $0x2,%rcx
add    %rcx,%rax
mov    %rax,-0x8(%rbp)
mov    -0xc(%rbp),%ecx
mov    -0x18(%rbp),%rax
movslq %ecx,%rcx
shl    $0x2,%rcx
add    %rcx,%rax
mov    %rax,-0x18(%rbp)
mov    -0xc(%rbp),%ecx
mov    -0x28(%rbp),%rax
movslq %ecx,%rcx
shl    $0x2,%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_float+0xc4>
mov    -0x18(%rbp),%rax
mov    %rax,%rcx
add    $0x4,%rcx
mov    %rcx,-0x18(%rbp)
movss  (%rax),%xmm1
movss  -0x1c(%rbp),%xmm0
mov    -0x28(%rbp),%rax
mov    %rax,%rcx
add    $0x4,%rcx
mov    %rcx,-0x28(%rbp)
movss  (%rax),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
mov    %rax,%rcx
add    $0x4,%rcx
mov    %rcx,-0x8(%rbp)
movss  %xmm0,(%rax)
mov    -0x2c(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x2c(%rbp)
jmp    1166 <update_float+0x66>
pop    %rbp
ret

x86 assembly (clang, -O2):
<update_float>:
push   %rbp
push   %r15
push   %r14
push   %r12
push   %rbx
sub    %esi,%edx
mov    %edx,%r10d
inc    %r10d
je     1317 <update_float+0x217>
movslq %esi,%r15
lea    (%r8,%r15,4),%r9
lea    (%rcx,%r15,4),%rsi
lea    (%rdi,%r15,4),%rax
cmp    $0x7,%edx
jae    11e6 <update_float+0xe6>
lea    -0x1(%r10),%ecx
mov    %r10d,%edx
and    $0x3,%edx
je     1167 <update_float+0x67>
xor    %edi,%edi
nopl   0x0(%rax)
movss  (%r9),%xmm1
mulss  %xmm0,%xmm1
add    $0x4,%r9
addss  (%rsi),%xmm1
add    $0x4,%rsi
movss  %xmm1,(%rax)
add    $0x4,%rax
add    $0x1,%edi
cmp    %edi,%edx
jne    1140 <update_float+0x40>
sub    %edi,%r10d
cmp    $0x3,%ecx
jb     1317 <update_float+0x217>
mov    %r10d,%ecx
xor    %edx,%edx
cs nopw 0x0(%rax,%rax,1)
nop
movss  (%r9,%rdx,4),%xmm1
mulss  %xmm0,%xmm1
addss  (%rsi,%rdx,4),%xmm1
movss  %xmm1,(%rax,%rdx,4)
movss  0x4(%r9,%rdx,4),%xmm1
mulss  %xmm0,%xmm1
addss  0x4(%rsi,%rdx,4),%xmm1
movss  %xmm1,0x4(%rax,%rdx,4)
movss  0x8(%r9,%rdx,4),%xmm1
mulss  %xmm0,%xmm1
addss  0x8(%rsi,%rdx,4),%xmm1
movss  %xmm1,0x8(%rax,%rdx,4)
movss  0xc(%r9,%rdx,4),%xmm1
mulss  %xmm0,%xmm1
addss  0xc(%rsi,%rdx,4),%xmm1
movss  %xmm1,0xc(%rax,%rdx,4)
add    $0x4,%rdx
cmp    %edx,%ecx
jne    1180 <update_float+0x80>
jmp    1317 <update_float+0x217>
lea    (%r15,%rdx,1),%rbx
lea    (%rdi,%rbx,4),%r11
add    $0x4,%r11
lea    (%rcx,%rbx,4),%r14
add    $0x4,%r14
lea    (%r8,%rbx,4),%rbx
add    $0x4,%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_float+0x2e>
and    %r11b,%bl
jne    112e <update_float+0x2e>
add    $0x1,%rdx
mov    %rdx,%r11
and    $0xfffffffffffffff8,%r11
movaps %xmm0,%xmm1
shufps $0x0,%xmm0,%xmm1
lea    -0x8(%r11),%rbx
mov    %rbx,%r14
shr    $0x3,%r14
add    $0x1,%r14
test   %rbx,%rbx
je     1334 <update_float+0x234>
mov    %r14,%rbx
and    $0xfffffffffffffffe,%rbx
lea    (%r8,%r15,4),%r8
add    $0x30,%r8
lea    (%rdi,%r15,4),%r12
add    $0x30,%r12
lea    (%rcx,%r15,4),%rdi
add    $0x30,%rdi
xor    %ecx,%ecx
nopw   0x0(%rax,%rax,1)
movups -0x30(%rdi,%rcx,4),%xmm2
movups -0x20(%rdi,%rcx,4),%xmm3
movups -0x30(%r8,%rcx,4),%xmm4
movups -0x20(%r8,%rcx,4),%xmm5
mulps  %xmm1,%xmm4
addps  %xmm2,%xmm4
mulps  %xmm1,%xmm5
addps  %xmm3,%xmm5
movups %xmm4,-0x30(%r12,%rcx,4)
movups %xmm5,-0x20(%r12,%rcx,4)
movups -0x10(%rdi,%rcx,4),%xmm2
movups (%rdi,%rcx,4),%xmm3
movups -0x10(%r8,%rcx,4),%xmm4
movups (%r8,%rcx,4),%xmm5
mulps  %xmm1,%xmm4
addps  %xmm2,%xmm4
mulps  %xmm1,%xmm5
addps  %xmm3,%xmm5
movups %xmm4,-0x10(%r12,%rcx,4)
movups %xmm5,(%r12,%rcx,4)
add    $0x10,%rcx
add    $0xfffffffffffffffe,%rbx
jne    1280 <update_float+0x180>
test   $0x1,%r14b
je     1312 <update_float+0x212>
movups (%rsi,%rcx,4),%xmm2
movups 0x10(%rsi,%rcx,4),%xmm3
movups (%r9,%rcx,4),%xmm4
movups 0x10(%r9,%rcx,4),%xmm5
mulps  %xmm1,%xmm4
addps  %xmm2,%xmm4
mulps  %xmm1,%xmm5
addps  %xmm3,%xmm5
movups %xmm4,(%rax,%rcx,4)
movups %xmm5,0x10(%rax,%rcx,4)
cmp    %r11,%rdx
jne    1320 <update_float+0x220>
pop    %rbx
pop    %r12
pop    %r14
pop    %r15
pop    %rbp
ret
sub    %r11d,%r10d
lea    (%r9,%r11,4),%r9
lea    (%rax,%r11,4),%rax
lea    (%rsi,%r11,4),%rsi
jmp    112e <update_float+0x2e>
xor    %ecx,%ecx
test   $0x1,%r14b
jne    12e9 <update_float+0x1e9>
jmp    1312 <update_float+0x212>

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