void complex_product(float *xy, float *x, float *y) {
  xy[0] = x[0] * y[0] - x[1] * y[1];
  xy[1] = x[0] * y[1] + x[1] * y[0];
}


x86 assembly (clang, -O0):
<complex_product>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    -0x10(%rbp),%rax
movss  (%rax),%xmm0
mov    -0x18(%rbp),%rax
movss  (%rax),%xmm2
mov    -0x10(%rbp),%rax
movss  0x4(%rax),%xmm1
mov    -0x18(%rbp),%rax
mulss  0x4(%rax),%xmm1
movd   %xmm1,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
movss  %xmm0,(%rax)
mov    -0x10(%rbp),%rax
movss  (%rax),%xmm0
mov    -0x18(%rbp),%rax
movss  0x4(%rax),%xmm2
mov    -0x10(%rbp),%rax
movss  0x4(%rax),%xmm1
mov    -0x18(%rbp),%rax
mulss  (%rax),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
movss  %xmm0,0x4(%rax)
pop    %rbp
ret

x86 assembly (clang, -O2):
<complex_product>:
movss  (%rsi),%xmm0
movss  0x4(%rdx),%xmm1
mulss  0x4(%rsi),%xmm1
mulss  (%rdx),%xmm0
subss  %xmm1,%xmm0
movss  %xmm0,(%rdi)
movss  (%rsi),%xmm0
movss  0x4(%rsi),%xmm1
mulss  (%rdx),%xmm1
mulss  0x4(%rdx),%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,0x4(%rdi)
ret

LLVM IR (-O0):
define void @complex_product(float* %0, float* %1, float* %2) {
%4 = alloca float*
%5 = alloca float*
%6 = alloca float*
store float* %0, float** %4
store float* %1, float** %5
store float* %2, float** %6
%7 = load float*, float** %5
%8 = getelementptr inbounds float, float* %7, i64 0
%9 = load float, float* %8
%10 = load float*, float** %6
%11 = getelementptr inbounds float, float* %10, i64 0
%12 = load float, float* %11
%13 = load float*, float** %5
%14 = getelementptr inbounds float, float* %13, i64 1
%15 = load float, float* %14
%16 = load float*, float** %6
%17 = getelementptr inbounds float, float* %16, i64 1
%18 = load float, float* %17
%19 = fmul float %15, %18
%20 = fneg float %19
%21 = call float @llvm.fmuladd.f32(float %9, float %12, float %20)
%22 = load float*, float** %4
%23 = getelementptr inbounds float, float* %22, i64 0
store float %21, float* %23
%24 = load float*, float** %5
%25 = getelementptr inbounds float, float* %24, i64 0
%26 = load float, float* %25
%27 = load float*, float** %6
%28 = getelementptr inbounds float, float* %27, i64 1
%29 = load float, float* %28
%30 = load float*, float** %5
%31 = getelementptr inbounds float, float* %30, i64 1
%32 = load float, float* %31
%33 = load float*, float** %6
%34 = getelementptr inbounds float, float* %33, i64 0
%35 = load float, float* %34
%36 = fmul float %32, %35
%37 = call float @llvm.fmuladd.f32(float %26, float %29, float %36)
%38 = load float*, float** %4
%39 = getelementptr inbounds float, float* %38, i64 1
store float %37, float* %39
ret void
}

LLVM IR (-O2):
define void @complex_product(float* nocapture writeonly %0, float* nocapture %1, float* nocapture %2) {
%4 = load float, float* %1
%5 = load float, float* %2
%6 = getelementptr inbounds float, float* %1, i64 1
%7 = load float, float* %6
%8 = getelementptr inbounds float, float* %2, i64 1
%9 = load float, float* %8
%10 = fneg float %7
%11 = fmul float %9, %10
%12 = tail call float @llvm.fmuladd.f32(float %4, float %5, float %11)
store float %12, float* %0
%13 = load float, float* %1
%14 = load float, float* %8
%15 = load float, float* %6
%16 = load float, float* %2
%17 = fmul float %15, %16
%18 = tail call float @llvm.fmuladd.f32(float %13, float %14, float %17)
%19 = getelementptr inbounds float, float* %0, i64 1
store float %18, float* %19
ret void
}