#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 atan(double);
double cosTyl(double, int);
double sinTyl(double, int);
void resp_phase(double *num, int lnum, double *den, int lden, double *res,
                int N) {
  double w;
  int m, i;
  double out_numRe[N + 1], old_out_r;
  double out_numIm[N + 1];
  double out_denRe[N + 1], out_denIm[N + 1];
  for (w = 0, i = 0; w <= 3.14159265358979323846;
       w += 3.14159265358979323846 / N, ++i) {
    out_numRe[i] = num[0];
    out_numIm[i] = 0;
    for (m = 1; m < lnum; ++m) {
      old_out_r = out_numRe[i];
      out_numRe[i] =
          cosTyl(w, 7) * out_numRe[i] - sinTyl(w, 7) * out_numIm[i] + num[m];
      out_numIm[i] = sinTyl(w, 7) * old_out_r + cosTyl(w, 7) * out_numIm[i];
    }
    out_denRe[i] = den[0];
    out_denIm[i] = 0;
    for (m = 1; m < lden; ++m) {
      old_out_r = out_denRe[i];
      out_denRe[i] =
          cosTyl(w, 7) * out_denRe[i] - sinTyl(w, 7) * out_denIm[i] + den[m];
      out_denIm[i] = sinTyl(w, 7) * old_out_r + cosTyl(w, 7) * out_denIm[i];
    }
    res[i] = atan(out_numIm[i] / out_numRe[i]);
    res[i] = res[i] - atan(out_denIm[i] / out_denRe[i]);
  }
}


x86 assembly (clang, -O0):
<resp_phase>:
push   %rbp
mov    %rsp,%rbp
sub    $0xe0,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9d,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %rsp,%rcx
mov    %rcx,-0x48(%rbp)
lea    0xf(,%rax,8),%rdx
and    $0xfffffffffffffff0,%rdx
mov    %rsp,%rcx
sub    %rdx,%rcx
mov    %rcx,-0x90(%rbp)
mov    %rcx,%rsp
mov    %rax,-0x50(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
lea    0xf(,%rax,8),%rdx
and    $0xfffffffffffffff0,%rdx
mov    %rsp,%rcx
sub    %rdx,%rcx
mov    %rcx,-0x88(%rbp)
mov    %rcx,%rsp
mov    %rax,-0x60(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
lea    0xf(,%rax,8),%rdx
and    $0xfffffffffffffff0,%rdx
mov    %rsp,%rcx
sub    %rdx,%rcx
mov    %rcx,-0x80(%rbp)
mov    %rcx,%rsp
mov    %rax,-0x68(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
lea    0xf(,%rax,8),%rdx
and    $0xfffffffffffffff0,%rdx
mov    %rsp,%rcx
sub    %rdx,%rcx
mov    %rcx,-0x78(%rbp)
mov    %rcx,%rsp
mov    %rax,-0x70(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x38(%rbp)
movl   $0x0,-0x40(%rbp)
movsd  0xdff(%rip),%xmm0
ucomisd -0x38(%rbp),%xmm0
jb     156c <resp_phase+0x43c>
mov    -0x88(%rbp),%rax
mov    -0x90(%rbp),%rcx
mov    -0x8(%rbp),%rdx
movsd  (%rdx),%xmm0
movslq -0x40(%rbp),%rdx
movsd  %xmm0,(%rcx,%rdx,8)
movslq -0x40(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movl   $0x1,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1371 <resp_phase+0x241>
mov    -0x90(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x38(%rbp),%xmm0
mov    $0x7,%edi
call   1050 <cosTyl@plt>
mov    -0x90(%rbp),%rax
movsd  %xmm0,-0xa8(%rbp)
ff
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0xb0(%rbp)
ff
movsd  -0x38(%rbp),%xmm0
mov    $0x7,%edi
call   1030 <sinTyl@plt>
movsd  -0xb0(%rbp),%xmm2
ff
mov    -0x90(%rbp),%rax
mov    -0x88(%rbp),%rcx
movaps %xmm0,%xmm1
movsd  -0xa8(%rbp),%xmm0
ff
movslq -0x40(%rbp),%rdx
mulsd  (%rcx,%rdx,8),%xmm1
movq   %xmm1,%rcx
movabs $0x8000000000000000,%rdx
xor    %rdx,%rcx
movq   %rcx,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x8(%rbp),%rcx
movslq -0x3c(%rbp),%rdx
addsd  (%rcx,%rdx,8),%xmm0
movslq -0x40(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x38(%rbp),%xmm0
mov    $0x7,%edi
call   1030 <sinTyl@plt>
movsd  %xmm0,-0x98(%rbp)
ff
movsd  -0x58(%rbp),%xmm0
movsd  %xmm0,-0xa0(%rbp)
ff
movsd  -0x38(%rbp),%xmm0
mov    $0x7,%edi
call   1050 <cosTyl@plt>
movsd  -0xa0(%rbp),%xmm2
ff
mov    -0x88(%rbp),%rax
movaps %xmm0,%xmm1
movsd  -0x98(%rbp),%xmm0
ff
movslq -0x40(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movslq -0x40(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    123e <resp_phase+0x10e>
mov    -0x78(%rbp),%rax
mov    -0x80(%rbp),%rcx
mov    -0x18(%rbp),%rdx
movsd  (%rdx),%xmm0
movslq -0x40(%rbp),%rdx
movsd  %xmm0,(%rcx,%rdx,8)
movslq -0x40(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movl   $0x1,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    14c1 <resp_phase+0x391>
mov    -0x80(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x38(%rbp),%xmm0
mov    $0x7,%edi
call   1050 <cosTyl@plt>
mov    -0x80(%rbp),%rax
movsd  %xmm0,-0xc8(%rbp)
ff
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0xd0(%rbp)
ff
movsd  -0x38(%rbp),%xmm0
mov    $0x7,%edi
call   1030 <sinTyl@plt>
movsd  -0xd0(%rbp),%xmm2
ff
mov    -0x80(%rbp),%rax
mov    -0x78(%rbp),%rcx
movaps %xmm0,%xmm1
movsd  -0xc8(%rbp),%xmm0
ff
movslq -0x40(%rbp),%rdx
mulsd  (%rcx,%rdx,8),%xmm1
movq   %xmm1,%rcx
movabs $0x8000000000000000,%rdx
xor    %rdx,%rcx
movq   %rcx,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rcx
movslq -0x3c(%rbp),%rdx
addsd  (%rcx,%rdx,8),%xmm0
movslq -0x40(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x38(%rbp),%xmm0
mov    $0x7,%edi
call   1030 <sinTyl@plt>
movsd  %xmm0,-0xb8(%rbp)
ff
movsd  -0x58(%rbp),%xmm0
movsd  %xmm0,-0xc0(%rbp)
ff
movsd  -0x38(%rbp),%xmm0
mov    $0x7,%edi
call   1050 <cosTyl@plt>
movsd  -0xc0(%rbp),%xmm2
ff
mov    -0x78(%rbp),%rax
movaps %xmm0,%xmm1
movsd  -0xb8(%rbp),%xmm0
ff
movslq -0x40(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movslq -0x40(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    139d <resp_phase+0x26d>
mov    -0x90(%rbp),%rax
mov    -0x88(%rbp),%rcx
movslq -0x40(%rbp),%rdx
movsd  (%rcx,%rdx,8),%xmm0
movslq -0x40(%rbp),%rcx
divsd  (%rax,%rcx,8),%xmm0
call   1040 <atan@plt>
mov    -0x78(%rbp),%rcx
mov    -0x80(%rbp),%rax
mov    -0x28(%rbp),%rdx
movslq -0x40(%rbp),%rsi
movsd  %xmm0,(%rdx,%rsi,8)
mov    -0x28(%rbp),%rdx
movslq -0x40(%rbp),%rsi
movsd  (%rdx,%rsi,8),%xmm0
movsd  %xmm0,-0xd8(%rbp)
ff
movslq -0x40(%rbp),%rdx
movsd  (%rcx,%rdx,8),%xmm0
movslq -0x40(%rbp),%rcx
divsd  (%rax,%rcx,8),%xmm0
call   1040 <atan@plt>
movaps %xmm0,%xmm1
movsd  -0xd8(%rbp),%xmm0
ff
subsd  %xmm1,%xmm0
mov    -0x28(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
cvtsi2sdl -0x2c(%rbp),%xmm1
movsd  0xab0(%rip),%xmm0
divsd  %xmm1,%xmm0
addsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    11f9 <resp_phase+0xc9>
mov    -0x48(%rbp),%rax
mov    %rax,%rsp
mov    %rbp,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<resp_phase>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x58,%rsp
mov    %r8,0x48(%rsp)
mov    %rdx,%r15
mov    %rdi,%rbp
cvtsi2sd %r9d,%xmm0
movsd  0xeaa(%rip),%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,0x40(%rsp)
mov    %esi,0x30(%rsp)
mov    %esi,%r13d
mov    %ecx,0x34(%rsp)
mov    %ecx,%r14d
xorpd  %xmm1,%xmm1
xor    %eax,%eax
jmp    11ed <resp_phase+0xbd>
cs nopw 0x0(%rax,%rax,1)
movapd %xmm3,%xmm0
divsd  %xmm0,%xmm4
movsd  %xmm4,0x20(%rsp)
movsd  0x28(%rsp),%xmm0
divsd  0x8(%rsp),%xmm0
call   1040 <atan@plt>
movsd  %xmm0,0x8(%rsp)
movsd  0x20(%rsp),%xmm0
call   1040 <atan@plt>
movsd  0x8(%rsp),%xmm1
subsd  %xmm0,%xmm1
mov    0x48(%rsp),%rax
mov    0x50(%rsp),%rcx
movsd  %xmm1,(%rax,%rcx,8)
mov    %rcx,%rax
movsd  (%rsp),%xmm1
addsd  0x40(%rsp),%xmm1
add    $0x1,%rax
movsd  0xe1d(%rip),%xmm0
ucomisd %xmm1,%xmm0
jb     13ac <resp_phase+0x27c>
movsd  0x0(%rbp),%xmm3
cmpl   $0x2,0x30(%rsp)
movsd  %xmm1,(%rsp)
mov    %rax,0x50(%rsp)
jl     12c0 <resp_phase+0x190>
mov    $0x1,%r12d
xorps  %xmm4,%xmm4
movsd  %xmm4,0x18(%rsp)
movsd  %xmm3,0x10(%rsp)
movapd %xmm1,%xmm0
mov    $0x7,%edi
call   1050 <cosTyl@plt>
movsd  %xmm0,0x8(%rsp)
movsd  (%rsp),%xmm0
mov    $0x7,%edi
call   1030 <sinTyl@plt>
mulsd  0x18(%rsp),%xmm0
movsd  0x8(%rsp),%xmm1
mulsd  0x10(%rsp),%xmm1
subsd  %xmm0,%xmm1
addsd  0x0(%rbp,%r12,8),%xmm1
movsd  %xmm1,0x8(%rsp)
movsd  (%rsp),%xmm0
mov    $0x7,%edi
call   1030 <sinTyl@plt>
movsd  %xmm0,0x28(%rsp)
movsd  (%rsp),%xmm0
mov    $0x7,%edi
call   1050 <cosTyl@plt>
movsd  0x28(%rsp),%xmm2
movsd  (%rsp),%xmm1
mulsd  0x18(%rsp),%xmm0
mulsd  0x10(%rsp),%xmm2
addsd  %xmm0,%xmm2
movsd  0x8(%rsp),%xmm0
add    $0x1,%r12
movapd %xmm2,%xmm4
movapd %xmm0,%xmm3
cmp    %r12,%r13
jne    1210 <resp_phase+0xe0>
jmp    12c8 <resp_phase+0x198>
xchg   %ax,%ax
movapd %xmm3,%xmm0
xorpd  %xmm2,%xmm2
movsd  (%r15),%xmm3
xorpd  %xmm4,%xmm4
cmpl   $0x2,0x34(%rsp)
movsd  %xmm0,0x8(%rsp)
movsd  %xmm2,0x28(%rsp)
jl     1180 <resp_phase+0x50>
mov    $0x1,%ebx
xorpd  %xmm2,%xmm2
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movsd  %xmm2,0x38(%rsp)
movsd  %xmm3,0x18(%rsp)
movsd  (%rsp),%xmm0
mov    $0x7,%edi
call   1050 <cosTyl@plt>
movsd  %xmm0,0x10(%rsp)
movsd  (%rsp),%xmm0
mov    $0x7,%edi
call   1030 <sinTyl@plt>
mulsd  0x38(%rsp),%xmm0
movsd  0x10(%rsp),%xmm1
mulsd  0x18(%rsp),%xmm1
subsd  %xmm0,%xmm1
addsd  (%r15,%rbx,8),%xmm1
movsd  %xmm1,0x10(%rsp)
movsd  (%rsp),%xmm0
mov    $0x7,%edi
call   1030 <sinTyl@plt>
movsd  %xmm0,0x20(%rsp)
movsd  (%rsp),%xmm0
mov    $0x7,%edi
call   1050 <cosTyl@plt>
movsd  0x20(%rsp),%xmm4
mulsd  0x38(%rsp),%xmm0
mulsd  0x18(%rsp),%xmm4
addsd  %xmm0,%xmm4
movsd  0x10(%rsp),%xmm0
add    $0x1,%rbx
movapd %xmm4,%xmm2
movapd %xmm0,%xmm3
cmp    %rbx,%r14
jne    1300 <resp_phase+0x1d0>
jmp    1184 <resp_phase+0x54>
add    $0x58,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define void @resp_phase(double* nocapture %0, i32 %1, double* nocapture %2, i32 %3, double* nocapture writeonly %4, i32 %5) {
%7 = icmp sgt i32 %1, 1
%8 = icmp sgt i32 %3, 1
%9 = sitofp i32 %5 to double
%10 = fdiv double 0x400921FB54442D18, %9
%11 = zext i32 %1 to i64
%12 = zext i32 %3 to i64
br label %13
13: ; preds = %6, %57
%14 = phi i64 [ 0, %6 ], [ %67, %57 ]
%15 = phi double [ 0.000000e+00, %6 ], [ %66, %57 ]
%16 = load double, double* %0
br i1 %7, label %17, label %35
17: ; preds = %13, %17
%18 = phi double [ %32, %17 ], [ 0.000000e+00, %13 ]
%19 = phi i64 [ %33, %17 ], [ 1, %13 ]
%20 = phi double [ %28, %17 ], [ %16, %13 ]
%21 = tail call double @cosTyl(double %15, i32 7)
%22 = tail call double @sinTyl(double %15, i32 7)
%23 = fneg double %22
%24 = fmul double %18, %23
%25 = tail call double @llvm.fmuladd.f64(double %21, double %20, double %24)
%26 = getelementptr inbounds double, double* %0, i64 %19
%27 = load double, double* %26
%28 = fadd double %27, %25
%29 = tail call double @sinTyl(double %15, i32 7)
%30 = tail call double @cosTyl(double %15, i32 7)
%31 = fmul double %18, %30
%32 = tail call double @llvm.fmuladd.f64(double %29, double %20, double %31)
%33 = add nuw nsw i64 %19, 1
%34 = icmp eq i64 %33, %11
br i1 %34, label %35, label %17
35: ; preds = %17, %13
%36 = phi double [ %16, %13 ], [ %28, %17 ]
%37 = phi double [ 0.000000e+00, %13 ], [ %32, %17 ]
%38 = load double, double* %2
br i1 %8, label %39, label %57
39: ; preds = %35, %39
%40 = phi double [ %54, %39 ], [ 0.000000e+00, %35 ]
%41 = phi i64 [ %55, %39 ], [ 1, %35 ]
%42 = phi double [ %50, %39 ], [ %38, %35 ]
%43 = tail call double @cosTyl(double %15, i32 7)
%44 = tail call double @sinTyl(double %15, i32 7)
%45 = fneg double %44
%46 = fmul double %40, %45
%47 = tail call double @llvm.fmuladd.f64(double %43, double %42, double %46)
%48 = getelementptr inbounds double, double* %2, i64 %41
%49 = load double, double* %48
%50 = fadd double %49, %47
%51 = tail call double @sinTyl(double %15, i32 7)
%52 = tail call double @cosTyl(double %15, i32 7)
%53 = fmul double %40, %52
%54 = tail call double @llvm.fmuladd.f64(double %51, double %42, double %53)
%55 = add nuw nsw i64 %41, 1
%56 = icmp eq i64 %55, %12
br i1 %56, label %57, label %39
57: ; preds = %39, %35
%58 = phi double [ %38, %35 ], [ %50, %39 ]
%59 = phi double [ 0.000000e+00, %35 ], [ %54, %39 ]
%60 = fdiv double %37, %36
%61 = tail call double @atan(double %60)
%62 = getelementptr inbounds double, double* %4, i64 %14
%63 = fdiv double %59, %58
%64 = tail call double @atan(double %63)
%65 = fsub double %61, %64
store double %65, double* %62
%66 = fadd double %10, %15
%67 = add nuw i64 %14, 1
%68 = fcmp ugt double %66, 0x400921FB54442D18
br i1 %68, label %69, label %13
69: ; preds = %57
ret void
}