#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 */
double atan2(double, double);
double log(double);
double sqrt(double);
void gp_pt_alp(double *to, const double *p, int cnt, int siz, int rem) {
  double x, y;
  int i;
  if (siz == 1) {
    for (i = 0; i < cnt; i++, to += 8, p += 2)
      x = *p, y = p[1], to[3] = log(to[0] = sqrt(x * x + y * y)),
      to[6] = atan2(y, x);
    return;
  }
  double z, xa, ya, za, zm, zM, sz_1 = 1.0 / (double)siz;
  int j, i1 = cnt + !!rem;
  for (i = 0; i < i1; i++, to += 8, p++) {
    if (i == cnt)
      siz = rem, sz_1 = 1.0 / (double)(rem);
    for (j = 1, xa = *p, ya = *++p, zm = zM = za = xa * xa + ya * ya; j < siz;
         j++) {
      x = *++p, y = *++p, xa += x, ya += y;
      za += (z = x * x + y * y);
      if (z > zM)
        zM = z;
      else if (z < zm)
        zm = z;
    }
    to[3] = log(to[0] = sqrt(za * sz_1));
    to[4] = log(to[1] = sqrt(zm));
    to[6] = atan2(ya, xa);
    to[5] = log(to[2] = sqrt(zM));
  }
}


x86 assembly (clang, -O0):
<gp_pt_alp>:
push   %rbp
mov    %rsp,%rbp
sub    $0x80,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
mov    %r8d,-0x1c(%rbp)
cmpl   $0x1,-0x18(%rbp)
jne    11ff <gp_pt_alp+0xcf>
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    11fa <gp_pt_alp+0xca>
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0x10(%rbp),%rax
movsd  0x8(%rax),%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x28(%rbp),%xmm0
movsd  -0x28(%rbp),%xmm2
movsd  -0x30(%rbp),%xmm1
mulsd  -0x30(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
call   1050 <sqrt@plt>
mov    -0x8(%rbp),%rax
movsd  %xmm0,(%rax)
call   1040 <log@plt>
mov    -0x8(%rbp),%rax
movsd  %xmm0,0x18(%rax)
movsd  -0x30(%rbp),%xmm0
movsd  -0x28(%rbp),%xmm1
call   1030 <atan2@plt>
mov    -0x8(%rbp),%rax
movsd  %xmm0,0x30(%rax)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
mov    -0x8(%rbp),%rax
add    $0x40,%rax
mov    %rax,-0x8(%rbp)
mov    -0x10(%rbp),%rax
add    $0x10,%rax
mov    %rax,-0x10(%rbp)
jmp    115e <gp_pt_alp+0x2e>
jmp    1444 <gp_pt_alp+0x314>
cvtsi2sdl -0x18(%rbp),%xmm1
movsd  0xdf4(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    -0x14(%rbp),%eax
cmpl   $0x0,-0x1c(%rbp)
setne  %cl
xor    $0xff,%cl
xor    $0xff,%cl
and    $0x1,%cl
movzbl %cl,%ecx
add    %ecx,%eax
mov    %eax,-0x78(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x78(%rbp),%eax
jge    1444 <gp_pt_alp+0x314>
mov    -0x34(%rbp),%eax
cmp    -0x14(%rbp),%eax
jne    126b <gp_pt_alp+0x13b>
mov    -0x1c(%rbp),%eax
mov    %eax,-0x18(%rbp)
cvtsi2sdl -0x1c(%rbp),%xmm1
movsd  0xd9e(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x70(%rbp)
movl   $0x1,-0x74(%rbp)
mov    -0x10(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x10(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x10(%rbp)
movsd  0x8(%rax),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x48(%rbp),%xmm0
movsd  -0x48(%rbp),%xmm2
movsd  -0x50(%rbp),%xmm1
mulsd  -0x50(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  %xmm0,-0x68(%rbp)
movsd  %xmm0,-0x60(%rbp)
mov    -0x74(%rbp),%eax
cmp    -0x18(%rbp),%eax
jge    139f <gp_pt_alp+0x26f>
mov    -0x10(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x10(%rbp)
movsd  0x8(%rax),%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0x10(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x10(%rbp)
movsd  0x8(%rax),%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x28(%rbp),%xmm0
addsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x30(%rbp),%xmm0
addsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x28(%rbp),%xmm0
movsd  -0x28(%rbp),%xmm2
movsd  -0x30(%rbp),%xmm1
mulsd  -0x30(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
addsd  -0x58(%rbp),%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x40(%rbp),%xmm0
ucomisd -0x68(%rbp),%xmm0
jbe    1369 <gp_pt_alp+0x239>
movsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x68(%rbp)
jmp    138c <gp_pt_alp+0x25c>
movsd  -0x40(%rbp),%xmm1
movsd  -0x60(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    1387 <gp_pt_alp+0x257>
movsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x60(%rbp)
jmp    138c <gp_pt_alp+0x25c>
jmp    1391 <gp_pt_alp+0x261>
mov    -0x74(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x74(%rbp)
jmp    12c3 <gp_pt_alp+0x193>
movsd  -0x58(%rbp),%xmm0
mulsd  -0x70(%rbp),%xmm0
call   1050 <sqrt@plt>
mov    -0x8(%rbp),%rax
movsd  %xmm0,(%rax)
call   1040 <log@plt>
mov    -0x8(%rbp),%rax
movsd  %xmm0,0x18(%rax)
movsd  -0x60(%rbp),%xmm0
call   1050 <sqrt@plt>
mov    -0x8(%rbp),%rax
movsd  %xmm0,0x8(%rax)
call   1040 <log@plt>
mov    -0x8(%rbp),%rax
movsd  %xmm0,0x20(%rax)
movsd  -0x50(%rbp),%xmm0
movsd  -0x48(%rbp),%xmm1
call   1030 <atan2@plt>
mov    -0x8(%rbp),%rax
movsd  %xmm0,0x30(%rax)
movsd  -0x68(%rbp),%xmm0
call   1050 <sqrt@plt>
mov    -0x8(%rbp),%rax
movsd  %xmm0,0x10(%rax)
call   1040 <log@plt>
mov    -0x8(%rbp),%rax
movsd  %xmm0,0x28(%rax)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
mov    -0x8(%rbp),%rax
add    $0x40,%rax
mov    %rax,-0x8(%rbp)
mov    -0x10(%rbp),%rax
add    $0x8,%rax
mov    %rax,-0x10(%rbp)
jmp    1237 <gp_pt_alp+0x107>
add    $0x80,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<gp_pt_alp>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x48,%rsp
mov    %edx,%r14d
mov    %rsi,%rbp
mov    %rdi,%r13
cmp    $0x1,%ecx
jne    11d8 <gp_pt_alp+0xa8>
test   %r14d,%r14d
jle    149b <gp_pt_alp+0x36b>
xor    %ebx,%ebx
jmp    119d <gp_pt_alp+0x6d>
nopl   (%rax)
call   1050 <sqrt@plt>
movsd  %xmm0,0x0(%r13,%rbx,4)
call   1040 <log@plt>
movsd  %xmm0,0x18(%r13,%rbx,4)
movsd  0x10(%rsp),%xmm0
movsd  (%rsp),%xmm1
call   1030 <atan2@plt>
movsd  %xmm0,0x30(%r13,%rbx,4)
add    $0x10,%rbx
add    $0xffffffff,%r14d
je     149b <gp_pt_alp+0x36b>
movsd  0x0(%rbp,%rbx,1),%xmm2
movsd  0x8(%rbp,%rbx,1),%xmm3
movapd %xmm3,%xmm1
mulsd  %xmm3,%xmm1
movapd %xmm2,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
ucomisd 0xe43(%rip),%xmm0
movsd  %xmm2,(%rsp)
movsd  %xmm3,0x10(%rsp)
jb     1160 <gp_pt_alp+0x30>
sqrtsd %xmm0,%xmm0
jmp    1165 <gp_pt_alp+0x35>
cmp    $0x1,%r8d
mov    %r14d,%ebx
sbb    $0xffffffff,%ebx
test   %ebx,%ebx
jle    149b <gp_pt_alp+0x36b>
mov    %ecx,%r12d
cvtsi2sd %ecx,%xmm1
movsd  0xe07(%rip),%xmm3
movapd %xmm3,%xmm0
cvtsi2sd %r8d,%xmm2
divsd  %xmm1,%xmm0
divsd  %xmm2,%xmm3
movsd  %xmm3,0x20(%rsp)
xor    %r15d,%r15d
mov    %r8d,0xc(%rsp)
jmp    1255 <gp_pt_alp+0x125>
nopw   0x0(%rax,%rax,1)
call   1050 <sqrt@plt>
movsd  %xmm0,0x10(%r13)
call   1040 <log@plt>
movsd  %xmm0,0x28(%r13)
add    $0x1,%r15d
add    $0x40,%r13
add    $0x8,%rbp
movsd  (%rsp),%xmm0
cmp    %ebx,%r15d
mov    0xc(%rsp),%r8d
je     149b <gp_pt_alp+0x36b>
movsd  0x20(%rsp),%xmm8
cmp    %r14d,%r15d
cmove  %r8d,%r12d
je     126a <gp_pt_alp+0x13a>
movapd %xmm0,%xmm8
lea    0x8(%rbp),%rcx
movupd 0x0(%rbp),%xmm7
movapd %xmm7,%xmm0
mulpd  %xmm7,%xmm0
unpckhpd %xmm0,%xmm0
movapd %xmm7,%xmm6
mulsd  %xmm7,%xmm6
addsd  %xmm0,%xmm6
cmp    $0x2,%r12d
jl     12c0 <gp_pt_alp+0x190>
lea    -0x1(%r12),%eax
cmp    $0x2,%r12d
jne    1380 <gp_pt_alp+0x250>
movapd %xmm6,%xmm9
movapd %xmm6,%xmm0
test   $0x1,%al
je     12cc <gp_pt_alp+0x19c>
jmp    144f <gp_pt_alp+0x31f>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
mov    %rcx,%rbp
movapd %xmm6,%xmm0
movapd %xmm6,%xmm9
mulsd  %xmm8,%xmm0
ucomisd 0xd2f(%rip),%xmm0
movsd  %xmm8,(%rsp)
movapd %xmm7,0x10(%rsp)
movapd %xmm6,0x30(%rsp)
movsd  %xmm9,0x28(%rsp)
jb     1300 <gp_pt_alp+0x1d0>
sqrtsd %xmm0,%xmm0
jmp    1305 <gp_pt_alp+0x1d5>
nopw   0x0(%rax,%rax,1)
call   1050 <sqrt@plt>
movsd  %xmm0,0x0(%r13)
call   1040 <log@plt>
movsd  %xmm0,0x18(%r13)
movsd  0x28(%rsp),%xmm0
ucomisd 0xce4(%rip),%xmm0
jb     1330 <gp_pt_alp+0x200>
sqrtsd %xmm0,%xmm0
jmp    1335 <gp_pt_alp+0x205>
nopl   0x0(%rax)
call   1050 <sqrt@plt>
movsd  %xmm0,0x8(%r13)
call   1040 <log@plt>
movsd  %xmm0,0x20(%r13)
movapd 0x10(%rsp),%xmm1
movapd %xmm1,%xmm0
unpckhpd %xmm1,%xmm0
call   1030 <atan2@plt>
movsd  %xmm0,0x30(%r13)
movapd 0x30(%rsp),%xmm0
ucomisd 0xc9b(%rip),%xmm0
jb     1220 <gp_pt_alp+0xf0>
sqrtsd %xmm0,%xmm0
jmp    1225 <gp_pt_alp+0xf5>
nopl   0x0(%rax)
mov    %eax,%ecx
and    $0xfffffffe,%ecx
add    $0x20,%rbp
movapd %xmm6,%xmm9
movapd %xmm6,%xmm0
jmp    13c1 <gp_pt_alp+0x291>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movapd %xmm4,%xmm6
addpd  %xmm1,%xmm7
addsd  %xmm2,%xmm0
addpd  %xmm3,%xmm7
addsd  %xmm4,%xmm0
add    $0x20,%rbp
add    $0xfffffffe,%ecx
je     1440 <gp_pt_alp+0x310>
movupd -0x10(%rbp),%xmm1
movapd %xmm1,%xmm3
mulpd  %xmm1,%xmm3
unpckhpd %xmm3,%xmm3
movapd %xmm1,%xmm2
mulsd  %xmm1,%xmm2
addsd  %xmm3,%xmm2
ucomisd %xmm6,%xmm2
jbe    13f0 <gp_pt_alp+0x2c0>
movapd %xmm2,%xmm6
jmp    13fc <gp_pt_alp+0x2cc>
nopw   0x0(%rax,%rax,1)
ucomisd %xmm2,%xmm9
jbe    13fc <gp_pt_alp+0x2cc>
movapd %xmm2,%xmm9
movupd 0x0(%rbp),%xmm3
movapd %xmm3,%xmm5
mulpd  %xmm3,%xmm5
unpckhpd %xmm5,%xmm5
movapd %xmm3,%xmm4
mulsd  %xmm3,%xmm4
addsd  %xmm5,%xmm4
ucomisd %xmm6,%xmm4
ja     13a0 <gp_pt_alp+0x270>
ucomisd %xmm4,%xmm9
jbe    13a4 <gp_pt_alp+0x274>
movapd %xmm4,%xmm9
jmp    13a4 <gp_pt_alp+0x274>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
add    $0xffffffffffffffe8,%rbp
mov    %rbp,%rcx
test   $0x1,%al
je     12cc <gp_pt_alp+0x19c>
movupd 0x8(%rcx),%xmm1
add    $0x10,%rcx
addpd  %xmm1,%xmm7
movapd %xmm1,%xmm2
mulpd  %xmm1,%xmm2
unpckhpd %xmm2,%xmm2
mulsd  %xmm1,%xmm1
addsd  %xmm2,%xmm1
addsd  %xmm1,%xmm0
mov    %rcx,%rbp
ucomisd %xmm6,%xmm1
jbe    1486 <gp_pt_alp+0x356>
movapd %xmm1,%xmm6
jmp    12cc <gp_pt_alp+0x19c>
ucomisd %xmm1,%xmm9
jbe    12cc <gp_pt_alp+0x19c>
movapd %xmm1,%xmm9
jmp    12cc <gp_pt_alp+0x19c>
add    $0x48,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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