#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 asin(double);
double sin(double);
double sqrt(double);
void gcpt(double *pxn, int kig1s, int kig1e, int kig2s, int kig2e, int knd,
          int kjd, double pgamma, int ki1, int kj1, int ki2, int kj2, int ki,
          int kj) {
  int id1, id2, id3, ioffset;
  double r1, r2, r3;
  double zbeta, zalpha, zchord, ztheta;
  id1 = kig1e - kig1s + 1;
  id2 = id1 * (kig2e - kig2s + 1);
  id3 = id2 * 3;
  ioffset = -(id1 + id2 + id3);
  r1 = (pxn[ki2 + id1 * kj2 + id2 * 1 + id3 * kjd + ioffset] -
        pxn[ki1 + id1 * kj1 + id2 * 1 + id3 * kjd + ioffset]);
  r2 = (pxn[ki2 + id1 * kj2 + id2 * 2 + id3 * kjd + ioffset] -
        pxn[ki1 + id1 * kj1 + id2 * 2 + id3 * kjd + ioffset]);
  r3 = (pxn[ki2 + id1 * kj2 + id2 * 3 + id3 * kjd + ioffset] -
        pxn[ki1 + id1 * kj1 + id2 * 3 + id3 * kjd + ioffset]);
  zchord = sqrt((r1 * r1) + (r2 * r2) + (r3 * r3));
  ztheta = 2.0 * asin(zchord * 0.5);
  zbeta = sin(pgamma * ztheta) / sin(ztheta);
  zalpha = sin((1.0 - pgamma) * ztheta) / sin(ztheta);
  pxn[ki + id1 * kj + id2 * 1 + id3 * kjd + ioffset] =
      zalpha * pxn[ki1 + id1 * kj1 + id2 * 1 + id3 * kjd + ioffset] +
      zbeta * pxn[ki2 + id1 * kj2 + id2 * 1 + id3 * kjd + ioffset];
  pxn[ki + id1 * kj + id2 * 2 + id3 * kjd + ioffset] =
      zalpha * pxn[ki1 + id1 * kj1 + id2 * 2 + id3 * kjd + ioffset] +
      zbeta * pxn[ki2 + id1 * kj2 + id2 * 2 + id3 * kjd + ioffset];
  pxn[ki + id1 * kj + id2 * 3 + id3 * kjd + ioffset] =
      zalpha * pxn[ki1 + id1 * kj1 + id2 * 3 + id3 * kjd + ioffset] +
      zbeta * pxn[ki2 + id1 * kj2 + id2 * 3 + id3 * kjd + ioffset];
  return;
}


x86 assembly (clang, -O0):
<gcpt>:
push   %rbp
mov    %rsp,%rbp
sub    $0x80,%rsp
mov    0x40(%rbp),%eax
mov    0x38(%rbp),%eax
mov    0x30(%rbp),%eax
mov    0x28(%rbp),%eax
mov    0x20(%rbp),%eax
mov    0x18(%rbp),%eax
mov    0x10(%rbp),%eax
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
mov    %ecx,-0x14(%rbp)
mov    %r8d,-0x18(%rbp)
mov    %r9d,-0x1c(%rbp)
movsd  %xmm0,-0x28(%rbp)
mov    -0x10(%rbp),%eax
sub    -0xc(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
mov    -0x18(%rbp),%ecx
sub    -0x14(%rbp),%ecx
add    $0x1,%ecx
imul   %ecx,%eax
mov    %eax,-0x30(%rbp)
imul   $0x3,-0x30(%rbp),%eax
mov    %eax,-0x34(%rbp)
mov    -0x2c(%rbp),%ecx
add    -0x30(%rbp),%ecx
add    -0x34(%rbp),%ecx
xor    %eax,%eax
sub    %ecx,%eax
mov    %eax,-0x38(%rbp)
mov    -0x8(%rbp),%rax
mov    0x28(%rbp),%ecx
mov    -0x2c(%rbp),%edx
imul   0x30(%rbp),%edx
add    %edx,%ecx
mov    -0x30(%rbp),%edx
shl    $0x0,%edx
add    %edx,%ecx
mov    -0x34(%rbp),%edx
imul   0x10(%rbp),%edx
add    %edx,%ecx
add    -0x38(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    0x18(%rbp),%ecx
mov    -0x2c(%rbp),%edx
imul   0x20(%rbp),%edx
add    %edx,%ecx
mov    -0x30(%rbp),%edx
shl    $0x0,%edx
add    %edx,%ecx
mov    -0x34(%rbp),%edx
imul   0x10(%rbp),%edx
add    %edx,%ecx
add    -0x38(%rbp),%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x8(%rbp),%rax
mov    0x28(%rbp),%ecx
mov    -0x2c(%rbp),%edx
imul   0x30(%rbp),%edx
add    %edx,%ecx
mov    -0x30(%rbp),%edx
shl    $0x1,%edx
add    %edx,%ecx
mov    -0x34(%rbp),%edx
imul   0x10(%rbp),%edx
add    %edx,%ecx
add    -0x38(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    0x18(%rbp),%ecx
mov    -0x2c(%rbp),%edx
imul   0x20(%rbp),%edx
add    %edx,%ecx
mov    -0x30(%rbp),%edx
shl    $0x1,%edx
add    %edx,%ecx
mov    -0x34(%rbp),%edx
imul   0x10(%rbp),%edx
add    %edx,%ecx
add    -0x38(%rbp),%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x8(%rbp),%rax
mov    0x28(%rbp),%ecx
mov    -0x2c(%rbp),%edx
imul   0x30(%rbp),%edx
add    %edx,%ecx
imul   $0x3,-0x30(%rbp),%edx
add    %edx,%ecx
mov    -0x34(%rbp),%edx
imul   0x10(%rbp),%edx
add    %edx,%ecx
add    -0x38(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    0x18(%rbp),%ecx
mov    -0x2c(%rbp),%edx
imul   0x20(%rbp),%edx
add    %edx,%ecx
imul   $0x3,-0x30(%rbp),%edx
add    %edx,%ecx
mov    -0x34(%rbp),%edx
imul   0x10(%rbp),%edx
add    %edx,%ecx
add    -0x38(%rbp),%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x40(%rbp),%xmm1
movsd  -0x40(%rbp),%xmm2
movsd  -0x48(%rbp),%xmm0
mulsd  -0x48(%rbp),%xmm0
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
movsd  -0x50(%rbp),%xmm0
movsd  -0x50(%rbp),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
call   1050 <sqrt@plt>
movsd  %xmm0,-0x68(%rbp)
movsd  0xd1e(%rip),%xmm0
mulsd  -0x68(%rbp),%xmm0
call   1040 <asin@plt>
movaps %xmm0,%xmm1
movsd  0xd01(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x28(%rbp),%xmm0
mulsd  -0x70(%rbp),%xmm0
call   1030 <sin@plt>
movsd  %xmm0,-0x80(%rbp)
movsd  -0x70(%rbp),%xmm0
call   1030 <sin@plt>
movaps %xmm0,%xmm1
movsd  -0x80(%rbp),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  0xcb9(%rip),%xmm0
subsd  -0x28(%rbp),%xmm0
mulsd  -0x70(%rbp),%xmm0
call   1030 <sin@plt>
movsd  %xmm0,-0x78(%rbp)
movsd  -0x70(%rbp),%xmm0
call   1030 <sin@plt>
movaps %xmm0,%xmm1
movsd  -0x78(%rbp),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  -0x60(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    0x18(%rbp),%ecx
mov    -0x2c(%rbp),%edx
imul   0x20(%rbp),%edx
add    %edx,%ecx
mov    -0x30(%rbp),%edx
shl    $0x0,%edx
add    %edx,%ecx
mov    -0x34(%rbp),%edx
imul   0x10(%rbp),%edx
add    %edx,%ecx
add    -0x38(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x58(%rbp),%xmm1
mov    -0x8(%rbp),%rax
mov    0x28(%rbp),%ecx
mov    -0x2c(%rbp),%edx
imul   0x30(%rbp),%edx
add    %edx,%ecx
mov    -0x30(%rbp),%edx
shl    $0x0,%edx
add    %edx,%ecx
mov    -0x34(%rbp),%edx
imul   0x10(%rbp),%edx
add    %edx,%ecx
add    -0x38(%rbp),%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
mov    0x38(%rbp),%ecx
mov    -0x2c(%rbp),%edx
imul   0x40(%rbp),%edx
add    %edx,%ecx
mov    -0x30(%rbp),%edx
shl    $0x0,%edx
add    %edx,%ecx
mov    -0x34(%rbp),%edx
imul   0x10(%rbp),%edx
add    %edx,%ecx
add    -0x38(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x60(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    0x18(%rbp),%ecx
mov    -0x2c(%rbp),%edx
imul   0x20(%rbp),%edx
add    %edx,%ecx
mov    -0x30(%rbp),%edx
shl    $0x1,%edx
add    %edx,%ecx
mov    -0x34(%rbp),%edx
imul   0x10(%rbp),%edx
add    %edx,%ecx
add    -0x38(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x58(%rbp),%xmm1
mov    -0x8(%rbp),%rax
mov    0x28(%rbp),%ecx
mov    -0x2c(%rbp),%edx
imul   0x30(%rbp),%edx
add    %edx,%ecx
mov    -0x30(%rbp),%edx
shl    $0x1,%edx
add    %edx,%ecx
mov    -0x34(%rbp),%edx
imul   0x10(%rbp),%edx
add    %edx,%ecx
add    -0x38(%rbp),%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
mov    0x38(%rbp),%ecx
mov    -0x2c(%rbp),%edx
imul   0x40(%rbp),%edx
add    %edx,%ecx
mov    -0x30(%rbp),%edx
shl    $0x1,%edx
add    %edx,%ecx
mov    -0x34(%rbp),%edx
imul   0x10(%rbp),%edx
add    %edx,%ecx
add    -0x38(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x60(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    0x18(%rbp),%ecx
mov    -0x2c(%rbp),%edx
imul   0x20(%rbp),%edx
add    %edx,%ecx
imul   $0x3,-0x30(%rbp),%edx
add    %edx,%ecx
mov    -0x34(%rbp),%edx
imul   0x10(%rbp),%edx
add    %edx,%ecx
add    -0x38(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x58(%rbp),%xmm1
mov    -0x8(%rbp),%rax
mov    0x28(%rbp),%ecx
mov    -0x2c(%rbp),%edx
imul   0x30(%rbp),%edx
add    %edx,%ecx
imul   $0x3,-0x30(%rbp),%edx
add    %edx,%ecx
mov    -0x34(%rbp),%edx
imul   0x10(%rbp),%edx
add    %edx,%ecx
add    -0x38(%rbp),%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
mov    0x38(%rbp),%ecx
mov    -0x2c(%rbp),%edx
imul   0x40(%rbp),%edx
add    %edx,%ecx
imul   $0x3,-0x30(%rbp),%edx
add    %edx,%ecx
mov    -0x34(%rbp),%edx
imul   0x10(%rbp),%edx
add    %edx,%ecx
add    -0x38(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
add    $0x80,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<gcpt>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x48,%rsp
movsd  %xmm0,0x8(%rsp)
mov    %edx,%ebx
mov    %rdi,%r14
sub    %esi,%ebx
mov    %ebx,%edx
not    %edx
add    $0x1,%ebx
sub    %ecx,%r8d
add    $0x1,%r8d
imul   %ebx,%r8d
lea    (%r8,%r8,2),%eax
lea    0x0(,%r8,4),%ecx
sub    %ecx,%edx
mov    0xa0(%rsp),%ecx
imul   %ebx,%ecx
add    0x98(%rsp),%ecx
mov    0x80(%rsp),%r15d
imul   %eax,%r15d
add    %edx,%r15d
lea    (%r15,%r8,2),%r13d
mov    %r8d,%ebp
add    %r15d,%ebp
lea    (%rcx,%rbp,1),%edx
movslq %edx,%r12
movsd  (%rdi,%r12,8),%xmm1
mov    0x90(%rsp),%edx
imul   %ebx,%edx
add    0x88(%rsp),%edx
lea    (%rdx,%rbp,1),%esi
movslq %esi,%rsi
mov    %rsi,0x20(%rsp)
subsd  (%rdi,%rsi,8),%xmm1
lea    0x0(%r13,%rcx,1),%esi
movslq %esi,%rsi
mov    %rsi,0x28(%rsp)
movsd  (%rdi,%rsi,8),%xmm2
lea    0x0(%r13,%rdx,1),%esi
movslq %esi,%rsi
mov    %rsi,0x30(%rsp)
subsd  (%rdi,%rsi,8),%xmm2
add    %eax,%r15d
add    %r15d,%ecx
movslq %ecx,%rax
mov    %rax,0x40(%rsp)
movsd  (%rdi,%rax,8),%xmm0
add    %r15d,%edx
movslq %edx,%rax
mov    %rax,0x38(%rsp)
subsd  (%rdi,%rax,8),%xmm0
mulsd  %xmm2,%xmm2
mulsd  %xmm1,%xmm1
addsd  %xmm2,%xmm1
mulsd  %xmm0,%xmm0
addsd  %xmm1,%xmm0
xorpd  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jb     122d <gcpt+0xfd>
sqrtsd %xmm0,%xmm0
jmp    1232 <gcpt+0x102>
call   1050 <sqrt@plt>
mulsd  0xdc6(%rip),%xmm0
call   1040 <asin@plt>
addsd  %xmm0,%xmm0
movsd  %xmm0,0x10(%rsp)
mulsd  0x8(%rsp),%xmm0
call   1030 <sin@plt>
movsd  %xmm0,0x18(%rsp)
movsd  0x10(%rsp),%xmm0
call   1030 <sin@plt>
movsd  0x18(%rsp),%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,0x18(%rsp)
movsd  0xd8b(%rip),%xmm0
subsd  0x8(%rsp),%xmm0
mulsd  0x10(%rsp),%xmm0
call   1030 <sin@plt>
movsd  %xmm0,0x8(%rsp)
movsd  0x10(%rsp),%xmm0
call   1030 <sin@plt>
movsd  0x8(%rsp),%xmm3
divsd  %xmm0,%xmm3
movsd  (%r14,%r12,8),%xmm0
movsd  0x18(%rsp),%xmm2
mulsd  %xmm2,%xmm0
mov    0x20(%rsp),%rax
movsd  (%r14,%rax,8),%xmm1
mulsd  %xmm3,%xmm1
addsd  %xmm0,%xmm1
imul   0xb0(%rsp),%ebx
add    0xa8(%rsp),%ebx
add    %ebx,%ebp
movslq %ebp,%rax
movsd  %xmm1,(%r14,%rax,8)
mov    0x28(%rsp),%rax
movsd  (%r14,%rax,8),%xmm0
mulsd  %xmm2,%xmm0
mov    0x30(%rsp),%rax
movsd  (%r14,%rax,8),%xmm1
mulsd  %xmm3,%xmm1
addsd  %xmm0,%xmm1
add    %ebx,%r13d
movslq %r13d,%rax
movsd  %xmm1,(%r14,%rax,8)
mov    0x40(%rsp),%rax
movapd %xmm2,%xmm0
mulsd  (%r14,%rax,8),%xmm0
mov    0x38(%rsp),%rax
movapd %xmm3,%xmm1
mulsd  (%r14,%rax,8),%xmm1
addsd  %xmm0,%xmm1
add    %ebx,%r15d
movslq %r15d,%rax
movsd  %xmm1,(%r14,%rax,8)
add    $0x48,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @gcpt(double* %0, i32 %1, i32 %2, i32 %3, i32 %4, i32 %5, i32 %6, double %7, i32 %8, i32 %9, i32 %10, i32 %11, i32 %12, i32 %13) {
%15 = alloca double*
%16 = alloca i32
%17 = alloca i32
%18 = alloca i32
%19 = alloca i32
%20 = alloca i32
%21 = alloca i32
%22 = alloca double
%23 = alloca i32
%24 = alloca i32
%25 = alloca i32
%26 = alloca i32
%27 = alloca i32
%28 = alloca i32
%29 = alloca i32
%30 = alloca i32
%31 = alloca i32
%32 = alloca i32
%33 = alloca double
%34 = alloca double
%35 = alloca double
%36 = alloca double
%37 = alloca double
%38 = alloca double
%39 = alloca double
store double* %0, double** %15
store i32 %1, i32* %16
store i32 %2, i32* %17
store i32 %3, i32* %18
store i32 %4, i32* %19
store i32 %5, i32* %20
store i32 %6, i32* %21
store double %7, double* %22
store i32 %8, i32* %23
store i32 %9, i32* %24
store i32 %10, i32* %25
store i32 %11, i32* %26
store i32 %12, i32* %27
store i32 %13, i32* %28
%40 = load i32, i32* %17
%41 = load i32, i32* %16
%42 = sub nsw i32 %40, %41
%43 = add nsw i32 %42, 1
store i32 %43, i32* %29
%44 = load i32, i32* %29
%45 = load i32, i32* %19
%46 = load i32, i32* %18
%47 = sub nsw i32 %45, %46
%48 = add nsw i32 %47, 1
%49 = mul nsw i32 %44, %48
store i32 %49, i32* %30
%50 = load i32, i32* %30
%51 = mul nsw i32 %50, 3
store i32 %51, i32* %31
%52 = load i32, i32* %29
%53 = load i32, i32* %30
%54 = add nsw i32 %52, %53
%55 = load i32, i32* %31
%56 = add nsw i32 %54, %55
%57 = sub nsw i32 0, %56
store i32 %57, i32* %32
%58 = load double*, double** %15
%59 = load i32, i32* %25
%60 = load i32, i32* %29
%61 = load i32, i32* %26
%62 = mul nsw i32 %60, %61
%63 = add nsw i32 %59, %62
%64 = load i32, i32* %30
%65 = mul nsw i32 %64, 1
%66 = add nsw i32 %63, %65
%67 = load i32, i32* %31
%68 = load i32, i32* %21
%69 = mul nsw i32 %67, %68
%70 = add nsw i32 %66, %69
%71 = load i32, i32* %32
%72 = add nsw i32 %70, %71
%73 = sext i32 %72 to i64
%74 = getelementptr inbounds double, double* %58, i64 %73
%75 = load double, double* %74
%76 = load double*, double** %15
%77 = load i32, i32* %23
%78 = load i32, i32* %29
%79 = load i32, i32* %24
%80 = mul nsw i32 %78, %79
%81 = add nsw i32 %77, %80
%82 = load i32, i32* %30
%83 = mul nsw i32 %82, 1
%84 = add nsw i32 %81, %83
%85 = load i32, i32* %31
%86 = load i32, i32* %21
%87 = mul nsw i32 %85, %86
%88 = add nsw i32 %84, %87
%89 = load i32, i32* %32
%90 = add nsw i32 %88, %89
%91 = sext i32 %90 to i64
%92 = getelementptr inbounds double, double* %76, i64 %91
%93 = load double, double* %92
%94 = fsub double %75, %93
store double %94, double* %33
%95 = load double*, double** %15
%96 = load i32, i32* %25
%97 = load i32, i32* %29
%98 = load i32, i32* %26
%99 = mul nsw i32 %97, %98
%100 = add nsw i32 %96, %99
%101 = load i32, i32* %30
%102 = mul nsw i32 %101, 2
%103 = add nsw i32 %100, %102
%104 = load i32, i32* %31
%105 = load i32, i32* %21
%106 = mul nsw i32 %104, %105
%107 = add nsw i32 %103, %106
%108 = load i32, i32* %32
%109 = add nsw i32 %107, %108
%110 = sext i32 %109 to i64
%111 = getelementptr inbounds double, double* %95, i64 %110
%112 = load double, double* %111
%113 = load double*, double** %15
%114 = load i32, i32* %23
%115 = load i32, i32* %29
%116 = load i32, i32* %24
%117 = mul nsw i32 %115, %116
%118 = add nsw i32 %114, %117
%119 = load i32, i32* %30
%120 = mul nsw i32 %119, 2
%121 = add nsw i32 %118, %120
%122 = load i32, i32* %31
%123 = load i32, i32* %21
%124 = mul nsw i32 %122, %123
%125 = add nsw i32 %121, %124
%126 = load i32, i32* %32
%127 = add nsw i32 %125, %126
%128 = sext i32 %127 to i64
%129 = getelementptr inbounds double, double* %113, i64 %128
%130 = load double, double* %129
%131 = fsub double %112, %130
store double %131, double* %34
%132 = load double*, double** %15
%133 = load i32, i32* %25
%134 = load i32, i32* %29
%135 = load i32, i32* %26
%136 = mul nsw i32 %134, %135
%137 = add nsw i32 %133, %136
%138 = load i32, i32* %30
%139 = mul nsw i32 %138, 3
%140 = add nsw i32 %137, %139
%141 = load i32, i32* %31
%142 = load i32, i32* %21
%143 = mul nsw i32 %141, %142
%144 = add nsw i32 %140, %143
%145 = load i32, i32* %32
%146 = add nsw i32 %144, %145
%147 = sext i32 %146 to i64
%148 = getelementptr inbounds double, double* %132, i64 %147
%149 = load double, double* %148
%150 = load double*, double** %15
%151 = load i32, i32* %23
%152 = load i32, i32* %29
%153 = load i32, i32* %24
%154 = mul nsw i32 %152, %153
%155 = add nsw i32 %151, %154
%156 = load i32, i32* %30
%157 = mul nsw i32 %156, 3
%158 = add nsw i32 %155, %157
%159 = load i32, i32* %31
%160 = load i32, i32* %21
%161 = mul nsw i32 %159, %160
%162 = add nsw i32 %158, %161
%163 = load i32, i32* %32
%164 = add nsw i32 %162, %163
%165 = sext i32 %164 to i64
%166 = getelementptr inbounds double, double* %150, i64 %165
%167 = load double, double* %166
%168 = fsub double %149, %167
store double %168, double* %35
%169 = load double, double* %33
%170 = load double, double* %33
%171 = load double, double* %34
%172 = load double, double* %34
%173 = fmul double %171, %172
%174 = call double @llvm.fmuladd.f64(double %169, double %170, double %173)
%175 = load double, double* %35
%176 = load double, double* %35
%177 = call double @llvm.fmuladd.f64(double %175, double %176, double %174)
%178 = call double @sqrt(double %177)
store double %178, double* %38
%179 = load double, double* %38
%180 = fmul double %179, 5.000000e-01
%181 = call double @asin(double %180)
%182 = fmul double 2.000000e+00, %181
store double %182, double* %39
%183 = load double, double* %22
%184 = load double, double* %39
%185 = fmul double %183, %184
%186 = call double @sin(double %185)
%187 = load double, double* %39
%188 = call double @sin(double %187)
%189 = fdiv double %186, %188
store double %189, double* %36
%190 = load double, double* %22
%191 = fsub double 1.000000e+00, %190
%192 = load double, double* %39
%193 = fmul double %191, %192
%194 = call double @sin(double %193)
%195 = load double, double* %39
%196 = call double @sin(double %195)
%197 = fdiv double %194, %196
store double %197, double* %37
%198 = load double, double* %37
%199 = load double*, double** %15
%200 = load i32, i32* %23
%201 = load i32, i32* %29
%202 = load i32, i32* %24
%203 = mul nsw i32 %201, %202
%204 = add nsw i32 %200, %203
%205 = load i32, i32* %30
%206 = mul nsw i32 %205, 1
%207 = add nsw i32 %204, %206
%208 = load i32, i32* %31
%209 = load i32, i32* %21
%210 = mul nsw i32 %208, %209
%211 = add nsw i32 %207, %210
%212 = load i32, i32* %32
%213 = add nsw i32 %211, %212
%214 = sext i32 %213 to i64
%215 = getelementptr inbounds double, double* %199, i64 %214
%216 = load double, double* %215
%217 = load double, double* %36
%218 = load double*, double** %15
%219 = load i32, i32* %25
%220 = load i32, i32* %29
%221 = load i32, i32* %26
%222 = mul nsw i32 %220, %221
%223 = add nsw i32 %219, %222
%224 = load i32, i32* %30
%225 = mul nsw i32 %224, 1
%226 = add nsw i32 %223, %225
%227 = load i32, i32* %31
%228 = load i32, i32* %21
%229 = mul nsw i32 %227, %228
%230 = add nsw i32 %226, %229
%231 = load i32, i32* %32
%232 = add nsw i32 %230, %231
%233 = sext i32 %232 to i64
%234 = getelementptr inbounds double, double* %218, i64 %233
%235 = load double, double* %234
%236 = fmul double %217, %235
%237 = call double @llvm.fmuladd.f64(double %198, double %216, double %236)
%238 = load double*, double** %15
%239 = load i32, i32* %27
%240 = load i32, i32* %29
%241 = load i32, i32* %28
%242 = mul nsw i32 %240, %241
%243 = add nsw i32 %239, %242
%244 = load i32, i32* %30
%245 = mul nsw i32 %244, 1
%246 = add nsw i32 %243, %245
%247 = load i32, i32* %31
%248 = load i32, i32* %21
%249 = mul nsw i32 %247, %248
%250 = add nsw i32 %246, %249
%251 = load i32, i32* %32
%252 = add nsw i32 %250, %251
%253 = sext i32 %252 to i64
%254 = getelementptr inbounds double, double* %238, i64 %253
store double %237, double* %254
%255 = load double, double* %37
%256 = load double*, double** %15
%257 = load i32, i32* %23
%258 = load i32, i32* %29
%259 = load i32, i32* %24
%260 = mul nsw i32 %258, %259
%261 = add nsw i32 %257, %260
%262 = load i32, i32* %30
%263 = mul nsw i32 %262, 2
%264 = add nsw i32 %261, %263
%265 = load i32, i32* %31
%266 = load i32, i32* %21
%267 = mul nsw i32 %265, %266
%268 = add nsw i32 %264, %267
%269 = load i32, i32* %32
%270 = add nsw i32 %268, %269
%271 = sext i32 %270 to i64
%272 = getelementptr inbounds double, double* %256, i64 %271
%273 = load double, double* %272
%274 = load double, double* %36
%275 = load double*, double** %15
%276 = load i32, i32* %25
%277 = load i32, i32* %29
%278 = load i32, i32* %26
%279 = mul nsw i32 %277, %278
%280 = add nsw i32 %276, %279
%281 = load i32, i32* %30
%282 = mul nsw i32 %281, 2
%283 = add nsw i32 %280, %282
%284 = load i32, i32* %31
%285 = load i32, i32* %21
%286 = mul nsw i32 %284, %285
%287 = add nsw i32 %283, %286
%288 = load i32, i32* %32
%289 = add nsw i32 %287, %288
%290 = sext i32 %289 to i64
%291 = getelementptr inbounds double, double* %275, i64 %290
%292 = load double, double* %291
%293 = fmul double %274, %292
%294 = call double @llvm.fmuladd.f64(double %255, double %273, double %293)
%295 = load double*, double** %15
%296 = load i32, i32* %27
%297 = load i32, i32* %29
%298 = load i32, i32* %28
%299 = mul nsw i32 %297, %298
%300 = add nsw i32 %296, %299
%301 = load i32, i32* %30
%302 = mul nsw i32 %301, 2
%303 = add nsw i32 %300, %302
%304 = load i32, i32* %31
%305 = load i32, i32* %21
%306 = mul nsw i32 %304, %305
%307 = add nsw i32 %303, %306
%308 = load i32, i32* %32
%309 = add nsw i32 %307, %308
%310 = sext i32 %309 to i64
%311 = getelementptr inbounds double, double* %295, i64 %310
store double %294, double* %311
%312 = load double, double* %37
%313 = load double*, double** %15
%314 = load i32, i32* %23
%315 = load i32, i32* %29
%316 = load i32, i32* %24
%317 = mul nsw i32 %315, %316
%318 = add nsw i32 %314, %317
%319 = load i32, i32* %30
%320 = mul nsw i32 %319, 3
%321 = add nsw i32 %318, %320
%322 = load i32, i32* %31
%323 = load i32, i32* %21
%324 = mul nsw i32 %322, %323
%325 = add nsw i32 %321, %324
%326 = load i32, i32* %32
%327 = add nsw i32 %325, %326
%328 = sext i32 %327 to i64
%329 = getelementptr inbounds double, double* %313, i64 %328
%330 = load double, double* %329
%331 = load double, double* %36
%332 = load double*, double** %15
%333 = load i32, i32* %25
%334 = load i32, i32* %29
%335 = load i32, i32* %26
%336 = mul nsw i32 %334, %335
%337 = add nsw i32 %333, %336
%338 = load i32, i32* %30
%339 = mul nsw i32 %338, 3
%340 = add nsw i32 %337, %339
%341 = load i32, i32* %31
%342 = load i32, i32* %21
%343 = mul nsw i32 %341, %342
%344 = add nsw i32 %340, %343
%345 = load i32, i32* %32
%346 = add nsw i32 %344, %345
%347 = sext i32 %346 to i64
%348 = getelementptr inbounds double, double* %332, i64 %347
%349 = load double, double* %348
%350 = fmul double %331, %349
%351 = call double @llvm.fmuladd.f64(double %312, double %330, double %350)
%352 = load double*, double** %15
%353 = load i32, i32* %27
%354 = load i32, i32* %29
%355 = load i32, i32* %28
%356 = mul nsw i32 %354, %355
%357 = add nsw i32 %353, %356
%358 = load i32, i32* %30
%359 = mul nsw i32 %358, 3
%360 = add nsw i32 %357, %359
%361 = load i32, i32* %31
%362 = load i32, i32* %21
%363 = mul nsw i32 %361, %362
%364 = add nsw i32 %360, %363
%365 = load i32, i32* %32
%366 = add nsw i32 %364, %365
%367 = sext i32 %366 to i64
%368 = getelementptr inbounds double, double* %352, i64 %367
store double %351, double* %368
ret void
}

LLVM IR (-O2):
define void @gcpt(double* nocapture %0, i32 %1, i32 %2, i32 %3, i32 %4, i32 %5, i32 %6, double %7, i32 %8, i32 %9, i32 %10, i32 %11, i32 %12, i32 %13) {
%15 = sub nsw i32 %2, %1
%16 = xor i32 %15, -1
%17 = add nsw i32 %15, 1
%18 = sub i32 1, %3
%19 = add i32 %18, %4
%20 = mul nsw i32 %19, %17
%21 = mul nsw i32 %20, 3
%22 = mul i32 %20, -4
%23 = add i32 %22, %16
%24 = mul nsw i32 %17, %11
%25 = add nsw i32 %24, %10
%26 = mul nsw i32 %21, %6
%27 = add i32 %23, %26
%28 = add i32 %27, %20
%29 = add i32 %28, %25
%30 = sext i32 %29 to i64
%31 = getelementptr inbounds double, double* %0, i64 %30
%32 = load double, double* %31
%33 = mul nsw i32 %17, %9
%34 = add nsw i32 %33, %8
%35 = add i32 %28, %34
%36 = sext i32 %35 to i64
%37 = getelementptr inbounds double, double* %0, i64 %36
%38 = load double, double* %37
%39 = fsub double %32, %38
%40 = shl nsw i32 %20, 1
%41 = add i32 %27, %40
%42 = add i32 %41, %25
%43 = sext i32 %42 to i64
%44 = getelementptr inbounds double, double* %0, i64 %43
%45 = load double, double* %44
%46 = add i32 %41, %34
%47 = sext i32 %46 to i64
%48 = getelementptr inbounds double, double* %0, i64 %47
%49 = load double, double* %48
%50 = fsub double %45, %49
%51 = add i32 %27, %21
%52 = add i32 %51, %25
%53 = sext i32 %52 to i64
%54 = getelementptr inbounds double, double* %0, i64 %53
%55 = load double, double* %54
%56 = add i32 %51, %34
%57 = sext i32 %56 to i64
%58 = getelementptr inbounds double, double* %0, i64 %57
%59 = load double, double* %58
%60 = fsub double %55, %59
%61 = fmul double %50, %50
%62 = tail call double @llvm.fmuladd.f64(double %39, double %39, double %61)
%63 = tail call double @llvm.fmuladd.f64(double %60, double %60, double %62)
%64 = tail call double @sqrt(double %63)
%65 = fmul double %64, 5.000000e-01
%66 = tail call double @asin(double %65)
%67 = fmul double %66, 2.000000e+00
%68 = fmul double %67, %7
%69 = tail call double @sin(double %68)
%70 = tail call double @sin(double %67)
%71 = fdiv double %69, %70
%72 = fsub double 1.000000e+00, %7
%73 = fmul double %72, %67
%74 = tail call double @sin(double %73)
%75 = tail call double @sin(double %67)
%76 = fdiv double %74, %75
%77 = load double, double* %37
%78 = load double, double* %31
%79 = fmul double %71, %78
%80 = tail call double @llvm.fmuladd.f64(double %76, double %77, double %79)
%81 = mul nsw i32 %17, %13
%82 = add nsw i32 %81, %12
%83 = add i32 %28, %82
%84 = sext i32 %83 to i64
%85 = getelementptr inbounds double, double* %0, i64 %84
store double %80, double* %85
%86 = load double, double* %48
%87 = load double, double* %44
%88 = fmul double %71, %87
%89 = tail call double @llvm.fmuladd.f64(double %76, double %86, double %88)
%90 = add i32 %41, %82
%91 = sext i32 %90 to i64
%92 = getelementptr inbounds double, double* %0, i64 %91
store double %89, double* %92
%93 = load double, double* %58
%94 = load double, double* %54
%95 = fmul double %71, %94
%96 = tail call double @llvm.fmuladd.f64(double %76, double %93, double %95)
%97 = add i32 %51, %82
%98 = sext i32 %97 to i64
%99 = getelementptr inbounds double, double* %0, i64 %98
store double %96, double* %99
ret void
}