#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
double log10(double);
double sqrt(double);
extern double red_diff_uyvy(const unsigned char *video) {
  int sumSqDiff = 0;
  unsigned char red[] = {0x5f, 0x4b, 0xe6, 0x4b};
  sumSqDiff += (red[0] - video[0]) * (red[0] - video[0]);
  sumSqDiff += (red[1] - video[1]) * (red[1] - video[1]);
  sumSqDiff += (red[2] - video[2]) * (red[2] - video[2]);
  sumSqDiff += (red[3] - video[3]) * (red[3] - video[3]);
  if (sumSqDiff == 0) {
    return 20.0 * log10(255.0 / sqrt(1.0 / 4));
  }
  return 20.0 * log10(255.0 / sqrt((double)sumSqDiff / 4));
}


x86 assembly (clang, -O0):
<red_diff_uyvy>:
push   %rbp
mov    %rsp,%rbp
sub    $0x20,%rsp
mov    %rdi,-0x10(%rbp)
movl   $0x0,-0x14(%rbp)
mov    0xee7(%rip),%eax
mov    %eax,-0x18(%rbp)
movzbl -0x18(%rbp),%eax
mov    -0x10(%rbp),%rcx
movzbl (%rcx),%ecx
sub    %ecx,%eax
movzbl -0x18(%rbp),%ecx
mov    -0x10(%rbp),%rdx
movzbl (%rdx),%edx
sub    %edx,%ecx
imul   %ecx,%eax
add    -0x14(%rbp),%eax
mov    %eax,-0x14(%rbp)
movzbl -0x17(%rbp),%eax
mov    -0x10(%rbp),%rcx
movzbl 0x1(%rcx),%ecx
sub    %ecx,%eax
movzbl -0x17(%rbp),%ecx
mov    -0x10(%rbp),%rdx
movzbl 0x1(%rdx),%edx
sub    %edx,%ecx
imul   %ecx,%eax
add    -0x14(%rbp),%eax
mov    %eax,-0x14(%rbp)
movzbl -0x16(%rbp),%eax
mov    -0x10(%rbp),%rcx
movzbl 0x2(%rcx),%ecx
sub    %ecx,%eax
movzbl -0x16(%rbp),%ecx
mov    -0x10(%rbp),%rdx
movzbl 0x2(%rdx),%edx
sub    %edx,%ecx
imul   %ecx,%eax
add    -0x14(%rbp),%eax
mov    %eax,-0x14(%rbp)
movzbl -0x15(%rbp),%eax
mov    -0x10(%rbp),%rcx
movzbl 0x3(%rcx),%ecx
sub    %ecx,%eax
movzbl -0x15(%rbp),%ecx
mov    -0x10(%rbp),%rdx
movzbl 0x3(%rdx),%edx
sub    %edx,%ecx
imul   %ecx,%eax
add    -0x14(%rbp),%eax
mov    %eax,-0x14(%rbp)
cmpl   $0x0,-0x14(%rbp)
jne    1212 <red_diff_uyvy+0xf2>
movsd  0xe38(%rip),%xmm0
call   1040 <sqrt@plt>
movaps %xmm0,%xmm1
movsd  0xe18(%rip),%xmm0
divsd  %xmm1,%xmm0
call   1030 <log10@plt>
movaps %xmm0,%xmm1
movsd  0xdfc(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    1250 <red_diff_uyvy+0x130>
cvtsi2sdl -0x14(%rbp),%xmm0
movsd  0xdf1(%rip),%xmm1
divsd  %xmm1,%xmm0
call   1040 <sqrt@plt>
movaps %xmm0,%xmm1
movsd  0xdd5(%rip),%xmm0
divsd  %xmm1,%xmm0
call   1030 <log10@plt>
movaps %xmm0,%xmm1
movsd  0xdb9(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
add    $0x20,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<red_diff_uyvy>:
push   %rax
movd   (%rdi),%xmm0
pxor   %xmm1,%xmm1
punpcklbw %xmm1,%xmm0
punpcklwd %xmm1,%xmm0
movdqa 0xec7(%rip),%xmm1
psubd  %xmm0,%xmm1
packssdw %xmm1,%xmm1
movdqa %xmm1,%xmm0
pmulhw %xmm1,%xmm0
pmullw %xmm1,%xmm1
punpcklwd %xmm0,%xmm1
pshufd $0xee,%xmm1,%xmm0
punpcklqdq %xmm0,%xmm1
paddd  %xmm0,%xmm1
pshufd $0x55,%xmm1,%xmm0
paddd  %xmm1,%xmm0
movd   %xmm0,%eax
test   %eax,%eax
je     118e <red_diff_uyvy+0x6e>
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
mulsd  0xe9a(%rip),%xmm0
pxor   %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jb     1198 <red_diff_uyvy+0x78>
sqrtsd %xmm0,%xmm0
jmp    119d <red_diff_uyvy+0x7d>
movsd  0xe7a(%rip),%xmm0
pop    %rax
ret
call   1040 <sqrt@plt>
movsd  0xe7b(%rip),%xmm1
divsd  %xmm0,%xmm1
movapd %xmm1,%xmm0
call   1030 <log10@plt>
mulsd  0xe6e(%rip),%xmm0
pop    %rax
ret

LLVM IR (-O0):
@__const.red_diff_uyvy.red = [4 x i8] c"_K\E6K"
define double @red_diff_uyvy(i8* %0) {
%2 = alloca double
%3 = alloca i8*
%4 = alloca i32
%5 = alloca [4 x i8]
store i8* %0, i8** %3
store i32 0, i32* %4
%6 = bitcast [4 x i8]* %5 to i8*
%7 = getelementptr inbounds [4 x i8], [4 x i8]* %5, i64 0, i64 0
%8 = load i8, i8* %7
%9 = zext i8 %8 to i32
%10 = load i8*, i8** %3
%11 = getelementptr inbounds i8, i8* %10, i64 0
%12 = load i8, i8* %11
%13 = zext i8 %12 to i32
%14 = sub nsw i32 %9, %13
%15 = getelementptr inbounds [4 x i8], [4 x i8]* %5, i64 0, i64 0
%16 = load i8, i8* %15
%17 = zext i8 %16 to i32
%18 = load i8*, i8** %3
%19 = getelementptr inbounds i8, i8* %18, i64 0
%20 = load i8, i8* %19
%21 = zext i8 %20 to i32
%22 = sub nsw i32 %17, %21
%23 = mul nsw i32 %14, %22
%24 = load i32, i32* %4
%25 = add nsw i32 %24, %23
store i32 %25, i32* %4
%26 = getelementptr inbounds [4 x i8], [4 x i8]* %5, i64 0, i64 1
%27 = load i8, i8* %26
%28 = zext i8 %27 to i32
%29 = load i8*, i8** %3
%30 = getelementptr inbounds i8, i8* %29, i64 1
%31 = load i8, i8* %30
%32 = zext i8 %31 to i32
%33 = sub nsw i32 %28, %32
%34 = getelementptr inbounds [4 x i8], [4 x i8]* %5, i64 0, i64 1
%35 = load i8, i8* %34
%36 = zext i8 %35 to i32
%37 = load i8*, i8** %3
%38 = getelementptr inbounds i8, i8* %37, i64 1
%39 = load i8, i8* %38
%40 = zext i8 %39 to i32
%41 = sub nsw i32 %36, %40
%42 = mul nsw i32 %33, %41
%43 = load i32, i32* %4
%44 = add nsw i32 %43, %42
store i32 %44, i32* %4
%45 = getelementptr inbounds [4 x i8], [4 x i8]* %5, i64 0, i64 2
%46 = load i8, i8* %45
%47 = zext i8 %46 to i32
%48 = load i8*, i8** %3
%49 = getelementptr inbounds i8, i8* %48, i64 2
%50 = load i8, i8* %49
%51 = zext i8 %50 to i32
%52 = sub nsw i32 %47, %51
%53 = getelementptr inbounds [4 x i8], [4 x i8]* %5, i64 0, i64 2
%54 = load i8, i8* %53
%55 = zext i8 %54 to i32
%56 = load i8*, i8** %3
%57 = getelementptr inbounds i8, i8* %56, i64 2
%58 = load i8, i8* %57
%59 = zext i8 %58 to i32
%60 = sub nsw i32 %55, %59
%61 = mul nsw i32 %52, %60
%62 = load i32, i32* %4
%63 = add nsw i32 %62, %61
store i32 %63, i32* %4
%64 = getelementptr inbounds [4 x i8], [4 x i8]* %5, i64 0, i64 3
%65 = load i8, i8* %64
%66 = zext i8 %65 to i32
%67 = load i8*, i8** %3
%68 = getelementptr inbounds i8, i8* %67, i64 3
%69 = load i8, i8* %68
%70 = zext i8 %69 to i32
%71 = sub nsw i32 %66, %70
%72 = getelementptr inbounds [4 x i8], [4 x i8]* %5, i64 0, i64 3
%73 = load i8, i8* %72
%74 = zext i8 %73 to i32
%75 = load i8*, i8** %3
%76 = getelementptr inbounds i8, i8* %75, i64 3
%77 = load i8, i8* %76
%78 = zext i8 %77 to i32
%79 = sub nsw i32 %74, %78
%80 = mul nsw i32 %71, %79
%81 = load i32, i32* %4
%82 = add nsw i32 %81, %80
store i32 %82, i32* %4
%83 = load i32, i32* %4
%84 = icmp eq i32 %83, 0
br i1 %84, label %85, label %90
85: ; preds = %1
%86 = call double @sqrt(double 2.500000e-01)
%87 = fdiv double 2.550000e+02, %86
%88 = call double @log10(double %87)
%89 = fmul double 2.000000e+01, %88
store double %89, double* %2
br label %98
90: ; preds = %1
%91 = load i32, i32* %4
%92 = sitofp i32 %91 to double
%93 = fdiv double %92, 4.000000e+00
%94 = call double @sqrt(double %93)
%95 = fdiv double 2.550000e+02, %94
%96 = call double @log10(double %95)
%97 = fmul double 2.000000e+01, %96
store double %97, double* %2
br label %98
98: ; preds = %90, %85
%99 = load double, double* %2
ret double %99
}

LLVM IR (-O2):
define double @red_diff_uyvy(i8* nocapture %0) {
%2 = bitcast i8* %0 to <4 x i8>*
%3 = load <4 x i8>, <4 x i8>* %2
%4 = zext <4 x i8> %3 to <4 x i32>
%5 = sub nsw <4 x i32> <i32 95, i32 75, i32 230, i32 75>, %4
%6 = mul nsw <4 x i32> %5, %5
%7 = call i32 @llvm.vector.reduce.add.v4i32(<4 x i32> %6)
%8 = icmp eq i32 %7, 0
br i1 %8, label %16, label %9
9: ; preds = %1
%10 = sitofp i32 %7 to double
%11 = fmul double %10, 2.500000e-01
%12 = tail call double @sqrt(double %11)
%13 = fdiv double 2.550000e+02, %12
%14 = tail call double @log10(double %13)
%15 = fmul double %14, 2.000000e+01
br label %16
16: ; preds = %1, %9
%17 = phi double [ %15, %9 ], [ 0x404B136130CBA7EE, %1 ]
ret double %17
}