void vec_scaledxy_sub(double z[], const double x[], double x_scale,
                      const double y[], double y_scale, int n) {
  int i;
  for (i = 0; i < n; i++)
    z[i] = x[i] * x_scale - y[i] * y_scale;
}


x86 assembly (clang, -O0):
<vec_scaledxy_sub>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
movsd  %xmm0,-0x18(%rbp)
mov    %rdx,-0x20(%rbp)
movsd  %xmm1,-0x28(%rbp)
mov    %ecx,-0x2c(%rbp)
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x2c(%rbp),%eax
jge    118e <vec_scaledxy_sub+0x8e>
mov    -0x10(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x18(%rbp),%xmm2
mov    -0x20(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  -0x28(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1124 <vec_scaledxy_sub+0x24>
pop    %rbp
ret

x86 assembly (clang, -O2):
<vec_scaledxy_sub>:
test   %ecx,%ecx
jle    1232 <vec_scaledxy_sub+0x132>
mov    %ecx,%r8d
cmp    $0x4,%ecx
jae    1117 <vec_scaledxy_sub+0x17>
xor    %ecx,%ecx
jmp    11b7 <vec_scaledxy_sub+0xb7>
lea    (%rdi,%r8,8),%rcx
lea    (%rsi,%r8,8),%rax
cmp    %rdi,%rax
seta   %r10b
lea    (%rdx,%r8,8),%rax
cmp    %rsi,%rcx
seta   %r11b
cmp    %rdi,%rax
seta   %al
cmp    %rdx,%rcx
seta   %r9b
xor    %ecx,%ecx
test   %r11b,%r10b
jne    11b7 <vec_scaledxy_sub+0xb7>
and    %r9b,%al
jne    11b7 <vec_scaledxy_sub+0xb7>
mov    %r8d,%ecx
and    $0xfffffffc,%ecx
movapd %xmm1,%xmm2
unpcklpd %xmm1,%xmm2
movapd %xmm0,%xmm3
unpcklpd %xmm0,%xmm3
lea    0x0(,%r8,8),%r9
and    $0xffffffffffffffe0,%r9
xor    %eax,%eax
xchg   %ax,%ax
movupd (%rsi,%rax,1),%xmm4
movupd 0x10(%rsi,%rax,1),%xmm5
movupd (%rdx,%rax,1),%xmm6
movupd 0x10(%rdx,%rax,1),%xmm7
mulpd  %xmm2,%xmm6
mulpd  %xmm2,%xmm7
mulpd  %xmm3,%xmm4
subpd  %xmm6,%xmm4
mulpd  %xmm3,%xmm5
subpd  %xmm7,%xmm5
movupd %xmm4,(%rdi,%rax,1)
movupd %xmm5,0x10(%rdi,%rax,1)
add    $0x20,%rax
cmp    %rax,%r9
jne    1170 <vec_scaledxy_sub+0x70>
cmp    %r8,%rcx
je     1232 <vec_scaledxy_sub+0x132>
mov    %rcx,%rax
not    %rax
test   $0x1,%r8b
je     11e2 <vec_scaledxy_sub+0xe2>
movsd  (%rdx,%rcx,8),%xmm2
mulsd  %xmm1,%xmm2
movsd  (%rsi,%rcx,8),%xmm3
mulsd  %xmm0,%xmm3
subsd  %xmm2,%xmm3
movsd  %xmm3,(%rdi,%rcx,8)
or     $0x1,%rcx
add    %r8,%rax
je     1232 <vec_scaledxy_sub+0x132>
nopw   0x0(%rax,%rax,1)
movsd  (%rdx,%rcx,8),%xmm2
mulsd  %xmm1,%xmm2
movsd  (%rsi,%rcx,8),%xmm3
mulsd  %xmm0,%xmm3
subsd  %xmm2,%xmm3
movsd  %xmm3,(%rdi,%rcx,8)
movsd  0x8(%rdx,%rcx,8),%xmm2
mulsd  %xmm1,%xmm2
movsd  0x8(%rsi,%rcx,8),%xmm3
mulsd  %xmm0,%xmm3
subsd  %xmm2,%xmm3
movsd  %xmm3,0x8(%rdi,%rcx,8)
add    $0x2,%rcx
cmp    %rcx,%r8
jne    11f0 <vec_scaledxy_sub+0xf0>
ret

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

LLVM IR (-O2):
define void @vec_scaledxy_sub(double* nocapture writeonly %0, double* nocapture %1, double %2, double* nocapture %3, double %4, i32 %5) {
%7 = icmp sgt i32 %5, 0
br i1 %7, label %8, label %100
8: ; preds = %6
%9 = zext i32 %5 to i64
%10 = icmp ult i32 %5, 4
br i1 %10, label %60, label %11
11: ; preds = %8
%12 = getelementptr double, double* %0, i64 %9
%13 = getelementptr double, double* %1, i64 %9
%14 = getelementptr double, double* %3, i64 %9
%15 = icmp ugt double* %13, %0
%16 = icmp ugt double* %12, %1
%17 = and i1 %15, %16
%18 = icmp ugt double* %14, %0
%19 = icmp ugt double* %12, %3
%20 = and i1 %18, %19
%21 = or i1 %17, %20
br i1 %21, label %60, label %22
22: ; preds = %11
%23 = and i64 %9, 4294967292
%24 = insertelement <2 x double> poison, double %4, i64 0
%25 = shufflevector <2 x double> %24, <2 x double> poison, <2 x i32> zeroinitializer
%26 = insertelement <2 x double> poison, double %4, i64 0
%27 = shufflevector <2 x double> %26, <2 x double> poison, <2 x i32> zeroinitializer
%28 = insertelement <2 x double> poison, double %2, i64 0
%29 = shufflevector <2 x double> %28, <2 x double> poison, <2 x i32> zeroinitializer
%30 = insertelement <2 x double> poison, double %2, i64 0
%31 = shufflevector <2 x double> %30, <2 x double> poison, <2 x i32> zeroinitializer
br label %32
32: ; preds = %32, %22
%33 = phi i64 [ 0, %22 ], [ %56, %32 ]
%34 = getelementptr inbounds double, double* %1, i64 %33
%35 = bitcast double* %34 to <2 x double>*
%36 = load <2 x double>, <2 x double>* %35
%37 = getelementptr inbounds double, double* %34, i64 2
%38 = bitcast double* %37 to <2 x double>*
%39 = load <2 x double>, <2 x double>* %38
%40 = getelementptr inbounds double, double* %3, i64 %33
%41 = bitcast double* %40 to <2 x double>*
%42 = load <2 x double>, <2 x double>* %41
%43 = getelementptr inbounds double, double* %40, i64 2
%44 = bitcast double* %43 to <2 x double>*
%45 = load <2 x double>, <2 x double>* %44
%46 = fneg <2 x double> %42
%47 = fneg <2 x double> %45
%48 = fmul <2 x double> %25, %46
%49 = fmul <2 x double> %27, %47
%50 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %36, <2 x double> %29, <2 x double> %48)
%51 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %39, <2 x double> %31, <2 x double> %49)
%52 = getelementptr inbounds double, double* %0, i64 %33
%53 = bitcast double* %52 to <2 x double>*
store <2 x double> %50, <2 x double>* %53
%54 = getelementptr inbounds double, double* %52, i64 2
%55 = bitcast double* %54 to <2 x double>*
store <2 x double> %51, <2 x double>* %55
%56 = add nuw i64 %33, 4
%57 = icmp eq i64 %56, %23
br i1 %57, label %58, label %32
58: ; preds = %32
%59 = icmp eq i64 %23, %9
br i1 %59, label %100, label %60
60: ; preds = %11, %8, %58
%61 = phi i64 [ 0, %11 ], [ 0, %8 ], [ %23, %58 ]
%62 = xor i64 %61, -1
%63 = and i64 %9, 1
%64 = icmp eq i64 %63, 0
br i1 %64, label %75, label %65
65: ; preds = %60
%66 = getelementptr inbounds double, double* %1, i64 %61
%67 = load double, double* %66
%68 = getelementptr inbounds double, double* %3, i64 %61
%69 = load double, double* %68
%70 = fneg double %69
%71 = fmul double %70, %4
%72 = tail call double @llvm.fmuladd.f64(double %67, double %2, double %71)
%73 = getelementptr inbounds double, double* %0, i64 %61
store double %72, double* %73
%74 = or i64 %61, 1
br label %75
75: ; preds = %65, %60
%76 = phi i64 [ %61, %60 ], [ %74, %65 ]
%77 = sub nsw i64 0, %9
%78 = icmp eq i64 %62, %77
br i1 %78, label %100, label %79
79: ; preds = %75, %79
%80 = phi i64 [ %98, %79 ], [ %76, %75 ]
%81 = getelementptr inbounds double, double* %1, i64 %80
%82 = load double, double* %81
%83 = getelementptr inbounds double, double* %3, i64 %80
%84 = load double, double* %83
%85 = fneg double %84
%86 = fmul double %85, %4
%87 = tail call double @llvm.fmuladd.f64(double %82, double %2, double %86)
%88 = getelementptr inbounds double, double* %0, i64 %80
store double %87, double* %88
%89 = add nuw nsw i64 %80, 1
%90 = getelementptr inbounds double, double* %1, i64 %89
%91 = load double, double* %90
%92 = getelementptr inbounds double, double* %3, i64 %89
%93 = load double, double* %92
%94 = fneg double %93
%95 = fmul double %94, %4
%96 = tail call double @llvm.fmuladd.f64(double %91, double %2, double %95)
%97 = getelementptr inbounds double, double* %0, i64 %89
store double %96, double* %97
%98 = add nuw nsw i64 %80, 2
%99 = icmp eq i64 %98, %9
br i1 %99, label %100, label %79
100: ; preds = %75, %79, %58, %6
ret void
}