short int goertzel_output_power_v0(short int *delay, short int coefficient) {
  long product1, product2, product3;
  signed int goertzel_power;
  product1 = ((long)delay[0] * delay[0]);
  product2 = ((long)delay[1] * delay[1]);
  product3 = (((long)delay[0] * coefficient) >> 14);
  product3 = (product3 * delay[1]);
  goertzel_power = (short int)((product1 + product2 - product3) >> 14 - 5);
  delay[0] = 0;
  delay[1] = 0;
  return goertzel_power;
}


x86 assembly (clang, -O0):
<goertzel_output_power_v0>:
push   %rbp
mov    %rsp,%rbp
mov    %si,%ax
mov    %rdi,-0x8(%rbp)
mov    %ax,-0xa(%rbp)
mov    -0x8(%rbp),%rax
movswq (%rax),%rax
mov    -0x8(%rbp),%rcx
movswq (%rcx),%rcx
imul   %rcx,%rax
mov    %rax,-0x18(%rbp)
mov    -0x8(%rbp),%rax
movswq 0x2(%rax),%rax
mov    -0x8(%rbp),%rcx
movswq 0x2(%rcx),%rcx
imul   %rcx,%rax
mov    %rax,-0x20(%rbp)
mov    -0x8(%rbp),%rax
movswq (%rax),%rax
movswq -0xa(%rbp),%rcx
imul   %rcx,%rax
sar    $0xe,%rax
mov    %rax,-0x28(%rbp)
mov    -0x28(%rbp),%rax
mov    -0x8(%rbp),%rcx
movswq 0x2(%rcx),%rcx
imul   %rcx,%rax
mov    %rax,-0x28(%rbp)
mov    -0x18(%rbp),%rax
add    -0x20(%rbp),%rax
sub    -0x28(%rbp),%rax
sar    $0x9,%rax
cwtl
mov    %eax,-0x2c(%rbp)
mov    -0x8(%rbp),%rax
movw   $0x0,(%rax)
mov    -0x8(%rbp),%rax
movw   $0x0,0x2(%rax)
mov    -0x2c(%rbp),%eax
cwtl
pop    %rbp
ret

x86 assembly (clang, -O2):
<goertzel_output_power_v0>:
movswl (%rdi),%eax
cltq
movslq %esi,%rcx
imul   %rax,%rcx
mov    %eax,%edx
imul   %eax,%edx
movswl 0x2(%rdi),%esi
shr    $0xe,%rcx
mov    %esi,%eax
sub    %ecx,%eax
imul   %esi,%eax
add    %edx,%eax
shr    $0x9,%eax
movl   $0x0,(%rdi)
ret

LLVM IR (-O0):
define signext i16 @goertzel_output_power_v0(i16* %0, i16 signext %1) {
%3 = alloca i16*
%4 = alloca i16
%5 = alloca i64
%6 = alloca i64
%7 = alloca i64
%8 = alloca i32
store i16* %0, i16** %3
store i16 %1, i16* %4
%9 = load i16*, i16** %3
%10 = getelementptr inbounds i16, i16* %9, i64 0
%11 = load i16, i16* %10
%12 = sext i16 %11 to i64
%13 = load i16*, i16** %3
%14 = getelementptr inbounds i16, i16* %13, i64 0
%15 = load i16, i16* %14
%16 = sext i16 %15 to i64
%17 = mul nsw i64 %12, %16
store i64 %17, i64* %5
%18 = load i16*, i16** %3
%19 = getelementptr inbounds i16, i16* %18, i64 1
%20 = load i16, i16* %19
%21 = sext i16 %20 to i64
%22 = load i16*, i16** %3
%23 = getelementptr inbounds i16, i16* %22, i64 1
%24 = load i16, i16* %23
%25 = sext i16 %24 to i64
%26 = mul nsw i64 %21, %25
store i64 %26, i64* %6
%27 = load i16*, i16** %3
%28 = getelementptr inbounds i16, i16* %27, i64 0
%29 = load i16, i16* %28
%30 = sext i16 %29 to i64
%31 = load i16, i16* %4
%32 = sext i16 %31 to i64
%33 = mul nsw i64 %30, %32
%34 = ashr i64 %33, 14
store i64 %34, i64* %7
%35 = load i64, i64* %7
%36 = load i16*, i16** %3
%37 = getelementptr inbounds i16, i16* %36, i64 1
%38 = load i16, i16* %37
%39 = sext i16 %38 to i64
%40 = mul nsw i64 %35, %39
store i64 %40, i64* %7
%41 = load i64, i64* %5
%42 = load i64, i64* %6
%43 = add nsw i64 %41, %42
%44 = load i64, i64* %7
%45 = sub nsw i64 %43, %44
%46 = ashr i64 %45, 9
%47 = trunc i64 %46 to i16
%48 = sext i16 %47 to i32
store i32 %48, i32* %8
%49 = load i16*, i16** %3
%50 = getelementptr inbounds i16, i16* %49, i64 0
store i16 0, i16* %50
%51 = load i16*, i16** %3
%52 = getelementptr inbounds i16, i16* %51, i64 1
store i16 0, i16* %52
%53 = load i32, i32* %8
%54 = trunc i32 %53 to i16
ret i16 %54
}

LLVM IR (-O2):
define signext i16 @goertzel_output_power_v0(i16* nocapture %0, i16 signext %1) {
%3 = load i16, i16* %0
%4 = sext i16 %3 to i64
%5 = mul nsw i64 %4, %4
%6 = getelementptr inbounds i16, i16* %0, i64 1
%7 = load i16, i16* %6
%8 = sext i16 %7 to i64
%9 = sext i16 %1 to i64
%10 = mul nsw i64 %4, %9
%11 = ashr i64 %10, 14
%12 = sub nsw i64 %8, %11
%13 = mul nsw i64 %12, %8
%14 = add nsw i64 %13, %5
%15 = lshr i64 %14, 9
%16 = trunc i64 %15 to i16
store i16 0, i16* %0
store i16 0, i16* %6
ret i16 %16
}