void liquid_unwrap_phase2(float *_theta, unsigned int _n) {
  unsigned int i;
  float dphi = 0.0f;
  for (i = 1; i < _n; i++)
    dphi += _theta[i] - _theta[i - 1];
  dphi /= (float)(_n - 1);
  for (i = 1; i < _n; i++) {
    while ((_theta[i] - _theta[i - 1]) >
           (3.14159265358979323846264338327950288f) + dphi)
      _theta[i] -= 2 * (3.14159265358979323846264338327950288f);
    while ((_theta[i] - _theta[i - 1]) <
           -(3.14159265358979323846264338327950288f) + dphi)
      _theta[i] += 2 * (3.14159265358979323846264338327950288f);
  }
}


x86 assembly (clang, -O0):
<liquid_unwrap_phase2>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x14(%rbp)
movl   $0x1,-0x10(%rbp)
mov    -0x10(%rbp),%eax
cmp    -0xc(%rbp),%eax
jae    115b <liquid_unwrap_phase2+0x5b>
mov    -0x8(%rbp),%rax
mov    -0x10(%rbp),%ecx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x10(%rbp),%ecx
sub    $0x1,%ecx
mov    %ecx,%ecx
subss  (%rax,%rcx,4),%xmm0
addss  -0x14(%rbp),%xmm0
movss  %xmm0,-0x14(%rbp)
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    111a <liquid_unwrap_phase2+0x1a>
mov    -0xc(%rbp),%eax
add    $0xffffffff,%eax
cvtsi2ss %rax,%xmm1
movss  -0x14(%rbp),%xmm0
divss  %xmm1,%xmm0
movss  %xmm0,-0x14(%rbp)
movl   $0x1,-0x10(%rbp)
mov    -0x10(%rbp),%eax
cmp    -0xc(%rbp),%eax
jae    124a <liquid_unwrap_phase2+0x14a>
jmp    118c <liquid_unwrap_phase2+0x8c>
mov    -0x8(%rbp),%rax
mov    -0x10(%rbp),%ecx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x10(%rbp),%ecx
sub    $0x1,%ecx
mov    %ecx,%ecx
subss  (%rax,%rcx,4),%xmm0
movss  0xe4f(%rip),%xmm1
addss  -0x14(%rbp),%xmm1
ucomiss %xmm1,%xmm0
jbe    11e1 <liquid_unwrap_phase2+0xe1>
mov    -0x8(%rbp),%rax
mov    -0x10(%rbp),%ecx
movss  (%rax,%rcx,4),%xmm0
movss  0xe35(%rip),%xmm1
subss  %xmm1,%xmm0
movss  %xmm0,(%rax,%rcx,4)
jmp    118c <liquid_unwrap_phase2+0x8c>
jmp    11e6 <liquid_unwrap_phase2+0xe6>
mov    -0x8(%rbp),%rax
mov    -0x10(%rbp),%ecx
movss  (%rax,%rcx,4),%xmm1
mov    -0x8(%rbp),%rax
mov    -0x10(%rbp),%ecx
sub    $0x1,%ecx
mov    %ecx,%ecx
subss  (%rax,%rcx,4),%xmm1
movss  0xdf9(%rip),%xmm0
addss  -0x14(%rbp),%xmm0
ucomiss %xmm1,%xmm0
jbe    1237 <liquid_unwrap_phase2+0x137>
mov    -0x8(%rbp),%rax
mov    -0x10(%rbp),%ecx
movss  0xde0(%rip),%xmm0
addss  (%rax,%rcx,4),%xmm0
movss  %xmm0,(%rax,%rcx,4)
jmp    11e6 <liquid_unwrap_phase2+0xe6>
jmp    123c <liquid_unwrap_phase2+0x13c>
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    117b <liquid_unwrap_phase2+0x7b>
pop    %rbp
ret

x86 assembly (clang, -O2):
<liquid_unwrap_phase2>:
cmp    $0x2,%esi
jb     1277 <liquid_unwrap_phase2+0x177>
mov    %esi,%r8d
movss  (%rdi),%xmm1
lea    -0x1(%r8),%rcx
lea    -0x2(%r8),%rdx
mov    %ecx,%eax
and    $0x3,%eax
cmp    $0x3,%rdx
jae    1135 <liquid_unwrap_phase2+0x35>
xorps  %xmm0,%xmm0
mov    $0x1,%edx
test   %rax,%rax
jne    1190 <liquid_unwrap_phase2+0x90>
jmp    11bc <liquid_unwrap_phase2+0xbc>
and    $0xfffffffffffffffc,%rcx
xorps  %xmm0,%xmm0
xor    %edx,%edx
xchg   %ax,%ax
movss  0x4(%rdi,%rdx,4),%xmm2
movss  0x8(%rdi,%rdx,4),%xmm3
movss  0xc(%rdi,%rdx,4),%xmm4
movaps %xmm4,%xmm5
subss  %xmm3,%xmm5
subss  %xmm2,%xmm3
subss  %xmm1,%xmm2
addss  %xmm0,%xmm2
addss  %xmm2,%xmm3
addss  %xmm3,%xmm5
movss  0x10(%rdi,%rdx,4),%xmm1
movaps %xmm1,%xmm0
subss  %xmm4,%xmm0
addss  %xmm5,%xmm0
add    $0x4,%rdx
cmp    %rdx,%rcx
jne    1140 <liquid_unwrap_phase2+0x40>
add    $0x1,%rdx
test   %rax,%rax
je     11bc <liquid_unwrap_phase2+0xbc>
lea    (%rdi,%rdx,4),%rcx
xor    %edx,%edx
cs nopw 0x0(%rax,%rax,1)
movss  (%rcx,%rdx,4),%xmm2
movaps %xmm2,%xmm3
subss  %xmm1,%xmm3
addss  %xmm3,%xmm0
add    $0x1,%rdx
movaps %xmm2,%xmm1
cmp    %rdx,%rax
jne    11a0 <liquid_unwrap_phase2+0xa0>
cmp    $0x2,%esi
jb     1277 <liquid_unwrap_phase2+0x177>
add    $0xffffffff,%esi
xorps  %xmm1,%xmm1
cvtsi2ss %rsi,%xmm1
divss  %xmm1,%xmm0
movss  0xe24(%rip),%xmm1
addss  %xmm0,%xmm1
addss  0xe1c(%rip),%xmm0
movss  (%rdi),%xmm5
mov    $0x1,%ecx
movss  0xe0f(%rip),%xmm2
movss  0xe0b(%rip),%xmm3
jmp    121c <liquid_unwrap_phase2+0x11c>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
add    $0x1,%rcx
movaps %xmm4,%xmm5
cmp    %r8,%rcx
je     1277 <liquid_unwrap_phase2+0x177>
movss  (%rdi,%rcx,4),%xmm4
movaps %xmm4,%xmm6
subss  %xmm5,%xmm6
ucomiss %xmm1,%xmm6
jbe    1245 <liquid_unwrap_phase2+0x145>
nopl   (%rax)
addss  %xmm2,%xmm4
movaps %xmm4,%xmm6
subss  %xmm5,%xmm6
ucomiss %xmm1,%xmm6
ja     1230 <liquid_unwrap_phase2+0x130>
movss  %xmm4,(%rdi,%rcx,4)
movss  -0x4(%rdi,%rcx,4),%xmm5
movaps %xmm4,%xmm6
subss  %xmm5,%xmm6
ucomiss %xmm6,%xmm0
jbe    1210 <liquid_unwrap_phase2+0x110>
nopw   0x0(%rax,%rax,1)
addss  %xmm3,%xmm4
movaps %xmm4,%xmm6
subss  %xmm5,%xmm6
ucomiss %xmm6,%xmm0
ja     1260 <liquid_unwrap_phase2+0x160>
movss  %xmm4,(%rdi,%rcx,4)
jmp    1210 <liquid_unwrap_phase2+0x110>
ret

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

LLVM IR (-O2):
define void @liquid_unwrap_phase2(float* nocapture %0, i32 %1) {
%3 = icmp ugt i32 %1, 1
br i1 %3, label %4, label %98
4: ; preds = %2
%5 = zext i32 %1 to i64
%6 = load float, float* %0
%7 = add nsw i64 %5, -1
%8 = add nsw i64 %5, -2
%9 = and i64 %7, 3
%10 = icmp ult i64 %8, 3
br i1 %10, label %40, label %11
11: ; preds = %4
%12 = and i64 %7, -4
br label %13
13: ; preds = %13, %11
%14 = phi float [ %6, %11 ], [ %34, %13 ]
%15 = phi i64 [ 1, %11 ], [ %37, %13 ]
%16 = phi float [ 0.000000e+00, %11 ], [ %36, %13 ]
%17 = phi i64 [ 0, %11 ], [ %38, %13 ]
%18 = getelementptr inbounds float, float* %0, i64 %15
%19 = load float, float* %18
%20 = fsub float %19, %14
%21 = fadd float %16, %20
%22 = add nuw nsw i64 %15, 1
%23 = getelementptr inbounds float, float* %0, i64 %22
%24 = load float, float* %23
%25 = fsub float %24, %19
%26 = fadd float %21, %25
%27 = add nuw nsw i64 %15, 2
%28 = getelementptr inbounds float, float* %0, i64 %27
%29 = load float, float* %28
%30 = fsub float %29, %24
%31 = fadd float %26, %30
%32 = add nuw nsw i64 %15, 3
%33 = getelementptr inbounds float, float* %0, i64 %32
%34 = load float, float* %33
%35 = fsub float %34, %29
%36 = fadd float %31, %35
%37 = add nuw nsw i64 %15, 4
%38 = add i64 %17, 4
%39 = icmp eq i64 %38, %12
br i1 %39, label %40, label %13
40: ; preds = %13, %4
%41 = phi float [ undef, %4 ], [ %36, %13 ]
%42 = phi float [ %6, %4 ], [ %34, %13 ]
%43 = phi i64 [ 1, %4 ], [ %37, %13 ]
%44 = phi float [ 0.000000e+00, %4 ], [ %36, %13 ]
%45 = icmp eq i64 %9, 0
br i1 %45, label %58, label %46
46: ; preds = %40, %46
%47 = phi float [ %52, %46 ], [ %42, %40 ]
%48 = phi i64 [ %55, %46 ], [ %43, %40 ]
%49 = phi float [ %54, %46 ], [ %44, %40 ]
%50 = phi i64 [ %56, %46 ], [ 0, %40 ]
%51 = getelementptr inbounds float, float* %0, i64 %48
%52 = load float, float* %51
%53 = fsub float %52, %47
%54 = fadd float %49, %53
%55 = add nuw nsw i64 %48, 1
%56 = add i64 %50, 1
%57 = icmp eq i64 %56, %9
br i1 %57, label %58, label %46
58: ; preds = %46, %40
%59 = phi float [ %41, %40 ], [ %54, %46 ]
%60 = add i32 %1, -1
%61 = uitofp i32 %60 to float
%62 = fdiv float %59, %61
br i1 %3, label %63, label %98
63: ; preds = %58
%64 = fadd float %62, 0x400921FB60000000
%65 = fadd float %62, 0xC00921FB60000000
%66 = zext i32 %1 to i64
%67 = load float, float* %0
br label %68
68: ; preds = %63, %94
%69 = phi float [ %67, %63 ], [ %95, %94 ]
%70 = phi i64 [ 1, %63 ], [ %96, %94 ]
%71 = getelementptr inbounds float, float* %0, i64 %70
%72 = add nsw i64 %70, -1
%73 = getelementptr inbounds float, float* %0, i64 %72
%74 = load float, float* %71
%75 = fsub float %74, %69
%76 = fcmp ogt float %75, %64
br i1 %76, label %83, label %78
77: ; preds = %83
store float %85, float* %71
br label %78
78: ; preds = %77, %68
%79 = phi float [ %85, %77 ], [ %74, %68 ]
%80 = load float, float* %73
%81 = fsub float %79, %80
%82 = fcmp olt float %81, %65
br i1 %82, label %88, label %94
83: ; preds = %68, %83
%84 = phi float [ %85, %83 ], [ %74, %68 ]
%85 = fadd float %84, 0xC01921FB60000000
%86 = fsub float %85, %69
%87 = fcmp ogt float %86, %64
br i1 %87, label %83, label %77
88: ; preds = %78, %88
%89 = phi float [ %90, %88 ], [ %79, %78 ]
%90 = fadd float %89, 0x401921FB60000000
%91 = fsub float %90, %80
%92 = fcmp olt float %91, %65
br i1 %92, label %88, label %93
93: ; preds = %88
store float %90, float* %71
br label %94
94: ; preds = %93, %78
%95 = phi float [ %79, %78 ], [ %90, %93 ]
%96 = add nuw nsw i64 %70, 1
%97 = icmp eq i64 %96, %66
br i1 %97, label %98, label %68
98: ; preds = %94, %2, %58
ret void
}