#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 */
int fabs(float);
int /*<<< orphan*/ free(float *);
float *malloc(int);
void equalization_meanx(float *ref, float *modified, float *out, int w, int h,
                        int w2, int h2, int pd) {
  float *mean_ref = malloc(pd * sizeof(float));
  float *mean_modified = malloc(pd * sizeof(float));
  for (int l = 0; l < pd; l++) {
    mean_ref[l] = 0;
    for (int i = 0; i < w * h; i++)
      mean_ref[l] += ref[i + l * w * h];
    mean_ref[l] /= (float)(w * h);
    mean_modified[l] = 0;
    for (int i = 0; i < w2 * h2; i++)
      mean_modified[l] += modified[i + l * h2 * w2];
    mean_modified[l] /= (float)(w2 * h2);
  }
  float factor;
  for (int l = 0; l < pd; l++) {
    if (fabs(mean_modified[l]) < 1e-3)
      factor = 1;
    else
      factor = mean_ref[l] / mean_modified[l];
    for (int i = 0; i < w2 * h2; i++)
      out[i + l * w2 * h2] = modified[i + l * w2 * h2] * factor;
  }
  free(mean_ref);
  free(mean_modified);
}


x86 assembly (clang, -O0):
<equalization_meanx>:
push   %rbp
mov    %rsp,%rbp
sub    $0x50,%rsp
mov    0x18(%rbp),%eax
mov    0x10(%rbp),%eax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    %r8d,-0x20(%rbp)
mov    %r9d,-0x24(%rbp)
movslq 0x18(%rbp),%rax
shl    $0x2,%rax
mov    %eax,%edi
call   1050 <malloc@plt>
mov    %rax,-0x30(%rbp)
movslq 0x18(%rbp),%rax
shl    $0x2,%rax
mov    %eax,%edi
call   1050 <malloc@plt>
mov    %rax,-0x38(%rbp)
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    0x18(%rbp),%eax
jge    12a8 <equalization_meanx+0x178>
mov    -0x30(%rbp),%rax
movslq -0x3c(%rbp),%rcx
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax,%rcx,4)
movl   $0x0,-0x40(%rbp)
mov    -0x40(%rbp),%eax
mov    -0x1c(%rbp),%ecx
imul   -0x20(%rbp),%ecx
cmp    %ecx,%eax
jge    11f3 <equalization_meanx+0xc3>
mov    -0x8(%rbp),%rax
mov    -0x40(%rbp),%ecx
mov    -0x3c(%rbp),%edx
imul   -0x1c(%rbp),%edx
imul   -0x20(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x30(%rbp),%rax
movslq -0x3c(%rbp),%rcx
addss  (%rax,%rcx,4),%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    11a5 <equalization_meanx+0x75>
mov    -0x1c(%rbp),%eax
imul   -0x20(%rbp),%eax
cvtsi2ss %eax,%xmm1
mov    -0x30(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
divss  %xmm1,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x38(%rbp),%rax
movslq -0x3c(%rbp),%rcx
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax,%rcx,4)
movl   $0x0,-0x44(%rbp)
mov    -0x44(%rbp),%eax
mov    -0x24(%rbp),%ecx
imul   0x10(%rbp),%ecx
cmp    %ecx,%eax
jge    1279 <equalization_meanx+0x149>
mov    -0x10(%rbp),%rax
mov    -0x44(%rbp),%ecx
mov    -0x3c(%rbp),%edx
imul   0x10(%rbp),%edx
imul   -0x24(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x38(%rbp),%rax
movslq -0x3c(%rbp),%rcx
addss  (%rax,%rcx,4),%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    122b <equalization_meanx+0xfb>
mov    -0x24(%rbp),%eax
imul   0x10(%rbp),%eax
cvtsi2ss %eax,%xmm1
mov    -0x38(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
divss  %xmm1,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    1182 <equalization_meanx+0x52>
movl   $0x0,-0x4c(%rbp)
mov    -0x4c(%rbp),%eax
cmp    0x18(%rbp),%eax
jge    138b <equalization_meanx+0x25b>
mov    -0x38(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
call   1040 <fabs@plt>
cvtsi2sd %eax,%xmm1
movsd  0xd27(%rip),%xmm0
ucomisd %xmm1,%xmm0
jbe    12f5 <equalization_meanx+0x1c5>
movss  0xd1d(%rip),%xmm0
movss  %xmm0,-0x48(%rbp)
jmp    1314 <equalization_meanx+0x1e4>
mov    -0x30(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x38(%rbp),%rax
movslq -0x4c(%rbp),%rcx
divss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x48(%rbp)
movl   $0x0,-0x50(%rbp)
mov    -0x50(%rbp),%eax
mov    -0x24(%rbp),%ecx
imul   0x10(%rbp),%ecx
cmp    %ecx,%eax
jge    1378 <equalization_meanx+0x248>
mov    -0x10(%rbp),%rax
mov    -0x50(%rbp),%ecx
mov    -0x4c(%rbp),%edx
imul   -0x24(%rbp),%edx
imul   0x10(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mulss  -0x48(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x50(%rbp),%ecx
mov    -0x4c(%rbp),%edx
imul   -0x24(%rbp),%edx
imul   0x10(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x50(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x50(%rbp)
jmp    131b <equalization_meanx+0x1eb>
jmp    137d <equalization_meanx+0x24d>
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    12af <equalization_meanx+0x17f>
mov    -0x30(%rbp),%rdi
call   1030 <free@plt>
mov    -0x38(%rbp),%rdi
call   1030 <free@plt>
add    $0x50,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<equalization_meanx>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x78,%rsp
mov    %r9d,0x18(%rsp)
mov    %r8d,%r15d
mov    %ecx,0x8(%rsp)
mov    %rdx,0x30(%rsp)
mov    %rsi,0x10(%rsp)
mov    %rdi,%r12
mov    0xb8(%rsp),%ebx
mov    0xb0(%rsp),%r14d
lea    0x0(,%rbx,4),%ebp
mov    %ebp,%edi
call   1050 <malloc@plt>
mov    %rax,0x40(%rsp)
mov    %ebp,%edi
call   1050 <malloc@plt>
mov    %rax,%r13
test   %ebx,%ebx
jle    159e <equalization_meanx+0x46e>
imul   0x8(%rsp),%r15d
cvtsi2ss %r15d,%xmm0
imul   0x18(%rsp),%r14d
cvtsi2ss %r14d,%xmm1
mov    %rbx,0x20(%rsp)
mov    %r12,%rcx
mov    %ebx,%r12d
lea    -0x1(%r15),%rax
mov    %rax,0x18(%rsp)
lea    -0x1(%r14),%rax
mov    %rax,0x8(%rsp)
mov    %r15d,%edx
and    $0x3,%edx
mov    %r15d,%esi
and    $0xfffffffc,%esi
mov    %r14d,%edi
and    $0x3,%edi
mov    %r14d,%eax
and    $0xfffffffc,%eax
mov    %rcx,0x38(%rsp)
add    $0xc,%rcx
mov    %rcx,0x28(%rsp)
mov    0x10(%rsp),%rcx
add    $0xc,%rcx
mov    %rcx,0x48(%rsp)
xor    %r8d,%r8d
xor    %r10d,%r10d
xor    %ebx,%ebx
jmp    121e <equalization_meanx+0xee>
nopl   0x0(%rax)
divss  %xmm1,%xmm2
movss  %xmm2,0x0(%r13,%rbx,4)
add    $0x1,%rbx
add    %r15d,%r10d
add    %r14d,%r8d
cmp    %r12,%rbx
je     1333 <equalization_meanx+0x203>
xorps  %xmm2,%xmm2
xorps  %xmm3,%xmm3
test   %r15d,%r15d
jle    129e <equalization_meanx+0x16e>
movslq %r10d,%r11
cmpq   $0x3,0x18(%rsp)
jae    1240 <equalization_meanx+0x110>
xorps  %xmm3,%xmm3
xor    %r9d,%r9d
jmp    1274 <equalization_meanx+0x144>
nopl   0x0(%rax)
mov    0x28(%rsp),%rcx
lea    (%rcx,%r11,4),%rcx
xorps  %xmm3,%xmm3
xor    %r9d,%r9d
nop
addss  -0xc(%rcx,%r9,4),%xmm3
addss  -0x8(%rcx,%r9,4),%xmm3
addss  -0x4(%rcx,%r9,4),%xmm3
addss  (%rcx,%r9,4),%xmm3
add    $0x4,%r9
cmp    %r9,%rsi
jne    1250 <equalization_meanx+0x120>
test   %rdx,%rdx
je     129e <equalization_meanx+0x16e>
add    %r11,%r9
mov    0x38(%rsp),%rcx
lea    (%rcx,%r9,4),%rcx
xor    %ebp,%ebp
nopw   0x0(%rax,%rax,1)
addss  (%rcx,%rbp,4),%xmm3
add    $0x1,%rbp
cmp    %rbp,%rdx
jne    1290 <equalization_meanx+0x160>
divss  %xmm0,%xmm3
mov    0x40(%rsp),%rcx
movss  %xmm3,(%rcx,%rbx,4)
test   %r14d,%r14d
jle    1200 <equalization_meanx+0xd0>
movslq %r8d,%r11
cmpq   $0x3,0x8(%rsp)
jae    12d0 <equalization_meanx+0x1a0>
xorps  %xmm2,%xmm2
xor    %ecx,%ecx
jmp    1301 <equalization_meanx+0x1d1>
nopw   0x0(%rax,%rax,1)
mov    0x48(%rsp),%rcx
lea    (%rcx,%r11,4),%rbp
xorps  %xmm2,%xmm2
xor    %ecx,%ecx
xchg   %ax,%ax
addss  -0xc(%rbp,%rcx,4),%xmm2
addss  -0x8(%rbp,%rcx,4),%xmm2
addss  -0x4(%rbp,%rcx,4),%xmm2
addss  0x0(%rbp,%rcx,4),%xmm2
add    $0x4,%rcx
cmp    %rcx,%rax
jne    12e0 <equalization_meanx+0x1b0>
test   %rdi,%rdi
je     1200 <equalization_meanx+0xd0>
add    %r11,%rcx
mov    0x10(%rsp),%rbp
lea    0x0(,%rcx,4),%rcx
add    %rbp,%rcx
xor    %ebp,%ebp
nop
addss  (%rcx,%rbp,4),%xmm2
add    $0x1,%rbp
cmp    %rbp,%rdi
jne    1320 <equalization_meanx+0x1f0>
jmp    1200 <equalization_meanx+0xd0>
cmpl   $0x0,0x20(%rsp)
jle    159e <equalization_meanx+0x46e>
mov    0x30(%rsp),%rdx
lea    (%rdx,%r14,4),%rax
mov    0x10(%rsp),%rsi
lea    (%rsi,%r14,4),%rcx
mov    %r14d,%edi
and    $0xfffffff8,%edi
mov    %rdi,0x28(%rsp)
add    $0xfffffffffffffff8,%rdi
mov    %rdi,0x20(%rsp)
shr    $0x3,%rdi
add    $0x1,%rdi
cmp    %rdx,%rcx
seta   %cl
cmp    %rsi,%rax
seta   %al
and    %cl,%al
mov    %al,0x8(%rsp)
mov    %rdi,%rax
mov    %rdi,0x70(%rsp)
and    $0xfffffffffffffffe,%rdi
mov    %rdi,0x60(%rsp)
mov    %r14d,%ebp
and    $0x3,%ebp
lea    0x30(%rdx),%rax
mov    %rax,0x58(%rsp)
lea    0x30(%rsi),%rax
mov    %rax,0x50(%rsp)
lea    0xc(%rdx),%rax
mov    %rax,0x68(%rsp)
xor    %r15d,%r15d
xor    %ebx,%ebx
mov    %rbp,0x38(%rsp)
jmp    13d0 <equalization_meanx+0x2a0>
nopl   0x0(%rax)
add    $0x1,%rbx
add    %r14d,%r15d
cmp    %r12,%rbx
je     159e <equalization_meanx+0x46e>
movss  0x0(%r13,%rbx,4),%xmm0
movss  %xmm0,0x18(%rsp)
call   1040 <fabs@plt>
movss  0xc16(%rip),%xmm0
test   %eax,%eax
jle    13fe <equalization_meanx+0x2ce>
mov    0x40(%rsp),%rax
movss  (%rax,%rbx,4),%xmm0
divss  0x18(%rsp),%xmm0
test   %r14d,%r14d
jle    13c0 <equalization_meanx+0x290>
movslq %r15d,%rax
cmp    $0x8,%r14d
setb   %cl
or     0x8(%rsp),%cl
je     1420 <equalization_meanx+0x2f0>
xor    %ecx,%ecx
jmp    14d2 <equalization_meanx+0x3a2>
nopw   0x0(%rax,%rax,1)
movaps %xmm0,%xmm1
shufps $0x0,%xmm0,%xmm1
cmpq   $0x0,0x20(%rsp)
je     158c <equalization_meanx+0x45c>
mov    0x58(%rsp),%rcx
lea    (%rcx,%rax,4),%rdx
mov    0x50(%rsp),%rcx
lea    (%rcx,%rax,4),%rsi
mov    0x60(%rsp),%rdi
xor    %ecx,%ecx
nopl   0x0(%rax)
movups -0x30(%rsi,%rcx,4),%xmm2
movups -0x20(%rsi,%rcx,4),%xmm3
mulps  %xmm1,%xmm2
mulps  %xmm1,%xmm3
movups %xmm2,-0x30(%rdx,%rcx,4)
movups %xmm3,-0x20(%rdx,%rcx,4)
movups -0x10(%rsi,%rcx,4),%xmm2
movups (%rsi,%rcx,4),%xmm3
mulps  %xmm1,%xmm2
mulps  %xmm1,%xmm3
movups %xmm2,-0x10(%rdx,%rcx,4)
movups %xmm3,(%rdx,%rcx,4)
add    $0x10,%rcx
add    $0xfffffffffffffffe,%rdi
jne    1450 <equalization_meanx+0x320>
testb  $0x1,0x70(%rsp)
je     14c1 <equalization_meanx+0x391>
mov    %r14d,%edx
imul   %ebx,%edx
movslq %edx,%rdx
add    %rdx,%rcx
mov    0x10(%rsp),%rdx
movups (%rdx,%rcx,4),%xmm2
movups 0x10(%rdx,%rcx,4),%xmm3
mulps  %xmm1,%xmm2
mulps  %xmm1,%xmm3
mov    0x30(%rsp),%rdx
movups %xmm2,(%rdx,%rcx,4)
movups %xmm3,0x10(%rdx,%rcx,4)
mov    0x28(%rsp),%rdx
mov    %rdx,%rcx
cmp    %r14,%rdx
je     13c0 <equalization_meanx+0x290>
mov    %rcx,%rdx
not    %rdx
add    %r14,%rdx
test   %rbp,%rbp
je     1518 <equalization_meanx+0x3e8>
mov    0x30(%rsp),%rsi
lea    (%rsi,%rax,4),%rsi
mov    0x10(%rsp),%rdi
lea    (%rdi,%rax,4),%rdi
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movss  (%rdi,%rcx,4),%xmm1
mulss  %xmm0,%xmm1
movss  %xmm1,(%rsi,%rcx,4)
add    $0x1,%rcx
add    $0xffffffffffffffff,%rbp
jne    1500 <equalization_meanx+0x3d0>
cmp    $0x3,%rdx
mov    0x38(%rsp),%rbp
jb     13c0 <equalization_meanx+0x290>
mov    0x68(%rsp),%rdx
lea    (%rdx,%rax,4),%rdx
mov    0x48(%rsp),%rsi
lea    (%rsi,%rax,4),%rax
nopl   0x0(%rax)
movss  -0xc(%rax,%rcx,4),%xmm1
mulss  %xmm0,%xmm1
movss  %xmm1,-0xc(%rdx,%rcx,4)
movss  -0x8(%rax,%rcx,4),%xmm1
mulss  %xmm0,%xmm1
movss  %xmm1,-0x8(%rdx,%rcx,4)
movss  -0x4(%rax,%rcx,4),%xmm1
mulss  %xmm0,%xmm1
movss  %xmm1,-0x4(%rdx,%rcx,4)
movss  (%rax,%rcx,4),%xmm1
mulss  %xmm0,%xmm1
movss  %xmm1,(%rdx,%rcx,4)
add    $0x4,%rcx
cmp    %rcx,%r14
jne    1540 <equalization_meanx+0x410>
jmp    13c0 <equalization_meanx+0x290>
xor    %ecx,%ecx
testb  $0x1,0x70(%rsp)
jne    1493 <equalization_meanx+0x363>
jmp    14c1 <equalization_meanx+0x391>
mov    0x40(%rsp),%rdi
call   1030 <free@plt>
mov    %r13,%rdi
add    $0x78,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
jmp    1030 <free@plt>

LLVM IR (-O0):
define void @equalization_meanx(float* %0, float* %1, float* %2, i32 %3, i32 %4, i32 %5, i32 %6, i32 %7) {
%9 = alloca float*
%10 = alloca float*
%11 = alloca float*
%12 = alloca i32
%13 = alloca i32
%14 = alloca i32
%15 = alloca i32
%16 = alloca i32
%17 = alloca float*
%18 = alloca float*
%19 = alloca i32
%20 = alloca i32
%21 = alloca i32
%22 = alloca float
%23 = alloca i32
%24 = alloca i32
store float* %0, float** %9
store float* %1, float** %10
store float* %2, float** %11
store i32 %3, i32* %12
store i32 %4, i32* %13
store i32 %5, i32* %14
store i32 %6, i32* %15
store i32 %7, i32* %16
%25 = load i32, i32* %16
%26 = sext i32 %25 to i64
%27 = mul i64 %26, 4
%28 = trunc i64 %27 to i32
%29 = call float* @malloc(i32 %28)
store float* %29, float** %17
%30 = load i32, i32* %16
%31 = sext i32 %30 to i64
%32 = mul i64 %31, 4
%33 = trunc i64 %32 to i32
%34 = call float* @malloc(i32 %33)
store float* %34, float** %18
store i32 0, i32* %19
br label %35
35: ; preds = %124, %8
%36 = load i32, i32* %19
%37 = load i32, i32* %16
%38 = icmp slt i32 %36, %37
br i1 %38, label %39, label %127
39: ; preds = %35
%40 = load float*, float** %17
%41 = load i32, i32* %19
%42 = sext i32 %41 to i64
%43 = getelementptr inbounds float, float* %40, i64 %42
store float 0.000000e+00, float* %43
store i32 0, i32* %20
br label %44
44: ; preds = %68, %39
%45 = load i32, i32* %20
%46 = load i32, i32* %12
%47 = load i32, i32* %13
%48 = mul nsw i32 %46, %47
%49 = icmp slt i32 %45, %48
br i1 %49, label %50, label %71
50: ; preds = %44
%51 = load float*, float** %9
%52 = load i32, i32* %20
%53 = load i32, i32* %19
%54 = load i32, i32* %12
%55 = mul nsw i32 %53, %54
%56 = load i32, i32* %13
%57 = mul nsw i32 %55, %56
%58 = add nsw i32 %52, %57
%59 = sext i32 %58 to i64
%60 = getelementptr inbounds float, float* %51, i64 %59
%61 = load float, float* %60
%62 = load float*, float** %17
%63 = load i32, i32* %19
%64 = sext i32 %63 to i64
%65 = getelementptr inbounds float, float* %62, i64 %64
%66 = load float, float* %65
%67 = fadd float %66, %61
store float %67, float* %65
br label %68
68: ; preds = %50
%69 = load i32, i32* %20
%70 = add nsw i32 %69, 1
store i32 %70, i32* %20
br label %44
71: ; preds = %44
%72 = load i32, i32* %12
%73 = load i32, i32* %13
%74 = mul nsw i32 %72, %73
%75 = sitofp i32 %74 to float
%76 = load float*, float** %17
%77 = load i32, i32* %19
%78 = sext i32 %77 to i64
%79 = getelementptr inbounds float, float* %76, i64 %78
%80 = load float, float* %79
%81 = fdiv float %80, %75
store float %81, float* %79
%82 = load float*, float** %18
%83 = load i32, i32* %19
%84 = sext i32 %83 to i64
%85 = getelementptr inbounds float, float* %82, i64 %84
store float 0.000000e+00, float* %85
store i32 0, i32* %21
br label %86
86: ; preds = %110, %71
%87 = load i32, i32* %21
%88 = load i32, i32* %14
%89 = load i32, i32* %15
%90 = mul nsw i32 %88, %89
%91 = icmp slt i32 %87, %90
br i1 %91, label %92, label %113
92: ; preds = %86
%93 = load float*, float** %10
%94 = load i32, i32* %21
%95 = load i32, i32* %19
%96 = load i32, i32* %15
%97 = mul nsw i32 %95, %96
%98 = load i32, i32* %14
%99 = mul nsw i32 %97, %98
%100 = add nsw i32 %94, %99
%101 = sext i32 %100 to i64
%102 = getelementptr inbounds float, float* %93, i64 %101
%103 = load float, float* %102
%104 = load float*, float** %18
%105 = load i32, i32* %19
%106 = sext i32 %105 to i64
%107 = getelementptr inbounds float, float* %104, i64 %106
%108 = load float, float* %107
%109 = fadd float %108, %103
store float %109, float* %107
br label %110
110: ; preds = %92
%111 = load i32, i32* %21
%112 = add nsw i32 %111, 1
store i32 %112, i32* %21
br label %86
113: ; preds = %86
%114 = load i32, i32* %14
%115 = load i32, i32* %15
%116 = mul nsw i32 %114, %115
%117 = sitofp i32 %116 to float
%118 = load float*, float** %18
%119 = load i32, i32* %19
%120 = sext i32 %119 to i64
%121 = getelementptr inbounds float, float* %118, i64 %120
%122 = load float, float* %121
%123 = fdiv float %122, %117
store float %123, float* %121
br label %124
124: ; preds = %113
%125 = load i32, i32* %19
%126 = add nsw i32 %125, 1
store i32 %126, i32* %19
br label %35
127: ; preds = %35
store i32 0, i32* %23
br label %128
128: ; preds = %189, %127
%129 = load i32, i32* %23
%130 = load i32, i32* %16
%131 = icmp slt i32 %129, %130
br i1 %131, label %132, label %192
132: ; preds = %128
%133 = load float*, float** %18
%134 = load i32, i32* %23
%135 = sext i32 %134 to i64
%136 = getelementptr inbounds float, float* %133, i64 %135
%137 = load float, float* %136
%138 = call i32 @fabs(float %137)
%139 = sitofp i32 %138 to double
%140 = fcmp olt double %139, 1.000000e-03
br i1 %140, label %141, label %142
141: ; preds = %132
store float 1.000000e+00, float* %22
br label %154
142: ; preds = %132
%143 = load float*, float** %17
%144 = load i32, i32* %23
%145 = sext i32 %144 to i64
%146 = getelementptr inbounds float, float* %143, i64 %145
%147 = load float, float* %146
%148 = load float*, float** %18
%149 = load i32, i32* %23
%150 = sext i32 %149 to i64
%151 = getelementptr inbounds float, float* %148, i64 %150
%152 = load float, float* %151
%153 = fdiv float %147, %152
store float %153, float* %22
br label %154
154: ; preds = %142, %141
store i32 0, i32* %24
br label %155
155: ; preds = %185, %154
%156 = load i32, i32* %24
%157 = load i32, i32* %14
%158 = load i32, i32* %15
%159 = mul nsw i32 %157, %158
%160 = icmp slt i32 %156, %159
br i1 %160, label %161, label %188
161: ; preds = %155
%162 = load float*, float** %10
%163 = load i32, i32* %24
%164 = load i32, i32* %23
%165 = load i32, i32* %14
%166 = mul nsw i32 %164, %165
%167 = load i32, i32* %15
%168 = mul nsw i32 %166, %167
%169 = add nsw i32 %163, %168
%170 = sext i32 %169 to i64
%171 = getelementptr inbounds float, float* %162, i64 %170
%172 = load float, float* %171
%173 = load float, float* %22
%174 = fmul float %172, %173
%175 = load float*, float** %11
%176 = load i32, i32* %24
%177 = load i32, i32* %23
%178 = load i32, i32* %14
%179 = mul nsw i32 %177, %178
%180 = load i32, i32* %15
%181 = mul nsw i32 %179, %180
%182 = add nsw i32 %176, %181
%183 = sext i32 %182 to i64
%184 = getelementptr inbounds float, float* %175, i64 %183
store float %174, float* %184
br label %185
185: ; preds = %161
%186 = load i32, i32* %24
%187 = add nsw i32 %186, 1
store i32 %187, i32* %24
br label %155
188: ; preds = %155
br label %189
189: ; preds = %188
%190 = load i32, i32* %23
%191 = add nsw i32 %190, 1
store i32 %191, i32* %23
br label %128
192: ; preds = %128
%193 = load float*, float** %17
%194 = call i32 @free(float* %193)
%195 = load float*, float** %18
%196 = call i32 @free(float* %195)
ret void
}

LLVM IR (-O2):
define void @equalization_meanx(float* nocapture %0, float* nocapture %1, float* nocapture writeonly %2, i32 %3, i32 %4, i32 %5, i32 %6, i32 %7) {
%9 = shl i32 %7, 2
%10 = tail call float* @malloc(i32 %9)
%11 = tail call float* @malloc(i32 %9)
%12 = icmp sgt i32 %7, 0
br i1 %12, label %13, label %159
13: ; preds = %8
%14 = mul i32 %4, %3
%15 = icmp sgt i32 %14, 0
%16 = sitofp i32 %14 to float
%17 = mul i32 %6, %5
%18 = icmp sgt i32 %17, 0
%19 = sitofp i32 %17 to float
%20 = zext i32 %7 to i64
%21 = zext i32 %14 to i64
%22 = zext i32 %17 to i64
%23 = add nsw i64 %21, -1
%24 = add nsw i64 %22, -1
%25 = and i64 %21, 3
%26 = icmp ult i64 %23, 3
%27 = and i64 %21, 4294967292
%28 = icmp eq i64 %25, 0
%29 = and i64 %22, 3
%30 = icmp ult i64 %24, 3
%31 = and i64 %22, 4294967292
%32 = icmp eq i64 %29, 0
br label %57
33: ; preds = %128
br i1 %12, label %34, label %159
34: ; preds = %33
%35 = mul i32 %6, %5
%36 = icmp sgt i32 %35, 0
%37 = zext i32 %7 to i64
%38 = zext i32 %35 to i64
%39 = getelementptr float, float* %2, i64 %22
%40 = getelementptr float, float* %1, i64 %22
%41 = and i64 %22, 4294967288
%42 = add nsw i64 %41, -8
%43 = lshr exact i64 %42, 3
%44 = add nuw nsw i64 %43, 1
%45 = icmp ult i32 %17, 8
%46 = icmp ugt float* %40, %2
%47 = icmp ugt float* %39, %1
%48 = and i1 %46, %47
%49 = and i64 %22, 4294967288
%50 = and i64 %44, 1
%51 = icmp eq i64 %42, 0
%52 = and i64 %44, 4611686018427387902
%53 = icmp eq i64 %50, 0
%54 = icmp eq i64 %49, %22
%55 = and i64 %22, 3
%56 = icmp eq i64 %55, 0
br label %162
57: ; preds = %13, %128
%58 = phi i64 [ 0, %13 ], [ %131, %128 ]
%59 = getelementptr inbounds float, float* %10, i64 %58
br i1 %15, label %60, label %79
60: ; preds = %57
%61 = trunc i64 %58 to i32
%62 = mul i32 %14, %61
%63 = sext i32 %62 to i64
br i1 %26, label %64, label %87
64: ; preds = %87, %60
%65 = phi float [ undef, %60 ], [ %109, %87 ]
%66 = phi i64 [ 0, %60 ], [ %110, %87 ]
%67 = phi float [ 0.000000e+00, %60 ], [ %109, %87 ]
br i1 %28, label %79, label %68
68: ; preds = %64, %68
%69 = phi i64 [ %76, %68 ], [ %66, %64 ]
%70 = phi float [ %75, %68 ], [ %67, %64 ]
%71 = phi i64 [ %77, %68 ], [ 0, %64 ]
%72 = add nsw i64 %69, %63
%73 = getelementptr inbounds float, float* %0, i64 %72
%74 = load float, float* %73
%75 = fadd float %74, %70
%76 = add nuw nsw i64 %69, 1
%77 = add i64 %71, 1
%78 = icmp eq i64 %77, %25
br i1 %78, label %79, label %68
79: ; preds = %64, %68, %57
%80 = phi float [ 0.000000e+00, %57 ], [ %65, %64 ], [ %75, %68 ]
%81 = fdiv float %80, %16
store float %81, float* %59
%82 = getelementptr inbounds float, float* %11, i64 %58
br i1 %18, label %83, label %128
83: ; preds = %79
%84 = trunc i64 %58 to i32
%85 = mul i32 %17, %84
%86 = sext i32 %85 to i64
br i1 %30, label %113, label %133
87: ; preds = %60, %87
%88 = phi i64 [ %110, %87 ], [ 0, %60 ]
%89 = phi float [ %109, %87 ], [ 0.000000e+00, %60 ]
%90 = phi i64 [ %111, %87 ], [ 0, %60 ]
%91 = add nsw i64 %88, %63
%92 = getelementptr inbounds float, float* %0, i64 %91
%93 = load float, float* %92
%94 = fadd float %93, %89
%95 = or i64 %88, 1
%96 = add nsw i64 %95, %63
%97 = getelementptr inbounds float, float* %0, i64 %96
%98 = load float, float* %97
%99 = fadd float %98, %94
%100 = or i64 %88, 2
%101 = add nsw i64 %100, %63
%102 = getelementptr inbounds float, float* %0, i64 %101
%103 = load float, float* %102
%104 = fadd float %103, %99
%105 = or i64 %88, 3
%106 = add nsw i64 %105, %63
%107 = getelementptr inbounds float, float* %0, i64 %106
%108 = load float, float* %107
%109 = fadd float %108, %104
%110 = add nuw nsw i64 %88, 4
%111 = add i64 %90, 4
%112 = icmp eq i64 %111, %27
br i1 %112, label %64, label %87
113: ; preds = %133, %83
%114 = phi float [ undef, %83 ], [ %155, %133 ]
%115 = phi i64 [ 0, %83 ], [ %156, %133 ]
%116 = phi float [ 0.000000e+00, %83 ], [ %155, %133 ]
br i1 %32, label %128, label %117
117: ; preds = %113, %117
%118 = phi i64 [ %125, %117 ], [ %115, %113 ]
%119 = phi float [ %124, %117 ], [ %116, %113 ]
%120 = phi i64 [ %126, %117 ], [ 0, %113 ]
%121 = add nsw i64 %118, %86
%122 = getelementptr inbounds float, float* %1, i64 %121
%123 = load float, float* %122
%124 = fadd float %123, %119
%125 = add nuw nsw i64 %118, 1
%126 = add i64 %120, 1
%127 = icmp eq i64 %126, %29
br i1 %127, label %128, label %117
128: ; preds = %113, %117, %79
%129 = phi float [ 0.000000e+00, %79 ], [ %114, %113 ], [ %124, %117 ]
%130 = fdiv float %129, %19
store float %130, float* %82
%131 = add nuw nsw i64 %58, 1
%132 = icmp eq i64 %131, %20
br i1 %132, label %33, label %57
133: ; preds = %83, %133
%134 = phi i64 [ %156, %133 ], [ 0, %83 ]
%135 = phi float [ %155, %133 ], [ 0.000000e+00, %83 ]
%136 = phi i64 [ %157, %133 ], [ 0, %83 ]
%137 = add nsw i64 %134, %86
%138 = getelementptr inbounds float, float* %1, i64 %137
%139 = load float, float* %138
%140 = fadd float %139, %135
%141 = or i64 %134, 1
%142 = add nsw i64 %141, %86
%143 = getelementptr inbounds float, float* %1, i64 %142
%144 = load float, float* %143
%145 = fadd float %144, %140
%146 = or i64 %134, 2
%147 = add nsw i64 %146, %86
%148 = getelementptr inbounds float, float* %1, i64 %147
%149 = load float, float* %148
%150 = fadd float %149, %145
%151 = or i64 %134, 3
%152 = add nsw i64 %151, %86
%153 = getelementptr inbounds float, float* %1, i64 %152
%154 = load float, float* %153
%155 = fadd float %154, %150
%156 = add nuw nsw i64 %134, 4
%157 = add i64 %136, 4
%158 = icmp eq i64 %157, %31
br i1 %158, label %113, label %133
159: ; preds = %252, %8, %33
%160 = tail call i32 @free(float* %10)
%161 = tail call i32 @free(float* %11)
ret void
162: ; preds = %34, %252
%163 = phi i64 [ 0, %34 ], [ %253, %252 ]
%164 = getelementptr inbounds float, float* %11, i64 %163
%165 = load float, float* %164
%166 = tail call i32 @fabs(float %165)
%167 = icmp slt i32 %166, 1
br i1 %167, label %172, label %168
168: ; preds = %162
%169 = getelementptr inbounds float, float* %10, i64 %163
%170 = load float, float* %169
%171 = fdiv float %170, %165
br label %172
172: ; preds = %162, %168
%173 = phi float [ %171, %168 ], [ 1.000000e+00, %162 ]
br i1 %36, label %174, label %252
174: ; preds = %172
%175 = trunc i64 %163 to i32
%176 = mul i32 %35, %175
%177 = sext i32 %176 to i64
%178 = or i1 %45, %48
br i1 %178, label %234, label %179
179: ; preds = %174
%180 = insertelement <4 x float> poison, float %173, i64 0
%181 = shufflevector <4 x float> %180, <4 x float> poison, <4 x i32> zeroinitializer
%182 = insertelement <4 x float> poison, float %173, i64 0
%183 = shufflevector <4 x float> %182, <4 x float> poison, <4 x i32> zeroinitializer
br i1 %51, label %217, label %184
184: ; preds = %179, %184
%185 = phi i64 [ %214, %184 ], [ 0, %179 ]
%186 = phi i64 [ %215, %184 ], [ 0, %179 ]
%187 = add nsw i64 %185, %177
%188 = getelementptr inbounds float, float* %1, i64 %187
%189 = bitcast float* %188 to <4 x float>*
%190 = load <4 x float>, <4 x float>* %189
%191 = getelementptr inbounds float, float* %188, i64 4
%192 = bitcast float* %191 to <4 x float>*
%193 = load <4 x float>, <4 x float>* %192
%194 = fmul <4 x float> %181, %190
%195 = fmul <4 x float> %183, %193
%196 = getelementptr inbounds float, float* %2, i64 %187
%197 = bitcast float* %196 to <4 x float>*
store <4 x float> %194, <4 x float>* %197
%198 = getelementptr inbounds float, float* %196, i64 4
%199 = bitcast float* %198 to <4 x float>*
store <4 x float> %195, <4 x float>* %199
%200 = or i64 %185, 8
%201 = add nsw i64 %200, %177
%202 = getelementptr inbounds float, float* %1, i64 %201
%203 = bitcast float* %202 to <4 x float>*
%204 = load <4 x float>, <4 x float>* %203
%205 = getelementptr inbounds float, float* %202, i64 4
%206 = bitcast float* %205 to <4 x float>*
%207 = load <4 x float>, <4 x float>* %206
%208 = fmul <4 x float> %181, %204
%209 = fmul <4 x float> %183, %207
%210 = getelementptr inbounds float, float* %2, i64 %201
%211 = bitcast float* %210 to <4 x float>*
store <4 x float> %208, <4 x float>* %211
%212 = getelementptr inbounds float, float* %210, i64 4
%213 = bitcast float* %212 to <4 x float>*
store <4 x float> %209, <4 x float>* %213
%214 = add nuw i64 %185, 16
%215 = add i64 %186, 2
%216 = icmp eq i64 %215, %52
br i1 %216, label %217, label %184
217: ; preds = %184, %179
%218 = phi i64 [ 0, %179 ], [ %214, %184 ]
br i1 %53, label %233, label %219
219: ; preds = %217
%220 = add nsw i64 %218, %177
%221 = getelementptr inbounds float, float* %1, i64 %220
%222 = bitcast float* %221 to <4 x float>*
%223 = load <4 x float>, <4 x float>* %222
%224 = getelementptr inbounds float, float* %221, i64 4
%225 = bitcast float* %224 to <4 x float>*
%226 = load <4 x float>, <4 x float>* %225
%227 = fmul <4 x float> %181, %223
%228 = fmul <4 x float> %183, %226
%229 = getelementptr inbounds float, float* %2, i64 %220
%230 = bitcast float* %229 to <4 x float>*
store <4 x float> %227, <4 x float>* %230
%231 = getelementptr inbounds float, float* %229, i64 4
%232 = bitcast float* %231 to <4 x float>*
store <4 x float> %228, <4 x float>* %232
br label %233
233: ; preds = %217, %219
br i1 %54, label %252, label %234
234: ; preds = %174, %233
%235 = phi i64 [ 0, %174 ], [ %49, %233 ]
%236 = xor i64 %235, -1
%237 = add nsw i64 %236, %22
br i1 %56, label %249, label %238
238: ; preds = %234, %238
%239 = phi i64 [ %246, %238 ], [ %235, %234 ]
%240 = phi i64 [ %247, %238 ], [ 0, %234 ]
%241 = add nsw i64 %239, %177
%242 = getelementptr inbounds float, float* %1, i64 %241
%243 = load float, float* %242
%244 = fmul float %173, %243
%245 = getelementptr inbounds float, float* %2, i64 %241
store float %244, float* %245
%246 = add nuw nsw i64 %239, 1
%247 = add i64 %240, 1
%248 = icmp eq i64 %247, %55
br i1 %248, label %249, label %238
249: ; preds = %238, %234
%250 = phi i64 [ %235, %234 ], [ %246, %238 ]
%251 = icmp ult i64 %237, 3
br i1 %251, label %252, label %255
252: ; preds = %249, %255, %233, %172
%253 = add nuw nsw i64 %163, 1
%254 = icmp eq i64 %253, %37
br i1 %254, label %159, label %162
255: ; preds = %249, %255
%256 = phi i64 [ %280, %255 ], [ %250, %249 ]
%257 = add nsw i64 %256, %177
%258 = getelementptr inbounds float, float* %1, i64 %257
%259 = load float, float* %258
%260 = fmul float %173, %259
%261 = getelementptr inbounds float, float* %2, i64 %257
store float %260, float* %261
%262 = add nuw nsw i64 %256, 1
%263 = add nsw i64 %262, %177
%264 = getelementptr inbounds float, float* %1, i64 %263
%265 = load float, float* %264
%266 = fmul float %173, %265
%267 = getelementptr inbounds float, float* %2, i64 %263
store float %266, float* %267
%268 = add nuw nsw i64 %256, 2
%269 = add nsw i64 %268, %177
%270 = getelementptr inbounds float, float* %1, i64 %269
%271 = load float, float* %270
%272 = fmul float %173, %271
%273 = getelementptr inbounds float, float* %2, i64 %269
store float %272, float* %273
%274 = add nuw nsw i64 %256, 3
%275 = add nsw i64 %274, %177
%276 = getelementptr inbounds float, float* %1, i64 %275
%277 = load float, float* %276
%278 = fmul float %173, %277
%279 = getelementptr inbounds float, float* %2, i64 %275
store float %278, float* %279
%280 = add nuw nsw i64 %256, 4
%281 = icmp eq i64 %280, %38
br i1 %281, label %252, label %255
}