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


x86 assembly (clang, -O0):
<equalization_affine>:
push   %rbp
mov    %rsp,%rbp
sub    $0x70,%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)
movslq 0x18(%rbp),%rax
shl    $0x2,%rax
mov    %eax,%edi
call   1050 <malloc@plt>
mov    %rax,-0x40(%rbp)
movslq 0x18(%rbp),%rax
shl    $0x2,%rax
mov    %eax,%edi
call   1050 <malloc@plt>
mov    %rax,-0x48(%rbp)
movslq 0x18(%rbp),%rax
shl    $0x2,%rax
mov    %eax,%edi
call   1050 <malloc@plt>
mov    %rax,-0x50(%rbp)
movl   $0x0,-0x54(%rbp)
mov    -0x54(%rbp),%eax
cmp    0x18(%rbp),%eax
jge    12f1 <equalization_affine+0x1b1>
mov    -0x30(%rbp),%rax
movslq -0x54(%rbp),%rcx
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax,%rcx,4)
movl   $0x0,-0x58(%rbp)
mov    -0x58(%rbp),%eax
mov    -0x1c(%rbp),%ecx
imul   -0x20(%rbp),%ecx
cmp    %ecx,%eax
jge    123c <equalization_affine+0xfc>
mov    -0x8(%rbp),%rax
mov    -0x58(%rbp),%ecx
mov    -0x54(%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 -0x54(%rbp),%rcx
addss  (%rax,%rcx,4),%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x58(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x58(%rbp)
jmp    11ee <equalization_affine+0xae>
mov    -0x1c(%rbp),%eax
imul   -0x20(%rbp),%eax
cvtsi2ss %eax,%xmm1
mov    -0x30(%rbp),%rax
movslq -0x54(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
divss  %xmm1,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x38(%rbp),%rax
movslq -0x54(%rbp),%rcx
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax,%rcx,4)
movl   $0x0,-0x5c(%rbp)
mov    -0x5c(%rbp),%eax
mov    -0x24(%rbp),%ecx
imul   0x10(%rbp),%ecx
cmp    %ecx,%eax
jge    12c2 <equalization_affine+0x182>
mov    -0x10(%rbp),%rax
mov    -0x5c(%rbp),%ecx
mov    -0x54(%rbp),%edx
imul   -0x24(%rbp),%edx
imul   0x10(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x38(%rbp),%rax
movslq -0x54(%rbp),%rcx
addss  (%rax,%rcx,4),%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x5c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x5c(%rbp)
jmp    1274 <equalization_affine+0x134>
mov    -0x24(%rbp),%eax
imul   0x10(%rbp),%eax
cvtsi2ss %eax,%xmm1
mov    -0x38(%rbp),%rax
movslq -0x54(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
divss  %xmm1,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x54(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x54(%rbp)
jmp    11cb <equalization_affine+0x8b>
movl   $0x0,-0x64(%rbp)
mov    -0x64(%rbp),%eax
cmp    0x18(%rbp),%eax
jge    1567 <equalization_affine+0x427>
mov    -0x40(%rbp),%rax
movslq -0x64(%rbp),%rcx
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax,%rcx,4)
movl   $0x0,-0x68(%rbp)
mov    -0x68(%rbp),%eax
mov    -0x1c(%rbp),%ecx
imul   -0x20(%rbp),%ecx
cmp    %ecx,%eax
jge    138d <equalization_affine+0x24d>
mov    -0x8(%rbp),%rax
mov    -0x68(%rbp),%ecx
mov    -0x64(%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 -0x64(%rbp),%rcx
subss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x60(%rbp)
movss  -0x60(%rbp),%xmm0
movss  -0x60(%rbp),%xmm2
mov    -0x40(%rbp),%rax
movslq -0x64(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x68(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x68(%rbp)
jmp    131b <equalization_affine+0x1db>
mov    -0x1c(%rbp),%eax
imul   -0x20(%rbp),%eax
cvtsi2ss %eax,%xmm1
mov    -0x40(%rbp),%rax
movslq -0x64(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
divss  %xmm1,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x48(%rbp),%rax
movslq -0x64(%rbp),%rcx
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax,%rcx,4)
movl   $0x0,-0x6c(%rbp)
mov    -0x6c(%rbp),%eax
mov    -0x24(%rbp),%ecx
imul   0x10(%rbp),%ecx
cmp    %ecx,%eax
jge    1437 <equalization_affine+0x2f7>
mov    -0x10(%rbp),%rax
mov    -0x6c(%rbp),%ecx
mov    -0x64(%rbp),%edx
imul   -0x24(%rbp),%edx
imul   0x10(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x38(%rbp),%rax
movslq -0x64(%rbp),%rcx
subss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x60(%rbp)
movss  -0x60(%rbp),%xmm0
movss  -0x60(%rbp),%xmm2
mov    -0x48(%rbp),%rax
movslq -0x64(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x6c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x6c(%rbp)
jmp    13c5 <equalization_affine+0x285>
mov    -0x24(%rbp),%eax
imul   0x10(%rbp),%eax
cvtsi2ss %eax,%xmm1
mov    -0x48(%rbp),%rax
movslq -0x64(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
divss  %xmm1,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x48(%rbp),%rax
movslq -0x64(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
call   1040 <fabs@plt>
cvtsi2sd %eax,%xmm1
movsd  0xb8a(%rip),%xmm0
ucomisd %xmm1,%xmm0
jbe    149a <equalization_affine+0x35a>
mov    -0x50(%rbp),%rax
movslq -0x64(%rbp),%rcx
movss  0xb78(%rip),%xmm0
movss  %xmm0,(%rax,%rcx,4)
jmp    14c6 <equalization_affine+0x386>
mov    -0x40(%rbp),%rax
movslq -0x64(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x48(%rbp),%rax
movslq -0x64(%rbp),%rcx
divss  (%rax,%rcx,4),%xmm0
call   1060 <sqrt@plt>
mov    -0x50(%rbp),%rax
movslq -0x64(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
movl   $0x0,-0x70(%rbp)
mov    -0x70(%rbp),%eax
mov    -0x24(%rbp),%ecx
imul   0x10(%rbp),%ecx
cmp    %ecx,%eax
jge    1554 <equalization_affine+0x414>
mov    -0x50(%rbp),%rax
movslq -0x64(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x70(%rbp),%ecx
mov    -0x64(%rbp),%edx
imul   0x10(%rbp),%edx
imul   -0x24(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm2
mov    -0x38(%rbp),%rax
movslq -0x64(%rbp),%rcx
subss  (%rax,%rcx,4),%xmm2
mov    -0x30(%rbp),%rax
movslq -0x64(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
mov    -0x70(%rbp),%ecx
mov    -0x64(%rbp),%edx
imul   -0x24(%rbp),%edx
imul   0x10(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x70(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x70(%rbp)
jmp    14cd <equalization_affine+0x38d>
jmp    1559 <equalization_affine+0x419>
mov    -0x64(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x64(%rbp)
jmp    12f8 <equalization_affine+0x1b8>
mov    -0x30(%rbp),%rdi
call   1030 <free@plt>
mov    -0x38(%rbp),%rdi
call   1030 <free@plt>
mov    -0x40(%rbp),%rdi
call   1030 <free@plt>
mov    -0x48(%rbp),%rdi
call   1030 <free@plt>
mov    -0x50(%rbp),%rdi
call   1030 <free@plt>
add    $0x70,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<equalization_affine>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0xf8,%rsp
mov    %r9d,%r15d
mov    %r8d,%ebp
mov    %ecx,%r13d
mov    %rdx,0x28(%rsp)
mov    %rsi,0x10(%rsp)
mov    %rdi,0x48(%rsp)
mov    0x138(%rsp),%r14d
mov    0x130(%rsp),%r12d
lea    0x0(,%r14,4),%ebx
mov    %ebx,%edi
call   1050 <malloc@plt>
mov    %rax,0x18(%rsp)
mov    %ebx,%edi
call   1050 <malloc@plt>
mov    %rax,0x20(%rsp)
mov    %ebx,%edi
call   1050 <malloc@plt>
mov    %rax,0x70(%rsp)
mov    %ebx,%edi
call   1050 <malloc@plt>
mov    %rax,0x68(%rsp)
mov    %ebx,%edi
call   1050 <malloc@plt>
mov    %rax,0x60(%rsp)
test   %r14d,%r14d
jle    18d1 <equalization_affine+0x791>
mov    %rbp,%r11
imul   %r13d,%r11d
cvtsi2ss %r11d,%xmm3
imul   %r15d,%r12d
cvtsi2ss %r12d,%xmm8
mov    %r14d,%eax
mov    %rax,0x98(%rsp)
lea    -0x1(%r11),%rax
mov    %rax,0x90(%rsp)
lea    -0x1(%r12),%rax
mov    %rax,0x88(%rsp)
mov    %r11d,%eax
and    $0x3,%eax
mov    %r11d,%ecx
and    $0xfffffffc,%ecx
mov    %r12d,%edx
and    $0x3,%edx
mov    %r12d,%esi
and    $0xfffffffc,%esi
mov    0x48(%rsp),%r15
lea    0xc(%r15),%rdi
mov    %rdi,0x58(%rsp)
mov    0x10(%rsp),%rdi
add    $0xc,%rdi
mov    %rdi,0x50(%rsp)
xor    %r9d,%r9d
xor    %r10d,%r10d
xor    %ebp,%ebp
mov    0x18(%rsp),%r13
jmp    1277 <equalization_affine+0x137>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
divss  %xmm8,%xmm0
mov    0x20(%rsp),%rdi
movss  %xmm0,(%rdi,%rbp,4)
add    $0x1,%rbp
add    %r11d,%r10d
add    %r12d,%r9d
cmp    0x98(%rsp),%rbp
je     1383 <equalization_affine+0x243>
xorps  %xmm0,%xmm0
xorps  %xmm1,%xmm1
test   %r11d,%r11d
jle    12ee <equalization_affine+0x1ae>
movslq %r10d,%r8
cmpq   $0x3,0x90(%rsp)
jae    12a0 <equalization_affine+0x160>
xorps  %xmm1,%xmm1
xor    %edi,%edi
jmp    12d0 <equalization_affine+0x190>
nopw   0x0(%rax,%rax,1)
mov    0x58(%rsp),%rdi
lea    (%rdi,%r8,4),%rbx
xorps  %xmm1,%xmm1
xor    %edi,%edi
xchg   %ax,%ax
addss  -0xc(%rbx,%rdi,4),%xmm1
addss  -0x8(%rbx,%rdi,4),%xmm1
addss  -0x4(%rbx,%rdi,4),%xmm1
addss  (%rbx,%rdi,4),%xmm1
add    $0x4,%rdi
cmp    %rdi,%rcx
jne    12b0 <equalization_affine+0x170>
test   %rax,%rax
je     12ee <equalization_affine+0x1ae>
add    %r8,%rdi
lea    (%r15,%rdi,4),%rdi
xor    %ebx,%ebx
xchg   %ax,%ax
addss  (%rdi,%rbx,4),%xmm1
add    $0x1,%rbx
cmp    %rbx,%rax
jne    12e0 <equalization_affine+0x1a0>
divss  %xmm3,%xmm1
movss  %xmm1,0x0(%r13,%rbp,4)
test   %r12d,%r12d
jle    1250 <equalization_affine+0x110>
movslq %r9d,%r8
cmpq   $0x3,0x88(%rsp)
jae    1320 <equalization_affine+0x1e0>
xorps  %xmm0,%xmm0
xor    %edi,%edi
jmp    1350 <equalization_affine+0x210>
nopw   0x0(%rax,%rax,1)
mov    0x50(%rsp),%rdi
lea    (%rdi,%r8,4),%rbx
xorps  %xmm0,%xmm0
xor    %edi,%edi
xchg   %ax,%ax
addss  -0xc(%rbx,%rdi,4),%xmm0
addss  -0x8(%rbx,%rdi,4),%xmm0
addss  -0x4(%rbx,%rdi,4),%xmm0
addss  (%rbx,%rdi,4),%xmm0
add    $0x4,%rdi
cmp    %rdi,%rsi
jne    1330 <equalization_affine+0x1f0>
test   %rdx,%rdx
je     1250 <equalization_affine+0x110>
add    %r8,%rdi
mov    0x10(%rsp),%rbx
lea    (%rbx,%rdi,4),%rdi
xor    %ebx,%ebx
nopw   0x0(%rax,%rax,1)
addss  (%rdi,%rbx,4),%xmm0
add    $0x1,%rbx
cmp    %rbx,%rdx
jne    1370 <equalization_affine+0x230>
jmp    1250 <equalization_affine+0x110>
test   %r14d,%r14d
jle    18d1 <equalization_affine+0x791>
mov    0x10(%rsp),%rsi
lea    (%rsi,%r12,4),%rax
mov    %r12d,%ecx
and    $0xfffffff8,%ecx
mov    %rcx,0xe0(%rsp)
add    $0xfffffffffffffff8,%rcx
mov    %rcx,0xd8(%rsp)
mov    %rcx,%rdi
shr    $0x3,%rdi
add    $0x1,%rdi
mov    %r11d,%ebx
and    $0x3,%ebx
mov    %r11d,%ecx
and    $0xfffffffc,%ecx
mov    %rcx,0xf0(%rsp)
mov    %r12d,%r14d
and    $0x3,%r14d
mov    %r12d,%r15d
and    $0xfffffffc,%r15d
mov    0x28(%rsp),%rdx
cmp    %rdx,%rax
seta   %al
lea    (%rdx,%r12,4),%rcx
cmp    %rsi,%rcx
seta   %cl
and    %al,%cl
mov    %cl,0xf(%rsp)
mov    %rdi,%rax
mov    %rdi,0xd0(%rsp)
and    $0xfffffffffffffffe,%rdi
mov    %rdi,0xb0(%rsp)
mov    %r12,%rax
neg    %rax
mov    %rax,0xe8(%rsp)
lea    0x30(%rdx),%rax
mov    %rax,0xa8(%rsp)
lea    0x30(%rsi),%rax
mov    %rax,0xa0(%rsp)
lea    0x4(%rdx),%rax
mov    %rax,0xc0(%rsp)
lea    0x4(%rsi),%rax
mov    %rax,0xb8(%rsp)
xor    %ebp,%ebp
xor    %r8d,%r8d
xor    %r13d,%r13d
mov    %r11,0x78(%rsp)
movss  %xmm3,0x38(%rsp)
movss  %xmm8,0x34(%rsp)
mov    %r12,0xc8(%rsp)
jmp    148e <equalization_affine+0x34e>
add    $0x1,%r13
add    %r11d,%r8d
add    %r12d,%ebp
cmp    0x98(%rsp),%r13
movss  0x38(%rsp),%xmm3
je     18d1 <equalization_affine+0x791>
xorps  %xmm4,%xmm4
test   %r11d,%r11d
mov    0xf0(%rsp),%rsi
jle    156a <equalization_affine+0x42a>
movslq %r8d,%rax
mov    0x18(%rsp),%rcx
movss  (%rcx,%r13,4),%xmm0
cmpq   $0x3,0x90(%rsp)
jae    14d0 <equalization_affine+0x390>
xorps  %xmm4,%xmm4
xor    %ecx,%ecx
jmp    1530 <equalization_affine+0x3f0>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
mov    0x58(%rsp),%rcx
lea    (%rcx,%rax,4),%rdx
xorps  %xmm4,%xmm4
xor    %ecx,%ecx
xchg   %ax,%ax
movss  -0xc(%rdx,%rcx,4),%xmm1
movss  -0x8(%rdx,%rcx,4),%xmm2
subss  %xmm0,%xmm1
mulss  %xmm1,%xmm1
addss  %xmm4,%xmm1
subss  %xmm0,%xmm2
mulss  %xmm2,%xmm2
addss  %xmm1,%xmm2
movss  -0x4(%rdx,%rcx,4),%xmm1
subss  %xmm0,%xmm1
mulss  %xmm1,%xmm1
addss  %xmm2,%xmm1
movss  (%rdx,%rcx,4),%xmm4
subss  %xmm0,%xmm4
mulss  %xmm4,%xmm4
addss  %xmm1,%xmm4
add    $0x4,%rcx
cmp    %rcx,%rsi
jne    14e0 <equalization_affine+0x3a0>
test   %rbx,%rbx
je     156a <equalization_affine+0x42a>
add    %rax,%rcx
mov    0x48(%rsp),%rax
lea    (%rax,%rcx,4),%rax
xor    %ecx,%ecx
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movss  (%rax,%rcx,4),%xmm1
subss  %xmm0,%xmm1
mulss  %xmm1,%xmm1
addss  %xmm1,%xmm4
add    $0x1,%rcx
cmp    %rcx,%rbx
jne    1550 <equalization_affine+0x410>
movslq %ebp,%rdx
divss  %xmm3,%xmm4
mov    0x70(%rsp),%rax
movss  %xmm4,(%rax,%r13,4)
test   %r12d,%r12d
mov    %r8d,0x44(%rsp)
movss  %xmm4,0x40(%rsp)
jle    15b0 <equalization_affine+0x470>
mov    0x20(%rsp),%rax
movss  (%rax,%r13,4),%xmm3
cmpq   $0x3,0x88(%rsp)
jae    15c0 <equalization_affine+0x480>
xorps  %xmm0,%xmm0
xor    %eax,%eax
jmp    1620 <equalization_affine+0x4e0>
nopl   0x0(%rax)
xorps  %xmm0,%xmm0
jmp    165a <equalization_affine+0x51a>
nopl   0x0(%rax,%rax,1)
mov    0x50(%rsp),%rax
lea    (%rax,%rdx,4),%rcx
xorps  %xmm0,%xmm0
xor    %eax,%eax
xchg   %ax,%ax
movss  -0xc(%rcx,%rax,4),%xmm1
movss  -0x8(%rcx,%rax,4),%xmm2
subss  %xmm3,%xmm1
mulss  %xmm1,%xmm1
addss  %xmm0,%xmm1
subss  %xmm3,%xmm2
mulss  %xmm2,%xmm2
addss  %xmm1,%xmm2
movss  -0x4(%rcx,%rax,4),%xmm1
subss  %xmm3,%xmm1
mulss  %xmm1,%xmm1
addss  %xmm2,%xmm1
movss  (%rcx,%rax,4),%xmm0
subss  %xmm3,%xmm0
mulss  %xmm0,%xmm0
addss  %xmm1,%xmm0
add    $0x4,%rax
cmp    %rax,%r15
jne    15d0 <equalization_affine+0x490>
test   %r14,%r14
je     165a <equalization_affine+0x51a>
add    %rdx,%rax
mov    0x10(%rsp),%rcx
lea    (%rcx,%rax,4),%rax
xor    %ecx,%ecx
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movss  (%rax,%rcx,4),%xmm1
subss  %xmm3,%xmm1
mulss  %xmm1,%xmm1
addss  %xmm1,%xmm0
add    $0x1,%rcx
cmp    %rcx,%r14
jne    1640 <equalization_affine+0x500>
mov    %rdx,0x80(%rsp)
divss  %xmm8,%xmm0
movss  %xmm0,0x3c(%rsp)
mov    0x68(%rsp),%rax
movss  %xmm0,(%rax,%r13,4)
call   1040 <fabs@plt>
movss  0x97b(%rip),%xmm0
test   %eax,%eax
mov    0x78(%rsp),%r11
mov    0x44(%rsp),%r8d
jle    16ce <equalization_affine+0x58e>
movss  0x40(%rsp),%xmm0
divss  0x3c(%rsp),%xmm0
ucomiss 0x95e(%rip),%xmm0
jb     16b0 <equalization_affine+0x570>
sqrtss %xmm0,%xmm0
jmp    16ce <equalization_affine+0x58e>
xchg   %ax,%ax
mov    %r15,%r12
mov    %r8d,%r15d
call   1060 <sqrt@plt>
mov    %r15d,%r8d
mov    %r12,%r15
mov    0xc8(%rsp),%r12
mov    0x78(%rsp),%r11
mov    0x60(%rsp),%rax
movss  %xmm0,(%rax,%r13,4)
test   %r12d,%r12d
movss  0x34(%rsp),%xmm8
jle    1470 <equalization_affine+0x330>
cmp    $0x8,%r12d
setb   %cl
mov    %r12d,%eax
imul   %r13d,%eax
mov    0x20(%rsp),%rdx
movss  (%rdx,%r13,4),%xmm1
mov    0x18(%rsp),%rdx
movss  (%rdx,%r13,4),%xmm2
cltq
or     0xf(%rsp),%cl
je     1720 <equalization_affine+0x5e0>
xor    %edx,%edx
jmp    1813 <equalization_affine+0x6d3>
nopl   0x0(%rax)
movaps %xmm1,%xmm5
shufps $0x0,%xmm1,%xmm5
movaps %xmm0,%xmm3
shufps $0x0,%xmm0,%xmm3
movaps %xmm2,%xmm4
shufps $0x0,%xmm2,%xmm4
cmpq   $0x0,0xd8(%rsp)
je     18bc <equalization_affine+0x77c>
mov    0xa8(%rsp),%rcx
mov    0x80(%rsp),%rsi
lea    (%rcx,%rsi,4),%rdx
mov    0xa0(%rsp),%rcx
lea    (%rcx,%rsi,4),%rsi
mov    0xb0(%rsp),%rdi
xor    %ecx,%ecx
xchg   %ax,%ax
movups -0x30(%rsi,%rcx,4),%xmm6
movups -0x20(%rsi,%rcx,4),%xmm7
subps  %xmm5,%xmm6
subps  %xmm5,%xmm7
mulps  %xmm3,%xmm6
addps  %xmm4,%xmm6
mulps  %xmm3,%xmm7
addps  %xmm4,%xmm7
movups %xmm6,-0x30(%rdx,%rcx,4)
movups %xmm7,-0x20(%rdx,%rcx,4)
movups -0x10(%rsi,%rcx,4),%xmm6
movups (%rsi,%rcx,4),%xmm7
subps  %xmm5,%xmm6
subps  %xmm5,%xmm7
mulps  %xmm3,%xmm6
addps  %xmm4,%xmm6
mulps  %xmm3,%xmm7
addps  %xmm4,%xmm7
movups %xmm6,-0x10(%rdx,%rcx,4)
movups %xmm7,(%rdx,%rcx,4)
add    $0x10,%rcx
add    $0xfffffffffffffffe,%rdi
jne    1770 <equalization_affine+0x630>
testb  $0x1,0xd0(%rsp)
01
je     17ff <equalization_affine+0x6bf>
add    %rax,%rcx
mov    0x10(%rsp),%rdx
movups (%rdx,%rcx,4),%xmm6
movups 0x10(%rdx,%rcx,4),%xmm7
subps  %xmm5,%xmm6
subps  %xmm5,%xmm7
mulps  %xmm3,%xmm6
addps  %xmm4,%xmm6
mulps  %xmm3,%xmm7
addps  %xmm4,%xmm7
mov    0x28(%rsp),%rdx
movups %xmm6,(%rdx,%rcx,4)
movups %xmm7,0x10(%rdx,%rcx,4)
mov    0xe0(%rsp),%rcx
mov    %rcx,%rdx
cmp    %r12,%rcx
je     1470 <equalization_affine+0x330>
mov    %rdx,%rcx
test   $0x1,%r12b
je     1846 <equalization_affine+0x706>
add    %rdx,%rax
mov    0x10(%rsp),%rcx
movss  (%rcx,%rax,4),%xmm3
subss  %xmm1,%xmm3
mulss  %xmm0,%xmm3
addss  %xmm2,%xmm3
mov    0x28(%rsp),%rcx
movss  %xmm3,(%rcx,%rax,4)
mov    %rdx,%rcx
or     $0x1,%rcx
not    %rdx
cmp    0xe8(%rsp),%rdx
je     1470 <equalization_affine+0x330>
mov    0xc0(%rsp),%rax
mov    0x80(%rsp),%rsi
lea    (%rax,%rsi,4),%rax
mov    0xb8(%rsp),%rdx
lea    (%rdx,%rsi,4),%rdx
nopw   0x0(%rax,%rax,1)
movss  -0x4(%rdx,%rcx,4),%xmm3
subss  %xmm1,%xmm3
mulss  %xmm0,%xmm3
addss  %xmm2,%xmm3
movss  %xmm3,-0x4(%rax,%rcx,4)
movss  (%rdx,%rcx,4),%xmm3
subss  %xmm1,%xmm3
mulss  %xmm0,%xmm3
addss  %xmm2,%xmm3
movss  %xmm3,(%rax,%rcx,4)
add    $0x2,%rcx
cmp    %rcx,%r12
jne    1880 <equalization_affine+0x740>
jmp    1470 <equalization_affine+0x330>
xor    %ecx,%ecx
testb  $0x1,0xd0(%rsp)
01
jne    17ce <equalization_affine+0x68e>
jmp    17ff <equalization_affine+0x6bf>
mov    0x18(%rsp),%rdi
call   1030 <free@plt>
mov    0x20(%rsp),%rdi
call   1030 <free@plt>
mov    0x70(%rsp),%rdi
call   1030 <free@plt>
mov    0x68(%rsp),%rdi
call   1030 <free@plt>
mov    0x60(%rsp),%rdi
add    $0xf8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
jmp    1030 <free@plt>

LLVM IR (-O0):
define void @equalization_affine(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 float*
%20 = alloca float*
%21 = alloca float*
%22 = alloca i32
%23 = alloca i32
%24 = alloca i32
%25 = alloca float
%26 = alloca i32
%27 = alloca i32
%28 = alloca i32
%29 = 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
%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** %17
%35 = load i32, i32* %16
%36 = sext i32 %35 to i64
%37 = mul i64 %36, 4
%38 = trunc i64 %37 to i32
%39 = call float* @malloc(i32 %38)
store float* %39, float** %18
%40 = load i32, i32* %16
%41 = sext i32 %40 to i64
%42 = mul i64 %41, 4
%43 = trunc i64 %42 to i32
%44 = call float* @malloc(i32 %43)
store float* %44, float** %19
%45 = load i32, i32* %16
%46 = sext i32 %45 to i64
%47 = mul i64 %46, 4
%48 = trunc i64 %47 to i32
%49 = call float* @malloc(i32 %48)
store float* %49, float** %20
%50 = load i32, i32* %16
%51 = sext i32 %50 to i64
%52 = mul i64 %51, 4
%53 = trunc i64 %52 to i32
%54 = call float* @malloc(i32 %53)
store float* %54, float** %21
store i32 0, i32* %22
br label %55
55: ; preds = %144, %8
%56 = load i32, i32* %22
%57 = load i32, i32* %16
%58 = icmp slt i32 %56, %57
br i1 %58, label %59, label %147
59: ; preds = %55
%60 = load float*, float** %17
%61 = load i32, i32* %22
%62 = sext i32 %61 to i64
%63 = getelementptr inbounds float, float* %60, i64 %62
store float 0.000000e+00, float* %63
store i32 0, i32* %23
br label %64
64: ; preds = %88, %59
%65 = load i32, i32* %23
%66 = load i32, i32* %12
%67 = load i32, i32* %13
%68 = mul nsw i32 %66, %67
%69 = icmp slt i32 %65, %68
br i1 %69, label %70, label %91
70: ; preds = %64
%71 = load float*, float** %9
%72 = load i32, i32* %23
%73 = load i32, i32* %22
%74 = load i32, i32* %12
%75 = mul nsw i32 %73, %74
%76 = load i32, i32* %13
%77 = mul nsw i32 %75, %76
%78 = add nsw i32 %72, %77
%79 = sext i32 %78 to i64
%80 = getelementptr inbounds float, float* %71, i64 %79
%81 = load float, float* %80
%82 = load float*, float** %17
%83 = load i32, i32* %22
%84 = sext i32 %83 to i64
%85 = getelementptr inbounds float, float* %82, i64 %84
%86 = load float, float* %85
%87 = fadd float %86, %81
store float %87, float* %85
br label %88
88: ; preds = %70
%89 = load i32, i32* %23
%90 = add nsw i32 %89, 1
store i32 %90, i32* %23
br label %64
91: ; preds = %64
%92 = load i32, i32* %12
%93 = load i32, i32* %13
%94 = mul nsw i32 %92, %93
%95 = sitofp i32 %94 to float
%96 = load float*, float** %17
%97 = load i32, i32* %22
%98 = sext i32 %97 to i64
%99 = getelementptr inbounds float, float* %96, i64 %98
%100 = load float, float* %99
%101 = fdiv float %100, %95
store float %101, float* %99
%102 = load float*, float** %18
%103 = load i32, i32* %22
%104 = sext i32 %103 to i64
%105 = getelementptr inbounds float, float* %102, i64 %104
store float 0.000000e+00, float* %105
store i32 0, i32* %24
br label %106
106: ; preds = %130, %91
%107 = load i32, i32* %24
%108 = load i32, i32* %14
%109 = load i32, i32* %15
%110 = mul nsw i32 %108, %109
%111 = icmp slt i32 %107, %110
br i1 %111, label %112, label %133
112: ; preds = %106
%113 = load float*, float** %10
%114 = load i32, i32* %24
%115 = load i32, i32* %22
%116 = load i32, i32* %14
%117 = mul nsw i32 %115, %116
%118 = load i32, i32* %15
%119 = mul nsw i32 %117, %118
%120 = add nsw i32 %114, %119
%121 = sext i32 %120 to i64
%122 = getelementptr inbounds float, float* %113, i64 %121
%123 = load float, float* %122
%124 = load float*, float** %18
%125 = load i32, i32* %22
%126 = sext i32 %125 to i64
%127 = getelementptr inbounds float, float* %124, i64 %126
%128 = load float, float* %127
%129 = fadd float %128, %123
store float %129, float* %127
br label %130
130: ; preds = %112
%131 = load i32, i32* %24
%132 = add nsw i32 %131, 1
store i32 %132, i32* %24
br label %106
133: ; preds = %106
%134 = load i32, i32* %14
%135 = load i32, i32* %15
%136 = mul nsw i32 %134, %135
%137 = sitofp i32 %136 to float
%138 = load float*, float** %18
%139 = load i32, i32* %22
%140 = sext i32 %139 to i64
%141 = getelementptr inbounds float, float* %138, i64 %140
%142 = load float, float* %141
%143 = fdiv float %142, %137
store float %143, float* %141
br label %144
144: ; preds = %133
%145 = load i32, i32* %22
%146 = add nsw i32 %145, 1
store i32 %146, i32* %22
br label %55
147: ; preds = %55
store i32 0, i32* %26
br label %148
148: ; preds = %333, %147
%149 = load i32, i32* %26
%150 = load i32, i32* %16
%151 = icmp slt i32 %149, %150
br i1 %151, label %152, label %336
152: ; preds = %148
%153 = load float*, float** %19
%154 = load i32, i32* %26
%155 = sext i32 %154 to i64
%156 = getelementptr inbounds float, float* %153, i64 %155
store float 0.000000e+00, float* %156
store i32 0, i32* %27
br label %157
157: ; preds = %189, %152
%158 = load i32, i32* %27
%159 = load i32, i32* %12
%160 = load i32, i32* %13
%161 = mul nsw i32 %159, %160
%162 = icmp slt i32 %158, %161
br i1 %162, label %163, label %192
163: ; preds = %157
%164 = load float*, float** %9
%165 = load i32, i32* %27
%166 = load i32, i32* %26
%167 = load i32, i32* %12
%168 = mul nsw i32 %166, %167
%169 = load i32, i32* %13
%170 = mul nsw i32 %168, %169
%171 = add nsw i32 %165, %170
%172 = sext i32 %171 to i64
%173 = getelementptr inbounds float, float* %164, i64 %172
%174 = load float, float* %173
%175 = load float*, float** %17
%176 = load i32, i32* %26
%177 = sext i32 %176 to i64
%178 = getelementptr inbounds float, float* %175, i64 %177
%179 = load float, float* %178
%180 = fsub float %174, %179
store float %180, float* %25
%181 = load float, float* %25
%182 = load float, float* %25
%183 = load float*, float** %19
%184 = load i32, i32* %26
%185 = sext i32 %184 to i64
%186 = getelementptr inbounds float, float* %183, i64 %185
%187 = load float, float* %186
%188 = call float @llvm.fmuladd.f32(float %181, float %182, float %187)
store float %188, float* %186
br label %189
189: ; preds = %163
%190 = load i32, i32* %27
%191 = add nsw i32 %190, 1
store i32 %191, i32* %27
br label %157
192: ; preds = %157
%193 = load i32, i32* %12
%194 = load i32, i32* %13
%195 = mul nsw i32 %193, %194
%196 = sitofp i32 %195 to float
%197 = load float*, float** %19
%198 = load i32, i32* %26
%199 = sext i32 %198 to i64
%200 = getelementptr inbounds float, float* %197, i64 %199
%201 = load float, float* %200
%202 = fdiv float %201, %196
store float %202, float* %200
%203 = load float*, float** %20
%204 = load i32, i32* %26
%205 = sext i32 %204 to i64
%206 = getelementptr inbounds float, float* %203, i64 %205
store float 0.000000e+00, float* %206
store i32 0, i32* %28
br label %207
207: ; preds = %239, %192
%208 = load i32, i32* %28
%209 = load i32, i32* %14
%210 = load i32, i32* %15
%211 = mul nsw i32 %209, %210
%212 = icmp slt i32 %208, %211
br i1 %212, label %213, label %242
213: ; preds = %207
%214 = load float*, float** %10
%215 = load i32, i32* %28
%216 = load i32, i32* %26
%217 = load i32, i32* %14
%218 = mul nsw i32 %216, %217
%219 = load i32, i32* %15
%220 = mul nsw i32 %218, %219
%221 = add nsw i32 %215, %220
%222 = sext i32 %221 to i64
%223 = getelementptr inbounds float, float* %214, i64 %222
%224 = load float, float* %223
%225 = load float*, float** %18
%226 = load i32, i32* %26
%227 = sext i32 %226 to i64
%228 = getelementptr inbounds float, float* %225, i64 %227
%229 = load float, float* %228
%230 = fsub float %224, %229
store float %230, float* %25
%231 = load float, float* %25
%232 = load float, float* %25
%233 = load float*, float** %20
%234 = load i32, i32* %26
%235 = sext i32 %234 to i64
%236 = getelementptr inbounds float, float* %233, i64 %235
%237 = load float, float* %236
%238 = call float @llvm.fmuladd.f32(float %231, float %232, float %237)
store float %238, float* %236
br label %239
239: ; preds = %213
%240 = load i32, i32* %28
%241 = add nsw i32 %240, 1
store i32 %241, i32* %28
br label %207
242: ; preds = %207
%243 = load i32, i32* %14
%244 = load i32, i32* %15
%245 = mul nsw i32 %243, %244
%246 = sitofp i32 %245 to float
%247 = load float*, float** %20
%248 = load i32, i32* %26
%249 = sext i32 %248 to i64
%250 = getelementptr inbounds float, float* %247, i64 %249
%251 = load float, float* %250
%252 = fdiv float %251, %246
store float %252, float* %250
%253 = load float*, float** %20
%254 = load i32, i32* %26
%255 = sext i32 %254 to i64
%256 = getelementptr inbounds float, float* %253, i64 %255
%257 = load float, float* %256
%258 = call i32 @fabs(float %257)
%259 = sitofp i32 %258 to double
%260 = fcmp olt double %259, 1.000000e-03
br i1 %260, label %261, label %266
261: ; preds = %242
%262 = load float*, float** %21
%263 = load i32, i32* %26
%264 = sext i32 %263 to i64
%265 = getelementptr inbounds float, float* %262, i64 %264
store float 1.000000e+00, float* %265
br label %283
266: ; preds = %242
%267 = load float*, float** %19
%268 = load i32, i32* %26
%269 = sext i32 %268 to i64
%270 = getelementptr inbounds float, float* %267, i64 %269
%271 = load float, float* %270
%272 = load float*, float** %20
%273 = load i32, i32* %26
%274 = sext i32 %273 to i64
%275 = getelementptr inbounds float, float* %272, i64 %274
%276 = load float, float* %275
%277 = fdiv float %271, %276
%278 = call float @sqrt(float %277)
%279 = load float*, float** %21
%280 = load i32, i32* %26
%281 = sext i32 %280 to i64
%282 = getelementptr inbounds float, float* %279, i64 %281
store float %278, float* %282
br label %283
283: ; preds = %266, %261
store i32 0, i32* %29
br label %284
284: ; preds = %329, %283
%285 = load i32, i32* %29
%286 = load i32, i32* %14
%287 = load i32, i32* %15
%288 = mul nsw i32 %286, %287
%289 = icmp slt i32 %285, %288
br i1 %289, label %290, label %332
290: ; preds = %284
%291 = load float*, float** %21
%292 = load i32, i32* %26
%293 = sext i32 %292 to i64
%294 = getelementptr inbounds float, float* %291, i64 %293
%295 = load float, float* %294
%296 = load float*, float** %10
%297 = load i32, i32* %29
%298 = load i32, i32* %26
%299 = load i32, i32* %15
%300 = mul nsw i32 %298, %299
%301 = load i32, i32* %14
%302 = mul nsw i32 %300, %301
%303 = add nsw i32 %297, %302
%304 = sext i32 %303 to i64
%305 = getelementptr inbounds float, float* %296, i64 %304
%306 = load float, float* %305
%307 = load float*, float** %18
%308 = load i32, i32* %26
%309 = sext i32 %308 to i64
%310 = getelementptr inbounds float, float* %307, i64 %309
%311 = load float, float* %310
%312 = fsub float %306, %311
%313 = load float*, float** %17
%314 = load i32, i32* %26
%315 = sext i32 %314 to i64
%316 = getelementptr inbounds float, float* %313, i64 %315
%317 = load float, float* %316
%318 = call float @llvm.fmuladd.f32(float %295, float %312, float %317)
%319 = load float*, float** %11
%320 = load i32, i32* %29
%321 = load i32, i32* %26
%322 = load i32, i32* %14
%323 = mul nsw i32 %321, %322
%324 = load i32, i32* %15
%325 = mul nsw i32 %323, %324
%326 = add nsw i32 %320, %325
%327 = sext i32 %326 to i64
%328 = getelementptr inbounds float, float* %319, i64 %327
store float %318, float* %328
br label %329
329: ; preds = %290
%330 = load i32, i32* %29
%331 = add nsw i32 %330, 1
store i32 %331, i32* %29
br label %284
332: ; preds = %284
br label %333
333: ; preds = %332
%334 = load i32, i32* %26
%335 = add nsw i32 %334, 1
store i32 %335, i32* %26
br label %148
336: ; preds = %148
%337 = load float*, float** %17
%338 = call i32 @free(float* %337)
%339 = load float*, float** %18
%340 = call i32 @free(float* %339)
%341 = load float*, float** %19
%342 = call i32 @free(float* %341)
%343 = load float*, float** %20
%344 = call i32 @free(float* %343)
%345 = load float*, float** %21
%346 = call i32 @free(float* %345)
ret void
}

LLVM IR (-O2):
define void @equalization_affine(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 = tail call float* @malloc(i32 %9)
%13 = tail call float* @malloc(i32 %9)
%14 = tail call float* @malloc(i32 %9)
%15 = icmp sgt i32 %7, 0
br i1 %15, label %16, label %175
16: ; preds = %8
%17 = mul i32 %4, %3
%18 = icmp sgt i32 %17, 0
%19 = sitofp i32 %17 to float
%20 = mul i32 %6, %5
%21 = icmp sgt i32 %20, 0
%22 = sitofp i32 %20 to float
%23 = zext i32 %7 to i64
%24 = zext i32 %17 to i64
%25 = zext i32 %20 to i64
%26 = add nsw i64 %24, -1
%27 = add nsw i64 %25, -1
%28 = and i64 %24, 3
%29 = icmp ult i64 %26, 3
%30 = and i64 %24, 4294967292
%31 = icmp eq i64 %28, 0
%32 = and i64 %25, 3
%33 = icmp ult i64 %27, 3
%34 = and i64 %25, 4294967292
%35 = icmp eq i64 %32, 0
br label %73
36: ; preds = %144
br i1 %15, label %37, label %175
37: ; preds = %36
%38 = mul i32 %4, %3
%39 = icmp sgt i32 %38, 0
%40 = sitofp i32 %38 to float
%41 = mul i32 %6, %5
%42 = icmp sgt i32 %41, 0
%43 = sitofp i32 %41 to float
%44 = zext i32 %7 to i64
%45 = zext i32 %41 to i64
%46 = getelementptr float, float* %2, i64 %25
%47 = getelementptr float, float* %1, i64 %25
%48 = and i64 %25, 4294967288
%49 = add nsw i64 %48, -8
%50 = lshr exact i64 %49, 3
%51 = add nuw nsw i64 %50, 1
%52 = and i64 %24, 3
%53 = icmp ult i64 %26, 3
%54 = and i64 %24, 4294967292
%55 = icmp eq i64 %52, 0
%56 = and i64 %25, 3
%57 = icmp ult i64 %27, 3
%58 = and i64 %25, 4294967292
%59 = icmp eq i64 %56, 0
%60 = icmp ult i32 %20, 8
%61 = icmp ugt float* %47, %2
%62 = icmp ugt float* %46, %1
%63 = and i1 %61, %62
%64 = and i64 %25, 4294967288
%65 = and i64 %51, 1
%66 = icmp eq i64 %49, 0
%67 = and i64 %51, 4611686018427387902
%68 = icmp eq i64 %65, 0
%69 = icmp eq i64 %64, %25
%70 = and i64 %25, 1
%71 = icmp eq i64 %70, 0
%72 = sub nsw i64 0, %25
br label %181
73: ; preds = %16, %144
%74 = phi i64 [ 0, %16 ], [ %147, %144 ]
%75 = getelementptr inbounds float, float* %10, i64 %74
br i1 %18, label %76, label %95
76: ; preds = %73
%77 = trunc i64 %74 to i32
%78 = mul i32 %17, %77
%79 = sext i32 %78 to i64
br i1 %29, label %80, label %103
80: ; preds = %103, %76
%81 = phi float [ undef, %76 ], [ %125, %103 ]
%82 = phi i64 [ 0, %76 ], [ %126, %103 ]
%83 = phi float [ 0.000000e+00, %76 ], [ %125, %103 ]
br i1 %31, label %95, label %84
84: ; preds = %80, %84
%85 = phi i64 [ %92, %84 ], [ %82, %80 ]
%86 = phi float [ %91, %84 ], [ %83, %80 ]
%87 = phi i64 [ %93, %84 ], [ 0, %80 ]
%88 = add nsw i64 %85, %79
%89 = getelementptr inbounds float, float* %0, i64 %88
%90 = load float, float* %89
%91 = fadd float %90, %86
%92 = add nuw nsw i64 %85, 1
%93 = add i64 %87, 1
%94 = icmp eq i64 %93, %28
br i1 %94, label %95, label %84
95: ; preds = %80, %84, %73
%96 = phi float [ 0.000000e+00, %73 ], [ %81, %80 ], [ %91, %84 ]
%97 = fdiv float %96, %19
store float %97, float* %75
%98 = getelementptr inbounds float, float* %11, i64 %74
br i1 %21, label %99, label %144
99: ; preds = %95
%100 = trunc i64 %74 to i32
%101 = mul i32 %20, %100
%102 = sext i32 %101 to i64
br i1 %33, label %129, label %149
103: ; preds = %76, %103
%104 = phi i64 [ %126, %103 ], [ 0, %76 ]
%105 = phi float [ %125, %103 ], [ 0.000000e+00, %76 ]
%106 = phi i64 [ %127, %103 ], [ 0, %76 ]
%107 = add nsw i64 %104, %79
%108 = getelementptr inbounds float, float* %0, i64 %107
%109 = load float, float* %108
%110 = fadd float %109, %105
%111 = or i64 %104, 1
%112 = add nsw i64 %111, %79
%113 = getelementptr inbounds float, float* %0, i64 %112
%114 = load float, float* %113
%115 = fadd float %114, %110
%116 = or i64 %104, 2
%117 = add nsw i64 %116, %79
%118 = getelementptr inbounds float, float* %0, i64 %117
%119 = load float, float* %118
%120 = fadd float %119, %115
%121 = or i64 %104, 3
%122 = add nsw i64 %121, %79
%123 = getelementptr inbounds float, float* %0, i64 %122
%124 = load float, float* %123
%125 = fadd float %124, %120
%126 = add nuw nsw i64 %104, 4
%127 = add i64 %106, 4
%128 = icmp eq i64 %127, %30
br i1 %128, label %80, label %103
129: ; preds = %149, %99
%130 = phi float [ undef, %99 ], [ %171, %149 ]
%131 = phi i64 [ 0, %99 ], [ %172, %149 ]
%132 = phi float [ 0.000000e+00, %99 ], [ %171, %149 ]
br i1 %35, label %144, label %133
133: ; preds = %129, %133
%134 = phi i64 [ %141, %133 ], [ %131, %129 ]
%135 = phi float [ %140, %133 ], [ %132, %129 ]
%136 = phi i64 [ %142, %133 ], [ 0, %129 ]
%137 = add nsw i64 %134, %102
%138 = getelementptr inbounds float, float* %1, i64 %137
%139 = load float, float* %138
%140 = fadd float %139, %135
%141 = add nuw nsw i64 %134, 1
%142 = add i64 %136, 1
%143 = icmp eq i64 %142, %32
br i1 %143, label %144, label %133
144: ; preds = %129, %133, %95
%145 = phi float [ 0.000000e+00, %95 ], [ %130, %129 ], [ %140, %133 ]
%146 = fdiv float %145, %22
store float %146, float* %98
%147 = add nuw nsw i64 %74, 1
%148 = icmp eq i64 %147, %23
br i1 %148, label %36, label %73
149: ; preds = %99, %149
%150 = phi i64 [ %172, %149 ], [ 0, %99 ]
%151 = phi float [ %171, %149 ], [ 0.000000e+00, %99 ]
%152 = phi i64 [ %173, %149 ], [ 0, %99 ]
%153 = add nsw i64 %150, %102
%154 = getelementptr inbounds float, float* %1, i64 %153
%155 = load float, float* %154
%156 = fadd float %155, %151
%157 = or i64 %150, 1
%158 = add nsw i64 %157, %102
%159 = getelementptr inbounds float, float* %1, i64 %158
%160 = load float, float* %159
%161 = fadd float %160, %156
%162 = or i64 %150, 2
%163 = add nsw i64 %162, %102
%164 = getelementptr inbounds float, float* %1, i64 %163
%165 = load float, float* %164
%166 = fadd float %165, %161
%167 = or i64 %150, 3
%168 = add nsw i64 %167, %102
%169 = getelementptr inbounds float, float* %1, i64 %168
%170 = load float, float* %169
%171 = fadd float %170, %166
%172 = add nuw nsw i64 %150, 4
%173 = add i64 %152, 4
%174 = icmp eq i64 %173, %34
br i1 %174, label %129, label %149
175: ; preds = %395, %8, %36
%176 = tail call i32 @free(float* %10)
%177 = tail call i32 @free(float* %11)
%178 = tail call i32 @free(float* %12)
%179 = tail call i32 @free(float* %13)
%180 = tail call i32 @free(float* %14)
ret void
181: ; preds = %37, %395
%182 = phi i64 [ 0, %37 ], [ %396, %395 ]
%183 = getelementptr inbounds float, float* %12, i64 %182
br i1 %39, label %184, label %206
184: ; preds = %181
%185 = trunc i64 %182 to i32
%186 = mul i32 %38, %185
%187 = getelementptr inbounds float, float* %10, i64 %182
%188 = load float, float* %187
%189 = sext i32 %186 to i64
br i1 %53, label %190, label %216
190: ; preds = %216, %184
%191 = phi float [ undef, %184 ], [ %242, %216 ]
%192 = phi i64 [ 0, %184 ], [ %243, %216 ]
%193 = phi float [ 0.000000e+00, %184 ], [ %242, %216 ]
br i1 %55, label %206, label %194
194: ; preds = %190, %194
%195 = phi i64 [ %203, %194 ], [ %192, %190 ]
%196 = phi float [ %202, %194 ], [ %193, %190 ]
%197 = phi i64 [ %204, %194 ], [ 0, %190 ]
%198 = add nsw i64 %195, %189
%199 = getelementptr inbounds float, float* %0, i64 %198
%200 = load float, float* %199
%201 = fsub float %200, %188
%202 = tail call float @llvm.fmuladd.f32(float %201, float %201, float %196)
%203 = add nuw nsw i64 %195, 1
%204 = add i64 %197, 1
%205 = icmp eq i64 %204, %52
br i1 %205, label %206, label %194
206: ; preds = %190, %194, %181
%207 = phi float [ 0.000000e+00, %181 ], [ %191, %190 ], [ %202, %194 ]
%208 = fdiv float %207, %40
store float %208, float* %183
%209 = getelementptr inbounds float, float* %13, i64 %182
br i1 %42, label %210, label %262
210: ; preds = %206
%211 = trunc i64 %182 to i32
%212 = mul i32 %41, %211
%213 = getelementptr inbounds float, float* %11, i64 %182
%214 = load float, float* %213
%215 = sext i32 %212 to i64
br i1 %57, label %246, label %267
216: ; preds = %184, %216
%217 = phi i64 [ %243, %216 ], [ 0, %184 ]
%218 = phi float [ %242, %216 ], [ 0.000000e+00, %184 ]
%219 = phi i64 [ %244, %216 ], [ 0, %184 ]
%220 = add nsw i64 %217, %189
%221 = getelementptr inbounds float, float* %0, i64 %220
%222 = load float, float* %221
%223 = fsub float %222, %188
%224 = tail call float @llvm.fmuladd.f32(float %223, float %223, float %218)
%225 = or i64 %217, 1
%226 = add nsw i64 %225, %189
%227 = getelementptr inbounds float, float* %0, i64 %226
%228 = load float, float* %227
%229 = fsub float %228, %188
%230 = tail call float @llvm.fmuladd.f32(float %229, float %229, float %224)
%231 = or i64 %217, 2
%232 = add nsw i64 %231, %189
%233 = getelementptr inbounds float, float* %0, i64 %232
%234 = load float, float* %233
%235 = fsub float %234, %188
%236 = tail call float @llvm.fmuladd.f32(float %235, float %235, float %230)
%237 = or i64 %217, 3
%238 = add nsw i64 %237, %189
%239 = getelementptr inbounds float, float* %0, i64 %238
%240 = load float, float* %239
%241 = fsub float %240, %188
%242 = tail call float @llvm.fmuladd.f32(float %241, float %241, float %236)
%243 = add nuw nsw i64 %217, 4
%244 = add i64 %219, 4
%245 = icmp eq i64 %244, %54
br i1 %245, label %190, label %216
246: ; preds = %267, %210
%247 = phi float [ undef, %210 ], [ %293, %267 ]
%248 = phi i64 [ 0, %210 ], [ %294, %267 ]
%249 = phi float [ 0.000000e+00, %210 ], [ %293, %267 ]
br i1 %59, label %262, label %250
250: ; preds = %246, %250
%251 = phi i64 [ %259, %250 ], [ %248, %246 ]
%252 = phi float [ %258, %250 ], [ %249, %246 ]
%253 = phi i64 [ %260, %250 ], [ 0, %246 ]
%254 = add nsw i64 %251, %215
%255 = getelementptr inbounds float, float* %1, i64 %254
%256 = load float, float* %255
%257 = fsub float %256, %214
%258 = tail call float @llvm.fmuladd.f32(float %257, float %257, float %252)
%259 = add nuw nsw i64 %251, 1
%260 = add i64 %253, 1
%261 = icmp eq i64 %260, %56
br i1 %261, label %262, label %250
262: ; preds = %246, %250, %206
%263 = phi float [ 0.000000e+00, %206 ], [ %247, %246 ], [ %258, %250 ]
%264 = fdiv float %263, %43
store float %264, float* %209
%265 = tail call i32 @fabs(float %264)
%266 = icmp slt i32 %265, 1
br i1 %266, label %300, label %297
267: ; preds = %210, %267
%268 = phi i64 [ %294, %267 ], [ 0, %210 ]
%269 = phi float [ %293, %267 ], [ 0.000000e+00, %210 ]
%270 = phi i64 [ %295, %267 ], [ 0, %210 ]
%271 = add nsw i64 %268, %215
%272 = getelementptr inbounds float, float* %1, i64 %271
%273 = load float, float* %272
%274 = fsub float %273, %214
%275 = tail call float @llvm.fmuladd.f32(float %274, float %274, float %269)
%276 = or i64 %268, 1
%277 = add nsw i64 %276, %215
%278 = getelementptr inbounds float, float* %1, i64 %277
%279 = load float, float* %278
%280 = fsub float %279, %214
%281 = tail call float @llvm.fmuladd.f32(float %280, float %280, float %275)
%282 = or i64 %268, 2
%283 = add nsw i64 %282, %215
%284 = getelementptr inbounds float, float* %1, i64 %283
%285 = load float, float* %284
%286 = fsub float %285, %214
%287 = tail call float @llvm.fmuladd.f32(float %286, float %286, float %281)
%288 = or i64 %268, 3
%289 = add nsw i64 %288, %215
%290 = getelementptr inbounds float, float* %1, i64 %289
%291 = load float, float* %290
%292 = fsub float %291, %214
%293 = tail call float @llvm.fmuladd.f32(float %292, float %292, float %287)
%294 = add nuw nsw i64 %268, 4
%295 = add i64 %270, 4
%296 = icmp eq i64 %295, %58
br i1 %296, label %246, label %267
297: ; preds = %262
%298 = fdiv float %208, %264
%299 = tail call float @sqrt(float %298)
br label %300
300: ; preds = %262, %297
%301 = phi float [ %299, %297 ], [ 1.000000e+00, %262 ]
%302 = getelementptr inbounds float, float* %14, i64 %182
store float %301, float* %302
br i1 %42, label %303, label %395
303: ; preds = %300
%304 = trunc i64 %182 to i32
%305 = mul i32 %41, %304
%306 = getelementptr inbounds float, float* %11, i64 %182
%307 = load float, float* %306
%308 = getelementptr inbounds float, float* %10, i64 %182
%309 = load float, float* %308
%310 = sext i32 %305 to i64
%311 = or i1 %60, %63
br i1 %311, label %381, label %312
312: ; preds = %303
%313 = insertelement <4 x float> poison, float %307, i64 0
%314 = shufflevector <4 x float> %313, <4 x float> poison, <4 x i32> zeroinitializer
%315 = insertelement <4 x float> poison, float %307, i64 0
%316 = shufflevector <4 x float> %315, <4 x float> poison, <4 x i32> zeroinitializer
%317 = insertelement <4 x float> poison, float %301, i64 0
%318 = shufflevector <4 x float> %317, <4 x float> poison, <4 x i32> zeroinitializer
%319 = insertelement <4 x float> poison, float %309, i64 0
%320 = shufflevector <4 x float> %319, <4 x float> poison, <4 x i32> zeroinitializer
%321 = insertelement <4 x float> poison, float %301, i64 0
%322 = shufflevector <4 x float> %321, <4 x float> poison, <4 x i32> zeroinitializer
%323 = insertelement <4 x float> poison, float %309, i64 0
%324 = shufflevector <4 x float> %323, <4 x float> poison, <4 x i32> zeroinitializer
br i1 %66, label %362, label %325
325: ; preds = %312, %325
%326 = phi i64 [ %359, %325 ], [ 0, %312 ]
%327 = phi i64 [ %360, %325 ], [ 0, %312 ]
%328 = add nsw i64 %326, %310
%329 = getelementptr inbounds float, float* %1, i64 %328
%330 = bitcast float* %329 to <4 x float>*
%331 = load <4 x float>, <4 x float>* %330
%332 = getelementptr inbounds float, float* %329, i64 4
%333 = bitcast float* %332 to <4 x float>*
%334 = load <4 x float>, <4 x float>* %333
%335 = fsub <4 x float> %331, %314
%336 = fsub <4 x float> %334, %316
%337 = call <4 x float> @llvm.fmuladd.v4f32(<4 x float> %318, <4 x float> %335, <4 x float> %320)
%338 = call <4 x float> @llvm.fmuladd.v4f32(<4 x float> %322, <4 x float> %336, <4 x float> %324)
%339 = getelementptr inbounds float, float* %2, i64 %328
%340 = bitcast float* %339 to <4 x float>*
store <4 x float> %337, <4 x float>* %340
%341 = getelementptr inbounds float, float* %339, i64 4
%342 = bitcast float* %341 to <4 x float>*
store <4 x float> %338, <4 x float>* %342
%343 = or i64 %326, 8
%344 = add nsw i64 %343, %310
%345 = getelementptr inbounds float, float* %1, i64 %344
%346 = bitcast float* %345 to <4 x float>*
%347 = load <4 x float>, <4 x float>* %346
%348 = getelementptr inbounds float, float* %345, i64 4
%349 = bitcast float* %348 to <4 x float>*
%350 = load <4 x float>, <4 x float>* %349
%351 = fsub <4 x float> %347, %314
%352 = fsub <4 x float> %350, %316
%353 = call <4 x float> @llvm.fmuladd.v4f32(<4 x float> %318, <4 x float> %351, <4 x float> %320)
%354 = call <4 x float> @llvm.fmuladd.v4f32(<4 x float> %322, <4 x float> %352, <4 x float> %324)
%355 = getelementptr inbounds float, float* %2, i64 %344
%356 = bitcast float* %355 to <4 x float>*
store <4 x float> %353, <4 x float>* %356
%357 = getelementptr inbounds float, float* %355, i64 4
%358 = bitcast float* %357 to <4 x float>*
store <4 x float> %354, <4 x float>* %358
%359 = add nuw i64 %326, 16
%360 = add i64 %327, 2
%361 = icmp eq i64 %360, %67
br i1 %361, label %362, label %325
362: ; preds = %325, %312
%363 = phi i64 [ 0, %312 ], [ %359, %325 ]
br i1 %68, label %380, label %364
364: ; preds = %362
%365 = add nsw i64 %363, %310
%366 = getelementptr inbounds float, float* %1, i64 %365
%367 = bitcast float* %366 to <4 x float>*
%368 = load <4 x float>, <4 x float>* %367
%369 = getelementptr inbounds float, float* %366, i64 4
%370 = bitcast float* %369 to <4 x float>*
%371 = load <4 x float>, <4 x float>* %370
%372 = fsub <4 x float> %368, %314
%373 = fsub <4 x float> %371, %316
%374 = call <4 x float> @llvm.fmuladd.v4f32(<4 x float> %318, <4 x float> %372, <4 x float> %320)
%375 = call <4 x float> @llvm.fmuladd.v4f32(<4 x float> %322, <4 x float> %373, <4 x float> %324)
%376 = getelementptr inbounds float, float* %2, i64 %365
%377 = bitcast float* %376 to <4 x float>*
store <4 x float> %374, <4 x float>* %377
%378 = getelementptr inbounds float, float* %376, i64 4
%379 = bitcast float* %378 to <4 x float>*
store <4 x float> %375, <4 x float>* %379
br label %380
380: ; preds = %362, %364
br i1 %69, label %395, label %381
381: ; preds = %303, %380
%382 = phi i64 [ 0, %303 ], [ %64, %380 ]
%383 = xor i64 %382, -1
br i1 %71, label %392, label %384
384: ; preds = %381
%385 = add nsw i64 %382, %310
%386 = getelementptr inbounds float, float* %1, i64 %385
%387 = load float, float* %386
%388 = fsub float %387, %307
%389 = tail call float @llvm.fmuladd.f32(float %301, float %388, float %309)
%390 = getelementptr inbounds float, float* %2, i64 %385
store float %389, float* %390
%391 = or i64 %382, 1
br label %392
392: ; preds = %384, %381
%393 = phi i64 [ %382, %381 ], [ %391, %384 ]
%394 = icmp eq i64 %383, %72
br i1 %394, label %395, label %398
395: ; preds = %392, %398, %380, %300
%396 = add nuw nsw i64 %182, 1
%397 = icmp eq i64 %396, %44
br i1 %397, label %175, label %181
398: ; preds = %392, %398
%399 = phi i64 [ %413, %398 ], [ %393, %392 ]
%400 = add nsw i64 %399, %310
%401 = getelementptr inbounds float, float* %1, i64 %400
%402 = load float, float* %401
%403 = fsub float %402, %307
%404 = tail call float @llvm.fmuladd.f32(float %301, float %403, float %309)
%405 = getelementptr inbounds float, float* %2, i64 %400
store float %404, float* %405
%406 = add nuw nsw i64 %399, 1
%407 = add nsw i64 %406, %310
%408 = getelementptr inbounds float, float* %1, i64 %407
%409 = load float, float* %408
%410 = fsub float %409, %307
%411 = tail call float @llvm.fmuladd.f32(float %301, float %410, float %309)
%412 = getelementptr inbounds float, float* %2, i64 %407
store float %411, float* %412
%413 = add nuw nsw i64 %399, 2
%414 = icmp eq i64 %413, %45
br i1 %414, label %395, label %398
}