#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 */
typedef float uint32_t;

/* Variables and functions */
void build_channel_basecurve(int width_jpeg, int height_jpeg, float *buf_jpeg,
                             int offx_raw, int offy_raw, int width_raw,
                             float *buf_raw, int ch, float *curve,
                             uint32_t *cnt) {
  for (int j = 0; j < height_jpeg; j++) {
    for (int i = 0; i < width_jpeg; i++) {
      const int ri = offx_raw + i;
      const int rj = offy_raw + j;
      float jpegVal = buf_jpeg[3 * (width_jpeg * j + i) + ch];
      float rawVal = buf_raw[3 * (width_raw * rj + ri) + ch];
      size_t raw = (size_t)((rawVal * (float)(0x10000 - 1)) + 0.5f);
      curve[raw] = (curve[raw] * cnt[raw] + jpegVal) / (cnt[raw] + 1.0f);
      cnt[raw]++;
    }
  }
}


x86 assembly (clang, -O0):
<build_channel_basecurve>:
push   %rbp
mov    %rsp,%rbp
mov    0x28(%rbp),%rax
mov    0x20(%rbp),%rax
mov    0x18(%rbp),%eax
mov    0x10(%rbp),%rax
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %rdx,-0x10(%rbp)
mov    %ecx,-0x14(%rbp)
mov    %r8d,-0x18(%rbp)
mov    %r9d,-0x1c(%rbp)
movl   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    1285 <build_channel_basecurve+0x185>
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    1272 <build_channel_basecurve+0x172>
mov    -0x14(%rbp),%eax
mov    -0x24(%rbp),%ecx
add    %ecx,%eax
mov    %eax,-0x28(%rbp)
mov    -0x18(%rbp),%eax
mov    -0x20(%rbp),%ecx
add    %ecx,%eax
mov    %eax,-0x2c(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x4(%rbp),%edx
mov    -0x20(%rbp),%ecx
imul   %ecx,%edx
mov    -0x24(%rbp),%ecx
add    %ecx,%edx
mov    %edx,%ecx
lea    (%rcx,%rcx,2),%ecx
mov    0x18(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x30(%rbp)
mov    0x10(%rbp),%rax
mov    -0x1c(%rbp),%edx
mov    -0x2c(%rbp),%ecx
imul   %ecx,%edx
mov    -0x28(%rbp),%ecx
add    %ecx,%edx
mov    %edx,%ecx
lea    (%rcx,%rcx,2),%ecx
mov    0x18(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x34(%rbp)
movss  -0x34(%rbp),%xmm0
movss  0xe41(%rip),%xmm1
mulss  %xmm1,%xmm0
movss  0xe39(%rip),%xmm1
addss  %xmm1,%xmm0
movaps %xmm0,%xmm1
cvttss2si %xmm1,%rax
mov    %rax,%rdx
sar    $0x3f,%rdx
movss  0xe22(%rip),%xmm1
subss  %xmm1,%xmm0
cvttss2si %xmm0,%rcx
and    %rdx,%rcx
or     %rcx,%rax
mov    %rax,-0x40(%rbp)
mov    0x20(%rbp),%rax
mov    -0x40(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    0x28(%rbp),%rax
mov    -0x40(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
movss  -0x30(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    0x28(%rbp),%rax
mov    -0x40(%rbp),%rcx
movss  0xdcc(%rip),%xmm1
addss  (%rax,%rcx,4),%xmm1
divss  %xmm1,%xmm0
mov    0x20(%rbp),%rax
mov    -0x40(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    0x28(%rbp),%rax
mov    -0x40(%rbp),%rcx
movss  0xda6(%rip),%xmm0
addss  (%rax,%rcx,4),%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1142 <build_channel_basecurve+0x42>
jmp    1277 <build_channel_basecurve+0x177>
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    112f <build_channel_basecurve+0x2f>
pop    %rbp
ret

x86 assembly (clang, -O2):
<build_channel_basecurve>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
test   %esi,%esi
jle    120f <build_channel_basecurve+0x10f>
mov    0x50(%rsp),%rax
mov    0x48(%rsp),%r10
mov    0x40(%rsp),%ebx
mov    0x38(%rsp),%r11
movslq %ebx,%r12
movslq %edi,%rbx
mov    %esi,%r14d
mov    %ebx,%esi
lea    (%rdx,%r12,4),%r13
shl    $0x2,%rbx
lea    (%rbx,%rbx,2),%r15
imul   %r9d,%r8d
lea    (%r8,%r8,2),%ebx
add    %ebx,%r12d
lea    (%rcx,%rcx,2),%ecx
add    %r12d,%ecx
lea    (%r9,%r9,2),%r8d
lea    (%rsi,%rsi,2),%r12
xor    %r9d,%r9d
movss  0xea0(%rip),%xmm0
movss  0xe9c(%rip),%xmm1
movss  0xe98(%rip),%xmm2
movss  0xe94(%rip),%xmm3
jmp    1193 <build_channel_basecurve+0x93>
nopw   0x0(%rax,%rax,1)
add    $0x1,%r9
add    %r15,%r13
add    %r8,%rcx
cmp    %r14,%r9
je     120f <build_channel_basecurve+0x10f>
test   %edi,%edi
jle    1180 <build_channel_basecurve+0x80>
xor    %esi,%esi
nopl   0x0(%rax)
lea    (%rcx,%rsi,1),%ebp
movslq %ebp,%rbp
movss  (%r11,%rbp,4),%xmm4
mulss  %xmm0,%xmm4
addss  %xmm1,%xmm4
cvttss2si %xmm4,%rbx
mov    %rbx,%rbp
subss  %xmm2,%xmm4
cvttss2si %xmm4,%rdx
sar    $0x3f,%rbp
and    %rbp,%rdx
or     %rbx,%rdx
movss  (%rax,%rdx,4),%xmm4
movss  (%r10,%rdx,4),%xmm5
mulss  %xmm4,%xmm5
addss  0x0(%r13,%rsi,4),%xmm5
addss  %xmm3,%xmm4
divss  %xmm4,%xmm5
movss  %xmm5,(%r10,%rdx,4)
movss  (%rax,%rdx,4),%xmm4
addss  %xmm3,%xmm4
movss  %xmm4,(%rax,%rdx,4)
add    $0x3,%rsi
cmp    %rsi,%r12
jne    11a0 <build_channel_basecurve+0xa0>
jmp    1180 <build_channel_basecurve+0x80>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define void @build_channel_basecurve(i32 %0, i32 %1, float* nocapture %2, i32 %3, i32 %4, i32 %5, float* nocapture %6, i32 %7, float* nocapture %8, float* nocapture %9) {
%11 = icmp sgt i32 %1, 0
br i1 %11, label %12, label %26
12: ; preds = %10
%13 = icmp sgt i32 %0, 0
%14 = sext i32 %7 to i64
%15 = sext i32 %0 to i64
%16 = zext i32 %1 to i64
%17 = zext i32 %0 to i64
br label %18
18: ; preds = %12, %27
%19 = phi i64 [ 0, %12 ], [ %28, %27 ]
br i1 %13, label %20, label %27
20: ; preds = %18
%21 = mul nsw i64 %19, %15
%22 = trunc i64 %19 to i32
%23 = add i32 %22, %4
%24 = mul i32 %23, %5
%25 = add i32 %24, %3
br label %30
26: ; preds = %27, %10
ret void
27: ; preds = %30, %18
%28 = add nuw nsw i64 %19, 1
%29 = icmp eq i64 %28, %16
br i1 %29, label %26, label %18
30: ; preds = %20, %30
%31 = phi i64 [ 0, %20 ], [ %55, %30 ]
%32 = add nsw i64 %31, %21
%33 = mul nsw i64 %32, 3
%34 = add nsw i64 %33, %14
%35 = getelementptr inbounds float, float* %2, i64 %34
%36 = load float, float* %35
%37 = trunc i64 %31 to i32
%38 = add i32 %25, %37
%39 = mul nsw i32 %38, 3
%40 = add nsw i32 %39, %7
%41 = sext i32 %40 to i64
%42 = getelementptr inbounds float, float* %6, i64 %41
%43 = load float, float* %42
%44 = tail call float @llvm.fmuladd.f32(float %43, float 6.553500e+04, float 5.000000e-01)
%45 = fptoui float %44 to i64
%46 = getelementptr inbounds float, float* %8, i64 %45
%47 = load float, float* %46
%48 = getelementptr inbounds float, float* %9, i64 %45
%49 = load float, float* %48
%50 = tail call float @llvm.fmuladd.f32(float %47, float %49, float %36)
%51 = fadd float %49, 1.000000e+00
%52 = fdiv float %50, %51
store float %52, float* %46
%53 = load float, float* %48
%54 = fadd float %53, 1.000000e+00
store float %54, float* %48
%55 = add nuw nsw i64 %31, 1
%56 = icmp eq i64 %55, %17
br i1 %56, label %27, label %30
}