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


x86 assembly (clang, -O0):
<scadd_mixed>:
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    $0x2,%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     119c <scadd_mixed+0x9c>
movsd  -0x18(%rbp),%xmm0
mov    -0x20(%rbp),%rax
mov    %rax,%rcx
add    $0x4,%rcx
mov    %rcx,-0x20(%rbp)
movss  (%rax),%xmm1
cvtss2sd %xmm1,%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_mixed+0x49>
pop    %rbp
ret

x86 assembly (clang, -O2):
<scadd_mixed>:
push   %r14
push   %rbx
sub    %esi,%edx
mov    %edx,%r8d
inc    %r8d
je     1250 <scadd_mixed+0x150>
movslq %esi,%r14
lea    (%rcx,%r14,4),%r11
lea    (%rdi,%r14,8),%rsi
cmp    $0x3,%edx
jb     1222 <scadd_mixed+0x122>
add    $0x1,%rdx
mov    %rdx,%r9
and    $0xfffffffffffffffc,%r9
movapd %xmm0,%xmm1
unpcklpd %xmm0,%xmm1
lea    -0x4(%r9),%rax
mov    %rax,%r10
shr    $0x2,%r10
add    $0x1,%r10
test   %rax,%rax
je     1254 <scadd_mixed+0x154>
mov    %r10,%rax
and    $0xfffffffffffffffe,%rax
lea    (%rcx,%r14,4),%rbx
add    $0x18,%rbx
lea    (%rdi,%r14,8),%rdi
add    $0x30,%rdi
xor    %ecx,%ecx
nopl   0x0(%rax)
cvtps2pd -0x18(%rbx,%rcx,4),%xmm2
cvtps2pd -0x10(%rbx,%rcx,4),%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)
cvtps2pd -0x8(%rbx,%rcx,4),%xmm2
cvtps2pd (%rbx,%rcx,4),%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,%rax
jne    1170 <scadd_mixed+0x70>
test   $0x1,%r10b
je     1212 <scadd_mixed+0x112>
cvtps2pd (%r11,%rcx,4),%xmm2
cvtps2pd 0x8(%r11,%rcx,4),%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     1250 <scadd_mixed+0x150>
sub    %r9d,%r8d
lea    (%r11,%r9,4),%r11
lea    (%rsi,%r9,8),%rsi
mov    %r8d,%eax
xor    %ecx,%ecx
nopw   0x0(%rax,%rax,1)
movss  (%r11,%rcx,4),%xmm1
cvtss2sd %xmm1,%xmm1
mulsd  %xmm0,%xmm1
addsd  (%rsi,%rcx,8),%xmm1
movsd  %xmm1,(%rsi,%rcx,8)
add    $0x1,%rcx
cmp    %ecx,%eax
jne    1230 <scadd_mixed+0x130>
pop    %rbx
pop    %r14
ret
xor    %ecx,%ecx
test   $0x1,%r10b
jne    11e1 <scadd_mixed+0xe1>
jmp    1212 <scadd_mixed+0x112>

LLVM IR (-O0):
define void @scadd_mixed(double* %0, i32 %1, i32 %2, double %3, float* %4) {
%6 = alloca double*
%7 = alloca i32
%8 = alloca i32
%9 = alloca double
%10 = alloca float*
%11 = alloca i32
store double* %0, double** %6
store i32 %1, i32* %7
store i32 %2, i32* %8
store double %3, double* %9
store float* %4, float** %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 float*, float** %10
%17 = load i32, i32* %7
%18 = sext i32 %17 to i64
%19 = getelementptr inbounds float, float* %16, i64 %18
store float* %19, float** %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 = %37, %5
%25 = load i32, i32* %11
%26 = icmp ne i32 %25, 0
br i1 %26, label %27, label %40
27: ; preds = %24
%28 = load double, double* %9
%29 = load float*, float** %10
%30 = getelementptr inbounds float, float* %29, i32 1
store float* %30, float** %10
%31 = load float, float* %29
%32 = fpext float %31 to double
%33 = load double*, double** %6
%34 = getelementptr inbounds double, double* %33, i32 1
store double* %34, double** %6
%35 = load double, double* %33
%36 = call double @llvm.fmuladd.f64(double %28, double %32, double %35)
store double %36, double* %33
br label %37
37: ; preds = %27
%38 = load i32, i32* %11
%39 = add nsw i32 %38, -1
store i32 %39, i32* %11
br label %24
40: ; preds = %24
ret void
}

LLVM IR (-O2):
define void @scadd_mixed(double* nocapture %0, i32 %1, i32 %2, double %3, float* 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 %117, label %10
10: ; preds = %5
%11 = getelementptr inbounds float, float* %4, i64 %6
%12 = getelementptr inbounds 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 %101, label %17
17: ; preds = %10
%18 = and i64 %15, 8589934588
%19 = trunc i64 %18 to i32
%20 = sub i32 %8, %19
%21 = getelementptr float, float* %11, i64 %18
%22 = getelementptr double, double* %12, i64 %18
%23 = insertelement <2 x double> poison, double %3, i64 0
%24 = shufflevector <2 x double> %23, <2 x double> poison, <2 x i32> zeroinitializer
%25 = insertelement <2 x double> poison, double %3, i64 0
%26 = shufflevector <2 x double> %25, <2 x double> poison, <2 x i32> zeroinitializer
%27 = add nsw i64 %18, -4
%28 = lshr exact i64 %27, 2
%29 = add nuw nsw i64 %28, 1
%30 = and i64 %29, 1
%31 = icmp eq i64 %27, 0
br i1 %31, label %77, label %32
32: ; preds = %17
%33 = and i64 %29, 9223372036854775806
br label %34
34: ; preds = %34, %32
%35 = phi i64 [ 0, %32 ], [ %74, %34 ]
%36 = phi i64 [ 0, %32 ], [ %75, %34 ]
%37 = getelementptr float, float* %11, i64 %35
%38 = getelementptr double, double* %12, i64 %35
%39 = bitcast float* %37 to <2 x float>*
%40 = load <2 x float>, <2 x float>* %39
%41 = getelementptr float, float* %37, i64 2
%42 = bitcast float* %41 to <2 x float>*
%43 = load <2 x float>, <2 x float>* %42
%44 = fpext <2 x float> %40 to <2 x double>
%45 = fpext <2 x float> %43 to <2 x double>
%46 = bitcast double* %38 to <2 x double>*
%47 = load <2 x double>, <2 x double>* %46
%48 = getelementptr double, double* %38, i64 2
%49 = bitcast double* %48 to <2 x double>*
%50 = load <2 x double>, <2 x double>* %49
%51 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %24, <2 x double> %44, <2 x double> %47)
%52 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %26, <2 x double> %45, <2 x double> %50)
%53 = bitcast double* %38 to <2 x double>*
store <2 x double> %51, <2 x double>* %53
%54 = bitcast double* %48 to <2 x double>*
store <2 x double> %52, <2 x double>* %54
%55 = or i64 %35, 4
%56 = getelementptr float, float* %11, i64 %55
%57 = getelementptr double, double* %12, i64 %55
%58 = bitcast float* %56 to <2 x float>*
%59 = load <2 x float>, <2 x float>* %58
%60 = getelementptr float, float* %56, i64 2
%61 = bitcast float* %60 to <2 x float>*
%62 = load <2 x float>, <2 x float>* %61
%63 = fpext <2 x float> %59 to <2 x double>
%64 = fpext <2 x float> %62 to <2 x double>
%65 = bitcast double* %57 to <2 x double>*
%66 = load <2 x double>, <2 x double>* %65
%67 = getelementptr double, double* %57, i64 2
%68 = bitcast double* %67 to <2 x double>*
%69 = load <2 x double>, <2 x double>* %68
%70 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %24, <2 x double> %63, <2 x double> %66)
%71 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %26, <2 x double> %64, <2 x double> %69)
%72 = bitcast double* %57 to <2 x double>*
store <2 x double> %70, <2 x double>* %72
%73 = bitcast double* %67 to <2 x double>*
store <2 x double> %71, <2 x double>* %73
%74 = add nuw i64 %35, 8
%75 = add i64 %36, 2
%76 = icmp eq i64 %75, %33
br i1 %76, label %77, label %34
77: ; preds = %34, %17
%78 = phi i64 [ 0, %17 ], [ %74, %34 ]
%79 = icmp eq i64 %30, 0
br i1 %79, label %99, label %80
80: ; preds = %77
%81 = getelementptr float, float* %11, i64 %78
%82 = getelementptr double, double* %12, i64 %78
%83 = bitcast float* %81 to <2 x float>*
%84 = load <2 x float>, <2 x float>* %83
%85 = getelementptr float, float* %81, i64 2
%86 = bitcast float* %85 to <2 x float>*
%87 = load <2 x float>, <2 x float>* %86
%88 = fpext <2 x float> %84 to <2 x double>
%89 = fpext <2 x float> %87 to <2 x double>
%90 = bitcast double* %82 to <2 x double>*
%91 = load <2 x double>, <2 x double>* %90
%92 = getelementptr double, double* %82, i64 2
%93 = bitcast double* %92 to <2 x double>*
%94 = load <2 x double>, <2 x double>* %93
%95 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %24, <2 x double> %88, <2 x double> %91)
%96 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %26, <2 x double> %89, <2 x double> %94)
%97 = bitcast double* %82 to <2 x double>*
store <2 x double> %95, <2 x double>* %97
%98 = bitcast double* %92 to <2 x double>*
store <2 x double> %96, <2 x double>* %98
br label %99
99: ; preds = %77, %80
%100 = icmp eq i64 %15, %18
br i1 %100, label %117, label %101
101: ; preds = %10, %99
%102 = phi i32 [ %8, %10 ], [ %20, %99 ]
%103 = phi float* [ %11, %10 ], [ %21, %99 ]
%104 = phi double* [ %12, %10 ], [ %22, %99 ]
br label %105
105: ; preds = %101, %105
%106 = phi i32 [ %115, %105 ], [ %102, %101 ]
%107 = phi float* [ %109, %105 ], [ %103, %101 ]
%108 = phi double* [ %112, %105 ], [ %104, %101 ]
%109 = getelementptr inbounds float, float* %107, i64 1
%110 = load float, float* %107
%111 = fpext float %110 to double
%112 = getelementptr inbounds double, double* %108, i64 1
%113 = load double, double* %108
%114 = tail call double @llvm.fmuladd.f64(double %3, double %111, double %113)
store double %114, double* %108
%115 = add nsw i32 %106, -1
%116 = icmp eq i32 %115, 0
br i1 %116, label %117, label %105
117: ; preds = %105, %99, %5
ret void
}