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


x86 assembly (clang, -O0):
<yscale_up_rgba>:
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    1305 <yscale_up_rgba+0x1d5>
movl   $0x0,-0x28(%rbp)
cmpl   $0x4,-0x28(%rbp)
jge    1215 <yscale_up_rgba+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,-0x40(%rbp,%rax,4)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1161 <yscale_up_rgba+0x31>
movss  -0x34(%rbp),%xmm0
call   1040 <clampf@plt>
movss  %xmm0,-0x2c(%rbp)
movl   $0x0,-0x28(%rbp)
cmpl   $0x3,-0x28(%rbp)
jge    12d3 <yscale_up_rgba+0x1a3>
movss  -0x2c(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomiss %xmm1,%xmm0
jne    1251 <yscale_up_rgba+0x121>
jp     1251 <yscale_up_rgba+0x121>
jmp    12a4 <yscale_up_rgba+0x174>
movss  -0x2c(%rbp),%xmm0
movss  0xda6(%rip),%xmm1
ucomiss %xmm1,%xmm0
jne    1272 <yscale_up_rgba+0x142>
jp     1272 <yscale_up_rgba+0x142>
jmp    12a4 <yscale_up_rgba+0x174>
movss  -0x2c(%rbp),%xmm1
movslq -0x28(%rbp),%rax
movss  -0x40(%rbp,%rax,4),%xmm0
divss  %xmm1,%xmm0
movss  %xmm0,-0x40(%rbp,%rax,4)
movslq -0x28(%rbp),%rax
movss  -0x40(%rbp,%rax,4),%xmm0
call   1040 <clampf@plt>
movslq -0x28(%rbp),%rax
movss  %xmm0,-0x40(%rbp,%rax,4)
movslq -0x28(%rbp),%rax
movss  -0x40(%rbp,%rax,4),%xmm0
call   1030 <linear_sample_to_srgb@plt>
mov    %al,%dl
mov    -0x20(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    -0x28(%rbp),%ecx
movslq %ecx,%rcx
mov    %dl,(%rax,%rcx,1)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    122b <yscale_up_rgba+0xfb>
movss  0xd25(%rip),%xmm0
mulss  -0x2c(%rbp),%xmm0
call   1050 <f2i@plt>
mov    %al,%dl
mov    -0x20(%rbp),%rax
mov    -0x24(%rbp),%ecx
add    $0x3,%ecx
movslq %ecx,%rcx
mov    %dl,(%rax,%rcx,1)
mov    -0x24(%rbp),%eax
add    $0x4,%eax
mov    %eax,-0x24(%rbp)
jmp    114e <yscale_up_rgba+0x1e>
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<yscale_up_rgba>:
push   %rbp
push   %r15
push   %r14
push   %r12
push   %rbx
sub    $0x10,%rsp
test   %esi,%esi
jle    131e <yscale_up_rgba+0x1ee>
mov    %rcx,%r14
mov    %rdx,%r12
mov    %esi,%r15d
mov    %rdi,%rbp
xor    %ebx,%ebx
jmp    11b3 <yscale_up_rgba+0x83>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movss  0x4(%rsp),%xmm0
call   1030 <linear_sample_to_srgb@plt>
mov    %al,(%r14,%rbx,1)
movss  0x8(%rsp),%xmm0
call   1030 <linear_sample_to_srgb@plt>
mov    %al,0x1(%r14,%rbx,1)
movss  0xc(%rsp),%xmm0
call   1030 <linear_sample_to_srgb@plt>
mov    %al,0x2(%r14,%rbx,1)
movss  (%rsp),%xmm0
mulss  0xe68(%rip),%xmm0
call   1050 <f2i@plt>
mov    %al,0x3(%r14,%rbx,1)
add    $0x4,%rbx
cmp    %r15d,%ebx
jge    131e <yscale_up_rgba+0x1ee>
movss  (%r12),%xmm1
movss  0x4(%r12),%xmm3
mov    0x0(%rbp),%rdx
mov    0x8(%rbp),%rsi
movss  0x8(%r12),%xmm2
mov    0x10(%rbp),%rcx
movss  0xc(%r12),%xmm0
mov    0x18(%rbp),%rax
movss  (%rsi,%rbx,4),%xmm4
mulss  %xmm3,%xmm4
movss  (%rdx,%rbx,4),%xmm5
mulss  %xmm1,%xmm5
addss  %xmm4,%xmm5
movss  (%rcx,%rbx,4),%xmm4
mulss  %xmm2,%xmm4
addss  %xmm5,%xmm4
movss  (%rax,%rbx,4),%xmm5
mulss  %xmm0,%xmm5
addss  %xmm4,%xmm5
movss  %xmm5,0x4(%rsp)
movss  0x4(%rsi,%rbx,4),%xmm4
mulss  %xmm3,%xmm4
movss  0x4(%rdx,%rbx,4),%xmm5
mulss  %xmm1,%xmm5
addss  %xmm4,%xmm5
movss  0x4(%rcx,%rbx,4),%xmm4
mulss  %xmm2,%xmm4
addss  %xmm5,%xmm4
movss  0x4(%rax,%rbx,4),%xmm5
mulss  %xmm0,%xmm5
addss  %xmm4,%xmm5
movss  %xmm5,0x8(%rsp)
movss  0x8(%rsi,%rbx,4),%xmm4
mulss  %xmm3,%xmm4
movss  0x8(%rdx,%rbx,4),%xmm5
mulss  %xmm1,%xmm5
addss  %xmm4,%xmm5
movss  0x8(%rcx,%rbx,4),%xmm4
mulss  %xmm2,%xmm4
addss  %xmm5,%xmm4
movss  0x8(%rax,%rbx,4),%xmm5
mulss  %xmm0,%xmm5
mulss  0xc(%rsi,%rbx,4),%xmm3
mulss  0xc(%rdx,%rbx,4),%xmm1
addss  %xmm4,%xmm5
movss  %xmm5,0xc(%rsp)
addss  %xmm3,%xmm1
mulss  0xc(%rcx,%rbx,4),%xmm2
addss  %xmm1,%xmm2
mulss  0xc(%rax,%rbx,4),%xmm0
addss  %xmm2,%xmm0
call   1040 <clampf@plt>
ucomiss 0xd50(%rip),%xmm0
movss  %xmm0,(%rsp)
jne    12c5 <yscale_up_rgba+0x195>
jnp    1160 <yscale_up_rgba+0x30>
movaps %xmm0,%xmm1
ucomiss 0xd31(%rip),%xmm0
jne    12d7 <yscale_up_rgba+0x1a7>
jnp    1160 <yscale_up_rgba+0x30>
movss  0x4(%rsp),%xmm0
divss  %xmm1,%xmm0
call   1040 <clampf@plt>
call   1030 <linear_sample_to_srgb@plt>
mov    %al,(%r14,%rbx,1)
movss  0x8(%rsp),%xmm0
divss  (%rsp),%xmm0
call   1040 <clampf@plt>
call   1030 <linear_sample_to_srgb@plt>
mov    %al,0x1(%r14,%rbx,1)
movss  0xc(%rsp),%xmm0
divss  (%rsp),%xmm0
call   1040 <clampf@plt>
jmp    1185 <yscale_up_rgba+0x55>
add    $0x10,%rsp
pop    %rbx
pop    %r12
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @yscale_up_rgba(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 [4 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 = %132, %4
%14 = load i32, i32* %9
%15 = load i32, i32* %6
%16 = icmp slt i32 %14, %15
br i1 %16, label %17, label %135
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, 4
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 [4 x float], [4 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 [4 x float], [4 x float]* %12, i64 0, i64 3
%82 = load float, float* %81
%83 = call float @clampf(float %82)
store float %83, float* %11
store i32 0, i32* %10
br label %84
84: ; preds = %120, %80
%85 = load i32, i32* %10
%86 = icmp slt i32 %85, 3
br i1 %86, label %87, label %123
87: ; preds = %84
%88 = load float, float* %11
%89 = fcmp une float %88, 0.000000e+00
br i1 %89, label %90, label %108
90: ; preds = %87
%91 = load float, float* %11
%92 = fcmp une float %91, 1.000000e+00
br i1 %92, label %93, label %108
93: ; preds = %90
%94 = load float, float* %11
%95 = load i32, i32* %10
%96 = sext i32 %95 to i64
%97 = getelementptr inbounds [4 x float], [4 x float]* %12, i64 0, i64 %96
%98 = load float, float* %97
%99 = fdiv float %98, %94
store float %99, float* %97
%100 = load i32, i32* %10
%101 = sext i32 %100 to i64
%102 = getelementptr inbounds [4 x float], [4 x float]* %12, i64 0, i64 %101
%103 = load float, float* %102
%104 = call float @clampf(float %103)
%105 = load i32, i32* %10
%106 = sext i32 %105 to i64
%107 = getelementptr inbounds [4 x float], [4 x float]* %12, i64 0, i64 %106
store float %104, float* %107
br label %108
108: ; preds = %93, %90, %87
%109 = load i32, i32* %10
%110 = sext i32 %109 to i64
%111 = getelementptr inbounds [4 x float], [4 x float]* %12, i64 0, i64 %110
%112 = load float, float* %111
%113 = call zeroext i8 @linear_sample_to_srgb(float %112)
%114 = load i8*, i8** %8
%115 = load i32, i32* %9
%116 = load i32, i32* %10
%117 = add nsw i32 %115, %116
%118 = sext i32 %117 to i64
%119 = getelementptr inbounds i8, i8* %114, i64 %118
store i8 %113, i8* %119
br label %120
120: ; preds = %108
%121 = load i32, i32* %10
%122 = add nsw i32 %121, 1
store i32 %122, i32* %10
br label %84
123: ; preds = %84
%124 = load float, float* %11
%125 = fmul float %124, 2.550000e+02
%126 = call zeroext i8 @f2i(float %125)
%127 = load i8*, i8** %8
%128 = load i32, i32* %9
%129 = add nsw i32 %128, 3
%130 = sext i32 %129 to i64
%131 = getelementptr inbounds i8, i8* %127, i64 %130
store i8 %126, i8* %131
br label %132
132: ; preds = %123
%133 = load i32, i32* %9
%134 = add nsw i32 %133, 4
store i32 %134, i32* %9
br label %13
135: ; preds = %13
ret void
}

LLVM IR (-O2):
define void @yscale_up_rgba(float** nocapture %0, i32 %1, float* nocapture %2, i8* nocapture writeonly %3) {
%5 = icmp sgt i32 %1, 0
br i1 %5, label %6, label %105
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, %94
%14 = phi i64 [ 0, %6 ], [ %102, %94 ]
%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 = or i64 %14, 2
%49 = getelementptr inbounds float, float* %16, i64 %48
%50 = load float, float* %49
%51 = getelementptr inbounds float, float* %18, i64 %48
%52 = load float, float* %51
%53 = fmul float %17, %52
%54 = tail call float @llvm.fmuladd.f32(float %15, float %50, float %53)
%55 = getelementptr inbounds float, float* %20, i64 %48
%56 = load float, float* %55
%57 = tail call float @llvm.fmuladd.f32(float %19, float %56, float %54)
%58 = getelementptr inbounds float, float* %22, i64 %48
%59 = load float, float* %58
%60 = tail call float @llvm.fmuladd.f32(float %21, float %59, float %57)
%61 = or i64 %14, 3
%62 = getelementptr inbounds float, float* %16, i64 %61
%63 = load float, float* %62
%64 = getelementptr inbounds float, float* %18, i64 %61
%65 = load float, float* %64
%66 = fmul float %17, %65
%67 = tail call float @llvm.fmuladd.f32(float %15, float %63, float %66)
%68 = getelementptr inbounds float, float* %20, i64 %61
%69 = load float, float* %68
%70 = tail call float @llvm.fmuladd.f32(float %19, float %69, float %67)
%71 = getelementptr inbounds float, float* %22, i64 %61
%72 = load float, float* %71
%73 = tail call float @llvm.fmuladd.f32(float %21, float %72, float %70)
%74 = tail call float @clampf(float %73)
%75 = fcmp une float %74, 0.000000e+00
%76 = fcmp une float %74, 1.000000e+00
%77 = and i1 %75, %76
br i1 %77, label %83, label %78
78: ; preds = %13
%79 = tail call zeroext i8 @linear_sample_to_srgb(float %34)
%80 = getelementptr inbounds i8, i8* %3, i64 %14
store i8 %79, i8* %80
%81 = tail call zeroext i8 @linear_sample_to_srgb(float %47)
%82 = getelementptr inbounds i8, i8* %3, i64 %35
store i8 %81, i8* %82
br label %94
83: ; preds = %13
%84 = fdiv float %34, %74
%85 = tail call float @clampf(float %84)
%86 = tail call zeroext i8 @linear_sample_to_srgb(float %85)
%87 = getelementptr inbounds i8, i8* %3, i64 %14
store i8 %86, i8* %87
%88 = fdiv float %47, %74
%89 = tail call float @clampf(float %88)
%90 = tail call zeroext i8 @linear_sample_to_srgb(float %89)
%91 = getelementptr inbounds i8, i8* %3, i64 %35
store i8 %90, i8* %91
%92 = fdiv float %60, %74
%93 = tail call float @clampf(float %92)
br label %94
94: ; preds = %78, %83
%95 = phi float [ %93, %83 ], [ %60, %78 ]
%96 = tail call zeroext i8 @linear_sample_to_srgb(float %95)
%97 = getelementptr inbounds i8, i8* %3, i64 %48
store i8 %96, i8* %97
%98 = fmul float %74, 2.550000e+02
%99 = tail call zeroext i8 @f2i(float %98)
%100 = or i64 %14, 3
%101 = getelementptr inbounds i8, i8* %3, i64 %100
store i8 %99, i8* %101
%102 = add nuw i64 %14, 4
%103 = trunc i64 %102 to i32
%104 = icmp slt i32 %103, %1
br i1 %104, label %13, label %105
105: ; preds = %94, %4
ret void
}