#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 */
int circlestarintersection(double, double, double, double *, double *);
double sqrt(double);
int ellipsestarintersection(double r, double dcen, double a, double b,
                            double v[2], double w[2]) {
  int k;
  double wp, wn, t0, t1, t2, aobsq;
  if (a == b)
    return circlestarintersection(r, dcen, a, v, w);
  aobsq = a * a / (b * b);
  t1 = 4.0 * dcen * dcen - 4.0 * (1.0 - aobsq) * (a * a + dcen * dcen - r * r);
  if (t1 < 0)
    return 0;
  t1 = sqrt(t1);
  t0 = (-2.0) * dcen;
  t2 = 2.0 * (1.0 - aobsq);
  wp = (t0 + t1) / t2;
  wn = (t0 - t1) / t2;
  k = 0;
  if (wn > (-b) && wn < b) {
    w[k] = wn;
    v[k] = sqrt(a * a - aobsq * wn * wn);
    k++;
  }
  if (wp > (-b) && wp < b) {
    w[k] = wp;
    v[k] = sqrt(a * a - aobsq * wp * wp);
    k++;
  }
  if (k == 2 && w[0] > w[1]) {
    t0 = v[0];
    t1 = w[0];
    v[0] = v[1];
    w[0] = w[1];
    v[1] = t0;
    w[1] = t1;
  }
  return k;
}


x86 assembly (clang, -O0):
<ellipsestarintersection>:
push   %rbp
mov    %rsp,%rbp
sub    $0x70,%rsp
movsd  %xmm0,-0x10(%rbp)
movsd  %xmm1,-0x18(%rbp)
movsd  %xmm2,-0x20(%rbp)
movsd  %xmm3,-0x28(%rbp)
mov    %rdi,-0x30(%rbp)
mov    %rsi,-0x38(%rbp)
movsd  -0x20(%rbp),%xmm0
ucomisd -0x28(%rbp),%xmm0
jne    117e <ellipsestarintersection+0x5e>
jp     117e <ellipsestarintersection+0x5e>
movsd  -0x10(%rbp),%xmm0
movsd  -0x18(%rbp),%xmm1
movsd  -0x20(%rbp),%xmm2
mov    -0x30(%rbp),%rdi
mov    -0x38(%rbp),%rsi
call   1030 <circlestarintersection@plt>
mov    %eax,-0x4(%rbp)
jmp    1482 <ellipsestarintersection+0x362>
movsd  -0x20(%rbp),%xmm0
mulsd  -0x20(%rbp),%xmm0
movsd  -0x28(%rbp),%xmm1
mulsd  -0x28(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  0xe5d(%rip),%xmm0
mulsd  -0x18(%rbp),%xmm0
movsd  -0x18(%rbp),%xmm2
movsd  0xe53(%rip),%xmm3
subsd  -0x70(%rbp),%xmm3
movsd  0xe3e(%rip),%xmm1
mulsd  %xmm3,%xmm1
movsd  -0x20(%rbp),%xmm4
movsd  -0x20(%rbp),%xmm5
movsd  -0x18(%rbp),%xmm3
mulsd  -0x18(%rbp),%xmm3
mulsd  %xmm5,%xmm4
addsd  %xmm3,%xmm4
movsd  -0x10(%rbp),%xmm3
movsd  -0x10(%rbp),%xmm5
movq   %xmm3,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm3
mulsd  %xmm5,%xmm3
addsd  %xmm4,%xmm3
mulsd  %xmm3,%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp)
xorps  %xmm0,%xmm0
ucomisd -0x60(%rbp),%xmm0
jbe    124d <ellipsestarintersection+0x12d>
movl   $0x0,-0x4(%rbp)
jmp    1482 <ellipsestarintersection+0x362>
movsd  -0x60(%rbp),%xmm0
call   1040 <sqrt@plt>
movsd  %xmm0,-0x60(%rbp)
movsd  0xdb4(%rip),%xmm0
mulsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  0xd92(%rip),%xmm1
subsd  -0x70(%rbp),%xmm1
movsd  0xd8d(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  -0x58(%rbp),%xmm0
addsd  -0x60(%rbp),%xmm0
divsd  -0x68(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x58(%rbp),%xmm0
subsd  -0x60(%rbp),%xmm0
divsd  -0x68(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
movl   $0x0,-0x3c(%rbp)
movsd  -0x50(%rbp),%xmm0
movsd  -0x28(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
ucomisd %xmm1,%xmm0
jbe    135f <ellipsestarintersection+0x23f>
movsd  -0x50(%rbp),%xmm1
movsd  -0x28(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    135f <ellipsestarintersection+0x23f>
movsd  -0x50(%rbp),%xmm0
mov    -0x38(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x20(%rbp),%xmm0
movsd  -0x20(%rbp),%xmm2
movsd  -0x70(%rbp),%xmm1
mulsd  -0x50(%rbp),%xmm1
mulsd  -0x50(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
call   1040 <sqrt@plt>
mov    -0x30(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
movsd  -0x48(%rbp),%xmm0
movsd  -0x28(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
ucomisd %xmm1,%xmm0
jbe    1403 <ellipsestarintersection+0x2e3>
movsd  -0x48(%rbp),%xmm1
movsd  -0x28(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    1403 <ellipsestarintersection+0x2e3>
movsd  -0x48(%rbp),%xmm0
mov    -0x38(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x20(%rbp),%xmm0
movsd  -0x20(%rbp),%xmm2
movsd  -0x70(%rbp),%xmm1
mulsd  -0x48(%rbp),%xmm1
mulsd  -0x48(%rbp),%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
call   1040 <sqrt@plt>
mov    -0x30(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
cmpl   $0x2,-0x3c(%rbp)
jne    147c <ellipsestarintersection+0x35c>
mov    -0x38(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x38(%rbp),%rax
ucomisd 0x8(%rax),%xmm0
jbe    147c <ellipsestarintersection+0x35c>
mov    -0x30(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x38(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    -0x30(%rbp),%rax
movsd  0x8(%rax),%xmm0
mov    -0x30(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x38(%rbp),%rax
movsd  0x8(%rax),%xmm0
mov    -0x38(%rbp),%rax
movsd  %xmm0,(%rax)
movsd  -0x58(%rbp),%xmm0
mov    -0x30(%rbp),%rax
movsd  %xmm0,0x8(%rax)
movsd  -0x60(%rbp),%xmm0
mov    -0x38(%rbp),%rax
movsd  %xmm0,0x8(%rax)
mov    -0x3c(%rbp),%eax
mov    %eax,-0x4(%rbp)
mov    -0x4(%rbp),%eax
add    $0x70,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<ellipsestarintersection>:
push   %rbp
push   %r15
push   %r14
push   %rbx
sub    $0x48,%rsp
mov    %rsi,%rbx
mov    %rdi,%r14
movapd %xmm0,%xmm6
ucomisd %xmm3,%xmm2
jne    1155 <ellipsestarintersection+0x35>
jp     1155 <ellipsestarintersection+0x35>
movapd %xmm6,%xmm0
mov    %r14,%rdi
mov    %rbx,%rsi
add    $0x48,%rsp
pop    %rbx
pop    %r14
pop    %r15
pop    %rbp
jmp    1030 <circlestarintersection@plt>
mulsd  %xmm2,%xmm2
movapd %xmm3,%xmm0
mulsd  %xmm3,%xmm0
movapd %xmm2,%xmm5
divsd  %xmm0,%xmm5
movsd  0xe8f(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  0xe8b(%rip),%xmm7
subsd  %xmm5,%xmm7
movapd %xmm1,%xmm4
mulsd  %xmm1,%xmm4
addsd  %xmm2,%xmm4
mulsd  %xmm6,%xmm6
subsd  %xmm6,%xmm4
movsd  0xe73(%rip),%xmm6
mulsd  %xmm7,%xmm6
mulsd  %xmm4,%xmm6
mulsd  %xmm1,%xmm0
addsd  %xmm6,%xmm0
xor    %eax,%eax
xorpd  %xmm4,%xmm4
ucomisd %xmm0,%xmm4
ja     134a <ellipsestarintersection+0x22a>
ucomisd %xmm4,%xmm0
jb     11c9 <ellipsestarintersection+0xa9>
sqrtsd %xmm0,%xmm0
jmp    120a <ellipsestarintersection+0xea>
movsd  %xmm2,0x18(%rsp)
movapd %xmm3,0x30(%rsp)
movsd  %xmm5,0x10(%rsp)
movsd  %xmm1,0x8(%rsp)
movsd  %xmm7,0x20(%rsp)
call   1040 <sqrt@plt>
movsd  0x20(%rsp),%xmm7
movsd  0x8(%rsp),%xmm1
movsd  0x10(%rsp),%xmm5
movapd 0x30(%rsp),%xmm3
movsd  0x18(%rsp),%xmm2
mulsd  0xe06(%rip),%xmm1
addsd  %xmm7,%xmm7
movapd %xmm1,%xmm6
addsd  %xmm0,%xmm6
subsd  %xmm0,%xmm1
divsd  %xmm7,%xmm1
movapd 0xdf1(%rip),%xmm8
xorpd  %xmm3,%xmm8
ucomisd %xmm8,%xmm1
seta   %al
ucomisd %xmm1,%xmm3
seta   %bpl
divsd  %xmm7,%xmm6
and    %al,%bpl
xor    %eax,%eax
cmp    $0x1,%bpl
jne    12c8 <ellipsestarintersection+0x1a8>
movsd  %xmm1,(%rbx)
movapd %xmm1,%xmm4
mulsd  %xmm5,%xmm4
mulsd  %xmm1,%xmm4
movapd %xmm2,%xmm0
subsd  %xmm4,%xmm0
xorpd  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jb     127b <ellipsestarintersection+0x15b>
sqrtsd %xmm0,%xmm0
jmp    12be <ellipsestarintersection+0x19e>
movsd  %xmm2,0x18(%rsp)
movapd %xmm3,0x30(%rsp)
movsd  %xmm5,0x10(%rsp)
movsd  %xmm6,0x8(%rsp)
movapd %xmm8,0x20(%rsp)
call   1040 <sqrt@plt>
movapd 0x20(%rsp),%xmm8
movsd  0x8(%rsp),%xmm6
movsd  0x10(%rsp),%xmm5
movapd 0x30(%rsp),%xmm3
movsd  0x18(%rsp),%xmm2
movsd  %xmm0,(%r14)
mov    $0x1,%eax
ucomisd %xmm8,%xmm6
jbe    134a <ellipsestarintersection+0x22a>
ucomisd %xmm6,%xmm3
jbe    134a <ellipsestarintersection+0x22a>
mov    %eax,%r15d
movsd  %xmm6,(%rbx,%r15,8)
mulsd  %xmm6,%xmm5
mulsd  %xmm6,%xmm5
subsd  %xmm5,%xmm2
xorpd  %xmm0,%xmm0
ucomisd %xmm0,%xmm2
jb     12fd <ellipsestarintersection+0x1dd>
xorps  %xmm0,%xmm0
sqrtsd %xmm2,%xmm0
jmp    1306 <ellipsestarintersection+0x1e6>
movapd %xmm2,%xmm0
call   1040 <sqrt@plt>
movsd  %xmm0,(%r14,%r15,8)
mov    $0x1,%eax
test   %bpl,%bpl
je     134a <ellipsestarintersection+0x22a>
movsd  (%rbx),%xmm0
mov    $0x2,%eax
ucomisd 0x8(%rbx),%xmm0
jbe    134a <ellipsestarintersection+0x22a>
movsd  (%r14),%xmm1
movsd  0x8(%r14),%xmm2
movsd  %xmm2,(%r14)
movsd  0x8(%rbx),%xmm2
movsd  %xmm2,(%rbx)
movsd  %xmm1,0x8(%r14)
movsd  %xmm0,0x8(%rbx)
add    $0x48,%rsp
pop    %rbx
pop    %r14
pop    %r15
pop    %rbp
ret

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

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