#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 /*<<< orphan*/ M_PI;
int /*<<< orphan*/ addfloats(float *, float *, float *, int);
int /*<<< orphan*/ calculateMagnitude(float *, float *, float *, int);
int /*<<< orphan*/ calculatePhase(float *, float *, float *, int,
                                  int /*<<< orphan*/, double);
float fmax(float, float);
int /*<<< orphan*/ free(float *);
scalar_t__ malloc(int);
int /*<<< orphan*/ multiplyfloats(float *, float *, float *, int);
int /*<<< orphan*/ scalefloats(float *, int, double);
float sqrt(int);
int /*<<< orphan*/ unwrapPhase(float *, float *, int, int, int /*<<< orphan*/);
void calcSpiralDcf(float *gx, float *gy, float *kx, float *ky, int rolen,
                   float *denscomp) {
  float *kr, *gradientMagnitude;
  float *Gtwist;
  float *denscomp_it, *denscomp_is;
  float *kPhasexy, *kMagxy;
  int n;
  kr = (float *)malloc(rolen * sizeof(float));
  gradientMagnitude = (float *)malloc(rolen * sizeof(float));
  for (n = 0; n < rolen; n++) {
    kr[n] = sqrt(kx[n] * kx[n] + ky[n] * ky[n]);
    gradientMagnitude[n] = sqrt(gx[n] * gx[n] + gy[n] * gy[n]);
  }
  Gtwist = (float *)malloc(rolen * sizeof(float));
  kPhasexy = (float *)malloc(rolen * sizeof(float));
  kMagxy = (float *)malloc(rolen * sizeof(float));
  calculatePhase(kx, ky, kPhasexy, rolen, 0, 1.0);
  calculateMagnitude(kx, ky, kMagxy, rolen);
  unwrapPhase(kPhasexy, Gtwist, rolen, rolen, M_PI);
  for (n = 0; n < rolen - 1; n++)
    Gtwist[n] =
        fmax((Gtwist[n + 1] - Gtwist[n]) / (kr[n + 1] - kr[n]) * kr[n], 0.0f);
  Gtwist[rolen - 1] = Gtwist[rolen - 2];
  denscomp_it = (float *)malloc(rolen * sizeof(float));
  for (n = 0; n < rolen; n++) {
    if (Gtwist[n] >= 0)
      denscomp_it[n] = kMagxy[n] / sqrt(1 + Gtwist[n] * Gtwist[n]);
    else
      denscomp_it[n] = 0;
  }
  denscomp_is = (float *)malloc(rolen * sizeof(float));
  addfloats(gradientMagnitude, &(gradientMagnitude[1]), denscomp_is, rolen - 1);
  denscomp_is[rolen - 1] = 2 * gradientMagnitude[rolen - 1];
  scalefloats(denscomp_is, rolen, 0.5);
  multiplyfloats(denscomp_is, denscomp_it, denscomp, rolen);
  free(kr);
  free(gradientMagnitude);
  free(Gtwist);
  free(denscomp_it);
  free(denscomp_is);
  free(kPhasexy);
  free(kMagxy);
  return;
}


x86 assembly (clang, -O0):
<calcSpiralDcf>:
push   %rbp
mov    %rsp,%rbp
sub    $0x70,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8d,-0x24(%rbp)
mov    %r9,-0x30(%rbp)
movslq -0x24(%rbp),%rax
shl    $0x2,%rax
mov    %eax,%edi
call   1070 <malloc@plt>
mov    %rax,-0x38(%rbp)
movslq -0x24(%rbp),%rax
shl    $0x2,%rax
mov    %eax,%edi
call   1070 <malloc@plt>
mov    %rax,-0x40(%rbp)
movl   $0x0,-0x6c(%rbp)
mov    -0x6c(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    12ab <calcSpiralDcf+0x10b>
mov    -0x18(%rbp),%rax
movslq -0x6c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x6c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
mov    -0x20(%rbp),%rax
movslq -0x6c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
mov    -0x20(%rbp),%rax
movslq -0x6c(%rbp),%rcx
mulss  (%rax,%rcx,4),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
cvttss2si %xmm0,%edi
call   10c0 <sqrt@plt>
mov    -0x38(%rbp),%rax
movslq -0x6c(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x8(%rbp),%rax
movslq -0x6c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x6c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
mov    -0x10(%rbp),%rax
movslq -0x6c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
mov    -0x10(%rbp),%rax
movslq -0x6c(%rbp),%rcx
mulss  (%rax,%rcx,4),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
cvttss2si %xmm0,%edi
call   10c0 <sqrt@plt>
mov    -0x40(%rbp),%rax
movslq -0x6c(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x6c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x6c(%rbp)
jmp    11ed <calcSpiralDcf+0x4d>
movslq -0x24(%rbp),%rax
shl    $0x2,%rax
mov    %eax,%edi
call   1070 <malloc@plt>
mov    %rax,-0x48(%rbp)
movslq -0x24(%rbp),%rax
shl    $0x2,%rax
mov    %eax,%edi
call   1070 <malloc@plt>
mov    %rax,-0x60(%rbp)
movslq -0x24(%rbp),%rax
shl    $0x2,%rax
mov    %eax,%edi
call   1070 <malloc@plt>
mov    %rax,-0x68(%rbp)
mov    -0x18(%rbp),%rdi
mov    -0x20(%rbp),%rsi
mov    -0x60(%rbp),%rdx
mov    -0x24(%rbp),%ecx
xor    %r8d,%r8d
movsd  0xd02(%rip),%xmm0
call   1060 <calculatePhase@plt>
mov    -0x18(%rbp),%rdi
mov    -0x20(%rbp),%rsi
mov    -0x68(%rbp),%rdx
mov    -0x24(%rbp),%ecx
call   10b0 <calculateMagnitude@plt>
mov    -0x60(%rbp),%rdi
mov    -0x48(%rbp),%rsi
mov    -0x24(%rbp),%edx
mov    -0x24(%rbp),%ecx
mov    0x2cbc(%rip),%rax
mov    (%rax),%r8d
call   10a0 <unwrapPhase@plt>
movl   $0x0,-0x6c(%rbp)
mov    -0x6c(%rbp),%eax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    13be <calcSpiralDcf+0x21e>
mov    -0x48(%rbp),%rax
mov    -0x6c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x48(%rbp),%rax
movslq -0x6c(%rbp),%rcx
subss  (%rax,%rcx,4),%xmm0
mov    -0x38(%rbp),%rax
mov    -0x6c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm1
mov    -0x38(%rbp),%rax
movslq -0x6c(%rbp),%rcx
subss  (%rax,%rcx,4),%xmm1
divss  %xmm1,%xmm0
mov    -0x38(%rbp),%rax
movslq -0x6c(%rbp),%rcx
mulss  (%rax,%rcx,4),%xmm0
xorps  %xmm1,%xmm1
call   1090 <fmax@plt>
mov    -0x48(%rbp),%rax
movslq -0x6c(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x6c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x6c(%rbp)
jmp    133b <calcSpiralDcf+0x19b>
mov    -0x48(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x2,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x48(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
movslq -0x24(%rbp),%rax
shl    $0x2,%rax
mov    %eax,%edi
call   1070 <malloc@plt>
mov    %rax,-0x50(%rbp)
movl   $0x0,-0x6c(%rbp)
mov    -0x6c(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    14a7 <calcSpiralDcf+0x307>
mov    -0x48(%rbp),%rax
movslq -0x6c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
xorps  %xmm1,%xmm1
ucomiss %xmm1,%xmm0
jb     1484 <calcSpiralDcf+0x2e4>
mov    -0x68(%rbp),%rax
movslq -0x6c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x70(%rbp)
mov    -0x48(%rbp),%rax
movslq -0x6c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x48(%rbp),%rax
movslq -0x6c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
mulss  %xmm1,%xmm0
movss  0xbbb(%rip),%xmm1
addss  %xmm1,%xmm0
cvttss2si %xmm0,%edi
call   10c0 <sqrt@plt>
movaps %xmm0,%xmm1
movss  -0x70(%rbp),%xmm0
divss  %xmm1,%xmm0
mov    -0x50(%rbp),%rax
movslq -0x6c(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
jmp    1494 <calcSpiralDcf+0x2f4>
mov    -0x50(%rbp),%rax
movslq -0x6c(%rbp),%rcx
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax,%rcx,4)
jmp    1499 <calcSpiralDcf+0x2f9>
mov    -0x6c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x6c(%rbp)
jmp    13fc <calcSpiralDcf+0x25c>
movslq -0x24(%rbp),%rax
shl    $0x2,%rax
mov    %eax,%edi
call   1070 <malloc@plt>
mov    %rax,-0x58(%rbp)
mov    -0x40(%rbp),%rdi
mov    -0x40(%rbp),%rsi
add    $0x4,%rsi
mov    -0x58(%rbp),%rdx
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
call   1080 <addfloats@plt>
mov    -0x40(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  0xb26(%rip),%xmm0
mulss  (%rax,%rcx,4),%xmm0
mov    -0x58(%rbp),%rax
mov    -0x24(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x58(%rbp),%rdi
mov    -0x24(%rbp),%esi
movsd  0xaf8(%rip),%xmm0
call   1040 <scalefloats@plt>
mov    -0x58(%rbp),%rdi
mov    -0x50(%rbp),%rsi
mov    -0x30(%rbp),%rdx
mov    -0x24(%rbp),%ecx
call   1050 <multiplyfloats@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>
mov    -0x58(%rbp),%rdi
call   1030 <free@plt>
mov    -0x60(%rbp),%rdi
call   1030 <free@plt>
mov    -0x68(%rbp),%rdi
call   1030 <free@plt>
add    $0x70,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<calcSpiralDcf>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x58,%rsp
mov    %r9,0x50(%rsp)
mov    %r8d,%r12d
mov    %rcx,0x10(%rsp)
mov    %rdx,%r15
mov    %rsi,%r13
mov    %rdi,%rbp
lea    0x0(,%r12,4),%ebx
mov    %ebx,%edi
call   1070 <malloc@plt>
mov    %rax,0x18(%rsp)
mov    %ebx,0xc(%rsp)
mov    %ebx,%edi
call   1070 <malloc@plt>
mov    %rax,%r14
mov    %r12,%rax
mov    %eax,%r12d
mov    %rax,0x28(%rsp)
test   %eax,%eax
jle    1250 <calcSpiralDcf+0xc0>
xor    %ebx,%ebx
nopw   0x0(%rax,%rax,1)
movss  (%r15,%rbx,4),%xmm0
mov    0x10(%rsp),%rax
movss  (%rax,%rbx,4),%xmm1
mulss  %xmm1,%xmm1
mulss  %xmm0,%xmm0
addss  %xmm1,%xmm0
cvttss2si %xmm0,%edi
call   10b0 <sqrt@plt>
mov    0x18(%rsp),%rax
movss  %xmm0,(%rax,%rbx,4)
movss  0x0(%rbp,%rbx,4),%xmm0
movss  0x0(%r13,%rbx,4),%xmm1
mulss  %xmm1,%xmm1
mulss  %xmm0,%xmm0
addss  %xmm1,%xmm0
cvttss2si %xmm0,%edi
call   10b0 <sqrt@plt>
movss  %xmm0,(%r14,%rbx,4)
add    $0x1,%rbx
cmp    %rbx,%r12
jne    11f0 <calcSpiralDcf+0x60>
mov    %r12,0x30(%rsp)
mov    %r14,0x48(%rsp)
mov    0xc(%rsp),%ebx
mov    %ebx,%edi
call   1070 <malloc@plt>
mov    %rax,%r13
mov    %ebx,%edi
call   1070 <malloc@plt>
mov    %rax,%rbp
mov    %ebx,%edi
call   1070 <malloc@plt>
mov    %rax,%r12
movsd  0xd7c(%rip),%xmm0
mov    %r15,%r14
mov    %r15,%rdi
mov    0x10(%rsp),%r15
mov    %r15,%rsi
mov    %rbp,%rdx
mov    0x28(%rsp),%rbx
mov    %ebx,%ecx
xor    %r8d,%r8d
call   1060 <calculatePhase@plt>
mov    %r14,%rdi
mov    %r15,%rsi
mov    %r12,%r15
mov    %r12,%rdx
mov    %ebx,%ecx
call   10a0 <calculateMagnitude@plt>
mov    0x2d2a(%rip),%rax
mov    (%rax),%r8d
mov    %rbp,0x40(%rsp)
mov    %rbp,%rdi
mov    %r13,%rsi
mov    %ebx,%edx
mov    %ebx,%ecx
call   1090 <unwrapPhase@plt>
lea    -0x1(%rbx),%edi
cmp    $0x2,%ebx
mov    0x18(%rsp),%rsi
jl     135f <calcSpiralDcf+0x1cf>
mov    %edi,%eax
movss  0x0(%r13),%xmm0
cmp    $0x4,%edi
jb     1314 <calcSpiralDcf+0x184>
lea    (%rsi,%rax,4),%rcx
add    $0x4,%rcx
cmp    %rcx,%r13
jae    149e <calcSpiralDcf+0x30e>
lea    0x4(,%rax,4),%rcx
add    %r13,%rcx
cmp    %rcx,%rsi
jae    149e <calcSpiralDcf+0x30e>
xor    %ecx,%ecx
xorps  %xmm1,%xmm1
nopl   0x0(%rax)
movss  0x4(%r13,%rcx,4),%xmm2
movaps %xmm2,%xmm3
subss  %xmm0,%xmm3
movss  (%rsi,%rcx,4),%xmm0
movss  0x4(%rsi,%rcx,4),%xmm4
subss  %xmm0,%xmm4
divss  %xmm4,%xmm3
mulss  %xmm0,%xmm3
maxss  %xmm1,%xmm3
movss  %xmm3,0x0(%r13,%rcx,4)
lea    0x1(%rcx),%rdx
movaps %xmm2,%xmm0
mov    %rdx,%rcx
cmp    %rdx,%rax
jne    1320 <calcSpiralDcf+0x190>
movslq %ebx,%rax
movss  -0x8(%r13,%rax,4),%xmm0
mov    %edi,0x24(%rsp)
movslq %edi,%rax
mov    %rax,0x38(%rsp)
movss  %xmm0,0x0(%r13,%rax,4)
mov    0xc(%rsp),%r14d
mov    %r14d,%edi
call   1070 <malloc@plt>
mov    %rax,%rbp
test   %ebx,%ebx
mov    0x30(%rsp),%rax
jle    13f0 <calcSpiralDcf+0x260>
xor    %ebx,%ebx
jmp    13af <calcSpiralDcf+0x21f>
nopl   0x0(%rax)
movss  %xmm1,0x0(%rbp,%rbx,4)
add    $0x1,%rbx
cmp    %rbx,%rax
je     13f0 <calcSpiralDcf+0x260>
movss  0x0(%r13,%rbx,4),%xmm0
xorps  %xmm1,%xmm1
ucomiss %xmm1,%xmm0
jb     13a0 <calcSpiralDcf+0x210>
movss  (%r15,%rbx,4),%xmm1
movss  %xmm1,0x10(%rsp)
mulss  %xmm0,%xmm0
addss  0xc3a(%rip),%xmm0
cvttss2si %xmm0,%edi
call   10b0 <sqrt@plt>
movss  0x10(%rsp),%xmm1
mov    0x30(%rsp),%rax
divss  %xmm0,%xmm1
jmp    13a0 <calcSpiralDcf+0x210>
mov    %r14d,%edi
call   1070 <malloc@plt>
mov    %rax,%r14
mov    0x48(%rsp),%r12
mov    %r12,%rsi
add    $0x4,%rsi
mov    %r12,%rdi
mov    %rax,%rdx
mov    0x24(%rsp),%ecx
call   1080 <addfloats@plt>
mov    0x38(%rsp),%rax
movss  (%r12,%rax,4),%xmm0
addss  %xmm0,%xmm0
movss  %xmm0,(%r14,%rax,4)
movsd  0xbd5(%rip),%xmm0
mov    %r14,%rdi
mov    0x28(%rsp),%rbx
mov    %ebx,%esi
call   1040 <scalefloats@plt>
mov    %r14,%rdi
mov    %rbp,%rsi
mov    0x50(%rsp),%rdx
mov    %ebx,%ecx
call   1050 <multiplyfloats@plt>
mov    0x18(%rsp),%rdi
call   1030 <free@plt>
mov    %r12,%rdi
call   1030 <free@plt>
mov    %r13,%rdi
call   1030 <free@plt>
mov    %rbp,%rdi
call   1030 <free@plt>
mov    %r14,%rdi
call   1030 <free@plt>
mov    0x40(%rsp),%rdi
call   1030 <free@plt>
mov    %r15,%rdi
add    $0x58,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
jmp    1030 <free@plt>
mov    %eax,%ecx
and    $0xfffffffc,%ecx
shufps $0x0,%xmm0,%xmm0
xor    %edx,%edx
nopl   0x0(%rax)
movaps %xmm0,%xmm1
movups 0x4(%r13,%rdx,4),%xmm0
shufps $0x3,%xmm0,%xmm1
shufps $0x98,%xmm0,%xmm1
movaps %xmm0,%xmm2
subps  %xmm1,%xmm2
movups (%rsi,%rdx,4),%xmm1
movups 0x4(%rsi,%rdx,4),%xmm3
subps  %xmm1,%xmm3
divps  %xmm3,%xmm2
mulps  %xmm1,%xmm2
xorps  %xmm1,%xmm1
maxps  %xmm2,%xmm1
cmpordps %xmm2,%xmm2
andps  %xmm1,%xmm2
movups %xmm2,0x0(%r13,%rdx,4)
add    $0x4,%rdx
cmp    %rdx,%rcx
jne    14b0 <calcSpiralDcf+0x320>
cmp    %rax,%rcx
je     135f <calcSpiralDcf+0x1cf>
shufps $0xff,%xmm0,%xmm0
jmp    1316 <calcSpiralDcf+0x186>

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

LLVM IR (-O2):
@M_PI = global i32 0
define void @calcSpiralDcf(float* nocapture %0, float* nocapture %1, float* %2, float* %3, i32 %4, float* %5) {
%7 = shl i32 %4, 2
%8 = tail call i64 @malloc(i32 %7)
%9 = inttoptr i64 %8 to float*
%10 = inttoptr i64 %8 to i8*
%11 = tail call i64 @malloc(i32 %7)
%12 = inttoptr i64 %11 to float*
%13 = icmp sgt i32 %4, 0
br i1 %13, label %14, label %38
14: ; preds = %6
%15 = zext i32 %4 to i64
br label %16
16: ; preds = %14, %16
%17 = phi i64 [ 0, %14 ], [ %36, %16 ]
%18 = getelementptr inbounds float, float* %2, i64 %17
%19 = load float, float* %18
%20 = getelementptr inbounds float, float* %3, i64 %17
%21 = load float, float* %20
%22 = fmul float %21, %21
%23 = tail call float @llvm.fmuladd.f32(float %19, float %19, float %22)
%24 = fptosi float %23 to i32
%25 = tail call float @sqrt(i32 %24)
%26 = getelementptr inbounds float, float* %9, i64 %17
store float %25, float* %26
%27 = getelementptr inbounds float, float* %0, i64 %17
%28 = load float, float* %27
%29 = getelementptr inbounds float, float* %1, i64 %17
%30 = load float, float* %29
%31 = fmul float %30, %30
%32 = tail call float @llvm.fmuladd.f32(float %28, float %28, float %31)
%33 = fptosi float %32 to i32
%34 = tail call float @sqrt(i32 %33)
%35 = getelementptr inbounds float, float* %12, i64 %17
store float %34, float* %35
%36 = add nuw nsw i64 %17, 1
%37 = icmp eq i64 %36, %15
br i1 %37, label %38, label %16
38: ; preds = %16, %6
%39 = tail call i64 @malloc(i32 %7)
%40 = inttoptr i64 %39 to float*
%41 = inttoptr i64 %39 to i8*
%42 = tail call i64 @malloc(i32 %7)
%43 = inttoptr i64 %42 to float*
%44 = tail call i64 @malloc(i32 %7)
%45 = inttoptr i64 %44 to float*
%46 = tail call i32 @calculatePhase(float* %2, float* %3, float* %43, i32 %4, i32 0, double 1.000000e+00)
%47 = tail call i32 @calculateMagnitude(float* %2, float* %3, float* %45, i32 %4)
%48 = load i32, i32* @M_PI
%49 = tail call i32 @unwrapPhase(float* %43, float* %40, i32 %4, i32 %4, i32 %48)
%50 = add i32 %4, -1
%51 = icmp sgt i32 %4, 1
br i1 %51, label %52, label %114
52: ; preds = %38
%53 = zext i32 %50 to i64
%54 = load float, float* %40
%55 = icmp ult i32 %50, 4
br i1 %55, label %94, label %56
56: ; preds = %52
%57 = add nuw nsw i64 %53, 1
%58 = getelementptr float, float* %40, i64 %57
%59 = bitcast float* %58 to i8*
%60 = getelementptr float, float* %9, i64 %57
%61 = bitcast float* %60 to i8*
%62 = icmp ult i8* %41, %61
%63 = icmp ult i8* %10, %59
%64 = and i1 %62, %63
br i1 %64, label %94, label %65
65: ; preds = %56
%66 = and i64 %53, 4294967292
%67 = insertelement <4 x float> poison, float %54, i64 3
br label %68
68: ; preds = %68, %65
%69 = phi i64 [ 0, %65 ], [ %89, %68 ]
%70 = phi <4 x float> [ %67, %65 ], [ %74, %68 ]
%71 = or i64 %69, 1
%72 = getelementptr inbounds float, float* %40, i64 %71
%73 = bitcast float* %72 to <4 x float>*
%74 = load <4 x float>, <4 x float>* %73
%75 = shufflevector <4 x float> %70, <4 x float> %74, <4 x i32> <i32 3, i32 4, i32 5, i32 6>
%76 = getelementptr inbounds float, float* %40, i64 %69
%77 = fsub <4 x float> %74, %75
%78 = getelementptr inbounds float, float* %9, i64 %71
%79 = bitcast float* %78 to <4 x float>*
%80 = load <4 x float>, <4 x float>* %79
%81 = getelementptr inbounds float, float* %9, i64 %69
%82 = bitcast float* %81 to <4 x float>*
%83 = load <4 x float>, <4 x float>* %82
%84 = fsub <4 x float> %80, %83
%85 = fdiv <4 x float> %77, %84
%86 = fmul <4 x float> %83, %85
%87 = call nsz <4 x float> @llvm.maxnum.v4f32(<4 x float> %86, <4 x float> zeroinitializer)
%88 = bitcast float* %76 to <4 x float>*
store <4 x float> %87, <4 x float>* %88
%89 = add nuw i64 %69, 4
%90 = icmp eq i64 %89, %66
br i1 %90, label %91, label %68
91: ; preds = %68
%92 = icmp eq i64 %66, %53
%93 = extractelement <4 x float> %74, i64 3
br i1 %92, label %114, label %94
94: ; preds = %56, %52, %91
%95 = phi float [ %93, %91 ], [ %54, %52 ], [ %54, %56 ]
%96 = phi i64 [ %66, %91 ], [ 0, %52 ], [ 0, %56 ]
br label %97
97: ; preds = %94, %97
%98 = phi float [ %102, %97 ], [ %95, %94 ]
%99 = phi i64 [ %100, %97 ], [ %96, %94 ]
%100 = add nuw nsw i64 %99, 1
%101 = getelementptr inbounds float, float* %40, i64 %100
%102 = load float, float* %101
%103 = getelementptr inbounds float, float* %40, i64 %99
%104 = fsub float %102, %98
%105 = getelementptr inbounds float, float* %9, i64 %100
%106 = load float, float* %105
%107 = getelementptr inbounds float, float* %9, i64 %99
%108 = load float, float* %107
%109 = fsub float %106, %108
%110 = fdiv float %104, %109
%111 = fmul float %108, %110
%112 = tail call nsz float @llvm.maxnum.f32(float %111, float 0.000000e+00)
store float %112, float* %103
%113 = icmp eq i64 %100, %53
br i1 %113, label %114, label %97
114: ; preds = %97, %91, %38
%115 = add nsw i32 %4, -2
%116 = sext i32 %115 to i64
%117 = getelementptr inbounds float, float* %40, i64 %116
%118 = load float, float* %117
%119 = sext i32 %50 to i64
%120 = getelementptr inbounds float, float* %40, i64 %119
store float %118, float* %120
%121 = tail call i64 @malloc(i32 %7)
%122 = inttoptr i64 %121 to float*
br i1 %13, label %123, label %142
123: ; preds = %114
%124 = zext i32 %4 to i64
br label %125
125: ; preds = %123, %137
%126 = phi i64 [ 0, %123 ], [ %140, %137 ]
%127 = getelementptr inbounds float, float* %40, i64 %126
%128 = load float, float* %127
%129 = fcmp ult float %128, 0.000000e+00
br i1 %129, label %137, label %130
130: ; preds = %125
%131 = getelementptr inbounds float, float* %45, i64 %126
%132 = load float, float* %131
%133 = tail call float @llvm.fmuladd.f32(float %128, float %128, float 1.000000e+00)
%134 = fptosi float %133 to i32
%135 = tail call float @sqrt(i32 %134)
%136 = fdiv float %132, %135
br label %137
137: ; preds = %125, %130
%138 = phi float [ %136, %130 ], [ 0.000000e+00, %125 ]
%139 = getelementptr inbounds float, float* %122, i64 %126
store float %138, float* %139
%140 = add nuw nsw i64 %126, 1
%141 = icmp eq i64 %140, %124
br i1 %141, label %142, label %125
142: ; preds = %137, %114
%143 = tail call i64 @malloc(i32 %7)
%144 = inttoptr i64 %143 to float*
%145 = getelementptr inbounds float, float* %12, i64 1
%146 = tail call i32 @addfloats(float* %12, float* %145, float* %144, i32 %50)
%147 = getelementptr inbounds float, float* %12, i64 %119
%148 = load float, float* %147
%149 = fmul float %148, 2.000000e+00
%150 = getelementptr inbounds float, float* %144, i64 %119
store float %149, float* %150
%151 = tail call i32 @scalefloats(float* %144, i32 %4, double 5.000000e-01)
%152 = tail call i32 @multiplyfloats(float* %144, float* %122, float* %5, i32 %4)
%153 = tail call i32 @free(float* %9)
%154 = tail call i32 @free(float* %12)
%155 = tail call i32 @free(float* %40)
%156 = tail call i32 @free(float* %122)
%157 = tail call i32 @free(float* %144)
%158 = tail call i32 @free(float* %43)
%159 = tail call i32 @free(float* %45)
ret void
}