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


x86 assembly (clang, -O0):
<resp_mag>:
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,-0x98(%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,-0x90(%rbp)
mov    %rcx,%rsp
mov    %rax,-0x58(%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)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x38(%rbp)
movl   $0x0,-0x40(%rbp)
movsd  0xdfc(%rip),%xmm0
ucomisd -0x38(%rbp),%xmm0
jb     15a5 <resp_mag+0x475>
mov    -0x90(%rbp),%rax
mov    -0x98(%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    1374 <resp_mag+0x244>
mov    -0x98(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x38(%rbp),%xmm0
mov    $0x7,%edi
call   1050 <cosTyl@plt>
mov    -0x98(%rbp),%rax
movsd  %xmm0,-0xb0(%rbp)
ff
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0xb8(%rbp)
ff
movsd  -0x38(%rbp),%xmm0
mov    $0x7,%edi
call   1030 <sinTyl@plt>
movsd  -0xb8(%rbp),%xmm2
ff
mov    -0x98(%rbp),%rax
mov    -0x90(%rbp),%rcx
movaps %xmm0,%xmm1
movsd  -0xb0(%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,-0xa0(%rbp)
ff
movsd  -0x70(%rbp),%xmm0
movsd  %xmm0,-0xa8(%rbp)
ff
movsd  -0x38(%rbp),%xmm0
mov    $0x7,%edi
call   1050 <cosTyl@plt>
movsd  -0xa8(%rbp),%xmm2
ff
mov    -0x90(%rbp),%rax
movaps %xmm0,%xmm1
movsd  -0xa0(%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    1241 <resp_mag+0x111>
mov    -0x80(%rbp),%rax
mov    -0x88(%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    14d0 <resp_mag+0x3a0>
mov    -0x88(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x38(%rbp),%xmm0
mov    $0x7,%edi
call   1050 <cosTyl@plt>
mov    -0x88(%rbp),%rax
movsd  %xmm0,-0xd0(%rbp)
ff
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0xd8(%rbp)
ff
movsd  -0x38(%rbp),%xmm0
mov    $0x7,%edi
call   1030 <sinTyl@plt>
movsd  -0xd8(%rbp),%xmm2
ff
mov    -0x88(%rbp),%rax
mov    -0x80(%rbp),%rcx
movaps %xmm0,%xmm1
movsd  -0xd0(%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,-0xc0(%rbp)
ff
movsd  -0x70(%rbp),%xmm0
movsd  %xmm0,-0xc8(%rbp)
ff
movsd  -0x38(%rbp),%xmm0
mov    $0x7,%edi
call   1050 <cosTyl@plt>
movsd  -0xc8(%rbp),%xmm2
ff
mov    -0x80(%rbp),%rax
movaps %xmm0,%xmm1
movsd  -0xc0(%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    13a3 <resp_mag+0x273>
mov    -0x90(%rbp),%rax
mov    -0x98(%rbp),%rcx
movslq -0x40(%rbp),%rdx
movsd  (%rcx,%rdx,8),%xmm0
movslq -0x40(%rbp),%rdx
movsd  (%rcx,%rdx,8),%xmm2
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movslq -0x40(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
call   1040 <sqrt@plt>
mov    -0x88(%rbp),%rcx
mov    -0x80(%rbp),%rax
mov    -0x28(%rbp),%rdx
movslq -0x40(%rbp),%rsi
movsd  %xmm0,(%rdx,%rsi,8)
movslq -0x40(%rbp),%rdx
movsd  (%rcx,%rdx,8),%xmm0
movslq -0x40(%rbp),%rdx
movsd  (%rcx,%rdx,8),%xmm2
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movslq -0x40(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
call   1040 <sqrt@plt>
movsd  %xmm0,-0x78(%rbp)
mov    -0x28(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  -0x78(%rbp),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x40(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
cvtsi2sdl -0x2c(%rbp),%xmm1
movsd  0xa77(%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    11fc <resp_mag+0xcc>
mov    -0x48(%rbp),%rax
mov    %rax,%rsp
mov    %rbp,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<resp_mag>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x58,%rsp
cvtsi2sd %r9d,%xmm0
mov    %r8,0x38(%rsp)
mov    %rdx,%r15
mov    %rdi,%rbx
movsd  0xeaa(%rip),%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,0x50(%rsp)
mov    %esi,0x30(%rsp)
mov    %esi,%r13d
mov    %ecx,0x34(%rsp)
mov    %ecx,%r14d
xor    %ebp,%ebp
xorpd  %xmm2,%xmm2
jmp    11ba <resp_mag+0x8a>
cs nopw 0x0(%rax,%rax,1)
movapd %xmm3,%xmm0
call   1040 <sqrt@plt>
movsd  (%rsp),%xmm2
movsd  (%r12,%rbp,8),%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,(%r12,%rbp,8)
addsd  0x50(%rsp),%xmm2
add    $0x1,%rbp
movsd  0xe50(%rip),%xmm0
ucomisd %xmm2,%xmm0
jb     140e <resp_mag+0x2de>
movsd  (%rbx),%xmm0
cmpl   $0x2,0x30(%rsp)
movsd  %xmm2,(%rsp)
jl     1290 <resp_mag+0x160>
mov    $0x1,%r12d
xorps  %xmm3,%xmm3
nopw   0x0(%rax,%rax,1)
movsd  %xmm3,0x10(%rsp)
movsd  %xmm0,0x8(%rsp)
movapd %xmm2,%xmm0
mov    $0x7,%edi
call   1050 <cosTyl@plt>
movsd  %xmm0,0x18(%rsp)
movsd  (%rsp),%xmm0
mov    $0x7,%edi
call   1030 <sinTyl@plt>
mulsd  0x10(%rsp),%xmm0
movsd  0x18(%rsp),%xmm1
mulsd  0x8(%rsp),%xmm1
subsd  %xmm0,%xmm1
addsd  (%rbx,%r12,8),%xmm1
movsd  %xmm1,0x18(%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  0x18(%rsp),%xmm4
movsd  0x28(%rsp),%xmm1
movsd  (%rsp),%xmm2
mulsd  0x10(%rsp),%xmm0
mulsd  0x8(%rsp),%xmm1
addsd  %xmm0,%xmm1
add    $0x1,%r12
movapd %xmm1,%xmm3
movapd %xmm4,%xmm0
cmp    %r12,%r13
jne    11e0 <resp_mag+0xb0>
jmp    1298 <resp_mag+0x168>
nopl   (%rax)
xorpd  %xmm1,%xmm1
movapd %xmm0,%xmm4
movsd  %xmm4,0x18(%rsp)
movsd  (%r15),%xmm0
xorpd  %xmm4,%xmm4
cmpl   $0x2,0x34(%rsp)
mov    %rbp,0x40(%rsp)
movsd  %xmm1,0x28(%rsp)
jl     1380 <resp_mag+0x250>
mov    $0x1,%ebp
xorpd  %xmm1,%xmm1
mov    0x38(%rsp),%r12
nopl   0x0(%rax,%rax,1)
movsd  %xmm1,0x48(%rsp)
movsd  %xmm0,0x10(%rsp)
movapd %xmm2,%xmm0
mov    $0x7,%edi
call   1050 <cosTyl@plt>
movsd  %xmm0,0x8(%rsp)
movsd  (%rsp),%xmm0
mov    $0x7,%edi
call   1030 <sinTyl@plt>
mulsd  0x48(%rsp),%xmm0
movsd  0x8(%rsp),%xmm1
mulsd  0x10(%rsp),%xmm1
subsd  %xmm0,%xmm1
addsd  (%r15,%rbp,8),%xmm1
movsd  %xmm1,0x8(%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
movsd  0x8(%rsp),%xmm3
movsd  (%rsp),%xmm2
mulsd  0x48(%rsp),%xmm0
mulsd  0x10(%rsp),%xmm4
addsd  %xmm0,%xmm4
add    $0x1,%rbp
movapd %xmm4,%xmm1
movapd %xmm3,%xmm0
cmp    %rbp,%r14
jne    12d0 <resp_mag+0x1a0>
jmp    1389 <resp_mag+0x259>
nopl   (%rax)
movapd %xmm0,%xmm3
mov    0x38(%rsp),%r12
movsd  0x28(%rsp),%xmm1
mulsd  %xmm1,%xmm1
movsd  0x18(%rsp),%xmm0
mulsd  %xmm0,%xmm0
addsd  %xmm1,%xmm0
ucomisd 0xc5f(%rip),%xmm0
jb     13c0 <resp_mag+0x290>
sqrtsd %xmm0,%xmm0
jmp    13dd <resp_mag+0x2ad>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movsd  %xmm3,0x8(%rsp)
movsd  %xmm4,0x20(%rsp)
call   1040 <sqrt@plt>
movsd  0x20(%rsp),%xmm4
movsd  0x8(%rsp),%xmm3
mov    0x40(%rsp),%rbp
movsd  %xmm0,(%r12,%rbp,8)
mulsd  %xmm4,%xmm4
mulsd  %xmm3,%xmm3
addsd  %xmm4,%xmm3
ucomisd 0xc0c(%rip),%xmm3
jb     1180 <resp_mag+0x50>
xorps  %xmm0,%xmm0
sqrtsd %xmm3,%xmm0
jmp    1189 <resp_mag+0x59>
add    $0x58,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define void @resp_mag(double* nocapture %0, i32 %1, double* nocapture %2, i32 %3, double* nocapture %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 ], [ %70, %57 ]
%15 = phi double [ 0.000000e+00, %6 ], [ %69, %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 [ 0.000000e+00, %13 ], [ %32, %17 ]
%37 = phi double [ %16, %13 ], [ %28, %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 [ 0.000000e+00, %35 ], [ %54, %39 ]
%59 = phi double [ %38, %35 ], [ %50, %39 ]
%60 = fmul double %36, %36
%61 = tail call double @llvm.fmuladd.f64(double %37, double %37, double %60)
%62 = tail call double @sqrt(double %61)
%63 = getelementptr inbounds double, double* %4, i64 %14
store double %62, double* %63
%64 = fmul double %58, %58
%65 = tail call double @llvm.fmuladd.f64(double %59, double %59, double %64)
%66 = tail call double @sqrt(double %65)
%67 = load double, double* %63
%68 = fdiv double %67, %66
store double %68, double* %63
%69 = fadd double %10, %15
%70 = add nuw i64 %14, 1
%71 = fcmp ugt double %69, 0x400921FB54442D18
br i1 %71, label %72, label %13
72: ; preds = %57
ret void
}