#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 */
unsigned char clamp8(float);
float clampf(float);
unsigned char f2i(float);
void yscale_up_ga(float **in, int len, float *coeffs, unsigned char *out) {
  int i, j;
  float alpha, sums[2];
  for (i = 0; i < len; i += 2) {
    for (j = 0; j < 2; j++) {
      sums[j] = coeffs[0] * in[0][i + j] + coeffs[1] * in[1][i + j] +
                coeffs[2] * in[2][i + j] + coeffs[3] * in[3][i + j];
    }
    alpha = clampf(sums[1]);
    if (alpha != 0) {
      sums[0] /= alpha;
    }
    out[i] = clamp8(sums[0]);
    out[i + 1] = f2i(alpha * 255.0f);
  }
}


x86 assembly (clang, -O0):
<yscale_up_ga>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    129c <yscale_up_ga+0x16c>
movl   $0x0,-0x28(%rbp)
cmpl   $0x2,-0x28(%rbp)
jge    1215 <yscale_up_ga+0xe5>
mov    -0x18(%rbp),%rax
movss  (%rax),%xmm0
mov    -0x8(%rbp),%rax
mov    (%rax),%rax
mov    -0x24(%rbp),%ecx
add    -0x28(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm2
mov    -0x18(%rbp),%rax
movss  0x4(%rax),%xmm1
mov    -0x8(%rbp),%rax
mov    0x8(%rax),%rax
mov    -0x24(%rbp),%ecx
add    -0x28(%rbp),%ecx
movslq %ecx,%rcx
mulss  (%rax,%rcx,4),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
movss  0x8(%rax),%xmm1
mov    -0x8(%rbp),%rax
mov    0x10(%rax),%rax
mov    -0x24(%rbp),%ecx
add    -0x28(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm2
mulss  %xmm2,%xmm1
addss  %xmm0,%xmm1
mov    -0x18(%rbp),%rax
movss  0xc(%rax),%xmm0
mov    -0x8(%rbp),%rax
mov    0x18(%rax),%rax
mov    -0x24(%rbp),%ecx
add    -0x28(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movslq -0x28(%rbp),%rax
movss  %xmm0,-0x34(%rbp,%rax,4)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1161 <yscale_up_ga+0x31>
movss  -0x30(%rbp),%xmm0
call   1040 <clampf@plt>
movss  %xmm0,-0x2c(%rbp)
movss  -0x2c(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomiss %xmm1,%xmm0
jne    1240 <yscale_up_ga+0x110>
jp     1240 <yscale_up_ga+0x110>
jmp    1253 <yscale_up_ga+0x123>
movss  -0x2c(%rbp),%xmm1
movss  -0x34(%rbp),%xmm0
divss  %xmm1,%xmm0
movss  %xmm0,-0x34(%rbp)
movss  -0x34(%rbp),%xmm0
call   1030 <clamp8@plt>
mov    %al,%dl
mov    -0x20(%rbp),%rax
movslq -0x24(%rbp),%rcx
mov    %dl,(%rax,%rcx,1)
movss  0xd8e(%rip),%xmm0
mulss  -0x2c(%rbp),%xmm0
call   1050 <f2i@plt>
mov    %al,%dl
mov    -0x20(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
mov    %dl,(%rax,%rcx,1)
mov    -0x24(%rbp),%eax
add    $0x2,%eax
mov    %eax,-0x24(%rbp)
jmp    114e <yscale_up_ga+0x1e>
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<yscale_up_ga>:
push   %rbp
push   %r15
push   %r14
push   %r12
push   %rbx
sub    $0x20,%rsp
test   %esi,%esi
jle    1246 <yscale_up_ga+0x116>
mov    %rcx,%r14
mov    %rdx,%r12
mov    %esi,%r15d
mov    %rdi,%rbp
xor    %ebx,%ebx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movss  (%r12),%xmm1
movss  0x4(%r12),%xmm3
mov    0x0(%rbp),%rax
mov    0x8(%rbp),%rcx
movss  0x8(%r12),%xmm2
mov    0x10(%rbp),%rdx
movss  0xc(%r12),%xmm0
mov    0x18(%rbp),%rsi
movss  (%rcx,%rbx,4),%xmm4
mulss  %xmm3,%xmm4
movss  (%rax,%rbx,4),%xmm5
mulss  %xmm1,%xmm5
addss  %xmm4,%xmm5
movss  (%rdx,%rbx,4),%xmm4
mulss  %xmm2,%xmm4
addss  %xmm5,%xmm4
movss  (%rsi,%rbx,4),%xmm5
mulss  %xmm0,%xmm5
mulss  0x4(%rcx,%rbx,4),%xmm3
mulss  0x4(%rax,%rbx,4),%xmm1
addss  %xmm4,%xmm5
movss  %xmm5,0xc(%rsp)
addss  %xmm3,%xmm1
mulss  0x4(%rdx,%rbx,4),%xmm2
addss  %xmm1,%xmm2
mulss  0x4(%rsi,%rbx,4),%xmm0
addss  %xmm2,%xmm0
call   1040 <clampf@plt>
movaps %xmm0,0x10(%rsp)
movaps %xmm0,%xmm3
cmpneqss 0xe0d(%rip),%xmm3
movaps %xmm3,%xmm1
andps  %xmm0,%xmm1
movss  0xdf7(%rip),%xmm2
andnps %xmm2,%xmm3
orps   %xmm1,%xmm3
movss  0xc(%rsp),%xmm0
divss  %xmm3,%xmm0
call   1030 <clamp8@plt>
mov    %al,(%r14,%rbx,1)
movaps 0x10(%rsp),%xmm0
mulss  0xdd5(%rip),%xmm0
call   1050 <f2i@plt>
mov    %al,0x1(%r14,%rbx,1)
add    $0x2,%rbx
cmp    %r15d,%ebx
jl     1160 <yscale_up_ga+0x30>
add    $0x20,%rsp
pop    %rbx
pop    %r12
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @yscale_up_ga(float** %0, i32 %1, float* %2, i8* %3) {
%5 = alloca float**
%6 = alloca i32
%7 = alloca float*
%8 = alloca i8*
%9 = alloca i32
%10 = alloca i32
%11 = alloca float
%12 = alloca [2 x float]
store float** %0, float*** %5
store i32 %1, i32* %6
store float* %2, float** %7
store i8* %3, i8** %8
store i32 0, i32* %9
br label %13
13: ; preds = %107, %4
%14 = load i32, i32* %9
%15 = load i32, i32* %6
%16 = icmp slt i32 %14, %15
br i1 %16, label %17, label %110
17: ; preds = %13
store i32 0, i32* %10
br label %18
18: ; preds = %77, %17
%19 = load i32, i32* %10
%20 = icmp slt i32 %19, 2
br i1 %20, label %21, label %80
21: ; preds = %18
%22 = load float*, float** %7
%23 = getelementptr inbounds float, float* %22, i64 0
%24 = load float, float* %23
%25 = load float**, float*** %5
%26 = getelementptr inbounds float*, float** %25, i64 0
%27 = load float*, float** %26
%28 = load i32, i32* %9
%29 = load i32, i32* %10
%30 = add nsw i32 %28, %29
%31 = sext i32 %30 to i64
%32 = getelementptr inbounds float, float* %27, i64 %31
%33 = load float, float* %32
%34 = load float*, float** %7
%35 = getelementptr inbounds float, float* %34, i64 1
%36 = load float, float* %35
%37 = load float**, float*** %5
%38 = getelementptr inbounds float*, float** %37, i64 1
%39 = load float*, float** %38
%40 = load i32, i32* %9
%41 = load i32, i32* %10
%42 = add nsw i32 %40, %41
%43 = sext i32 %42 to i64
%44 = getelementptr inbounds float, float* %39, i64 %43
%45 = load float, float* %44
%46 = fmul float %36, %45
%47 = call float @llvm.fmuladd.f32(float %24, float %33, float %46)
%48 = load float*, float** %7
%49 = getelementptr inbounds float, float* %48, i64 2
%50 = load float, float* %49
%51 = load float**, float*** %5
%52 = getelementptr inbounds float*, float** %51, i64 2
%53 = load float*, float** %52
%54 = load i32, i32* %9
%55 = load i32, i32* %10
%56 = add nsw i32 %54, %55
%57 = sext i32 %56 to i64
%58 = getelementptr inbounds float, float* %53, i64 %57
%59 = load float, float* %58
%60 = call float @llvm.fmuladd.f32(float %50, float %59, float %47)
%61 = load float*, float** %7
%62 = getelementptr inbounds float, float* %61, i64 3
%63 = load float, float* %62
%64 = load float**, float*** %5
%65 = getelementptr inbounds float*, float** %64, i64 3
%66 = load float*, float** %65
%67 = load i32, i32* %9
%68 = load i32, i32* %10
%69 = add nsw i32 %67, %68
%70 = sext i32 %69 to i64
%71 = getelementptr inbounds float, float* %66, i64 %70
%72 = load float, float* %71
%73 = call float @llvm.fmuladd.f32(float %63, float %72, float %60)
%74 = load i32, i32* %10
%75 = sext i32 %74 to i64
%76 = getelementptr inbounds [2 x float], [2 x float]* %12, i64 0, i64 %75
store float %73, float* %76
br label %77
77: ; preds = %21
%78 = load i32, i32* %10
%79 = add nsw i32 %78, 1
store i32 %79, i32* %10
br label %18
80: ; preds = %18
%81 = getelementptr inbounds [2 x float], [2 x float]* %12, i64 0, i64 1
%82 = load float, float* %81
%83 = call float @clampf(float %82)
store float %83, float* %11
%84 = load float, float* %11
%85 = fcmp une float %84, 0.000000e+00
br i1 %85, label %86, label %91
86: ; preds = %80
%87 = load float, float* %11
%88 = getelementptr inbounds [2 x float], [2 x float]* %12, i64 0, i64 0
%89 = load float, float* %88
%90 = fdiv float %89, %87
store float %90, float* %88
br label %91
91: ; preds = %86, %80
%92 = getelementptr inbounds [2 x float], [2 x float]* %12, i64 0, i64 0
%93 = load float, float* %92
%94 = call zeroext i8 @clamp8(float %93)
%95 = load i8*, i8** %8
%96 = load i32, i32* %9
%97 = sext i32 %96 to i64
%98 = getelementptr inbounds i8, i8* %95, i64 %97
store i8 %94, i8* %98
%99 = load float, float* %11
%100 = fmul float %99, 2.550000e+02
%101 = call zeroext i8 @f2i(float %100)
%102 = load i8*, i8** %8
%103 = load i32, i32* %9
%104 = add nsw i32 %103, 1
%105 = sext i32 %104 to i64
%106 = getelementptr inbounds i8, i8* %102, i64 %105
store i8 %101, i8* %106
br label %107
107: ; preds = %91
%108 = load i32, i32* %9
%109 = add nsw i32 %108, 2
store i32 %109, i32* %9
br label %13
110: ; preds = %13
ret void
}

LLVM IR (-O2):
define void @yscale_up_ga(float** nocapture %0, i32 %1, float* nocapture %2, i8* nocapture writeonly %3) {
%5 = icmp sgt i32 %1, 0
br i1 %5, label %6, label %61
6: ; preds = %4
%7 = getelementptr inbounds float, float* %2, i64 1
%8 = getelementptr inbounds float*, float** %0, i64 1
%9 = getelementptr inbounds float, float* %2, i64 2
%10 = getelementptr inbounds float*, float** %0, i64 2
%11 = getelementptr inbounds float, float* %2, i64 3
%12 = getelementptr inbounds float*, float** %0, i64 3
br label %13
13: ; preds = %6, %13
%14 = phi i64 [ 0, %6 ], [ %58, %13 ]
%15 = load float, float* %2
%16 = load float*, float** %0
%17 = load float, float* %7
%18 = load float*, float** %8
%19 = load float, float* %9
%20 = load float*, float** %10
%21 = load float, float* %11
%22 = load float*, float** %12
%23 = getelementptr inbounds float, float* %16, i64 %14
%24 = load float, float* %23
%25 = getelementptr inbounds float, float* %18, i64 %14
%26 = load float, float* %25
%27 = fmul float %17, %26
%28 = tail call float @llvm.fmuladd.f32(float %15, float %24, float %27)
%29 = getelementptr inbounds float, float* %20, i64 %14
%30 = load float, float* %29
%31 = tail call float @llvm.fmuladd.f32(float %19, float %30, float %28)
%32 = getelementptr inbounds float, float* %22, i64 %14
%33 = load float, float* %32
%34 = tail call float @llvm.fmuladd.f32(float %21, float %33, float %31)
%35 = or i64 %14, 1
%36 = getelementptr inbounds float, float* %16, i64 %35
%37 = load float, float* %36
%38 = getelementptr inbounds float, float* %18, i64 %35
%39 = load float, float* %38
%40 = fmul float %17, %39
%41 = tail call float @llvm.fmuladd.f32(float %15, float %37, float %40)
%42 = getelementptr inbounds float, float* %20, i64 %35
%43 = load float, float* %42
%44 = tail call float @llvm.fmuladd.f32(float %19, float %43, float %41)
%45 = getelementptr inbounds float, float* %22, i64 %35
%46 = load float, float* %45
%47 = tail call float @llvm.fmuladd.f32(float %21, float %46, float %44)
%48 = tail call float @clampf(float %47)
%49 = fcmp une float %48, 0.000000e+00
%50 = select i1 %49, float %48, float 1.000000e+00
%51 = fdiv float %34, %50
%52 = tail call zeroext i8 @clamp8(float %51)
%53 = getelementptr inbounds i8, i8* %3, i64 %14
store i8 %52, i8* %53
%54 = fmul float %48, 2.550000e+02
%55 = tail call zeroext i8 @f2i(float %54)
%56 = or i64 %14, 1
%57 = getelementptr inbounds i8, i8* %3, i64 %56
store i8 %55, i8* %57
%58 = add nuw i64 %14, 2
%59 = trunc i64 %58 to i32
%60 = icmp slt i32 %59, %1
br i1 %60, label %13, label %61
61: ; preds = %13, %4
ret void
}