#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 fabs(double);
int /*<<< orphan*/ free_vector(double *, int, int);
int /*<<< orphan*/ nrerror(char *);
double *vector(int, int);
void ratint(double xa[], double ya[], int n, double x, double *y, double *dy) {
  int m, i, ns = 1;
  double w, t, hh, h, dd, *c, *d;
  c = vector(1, n);
  d = vector(1, n);
  hh = fabs(x - xa[1]);
  for (i = 1; i <= n; i++) {
    h = fabs(x - xa[i]);
    if (h == 0.0) {
      *y = ya[i];
      *dy = 0.0;
      {
        free_vector(d, 1, n);
        free_vector(c, 1, n);
        return;
      }
    } else if (h < hh) {
      ns = i;
      hh = h;
    }
    c[i] = ya[i];
    d[i] = ya[i] + 1.0e-25;
  }
  *y = ya[ns--];
  for (m = 1; m < n; m++) {
    for (i = 1; i <= n - m; i++) {
      w = c[i + 1] - d[i];
      h = xa[i + m] - x;
      t = (xa[i] - x) * d[i] / h;
      dd = t - c[i + 1];
      if (dd == 0.0)
        nrerror("Error in routine ratint");
      dd = w / dd;
      d[i] = c[i + 1] * dd;
      c[i] = t * dd;
    }
    *y += (*dy = (2 * ns < (n - m) ? c[ns + 1] : d[ns--]));
  }
  {
    free_vector(d, 1, n);
    free_vector(c, 1, n);
    return;
  }
}


x86 assembly (clang, -O0):
<ratint>:
push   %rbp
mov    %rsp,%rbp
sub    $0x80,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
movsd  %xmm0,-0x20(%rbp)
mov    %rcx,-0x28(%rbp)
mov    %r8,-0x30(%rbp)
movl   $0x1,-0x3c(%rbp)
mov    -0x14(%rbp),%esi
mov    $0x1,%edi
call   1050 <vector@plt>
mov    %rax,-0x70(%rbp)
mov    -0x14(%rbp),%esi
mov    $0x1,%edi
call   1050 <vector@plt>
mov    %rax,-0x78(%rbp)
movsd  -0x20(%rbp),%xmm0
mov    -0x8(%rbp),%rax
subsd  0x8(%rax),%xmm0
movaps 0xe6f(%rip),%xmm1
pand   %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
movl   $0x1,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x14(%rbp),%eax
jg     12a1 <ratint+0x171>
movsd  -0x20(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x38(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
movaps 0xe3a(%rip),%xmm1
pand   %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  -0x60(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    122e <ratint+0xfe>
jp     122e <ratint+0xfe>
mov    -0x10(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x30(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax)
mov    -0x78(%rbp),%rdi
mov    -0x14(%rbp),%edx
mov    $0x1,%esi
call   1040 <free_vector@plt>
mov    -0x70(%rbp),%rdi
mov    -0x14(%rbp),%edx
mov    $0x1,%esi
call   1040 <free_vector@plt>
jmp    1481 <ratint+0x351>
movsd  -0x60(%rbp),%xmm1
movsd  -0x58(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    1252 <ratint+0x122>
mov    -0x38(%rbp),%eax
mov    %eax,-0x3c(%rbp)
movsd  -0x60(%rbp),%xmm0
movsd  %xmm0,-0x58(%rbp)
jmp    1257 <ratint+0x127>
mov    -0x10(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x70(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  0xd8f(%rip),%xmm0
addsd  (%rax,%rcx,8),%xmm0
mov    -0x78(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    11a1 <ratint+0x71>
mov    -0x10(%rbp),%rax
mov    -0x3c(%rbp),%ecx
mov    %ecx,%edx
add    $0xffffffff,%edx
mov    %edx,-0x3c(%rbp)
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
movsd  %xmm0,(%rax)
movl   $0x1,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    145f <ratint+0x32f>
movl   $0x1,-0x38(%rbp)
mov    -0x38(%rbp),%eax
mov    -0x14(%rbp),%ecx
sub    -0x34(%rbp),%ecx
cmp    %ecx,%eax
jg     13ec <ratint+0x2bc>
mov    -0x70(%rbp),%rax
mov    -0x38(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x78(%rbp),%rax
movslq -0x38(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%ecx
add    -0x34(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x20(%rbp),%xmm0
mov    -0x78(%rbp),%rax
movslq -0x38(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
divsd  -0x60(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x50(%rbp),%xmm0
mov    -0x70(%rbp),%rax
mov    -0x38(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  -0x68(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1394 <ratint+0x264>
jp     1394 <ratint+0x264>
lea    0xc89(%rip),%rdi
call   1030 <nrerror@plt>
movsd  -0x48(%rbp),%xmm0
divsd  -0x68(%rbp),%xmm0
movsd  %xmm0,-0x68(%rbp)
mov    -0x70(%rbp),%rax
mov    -0x38(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mulsd  -0x68(%rbp),%xmm0
mov    -0x78(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x50(%rbp),%xmm0
mulsd  -0x68(%rbp),%xmm0
mov    -0x70(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    12da <ratint+0x1aa>
mov    -0x3c(%rbp),%eax
shl    $0x1,%eax
mov    -0x14(%rbp),%ecx
sub    -0x34(%rbp),%ecx
cmp    %ecx,%eax
jge    141c <ratint+0x2ec>
mov    -0x70(%rbp),%rax
mov    -0x3c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x80(%rbp)
jmp    1438 <ratint+0x308>
mov    -0x78(%rbp),%rax
mov    -0x3c(%rbp),%ecx
mov    %ecx,%edx
add    $0xffffffff,%edx
mov    %edx,-0x3c(%rbp)
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x80(%rbp)
movsd  -0x80(%rbp),%xmm0
mov    -0x30(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x28(%rbp),%rax
addsd  (%rax),%xmm0
movsd  %xmm0,(%rax)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    12c7 <ratint+0x197>
mov    -0x78(%rbp),%rdi
mov    -0x14(%rbp),%edx
mov    $0x1,%esi
call   1040 <free_vector@plt>
mov    -0x70(%rbp),%rdi
mov    -0x14(%rbp),%edx
mov    $0x1,%esi
call   1040 <free_vector@plt>
add    $0x80,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<ratint>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x68,%rsp
mov    %r8,0x18(%rsp)
mov    %rcx,%r14
movsd  %xmm0,0x20(%rsp)
mov    %edx,%r13d
mov    %rsi,%r12
mov    %rdi,%r15
mov    $0x1,%edi
mov    %edx,%esi
call   1050 <vector@plt>
mov    %rax,%rbx
mov    $0x1,%edi
mov    %r13d,%esi
call   1050 <vector@plt>
mov    %rax,%rbp
test   %r13d,%r13d
jle    1254 <ratint+0x124>
movsd  0x20(%rsp),%xmm6
movapd %xmm6,%xmm1
subsd  0x8(%r15),%xmm1
mov    $0x1,%eax
xorpd  %xmm0,%xmm0
ucomisd %xmm0,%xmm1
jne    11a2 <ratint+0x72>
jnp    123b <ratint+0x10b>
lea    0x1(%r13),%eax
andpd  0xe52(%rip),%xmm1
mov    $0x2,%ecx
sub    %rax,%rcx
mov    $0x1,%eax
mov    $0x1,%edx
mov    %rdx,0x8(%rsp)
movsd  0xe43(%rip),%xmm2
movapd 0xe2b(%rip),%xmm3
movapd %xmm1,%xmm4
nopl   0x0(%rax)
ucomisd %xmm1,%xmm4
mov    0x8(%rsp),%rdx
mov    %eax,%edi
cmovbe %edx,%edi
mov    %rdi,0x8(%rsp)
movsd  (%r12,%rax,8),%xmm5
movsd  %xmm5,(%rbx,%rax,8)
movsd  (%r12,%rax,8),%xmm5
addsd  %xmm2,%xmm5
movsd  %xmm5,0x0(%rbp,%rax,8)
lea    (%rcx,%rax,1),%rdx
cmp    $0x1,%rdx
je     128e <ratint+0x15e>
minsd  %xmm4,%xmm1
movapd %xmm1,%xmm4
movapd %xmm6,%xmm1
subsd  0x8(%r15,%rax,8),%xmm1
add    $0x1,%rax
ucomisd %xmm0,%xmm1
andpd  %xmm3,%xmm1
jne    11e0 <ratint+0xb0>
jp     11e0 <ratint+0xb0>
movsd  (%r12,%rax,8),%xmm0
movsd  %xmm0,(%r14)
mov    0x18(%rsp),%rax
movq   $0x0,(%rax)
jmp    1260 <ratint+0x130>
movsd  0x8(%r12),%xmm0
movsd  %xmm0,(%r14)
mov    %rbp,%rdi
mov    $0x1,%esi
mov    %r13d,%edx
call   1040 <free_vector@plt>
mov    %rbx,%rdi
mov    $0x1,%esi
mov    %r13d,%edx
add    $0x68,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
jmp    1040 <free_vector@plt>
movslq 0x8(%rsp),%rax
movsd  (%r12,%rax,8),%xmm0
mov    %r14,0x38(%rsp)
movsd  %xmm0,(%r14)
cmp    $0x2,%r13d
jl     1260 <ratint+0x130>
mov    0x8(%rsp),%rax
add    $0xffffffff,%eax
mov    %rax,0x8(%rsp)
mov    %r13d,%eax
mov    %rax,0x28(%rsp)
lea    0x10(%r15),%r14
mov    $0x1,%eax
mov    %r13d,%ecx
mov    %r13d,0x10(%rsp)
mov    %r13,0x30(%rsp)
jmp    1345 <ratint+0x215>
cs nopw 0x0(%rax,%rax,1)
mov    0x8(%rsp),%rdi
lea    (%rdi,%rdi,1),%eax
xor    %ecx,%ecx
cmp    0x14(%rsp),%eax
setge  %cl
movslq %edi,%rdi
lea    0x8(%rbx,%rdi,8),%rax
lea    0x0(%rbp,%rdi,8),%rdx
cmovl  %rax,%rdx
sub    %ecx,%edi
mov    %rdi,0x8(%rsp)
movsd  (%rdx),%xmm0
mov    0x18(%rsp),%rax
movsd  %xmm0,(%rax)
mov    0x38(%rsp),%rax
addsd  (%rax),%xmm0
movsd  %xmm0,(%rax)
mov    0x40(%rsp),%rax
add    $0x1,%rax
addl   $0xffffffff,0x10(%rsp)
add    $0x8,%r14
cmp    0x28(%rsp),%rax
mov    0x30(%rsp),%r13
je     1260 <ratint+0x130>
mov    %r13d,%ecx
mov    %rax,0x40(%rsp)
sub    %eax,%ecx
mov    %ecx,0x14(%rsp)
jle    12e0 <ratint+0x1b0>
mov    0x10(%rsp),%r12d
add    $0xffffffffffffffff,%r12
xor    %r13d,%r13d
jmp    139b <ratint+0x26b>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
subsd  %xmm3,%xmm1
divsd  %xmm4,%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,0x8(%rbp,%r13,8)
mulsd  %xmm1,%xmm2
movsd  %xmm2,0x8(%rbx,%r13,8)
add    $0x1,%r13
cmp    %r13,%r12
je     12e0 <ratint+0x1b0>
movsd  0x10(%rbx,%r13,8),%xmm1
movsd  0x8(%rbp,%r13,8),%xmm3
movsd  (%r14,%r13,8),%xmm0
subsd  %xmm6,%xmm0
movsd  0x8(%r15,%r13,8),%xmm2
subsd  %xmm6,%xmm2
mulsd  %xmm3,%xmm2
divsd  %xmm0,%xmm2
movapd %xmm2,%xmm4
subsd  %xmm1,%xmm4
ucomisd 0xc42(%rip),%xmm4
movapd %xmm1,%xmm0
jne    1370 <ratint+0x240>
jp     1370 <ratint+0x240>
lea    0xc3b(%rip),%rdi
movsd  %xmm1,0x60(%rsp)
movsd  %xmm2,0x58(%rsp)
movsd  %xmm3,0x50(%rsp)
movsd  %xmm4,0x48(%rsp)
call   1030 <nrerror@plt>
movsd  0x48(%rsp),%xmm4
movsd  0x50(%rsp),%xmm3
movsd  0x58(%rsp),%xmm2
movsd  0x60(%rsp),%xmm1
movsd  0x20(%rsp),%xmm6
movsd  0x10(%rbx,%r13,8),%xmm0
jmp    1370 <ratint+0x240>

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

LLVM IR (-O2):
@.str = [24 x i8] c"Error in routine ratint\00"
define void @ratint(double* nocapture %0, double* nocapture %1, i32 %2, double %3, double* nocapture %4, double* nocapture writeonly %5) {
%7 = tail call double* @vector(i32 1, i32 %2)
%8 = tail call double* @vector(i32 1, i32 %2)
%9 = getelementptr inbounds double, double* %0, i64 1
%10 = load double, double* %9
%11 = icmp slt i32 %2, 1
br i1 %11, label %12, label %15
12: ; preds = %6
%13 = getelementptr inbounds double, double* %1, i64 1
%14 = load double, double* %13
store double %14, double* %4
br label %117
15: ; preds = %6
%16 = add nuw i32 %2, 1
%17 = zext i32 %16 to i64
%18 = fsub double %3, %10
%19 = fcmp oeq double %18, 0.000000e+00
br i1 %19, label %26, label %20
20: ; preds = %15
%21 = tail call double @llvm.fabs.f64(double %18)
%22 = fsub double %3, %10
%23 = tail call double @llvm.fabs.f64(double %22)
br label %30
24: ; preds = %46
%25 = and i64 %44, 4294967295
br label %26
26: ; preds = %24, %15
%27 = phi i64 [ %25, %24 ], [ 1, %15 ]
%28 = getelementptr inbounds double, double* %1, i64 %27
%29 = load double, double* %28
store double %29, double* %4
store double 0.000000e+00, double* %5
br label %117
30: ; preds = %20, %46
%31 = phi double [ %21, %20 ], [ %51, %46 ]
%32 = phi i32 [ 1, %20 ], [ %37, %46 ]
%33 = phi double [ %23, %20 ], [ %47, %46 ]
%34 = phi i64 [ 1, %20 ], [ %44, %46 ]
%35 = fcmp olt double %31, %33
%36 = trunc i64 %34 to i32
%37 = select i1 %35, i32 %36, i32 %32
%38 = getelementptr inbounds double, double* %1, i64 %34
%39 = load double, double* %38
%40 = getelementptr inbounds double, double* %7, i64 %34
store double %39, double* %40
%41 = load double, double* %38
%42 = fadd double %41, 1.000000e-25
%43 = getelementptr inbounds double, double* %8, i64 %34
store double %42, double* %43
%44 = add nuw nsw i64 %34, 1
%45 = icmp eq i64 %44, %17
br i1 %45, label %53, label %46
46: ; preds = %30
%47 = select i1 %35, double %31, double %33
%48 = getelementptr inbounds double, double* %0, i64 %44
%49 = load double, double* %48
%50 = fsub double %3, %49
%51 = tail call double @llvm.fabs.f64(double %50)
%52 = fcmp oeq double %50, 0.000000e+00
br i1 %52, label %24, label %30
53: ; preds = %30
%54 = sext i32 %37 to i64
%55 = getelementptr inbounds double, double* %1, i64 %54
%56 = load double, double* %55
store double %56, double* %4
%57 = icmp sgt i32 %2, 1
br i1 %57, label %58, label %117
58: ; preds = %53
%59 = add nsw i32 %37, -1
%60 = zext i32 %2 to i64
br label %61
61: ; preds = %58, %99
%62 = phi i64 [ 1, %58 ], [ %114, %99 ]
%63 = phi i32 [ %2, %58 ], [ %115, %99 ]
%64 = phi i32 [ %59, %58 ], [ %109, %99 ]
%65 = trunc i64 %62 to i32
%66 = sub nsw i32 %2, %65
%67 = icmp slt i32 %66, 1
br i1 %67, label %99, label %68
68: ; preds = %61
%69 = zext i32 %63 to i64
br label %70
70: ; preds = %68, %92
%71 = phi i64 [ 1, %68 ], [ %72, %92 ]
%72 = add nuw nsw i64 %71, 1
%73 = getelementptr inbounds double, double* %7, i64 %72
%74 = load double, double* %73
%75 = getelementptr inbounds double, double* %8, i64 %71
%76 = load double, double* %75
%77 = fsub double %74, %76
%78 = add nuw nsw i64 %71, %62
%79 = getelementptr inbounds double, double* %0, i64 %78
%80 = load double, double* %79
%81 = fsub double %80, %3
%82 = getelementptr inbounds double, double* %0, i64 %71
%83 = load double, double* %82
%84 = fsub double %83, %3
%85 = fmul double %76, %84
%86 = fdiv double %85, %81
%87 = fsub double %86, %74
%88 = fcmp oeq double %87, 0.000000e+00
br i1 %88, label %89, label %92
89: ; preds = %70
%90 = tail call i32 @nrerror(i8* getelementptr inbounds ([24 x i8], [24 x i8]* @.str, i64 0, i64 0))
%91 = load double, double* %73
br label %92
92: ; preds = %89, %70
%93 = phi double [ %91, %89 ], [ %74, %70 ]
%94 = fdiv double %77, %87
%95 = fmul double %94, %93
store double %95, double* %75
%96 = fmul double %86, %94
%97 = getelementptr inbounds double, double* %7, i64 %71
store double %96, double* %97
%98 = icmp eq i64 %72, %69
br i1 %98, label %99, label %70
99: ; preds = %92, %61
%100 = shl nsw i32 %64, 1
%101 = icmp slt i32 %100, %66
%102 = add nsw i32 %64, 1
%103 = sext i32 %102 to i64
%104 = getelementptr inbounds double, double* %7, i64 %103
%105 = sext i32 %64 to i64
%106 = getelementptr inbounds double, double* %8, i64 %105
%107 = xor i1 %101, true
%108 = sext i1 %107 to i32
%109 = add nsw i32 %64, %108
%110 = select i1 %101, double* %104, double* %106
%111 = load double, double* %110
store double %111, double* %5
%112 = load double, double* %4
%113 = fadd double %111, %112
store double %113, double* %4
%114 = add nuw nsw i64 %62, 1
%115 = add i32 %63, -1
%116 = icmp eq i64 %114, %60
br i1 %116, label %117, label %61
117: ; preds = %99, %53, %12, %26
%118 = tail call i32 @free_vector(double* %8, i32 1, i32 %2)
%119 = tail call i32 @free_vector(double* %7, i32 1, i32 %2)
ret void
}