#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 M_PI;
void unwrapPhase(float *phasedata, float *uPhaseData, int numPoints, int frsize,
                 float tol) {
  int i;
  float nextPhase, currentPhase;
  float phaseDiff;
  static float offset = 0.0f;
  offset = -phasedata[0];
  uPhaseData[0] = 0;
  for (i = 0; i < (numPoints - 1); i++) {
    nextPhase = phasedata[i + 1];
    currentPhase = phasedata[i];
    phaseDiff = nextPhase - currentPhase;
    if ((i + 1) % frsize == 0)
      offset = -phasedata[i + 1];
    else {
      if (phaseDiff < -tol)
        offset += (float)(2 * M_PI);
      else if (phaseDiff > tol)
        offset -= (float)(2 * M_PI);
    }
    uPhaseData[i + 1] = phasedata[i + 1] + offset;
  }
  return;
}


x86 assembly (clang, -O0):
<unwrapPhase>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
movss  %xmm0,-0x1c(%rbp)
mov    -0x8(%rbp),%rax
movss  (%rax),%xmm0
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
movss  %xmm0,0x2ef0(%rip)
mov    -0x10(%rbp),%rax
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax)
movl   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%eax
mov    -0x14(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    128a <unwrapPhase+0x18a>
mov    -0x8(%rbp),%rax
mov    -0x20(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x24(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x20(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x28(%rbp)
movss  -0x24(%rbp),%xmm0
subss  -0x28(%rbp),%xmm0
movss  %xmm0,-0x2c(%rbp)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
cltd
idivl  -0x18(%rbp)
cmp    $0x0,%edx
jne    11ce <unwrapPhase+0xce>
mov    -0x8(%rbp),%rax
mov    -0x20(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
movss  %xmm0,0x2e5b(%rip)
jmp    1250 <unwrapPhase+0x150>
movss  -0x2c(%rbp),%xmm1
movss  -0x1c(%rbp),%xmm0
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
ucomiss %xmm1,%xmm0
jbe    1213 <unwrapPhase+0x113>
mov    0x2df3(%rip),%rax
mov    (%rax),%eax
shl    $0x1,%eax
cvtsi2ss %eax,%xmm0
addss  0x2e1e(%rip),%xmm0
movss  %xmm0,0x2e16(%rip)
jmp    124b <unwrapPhase+0x14b>
movss  -0x2c(%rbp),%xmm0
ucomiss -0x1c(%rbp),%xmm0
jbe    1246 <unwrapPhase+0x146>
mov    0x2dbf(%rip),%rax
mov    (%rax),%eax
shl    $0x1,%eax
cvtsi2ss %eax,%xmm1
movss  0x2dea(%rip),%xmm0
subss  %xmm1,%xmm0
movss  %xmm0,0x2dde(%rip)
jmp    124b <unwrapPhase+0x14b>
jmp    1250 <unwrapPhase+0x150>
mov    -0x8(%rbp),%rax
mov    -0x20(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
addss  0x2dba(%rip),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    1146 <unwrapPhase+0x46>
pop    %rbp
ret

x86 assembly (clang, -O2):
<unwrapPhase>:
movss  (%rdi),%xmm1
movl   $0x0,(%rsi)
cmp    $0x2,%edx
jl     11ab <unwrapPhase+0xab>
mov    %edx,%r8d
movaps 0xee3(%rip),%xmm2
xorps  %xmm2,%xmm1
mov    0x2ec1(%rip),%rax
mov    (%rax),%eax
add    %eax,%eax
cvtsi2ss %eax,%xmm3
add    $0xffffffff,%r8d
movaps %xmm0,%xmm4
xorps  %xmm2,%xmm4
neg    %r8
mov    $0x1,%r9d
jmp    1172 <unwrapPhase+0x72>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movaps %xmm5,%xmm1
xorps  %xmm2,%xmm1
addss  %xmm1,%xmm5
movss  %xmm5,(%rsi,%r9,4)
lea    (%r8,%r9,1),%rax
add    $0x1,%rax
add    $0x1,%r9
cmp    $0x1,%rax
je     11ab <unwrapPhase+0xab>
movss  (%rdi,%r9,4),%xmm5
mov    %r9d,%eax
cltd
idiv   %ecx
test   %edx,%edx
je     1150 <unwrapPhase+0x50>
movaps %xmm5,%xmm6
subss  -0x4(%rdi,%r9,4),%xmm6
ucomiss %xmm6,%xmm4
jbe    11a0 <unwrapPhase+0xa0>
addss  %xmm3,%xmm1
jmp    1156 <unwrapPhase+0x56>
nopw   0x0(%rax,%rax,1)
ucomiss %xmm0,%xmm6
jbe    1156 <unwrapPhase+0x56>
subss  %xmm3,%xmm1
jmp    1156 <unwrapPhase+0x56>
ret

LLVM IR (-O0):
@unwrapPhase.offset = internal global float 0.000000e+00
@M_PI = global i32 0
define void @unwrapPhase(float* %0, float* %1, i32 %2, i32 %3, float %4) {
%6 = alloca float*
%7 = alloca float*
%8 = alloca i32
%9 = alloca i32
%10 = alloca float
%11 = alloca i32
%12 = alloca float
%13 = alloca float
%14 = alloca float
store float* %0, float** %6
store float* %1, float** %7
store i32 %2, i32* %8
store i32 %3, i32* %9
store float %4, float* %10
%15 = load float*, float** %6
%16 = getelementptr inbounds float, float* %15, i64 0
%17 = load float, float* %16
%18 = fneg float %17
store float %18, float* @unwrapPhase.offset
%19 = load float*, float** %7
%20 = getelementptr inbounds float, float* %19, i64 0
store float 0.000000e+00, float* %20
store i32 0, i32* %11
br label %21
21: ; preds = %91, %5
%22 = load i32, i32* %11
%23 = load i32, i32* %8
%24 = sub nsw i32 %23, 1
%25 = icmp slt i32 %22, %24
br i1 %25, label %26, label %94
26: ; preds = %21
%27 = load float*, float** %6
%28 = load i32, i32* %11
%29 = add nsw i32 %28, 1
%30 = sext i32 %29 to i64
%31 = getelementptr inbounds float, float* %27, i64 %30
%32 = load float, float* %31
store float %32, float* %12
%33 = load float*, float** %6
%34 = load i32, i32* %11
%35 = sext i32 %34 to i64
%36 = getelementptr inbounds float, float* %33, i64 %35
%37 = load float, float* %36
store float %37, float* %13
%38 = load float, float* %12
%39 = load float, float* %13
%40 = fsub float %38, %39
store float %40, float* %14
%41 = load i32, i32* %11
%42 = add nsw i32 %41, 1
%43 = load i32, i32* %9
%44 = srem i32 %42, %43
%45 = icmp eq i32 %44, 0
br i1 %45, label %46, label %54
46: ; preds = %26
%47 = load float*, float** %6
%48 = load i32, i32* %11
%49 = add nsw i32 %48, 1
%50 = sext i32 %49 to i64
%51 = getelementptr inbounds float, float* %47, i64 %50
%52 = load float, float* %51
%53 = fneg float %52
store float %53, float* @unwrapPhase.offset
br label %77
54: ; preds = %26
%55 = load float, float* %14
%56 = load float, float* %10
%57 = fneg float %56
%58 = fcmp olt float %55, %57
br i1 %58, label %59, label %65
59: ; preds = %54
%60 = load i32, i32* @M_PI
%61 = mul nsw i32 2, %60
%62 = sitofp i32 %61 to float
%63 = load float, float* @unwrapPhase.offset
%64 = fadd float %63, %62
store float %64, float* @unwrapPhase.offset
br label %76
65: ; preds = %54
%66 = load float, float* %14
%67 = load float, float* %10
%68 = fcmp ogt float %66, %67
br i1 %68, label %69, label %75
69: ; preds = %65
%70 = load i32, i32* @M_PI
%71 = mul nsw i32 2, %70
%72 = sitofp i32 %71 to float
%73 = load float, float* @unwrapPhase.offset
%74 = fsub float %73, %72
store float %74, float* @unwrapPhase.offset
br label %75
75: ; preds = %69, %65
br label %76
76: ; preds = %75, %59
br label %77
77: ; preds = %76, %46
%78 = load float*, float** %6
%79 = load i32, i32* %11
%80 = add nsw i32 %79, 1
%81 = sext i32 %80 to i64
%82 = getelementptr inbounds float, float* %78, i64 %81
%83 = load float, float* %82
%84 = load float, float* @unwrapPhase.offset
%85 = fadd float %83, %84
%86 = load float*, float** %7
%87 = load i32, i32* %11
%88 = add nsw i32 %87, 1
%89 = sext i32 %88 to i64
%90 = getelementptr inbounds float, float* %86, i64 %89
store float %85, float* %90
br label %91
91: ; preds = %77
%92 = load i32, i32* %11
%93 = add nsw i32 %92, 1
store i32 %93, i32* %11
br label %21
94: ; preds = %21
ret void
}

LLVM IR (-O2):
@M_PI = global i32 0
define void @unwrapPhase(float* nocapture %0, float* nocapture writeonly %1, i32 %2, i32 %3, float %4) {
%6 = load float, float* %0
store float 0.000000e+00, float* %1
%7 = icmp sgt i32 %2, 1
br i1 %7, label %8, label %43
8: ; preds = %5
%9 = fneg float %6
%10 = add nsw i32 %2, -1
%11 = fneg float %4
%12 = load i32, i32* @M_PI
%13 = shl nsw i32 %12, 1
%14 = sitofp i32 %13 to float
%15 = zext i32 %10 to i64
br label %16
16: ; preds = %8, %38
%17 = phi i64 [ 0, %8 ], [ %19, %38 ]
%18 = phi float [ %9, %8 ], [ %39, %38 ]
%19 = add nuw nsw i64 %17, 1
%20 = getelementptr inbounds float, float* %0, i64 %19
%21 = load float, float* %20
%22 = getelementptr inbounds float, float* %0, i64 %17
%23 = load float, float* %22
%24 = fsub float %21, %23
%25 = trunc i64 %19 to i32
%26 = srem i32 %25, %3
%27 = icmp eq i32 %26, 0
br i1 %27, label %28, label %30
28: ; preds = %16
%29 = fneg float %21
br label %38
30: ; preds = %16
%31 = fcmp olt float %24, %11
br i1 %31, label %32, label %34
32: ; preds = %30
%33 = fadd float %18, %14
br label %38
34: ; preds = %30
%35 = fcmp ogt float %24, %4
br i1 %35, label %36, label %38
36: ; preds = %34
%37 = fsub float %18, %14
br label %38
38: ; preds = %32, %36, %28, %34
%39 = phi float [ %18, %34 ], [ %33, %32 ], [ %37, %36 ], [ %29, %28 ]
%40 = fadd float %21, %39
%41 = getelementptr inbounds float, float* %1, i64 %19
store float %40, float* %41
%42 = icmp eq i64 %19, %15
br i1 %42, label %43, label %16
43: ; preds = %38, %5
ret void
}