void hcconvolve(float *A, const float *B, const int N) {
  int R;
  int I;
  float T1;
  float T2;
  A[0] *= B[0];
  A[1] *= B[1];
  for (R = 2, I = 3; R < N; R += 2, I += 2) {
    T1 = A[R] * B[R];
    T2 = A[I] * B[I];
    A[I] = ((A[R] + A[I]) * (B[R] + B[I])) - (T1 + T2);
    A[R] = T1 - T2;
  }
}


x86 assembly (clang, -O0):
<hcconvolve>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    -0x10(%rbp),%rax
movss  (%rax),%xmm0
mov    -0x8(%rbp),%rax
mulss  (%rax),%xmm0
movss  %xmm0,(%rax)
mov    -0x10(%rbp),%rax
movss  0x4(%rax),%xmm0
mov    -0x8(%rbp),%rax
mulss  0x4(%rax),%xmm0
movss  %xmm0,0x4(%rax)
movl   $0x2,-0x18(%rbp)
movl   $0x3,-0x1c(%rbp)
mov    -0x18(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1220 <hcconvolve+0x120>
mov    -0x8(%rbp),%rax
movslq -0x18(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x18(%rbp),%rcx
mulss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x20(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
mulss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x24(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x18(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x1c(%rbp),%rcx
addss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x18(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
addss  (%rax,%rcx,4),%xmm2
movss  -0x20(%rbp),%xmm1
addss  -0x24(%rbp),%xmm1
movd   %xmm1,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
movss  -0x20(%rbp),%xmm0
subss  -0x24(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x18(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x18(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x18(%rbp)
mov    -0x1c(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x1c(%rbp)
jmp    1148 <hcconvolve+0x48>
pop    %rbp
ret

x86 assembly (clang, -O2):
<hcconvolve>:
movss  (%rsi),%xmm0
mulss  (%rdi),%xmm0
movss  %xmm0,(%rdi)
movss  0x4(%rsi),%xmm0
mulss  0x4(%rdi),%xmm0
movss  %xmm0,0x4(%rdi)
cmp    $0x3,%edx
jl     1271 <hcconvolve+0x171>
mov    %edx,%r10d
cmp    $0x5,%r10
mov    $0x4,%r11d
mov    $0x4,%r8d
cmovae %r10,%r8
add    $0xfffffffffffffffd,%r8
mov    $0x2,%edx
mov    $0x3,%ecx
cmp    $0x6,%r8
jb     1220 <hcconvolve+0x120>
lea    0x8(%rdi),%r9
cmp    $0x5,%r10
cmovae %r10,%r11
add    $0x1,%r11
and    $0xfffffffffffffffe,%r11
lea    (%rsi,%r11,4),%rax
cmp    %rax,%r9
jae    1181 <hcconvolve+0x81>
lea    (%rdi,%r11,4),%r9
lea    0x8(%rsi),%rax
cmp    %r9,%rax
jb     1220 <hcconvolve+0x120>
shr    %r8
add    $0x1,%r8
movabs $0x7ffffffffffffffc,%r9
ff ff 7f
and    %r8,%r9
lea    (%r9,%r9,1),%rcx
add    $0x3,%rcx
lea    (%r9,%r9,1),%rdx
add    $0x2,%rdx
xor    %eax,%eax
nopw   0x0(%rax,%rax,1)
movups 0x8(%rdi,%rax,8),%xmm0
movups 0x18(%rdi,%rax,8),%xmm1
movaps %xmm0,%xmm2
shufps $0x88,%xmm1,%xmm2
shufps $0xdd,%xmm1,%xmm0
movups 0x8(%rsi,%rax,8),%xmm1
movups 0x18(%rsi,%rax,8),%xmm3
movaps %xmm1,%xmm4
shufps $0x88,%xmm3,%xmm4
shufps $0xdd,%xmm3,%xmm1
movaps %xmm2,%xmm3
mulps  %xmm4,%xmm3
addps  %xmm0,%xmm2
mulps  %xmm1,%xmm0
addps  %xmm4,%xmm1
mulps  %xmm2,%xmm1
movaps %xmm3,%xmm2
addps  %xmm0,%xmm2
subps  %xmm2,%xmm1
subps  %xmm0,%xmm3
movaps %xmm3,%xmm0
unpckhps %xmm1,%xmm0
unpcklps %xmm1,%xmm3
movups %xmm3,0x8(%rdi,%rax,8)
movups %xmm0,0x18(%rdi,%rax,8)
add    $0x4,%rax
cmp    %rax,%r9
jne    11b0 <hcconvolve+0xb0>
cmp    %r9,%r8
je     1271 <hcconvolve+0x171>
nopl   0x0(%rax)
movss  (%rdi,%rdx,4),%xmm0
movss  (%rsi,%rdx,4),%xmm1
movaps %xmm0,%xmm2
mulss  %xmm1,%xmm2
movss  (%rdi,%rcx,4),%xmm3
movss  (%rsi,%rcx,4),%xmm4
addss  %xmm3,%xmm0
mulss  %xmm4,%xmm3
addss  %xmm4,%xmm1
mulss  %xmm0,%xmm1
movaps %xmm2,%xmm0
addss  %xmm3,%xmm0
subss  %xmm0,%xmm1
movss  %xmm1,(%rdi,%rcx,4)
subss  %xmm3,%xmm2
movss  %xmm2,(%rdi,%rdx,4)
add    $0x2,%rdx
add    $0x2,%rcx
cmp    %r10,%rdx
jb     1220 <hcconvolve+0x120>
ret

LLVM IR (-O0):
define void @hcconvolve(float* %0, float* %1, i32 %2) {
%4 = alloca float*
%5 = alloca float*
%6 = alloca i32
%7 = alloca i32
%8 = alloca i32
%9 = alloca float
%10 = alloca float
store float* %0, float** %4
store float* %1, float** %5
store i32 %2, i32* %6
%11 = load float*, float** %5
%12 = getelementptr inbounds float, float* %11, i64 0
%13 = load float, float* %12
%14 = load float*, float** %4
%15 = getelementptr inbounds float, float* %14, i64 0
%16 = load float, float* %15
%17 = fmul float %16, %13
store float %17, float* %15
%18 = load float*, float** %5
%19 = getelementptr inbounds float, float* %18, i64 1
%20 = load float, float* %19
%21 = load float*, float** %4
%22 = getelementptr inbounds float, float* %21, i64 1
%23 = load float, float* %22
%24 = fmul float %23, %20
store float %24, float* %22
store i32 2, i32* %7
store i32 3, i32* %8
br label %25
25: ; preds = %90, %3
%26 = load i32, i32* %7
%27 = load i32, i32* %6
%28 = icmp slt i32 %26, %27
br i1 %28, label %29, label %95
29: ; preds = %25
%30 = load float*, float** %4
%31 = load i32, i32* %7
%32 = sext i32 %31 to i64
%33 = getelementptr inbounds float, float* %30, i64 %32
%34 = load float, float* %33
%35 = load float*, float** %5
%36 = load i32, i32* %7
%37 = sext i32 %36 to i64
%38 = getelementptr inbounds float, float* %35, i64 %37
%39 = load float, float* %38
%40 = fmul float %34, %39
store float %40, float* %9
%41 = load float*, float** %4
%42 = load i32, i32* %8
%43 = sext i32 %42 to i64
%44 = getelementptr inbounds float, float* %41, i64 %43
%45 = load float, float* %44
%46 = load float*, float** %5
%47 = load i32, i32* %8
%48 = sext i32 %47 to i64
%49 = getelementptr inbounds float, float* %46, i64 %48
%50 = load float, float* %49
%51 = fmul float %45, %50
store float %51, float* %10
%52 = load float*, float** %4
%53 = load i32, i32* %7
%54 = sext i32 %53 to i64
%55 = getelementptr inbounds float, float* %52, i64 %54
%56 = load float, float* %55
%57 = load float*, float** %4
%58 = load i32, i32* %8
%59 = sext i32 %58 to i64
%60 = getelementptr inbounds float, float* %57, i64 %59
%61 = load float, float* %60
%62 = fadd float %56, %61
%63 = load float*, float** %5
%64 = load i32, i32* %7
%65 = sext i32 %64 to i64
%66 = getelementptr inbounds float, float* %63, i64 %65
%67 = load float, float* %66
%68 = load float*, float** %5
%69 = load i32, i32* %8
%70 = sext i32 %69 to i64
%71 = getelementptr inbounds float, float* %68, i64 %70
%72 = load float, float* %71
%73 = fadd float %67, %72
%74 = load float, float* %9
%75 = load float, float* %10
%76 = fadd float %74, %75
%77 = fneg float %76
%78 = call float @llvm.fmuladd.f32(float %62, float %73, float %77)
%79 = load float*, float** %4
%80 = load i32, i32* %8
%81 = sext i32 %80 to i64
%82 = getelementptr inbounds float, float* %79, i64 %81
store float %78, float* %82
%83 = load float, float* %9
%84 = load float, float* %10
%85 = fsub float %83, %84
%86 = load float*, float** %4
%87 = load i32, i32* %7
%88 = sext i32 %87 to i64
%89 = getelementptr inbounds float, float* %86, i64 %88
store float %85, float* %89
br label %90
90: ; preds = %29
%91 = load i32, i32* %7
%92 = add nsw i32 %91, 2
store i32 %92, i32* %7
%93 = load i32, i32* %8
%94 = add nsw i32 %93, 2
store i32 %94, i32* %8
br label %25
95: ; preds = %25
ret void
}

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