#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 */
scalar_t__ calloc(long, int);
int /*<<< orphan*/ free(float *);
float *gamma_transform(double, double, float *, long);
float median_filt(float *);
float *refineTransmap(float *input, long int width, long int height) {
  float *h = gamma_transform(0.1, 1.30, input, width * height);
  long int f, d, u, y;
  int p = 0;
  float buffer[9];
  float *med_f = (float *)calloc(width * height, sizeof(float));
  for (f = 0; f < height; f++) {
    for (d = 0; d < width; d++) {
      p = 0;
      for (u = (f - 1); u <= (f + 1); u++) {
        for (y = (d - 1); y <= (d + 1); y++) {
          if ((u < 0) || (u >= height) || (y < 0) || (y >= width))
            buffer[p] = 0;
          else
            buffer[p] = *(h + (u * width) + y);
          p = p + 1;
        }
      }
      *(med_f + (f * width) + d) = median_filt(buffer);
    }
  }
  for (u = 0; u < (width * height); u++) {
    if (((*(input + u)) - (*(med_f + u))) > 1)
      (*(input + u)) = 1;
    else if (((*(input + u)) - (*(med_f + u))) > 0)
      (*(input + u)) = (*(input + u)) - (*(med_f + u));
    else if (((*(input + u)) - (*(med_f + u))) < 0)
      (*(input + u)) = 0;
  }
  free(h);
  free(med_f);
  return input;
}


x86 assembly (clang, -O0):
<refineTransmap>:
push   %rbp
mov    %rsp,%rbp
sub    $0x80,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    -0x8(%rbp),%rdi
mov    -0x10(%rbp),%rsi
imul   -0x18(%rbp),%rsi
movsd  0xe94(%rip),%xmm0
movsd  0xe94(%rip),%xmm1
call   1030 <gamma_transform@plt>
mov    %rax,-0x20(%rbp)
movl   $0x0,-0x44(%rbp)
mov    -0x10(%rbp),%rdi
imul   -0x18(%rbp),%rdi
mov    $0x4,%esi
call   1060 <calloc@plt>
mov    %rax,-0x78(%rbp)
movq   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%rax
cmp    -0x18(%rbp),%rax
jge    12f8 <refineTransmap+0x1b8>
movq   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%rax
cmp    -0x10(%rbp),%rax
jge    12e2 <refineTransmap+0x1a2>
movl   $0x0,-0x44(%rbp)
mov    -0x28(%rbp),%rax
sub    $0x1,%rax
mov    %rax,-0x38(%rbp)
mov    -0x38(%rbp),%rax
mov    -0x28(%rbp),%rcx
add    $0x1,%rcx
cmp    %rcx,%rax
jg     12ab <refineTransmap+0x16b>
mov    -0x30(%rbp),%rax
sub    $0x1,%rax
mov    %rax,-0x40(%rbp)
mov    -0x40(%rbp),%rax
mov    -0x30(%rbp),%rcx
add    $0x1,%rcx
cmp    %rcx,%rax
jg     1295 <refineTransmap+0x155>
cmpq   $0x0,-0x38(%rbp)
jl     1242 <refineTransmap+0x102>
mov    -0x38(%rbp),%rax
cmp    -0x18(%rbp),%rax
jge    1242 <refineTransmap+0x102>
cmpq   $0x0,-0x40(%rbp)
jl     1242 <refineTransmap+0x102>
mov    -0x40(%rbp),%rax
cmp    -0x10(%rbp),%rax
jl     1254 <refineTransmap+0x114>
movslq -0x44(%rbp),%rax
xorps  %xmm0,%xmm0
movss  %xmm0,-0x70(%rbp,%rax,4)
jmp    127b <refineTransmap+0x13b>
mov    -0x20(%rbp),%rax
mov    -0x38(%rbp),%rcx
imul   -0x10(%rbp),%rcx
shl    $0x2,%rcx
add    %rcx,%rax
mov    -0x40(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movslq -0x44(%rbp),%rax
movss  %xmm0,-0x70(%rbp,%rax,4)
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
mov    -0x40(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x40(%rbp)
jmp    11fb <refineTransmap+0xbb>
jmp    129a <refineTransmap+0x15a>
mov    -0x38(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
jmp    11da <refineTransmap+0x9a>
lea    -0x70(%rbp),%rdi
call   1050 <median_filt@plt>
mov    -0x78(%rbp),%rax
mov    -0x28(%rbp),%rcx
imul   -0x10(%rbp),%rcx
shl    $0x2,%rcx
add    %rcx,%rax
mov    -0x30(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x30(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x30(%rbp)
jmp    11b9 <refineTransmap+0x79>
jmp    12e7 <refineTransmap+0x1a7>
mov    -0x28(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x28(%rbp)
jmp    11a3 <refineTransmap+0x63>
movq   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%rax
mov    -0x10(%rbp),%rcx
imul   -0x18(%rbp),%rcx
cmp    %rcx,%rax
jge    1403 <refineTransmap+0x2c3>
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x78(%rbp),%rax
mov    -0x38(%rbp),%rcx
subss  (%rax,%rcx,4),%xmm0
movss  0xcd8(%rip),%xmm1
ucomiss %xmm1,%xmm0
jbe    135b <refineTransmap+0x21b>
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  0xcbf(%rip),%xmm0
movss  %xmm0,(%rax,%rcx,4)
jmp    13ed <refineTransmap+0x2ad>
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x78(%rbp),%rax
mov    -0x38(%rbp),%rcx
subss  (%rax,%rcx,4),%xmm0
xorps  %xmm1,%xmm1
ucomiss %xmm1,%xmm0
jbe    13ad <refineTransmap+0x26d>
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x78(%rbp),%rax
mov    -0x38(%rbp),%rcx
subss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
jmp    13e8 <refineTransmap+0x2a8>
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
mov    -0x78(%rbp),%rax
mov    -0x38(%rbp),%rcx
subss  (%rax,%rcx,4),%xmm1
xorps  %xmm0,%xmm0
ucomiss %xmm1,%xmm0
jbe    13e3 <refineTransmap+0x2a3>
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax,%rcx,4)
jmp    13e8 <refineTransmap+0x2a8>
jmp    13ed <refineTransmap+0x2ad>
jmp    13f2 <refineTransmap+0x2b2>
mov    -0x38(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
jmp    1300 <refineTransmap+0x1c0>
mov    -0x20(%rbp),%rdi
call   1040 <free@plt>
mov    -0x78(%rbp),%rdi
call   1040 <free@plt>
mov    -0x8(%rbp),%rax
add    $0x80,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<refineTransmap>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x68,%rsp
mov    %rdx,%rbx
mov    %rsi,%rbp
mov    %rdx,%r14
imul   %rsi,%r14
movsd  0xe9d(%rip),%xmm0
movsd  0xe9d(%rip),%xmm1
mov    %rdi,0x48(%rsp)
mov    %r14,%rsi
call   1030 <gamma_transform@plt>
mov    %rax,0x38(%rsp)
mov    %r14,0x40(%rsp)
mov    %r14,%rdi
mov    $0x4,%esi
call   1060 <calloc@plt>
mov    %rax,%r14
test   %rbx,%rbx
jle    1394 <refineTransmap+0x254>
xor    %eax,%eax
mov    %rax,(%rsp)
mov    %r14,0x50(%rsp)
jmp    11c5 <refineTransmap+0x85>
nopl   0x0(%rax,%rax,1)
mov    %r15,%rax
mov    %r15,(%rsp)
mov    0x50(%rsp),%r14
cmp    %rbx,%r15
je     1394 <refineTransmap+0x254>
test   %rbp,%rbp
jle    1380 <refineTransmap+0x240>
mov    (%rsp),%rsi
lea    -0x1(%rsi),%rax
mov    %rsi,%rcx
imul   %rbp,%rcx
lea    (%r14,%rcx,4),%rdx
mov    %rdx,0x60(%rsp)
imul   %rbp,%rax
mov    0x38(%rsp),%rdx
lea    (%rdx,%rax,4),%rdi
lea    (%rdx,%rcx,4),%r14
lea    0x1(%rsi),%r15
mov    %r15,%rax
imul   %rbp,%rax
lea    (%rdx,%rax,4),%rax
mov    %rax,0x8(%rsp)
xor    %r13d,%r13d
mov    %rdi,0x58(%rsp)
jmp    1256 <refineTransmap+0x116>
cs nopw 0x0(%rax,%rax,1)
nop
movl   $0x0,0x2c(%rsp)
movl   $0x0,0x30(%rsp)
lea    0x10(%rsp),%rdi
call   1050 <median_filt@plt>
mov    0x60(%rsp),%rax
movss  %xmm0,(%rax,%r13,4)
mov    %r12,%r13
cmp    %rbp,%r12
mov    0x58(%rsp),%rdi
je     11b0 <refineTransmap+0x70>
cmpq   $0x0,(%rsp)
jle    1270 <refineTransmap+0x130>
test   %r13,%r13
je     1280 <refineTransmap+0x140>
movss  -0x4(%rdi,%r13,4),%xmm0
jmp    1283 <refineTransmap+0x143>
nopl   0x0(%rax,%rax,1)
movq   $0x0,0x10(%rsp)
lea    0x1(%r13),%r12
jmp    12c0 <refineTransmap+0x180>
nop
xorps  %xmm0,%xmm0
movss  %xmm0,0x10(%rsp)
movss  (%rdi,%r13,4),%xmm0
movss  %xmm0,0x14(%rsp)
lea    0x1(%r13),%r12
cmp    %rbp,%r12
jge    12c0 <refineTransmap+0x180>
movss  0x4(%rdi,%r13,4),%xmm0
movss  %xmm0,0x18(%rsp)
test   %r13,%r13
jne    12cd <refineTransmap+0x18d>
movl   $0x0,0x1c(%rsp)
jmp    12da <refineTransmap+0x19a>
nopw   0x0(%rax,%rax,1)
movl   $0x0,0x18(%rsp)
test   %r13,%r13
je     12b0 <refineTransmap+0x170>
movss  -0x4(%r14,%r13,4),%xmm0
movss  %xmm0,0x1c(%rsp)
movss  (%r14,%r13,4),%xmm0
movss  %xmm0,0x20(%rsp)
cmp    %rbp,%r12
jge    1300 <refineTransmap+0x1c0>
movss  (%r14,%r12,4),%xmm0
movss  %xmm0,0x24(%rsp)
cmp    %rbx,%r15
jl     130d <refineTransmap+0x1cd>
jmp    1330 <refineTransmap+0x1f0>
xchg   %ax,%ax
movl   $0x0,0x24(%rsp)
cmp    %rbx,%r15
jge    1330 <refineTransmap+0x1f0>
test   %r13,%r13
je     1330 <refineTransmap+0x1f0>
mov    0x8(%rsp),%rax
movss  -0x4(%rax,%r13,4),%xmm0
movss  %xmm0,0x28(%rsp)
cmp    %rbx,%r15
jge    1220 <refineTransmap+0xe0>
jmp    1341 <refineTransmap+0x201>
nop
movl   $0x0,0x28(%rsp)
cmp    %rbx,%r15
jge    1220 <refineTransmap+0xe0>
mov    0x8(%rsp),%rax
movss  (%rax,%r13,4),%xmm0
movss  %xmm0,0x2c(%rsp)
cmp    %rbp,%r12
jge    1228 <refineTransmap+0xe8>
mov    0x8(%rsp),%rax
movss  (%rax,%r12,4),%xmm0
movss  %xmm0,0x30(%rsp)
jmp    1230 <refineTransmap+0xf0>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    (%rsp),%r15
add    $0x1,%r15
mov    %r15,%rax
mov    %r15,(%rsp)
jmp    11bc <refineTransmap+0x7c>
mov    0x40(%rsp),%rdx
test   %rdx,%rdx
mov    0x48(%rsp),%rbx
jle    1473 <refineTransmap+0x333>
cmp    $0x1,%rdx
jne    13e8 <refineTransmap+0x2a8>
xor    %eax,%eax
test   $0x1,%dl
je     1473 <refineTransmap+0x333>
movss  (%rbx,%rax,4),%xmm1
subss  (%r14,%rax,4),%xmm1
movss  0xc45(%rip),%xmm0
ucomiss %xmm0,%xmm1
ja     146e <refineTransmap+0x32e>
xorps  %xmm0,%xmm0
ucomiss %xmm0,%xmm1
jbe    1469 <refineTransmap+0x329>
movaps %xmm1,%xmm0
jmp    146e <refineTransmap+0x32e>
mov    %rdx,%rcx
and    $0xfffffffffffffffe,%rcx
xor    %eax,%eax
movss  0xc17(%rip),%xmm0
xorps  %xmm1,%xmm1
jmp    1412 <refineTransmap+0x2d2>
xchg   %ax,%ax
movaps %xmm2,%xmm3
movss  %xmm3,0x4(%rbx,%rax,4)
add    $0x2,%rax
cmp    %rax,%rcx
je     13af <refineTransmap+0x26f>
movss  (%rbx,%rax,4),%xmm2
subss  (%r14,%rax,4),%xmm2
ucomiss %xmm0,%xmm2
movaps %xmm0,%xmm3
ja     1440 <refineTransmap+0x300>
ucomiss %xmm1,%xmm2
jbe    1430 <refineTransmap+0x2f0>
movaps %xmm2,%xmm3
jmp    1440 <refineTransmap+0x300>
nop
xorps  %xmm3,%xmm3
ucomiss %xmm2,%xmm3
jbe    1445 <refineTransmap+0x305>
nopl   0x0(%rax,%rax,1)
movss  %xmm3,(%rbx,%rax,4)
movss  0x4(%rbx,%rax,4),%xmm2
subss  0x4(%r14,%rax,4),%xmm2
ucomiss %xmm0,%xmm2
movaps %xmm0,%xmm3
ja     1403 <refineTransmap+0x2c3>
ucomiss %xmm1,%xmm2
ja     1400 <refineTransmap+0x2c0>
xorps  %xmm3,%xmm3
ucomiss %xmm2,%xmm3
ja     1403 <refineTransmap+0x2c3>
jmp    1409 <refineTransmap+0x2c9>
ucomiss %xmm1,%xmm0
jbe    1473 <refineTransmap+0x333>
movss  %xmm0,(%rbx,%rax,4)
mov    0x38(%rsp),%rdi
call   1040 <free@plt>
mov    %r14,%rdi
call   1040 <free@plt>
mov    %rbx,%rax
add    $0x68,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define float* @refineTransmap(float* returned %0, i64 %1, i64 %2) {
%4 = alloca [9 x float]
%5 = mul nsw i64 %2, %1
%6 = tail call float* @gamma_transform(double 1.000000e-01, double 1.300000e+00, float* %0, i64 %5)
%7 = bitcast [9 x float]* %4 to i8*
%8 = tail call i64 @calloc(i64 %5, i32 4)
%9 = inttoptr i64 %8 to float*
%10 = icmp sgt i64 %2, 0
br i1 %10, label %11, label %45
11: ; preds = %3
%12 = icmp sgt i64 %1, 0
%13 = getelementptr inbounds [9 x float], [9 x float]* %4, i64 0, i64 0
%14 = getelementptr inbounds [9 x float], [9 x float]* %4, i64 0, i64 1
%15 = getelementptr inbounds [9 x float], [9 x float]* %4, i64 0, i64 2
%16 = getelementptr inbounds [9 x float], [9 x float]* %4, i64 0, i64 1
%17 = getelementptr inbounds [9 x float], [9 x float]* %4, i64 0, i64 2
%18 = getelementptr inbounds [9 x float], [9 x float]* %4, i64 0, i64 3
%19 = getelementptr inbounds [9 x float], [9 x float]* %4, i64 0, i64 3
%20 = getelementptr inbounds [9 x float], [9 x float]* %4, i64 0, i64 4
%21 = getelementptr inbounds [9 x float], [9 x float]* %4, i64 0, i64 5
%22 = getelementptr inbounds [9 x float], [9 x float]* %4, i64 0, i64 5
%23 = getelementptr inbounds [9 x float], [9 x float]* %4, i64 0, i64 6
%24 = getelementptr inbounds [9 x float], [9 x float]* %4, i64 0, i64 6
%25 = getelementptr inbounds [9 x float], [9 x float]* %4, i64 0, i64 7
%26 = getelementptr inbounds [9 x float], [9 x float]* %4, i64 0, i64 7
%27 = getelementptr inbounds [9 x float], [9 x float]* %4, i64 0, i64 8
%28 = getelementptr inbounds [9 x float], [9 x float]* %4, i64 0, i64 8
br label %29
29: ; preds = %11, %109
%30 = phi i64 [ 0, %11 ], [ %110, %109 ]
br i1 %12, label %33, label %31
31: ; preds = %29
%32 = add nuw nsw i64 %30, 1
br label %109
33: ; preds = %29
%34 = add nsw i64 %30, -1
%35 = mul nsw i64 %30, %1
%36 = getelementptr inbounds float, float* %9, i64 %35
%37 = icmp slt i64 %30, 1
%38 = mul nsw i64 %34, %1
%39 = getelementptr inbounds float, float* %6, i64 %38
%40 = getelementptr inbounds float, float* %6, i64 %35
%41 = add nuw nsw i64 %30, 1
%42 = icmp slt i64 %41, %2
%43 = mul nsw i64 %41, %1
%44 = getelementptr inbounds float, float* %6, i64 %43
br label %52
45: ; preds = %109, %3
%46 = icmp sgt i64 %5, 0
br i1 %46, label %47, label %161
47: ; preds = %45
%48 = and i64 %5, 1
%49 = icmp eq i64 %5, 1
br i1 %49, label %145, label %50
50: ; preds = %47
%51 = and i64 %5, -2
br label %112
52: ; preds = %33, %105
%53 = phi i64 [ 0, %33 ], [ %74, %105 ]
%54 = add nsw i64 %53, -1
br i1 %37, label %60, label %55
55: ; preds = %52
%56 = icmp eq i64 %53, 0
br i1 %56, label %62, label %57
57: ; preds = %55
%58 = getelementptr inbounds float, float* %39, i64 %54
%59 = load float, float* %58
br label %62
60: ; preds = %52
store float 0.000000e+00, float* %13
store float 0.000000e+00, float* %16
%61 = add nuw nsw i64 %53, 1
br label %71
62: ; preds = %57, %55
%63 = phi float [ 0.000000e+00, %55 ], [ %59, %57 ]
store float %63, float* %13
%64 = getelementptr inbounds float, float* %39, i64 %53
%65 = load float, float* %64
store float %65, float* %14
%66 = add nuw nsw i64 %53, 1
%67 = icmp slt i64 %66, %1
br i1 %67, label %68, label %71
68: ; preds = %62
%69 = getelementptr inbounds float, float* %39, i64 %66
%70 = load float, float* %69
store float %70, float* %15
br label %73
71: ; preds = %60, %62
%72 = phi i64 [ %61, %60 ], [ %66, %62 ]
store float 0.000000e+00, float* %17
br label %73
73: ; preds = %68, %71
%74 = phi i64 [ %72, %71 ], [ %66, %68 ]
%75 = icmp eq i64 %53, 0
br i1 %75, label %79, label %76
76: ; preds = %73
%77 = getelementptr inbounds float, float* %40, i64 %54
%78 = load float, float* %77
store float %78, float* %18
br label %80
79: ; preds = %73
store float 0.000000e+00, float* %19
br label %80
80: ; preds = %79, %76
%81 = getelementptr inbounds float, float* %40, i64 %53
%82 = load float, float* %81
store float %82, float* %20
%83 = icmp slt i64 %74, %1
br i1 %83, label %84, label %87
84: ; preds = %80
%85 = getelementptr inbounds float, float* %40, i64 %74
%86 = load float, float* %85
store float %86, float* %22
br label %88
87: ; preds = %80
store float 0.000000e+00, float* %21
br label %88
88: ; preds = %84, %87
%89 = icmp ne i64 %53, 0
%90 = select i1 %42, i1 %89, i1 false
br i1 %90, label %91, label %94
91: ; preds = %88
%92 = getelementptr inbounds float, float* %44, i64 %54
%93 = load float, float* %92
store float %93, float* %24
br label %95
94: ; preds = %88
store float 0.000000e+00, float* %23
br label %95
95: ; preds = %91, %94
br i1 %42, label %97, label %96
96: ; preds = %95
store float 0.000000e+00, float* %25
br label %104
97: ; preds = %95
%98 = getelementptr inbounds float, float* %44, i64 %53
%99 = load float, float* %98
store float %99, float* %26
%100 = icmp slt i64 %74, %1
br i1 %100, label %101, label %104
101: ; preds = %97
%102 = getelementptr inbounds float, float* %44, i64 %74
%103 = load float, float* %102
store float %103, float* %28
br label %105
104: ; preds = %96, %97
store float 0.000000e+00, float* %27
br label %105
105: ; preds = %104, %101
%106 = call float @median_filt(float* %13)
%107 = getelementptr inbounds float, float* %36, i64 %53
store float %106, float* %107
%108 = icmp eq i64 %74, %1
br i1 %108, label %109, label %52
109: ; preds = %105, %31
%110 = phi i64 [ %32, %31 ], [ %41, %105 ]
%111 = icmp eq i64 %110, %2
br i1 %111, label %45, label %29
112: ; preds = %141, %50
%113 = phi i64 [ 0, %50 ], [ %142, %141 ]
%114 = phi i64 [ 0, %50 ], [ %143, %141 ]
%115 = getelementptr inbounds float, float* %0, i64 %113
%116 = load float, float* %115
%117 = getelementptr inbounds float, float* %9, i64 %113
%118 = load float, float* %117
%119 = fsub float %116, %118
%120 = fcmp ogt float %119, 1.000000e+00
br i1 %120, label %125, label %121
121: ; preds = %112
%122 = fcmp ogt float %119, 0.000000e+00
br i1 %122, label %125, label %123
123: ; preds = %121
%124 = fcmp olt float %119, 0.000000e+00
br i1 %124, label %125, label %127
125: ; preds = %123, %121, %112
%126 = phi float [ 1.000000e+00, %112 ], [ %119, %121 ], [ 0.000000e+00, %123 ]
store float %126, float* %115
br label %127
127: ; preds = %125, %123
%128 = or i64 %113, 1
%129 = getelementptr inbounds float, float* %0, i64 %128
%130 = load float, float* %129
%131 = getelementptr inbounds float, float* %9, i64 %128
%132 = load float, float* %131
%133 = fsub float %130, %132
%134 = fcmp ogt float %133, 1.000000e+00
br i1 %134, label %139, label %135
135: ; preds = %127
%136 = fcmp ogt float %133, 0.000000e+00
br i1 %136, label %139, label %137
137: ; preds = %135
%138 = fcmp olt float %133, 0.000000e+00
br i1 %138, label %139, label %141
139: ; preds = %137, %135, %127
%140 = phi float [ 1.000000e+00, %127 ], [ %133, %135 ], [ 0.000000e+00, %137 ]
store float %140, float* %129
br label %141
141: ; preds = %139, %137
%142 = add nuw nsw i64 %113, 2
%143 = add i64 %114, 2
%144 = icmp eq i64 %143, %51
br i1 %144, label %145, label %112
145: ; preds = %141, %47
%146 = phi i64 [ 0, %47 ], [ %142, %141 ]
%147 = icmp eq i64 %48, 0
br i1 %147, label %161, label %148
148: ; preds = %145
%149 = getelementptr inbounds float, float* %0, i64 %146
%150 = load float, float* %149
%151 = getelementptr inbounds float, float* %9, i64 %146
%152 = load float, float* %151
%153 = fsub float %150, %152
%154 = fcmp ogt float %153, 1.000000e+00
br i1 %154, label %159, label %155
155: ; preds = %148
%156 = fcmp ogt float %153, 0.000000e+00
br i1 %156, label %159, label %157
157: ; preds = %155
%158 = fcmp olt float %153, 0.000000e+00
br i1 %158, label %159, label %161
159: ; preds = %157, %155, %148
%160 = phi float [ 1.000000e+00, %148 ], [ %153, %155 ], [ 0.000000e+00, %157 ]
store float %160, float* %149
br label %161
161: ; preds = %145, %159, %157, %45
%162 = call i32 @free(float* %6)
%163 = call i32 @free(float* %9)
ret float* %0
}