#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 M_PI;
double bondi_rel_as(double, double);
double bondi_rel_solve(double, double, double, double, double, double, double);
double pow(double, double);
int /*<<< orphan*/ printf(char *, double, double, double, ...);
double sqrt(double);
int bondi_rel(double Mdot, double GM, double gam, double a0, double *r,
              double *rho, double *u, double *P, int N) {
  double as = bondi_rel_as(a0, gam);
  double us2 = as * as / (1 + 3 * as * as);
  double rs = GM / (2 * us2);
  double us = -sqrt(us2);
  double rhos = -Mdot / (4 * M_PI * rs * rs * us);
  double K = as * as / ((1 - as * as / (gam - 1)) * pow(rhos, gam - 1) * gam);
  printf("M: %.8lg rs: %.8lg us: %.8lg as: %.8lg\n", GM, rs, us, as);
  double rho0 =
      pow(a0 * a0 / (gam * K * (1.0 - a0 * a0 / (gam - 1))), 1.0 / (gam - 1));
  double h0 = 1.0 + gam / (gam - 1) * K * pow(rho0, gam - 1);
  printf("rho0: %.8lg h0: %.8lg k: %.8lg\n", rho0, h0, K);
  int i;
  for (i = 0; i < N; i++) {
    double ui = bondi_rel_solve(Mdot, GM, r[i], gam, h0, K, as);
    rho[i] = -Mdot / (4 * M_PI * r[i] * r[i] * ui);
    u[i] = ui;
    P[i] = K * pow(rho[i], gam);
  }
  return 0;
}


x86 assembly (clang, -O0):
<bondi_rel>:
push   %rbp
mov    %rsp,%rbp
sub    $0xc0,%rsp
movsd  %xmm0,-0x8(%rbp)
movsd  %xmm1,-0x10(%rbp)
movsd  %xmm2,-0x18(%rbp)
movsd  %xmm3,-0x20(%rbp)
mov    %rdi,-0x28(%rbp)
mov    %rsi,-0x30(%rbp)
mov    %rdx,-0x38(%rbp)
mov    %rcx,-0x40(%rbp)
mov    %r8d,-0x44(%rbp)
movsd  -0x20(%rbp),%xmm0
movsd  -0x18(%rbp),%xmm1
call   1040 <bondi_rel_as@plt>
movsd  %xmm0,-0x50(%rbp)
movsd  -0x50(%rbp),%xmm0
mulsd  -0x50(%rbp),%xmm0
movsd  0xe67(%rip),%xmm1
mulsd  -0x50(%rbp),%xmm1
movsd  -0x50(%rbp),%xmm2
mulsd  %xmm2,%xmm1
movsd  0xe41(%rip),%xmm2
addsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x10(%rbp),%xmm0
movsd  0xe2f(%rip),%xmm1
mulsd  -0x58(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  -0x58(%rbp),%xmm0
call   1070 <sqrt@plt>
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  -0x8(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    0x2db8(%rip),%rax
mov    (%rax),%eax
shl    $0x2,%eax
cvtsi2sd %eax,%xmm1
mulsd  -0x60(%rbp),%xmm1
mulsd  -0x60(%rbp),%xmm1
mulsd  -0x68(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x50(%rbp),%xmm0
mulsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0xa8(%rbp)
ff
movsd  -0x50(%rbp),%xmm1
mulsd  -0x50(%rbp),%xmm1
movsd  -0x18(%rbp),%xmm0
movsd  0xd86(%rip),%xmm2
subsd  %xmm2,%xmm0
divsd  %xmm0,%xmm1
movsd  0xd76(%rip),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,-0xb0(%rbp)
ff
movsd  -0x70(%rbp),%xmm0
movsd  -0x18(%rbp),%xmm1
movsd  0xd58(%rip),%xmm2
subsd  %xmm2,%xmm1
call   1060 <pow@plt>
movsd  -0xb0(%rbp),%xmm1
ff
movaps %xmm0,%xmm2
movsd  -0xa8(%rbp),%xmm0
ff
mulsd  %xmm2,%xmm1
mulsd  -0x18(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x78(%rbp)
movsd  -0x10(%rbp),%xmm0
movsd  -0x60(%rbp),%xmm1
movsd  -0x68(%rbp),%xmm2
movsd  -0x50(%rbp),%xmm3
lea    0xd27(%rip),%rdi
mov    $0x4,%al
call   1050 <printf@plt>
movsd  -0x20(%rbp),%xmm0
mulsd  -0x20(%rbp),%xmm0
movsd  -0x18(%rbp),%xmm1
mulsd  -0x78(%rbp),%xmm1
movsd  -0x20(%rbp),%xmm3
mulsd  -0x20(%rbp),%xmm3
movsd  -0x18(%rbp),%xmm2
movsd  0xcdd(%rip),%xmm4
subsd  %xmm4,%xmm2
divsd  %xmm2,%xmm3
movsd  0xccd(%rip),%xmm2
subsd  %xmm3,%xmm2
mulsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
movsd  -0x18(%rbp),%xmm2
movsd  0xcb4(%rip),%xmm1
subsd  %xmm1,%xmm2
movsd  0xca8(%rip),%xmm1
divsd  %xmm2,%xmm1
call   1060 <pow@plt>
movsd  %xmm0,-0x80(%rbp)
movsd  -0x18(%rbp),%xmm0
movsd  -0x18(%rbp),%xmm1
movsd  0xc88(%rip),%xmm2
subsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
mulsd  -0x78(%rbp),%xmm0
movsd  %xmm0,-0xa0(%rbp)
ff
movsd  -0x80(%rbp),%xmm0
movsd  -0x18(%rbp),%xmm1
movsd  0xc61(%rip),%xmm2
subsd  %xmm2,%xmm1
call   1060 <pow@plt>
movaps %xmm0,%xmm1
movsd  -0xa0(%rbp),%xmm0
ff
mulsd  %xmm1,%xmm0
movsd  0xc41(%rip),%xmm1
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
movsd  -0x80(%rbp),%xmm0
movsd  -0x88(%rbp),%xmm1
ff
movsd  -0x78(%rbp),%xmm2
lea    0xc5c(%rip),%rdi
mov    $0x3,%al
call   1050 <printf@plt>
movl   $0x0,-0x8c(%rbp)
mov    -0x8c(%rbp),%eax
cmp    -0x44(%rbp),%eax
jge    151c <bondi_rel+0x3cc>
movsd  -0x8(%rbp),%xmm0
movsd  -0x10(%rbp),%xmm1
mov    -0x28(%rbp),%rax
movslq -0x8c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x18(%rbp),%xmm3
movsd  -0x88(%rbp),%xmm4
ff
movsd  -0x78(%rbp),%xmm5
movsd  -0x50(%rbp),%xmm6
call   1030 <bondi_rel_solve@plt>
movsd  %xmm0,-0x98(%rbp)
ff
movsd  -0x8(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mov    0x2b83(%rip),%rax
mov    (%rax),%eax
shl    $0x2,%eax
cvtsi2sd %eax,%xmm1
mov    -0x28(%rbp),%rax
movslq -0x8c(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm1
mov    -0x28(%rbp),%rax
movslq -0x8c(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm1
mulsd  -0x98(%rbp),%xmm1
ff
divsd  %xmm1,%xmm0
mov    -0x30(%rbp),%rax
movslq -0x8c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x98(%rbp),%xmm0
ff
mov    -0x38(%rbp),%rax
movslq -0x8c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x78(%rbp),%xmm0
movsd  %xmm0,-0xb8(%rbp)
ff
mov    -0x30(%rbp),%rax
movslq -0x8c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x18(%rbp),%xmm1
call   1060 <pow@plt>
movaps %xmm0,%xmm1
movsd  -0xb8(%rbp),%xmm0
ff
mulsd  %xmm1,%xmm0
mov    -0x40(%rbp),%rax
movslq -0x8c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x8c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x8c(%rbp)
jmp    13f5 <bondi_rel+0x2a5>
xor    %eax,%eax
add    $0xc0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<bondi_rel>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0xa8,%rsp
mov    %r8d,%ebp
mov    %rcx,%r14
mov    %rdx,%r15
mov    %rsi,%r12
mov    %rdi,%rbx
movapd %xmm3,0x10(%rsp)
movsd  %xmm1,0x38(%rsp)
movaps %xmm0,0x70(%rsp)
movapd %xmm3,%xmm0
movsd  %xmm2,0x8(%rsp)
movaps %xmm2,%xmm1
call   1040 <bondi_rel_as@plt>
movsd  0xe65(%rip),%xmm1
mulsd  %xmm0,%xmm1
mulsd  %xmm0,%xmm1
addsd  0xe5d(%rip),%xmm1
movapd %xmm0,0x90(%rsp)
movapd 0x10(%rsp),%xmm3
unpcklpd %xmm0,%xmm3
mulpd  %xmm3,%xmm3
movapd %xmm3,0x10(%rsp)
unpckhpd %xmm3,%xmm3
movapd %xmm3,%xmm0
divsd  %xmm1,%xmm0
movapd %xmm0,%xmm1
addsd  %xmm0,%xmm1
movsd  0x38(%rsp),%xmm2
divsd  %xmm1,%xmm2
movsd  %xmm2,0x20(%rsp)
xorpd  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jb     11fc <bondi_rel+0xac>
sqrtsd %xmm0,%xmm0
jmp    120d <bondi_rel+0xbd>
movapd %xmm3,0x40(%rsp)
call   1070 <sqrt@plt>
movapd 0x40(%rsp),%xmm3
movapd 0xe0b(%rip),%xmm2
mov    0x2dcc(%rip),%rax
mov    (%rax),%eax
shl    $0x2,%eax
xorps  %xmm1,%xmm1
cvtsi2sd %eax,%xmm1
movsd  0x20(%rsp),%xmm4
mulsd  %xmm4,%xmm1
mulsd  %xmm4,%xmm1
mulsd  %xmm0,%xmm1
movapd %xmm2,0x40(%rsp)
xorpd  %xmm2,%xmm0
movapd %xmm0,0x60(%rsp)
movapd 0x70(%rsp),%xmm0
divsd  %xmm1,%xmm0
movsd  0xdb4(%rip),%xmm1
addsd  0x8(%rsp),%xmm1
movapd %xmm1,0x80(%rsp)
divsd  %xmm1,%xmm3
movsd  0xd91(%rip),%xmm2
subsd  %xmm3,%xmm2
movsd  %xmm2,0x50(%rsp)
call   1060 <pow@plt>
mulsd  0x50(%rsp),%xmm0
mulsd  0x8(%rsp),%xmm0
movapd %xmm0,0x50(%rsp)
lea    0xd91(%rip),%rdi
movsd  0x38(%rsp),%xmm0
movsd  0x20(%rsp),%xmm1
movaps 0x60(%rsp),%xmm2
movaps 0x90(%rsp),%xmm3
mov    $0x4,%al
call   1050 <printf@plt>
movapd 0x80(%rsp),%xmm5
movapd %xmm5,%xmm1
unpcklpd 0x50(%rsp),%xmm1
movapd 0x10(%rsp),%xmm0
movapd %xmm0,%xmm4
divpd  %xmm1,%xmm4
movsd  0xd20(%rip),%xmm3
movapd %xmm3,%xmm2
subsd  %xmm4,%xmm2
unpckhpd %xmm4,%xmm4
movapd %xmm4,0x20(%rsp)
movapd %xmm4,%xmm1
mulsd  0x8(%rsp),%xmm1
mulsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
movapd %xmm3,%xmm1
divsd  %xmm5,%xmm1
call   1060 <pow@plt>
movsd  %xmm0,0x10(%rsp)
movsd  0x8(%rsp),%xmm2
movapd 0x80(%rsp),%xmm1
divsd  %xmm1,%xmm2
mulsd  0x20(%rsp),%xmm2
movsd  %xmm2,0x60(%rsp)
call   1060 <pow@plt>
movapd %xmm0,%xmm1
mulsd  0x60(%rsp),%xmm1
addsd  0xcb3(%rip),%xmm1
lea    0xcfc(%rip),%rdi
movsd  0x10(%rsp),%xmm0
movsd  %xmm1,0x10(%rsp)
movapd 0x20(%rsp),%xmm2
mov    $0x3,%al
call   1050 <printf@plt>
movapd 0x70(%rsp),%xmm0
test   %ebp,%ebp
jle    143a <bondi_rel+0x2ea>
movapd 0x40(%rsp),%xmm1
xorpd  %xmm0,%xmm1
movapd %xmm1,0x40(%rsp)
mov    %ebp,%ebp
xor    %r13d,%r13d
nopl   0x0(%rax,%rax,1)
movsd  (%rbx,%r13,8),%xmm2
movapd 0x70(%rsp),%xmm0
movsd  0x38(%rsp),%xmm1
movsd  0x8(%rsp),%xmm3
movsd  0x10(%rsp),%xmm4
movapd 0x20(%rsp),%xmm5
movaps 0x90(%rsp),%xmm6
call   1030 <bondi_rel_solve@plt>
mov    0x2c10(%rip),%rax
mov    (%rax),%eax
shl    $0x2,%eax
xorps  %xmm1,%xmm1
cvtsi2sd %eax,%xmm1
movsd  (%rbx,%r13,8),%xmm2
mulsd  %xmm2,%xmm1
mulsd  %xmm2,%xmm1
mulsd  %xmm0,%xmm1
movapd 0x40(%rsp),%xmm2
divsd  %xmm1,%xmm2
movsd  %xmm2,(%r12,%r13,8)
movsd  %xmm0,(%r15,%r13,8)
movsd  (%r12,%r13,8),%xmm0
movsd  0x8(%rsp),%xmm1
call   1060 <pow@plt>
movapd 0x20(%rsp),%xmm2
mulsd  %xmm2,%xmm0
movsd  %xmm0,(%r14,%r13,8)
add    $0x1,%r13
cmp    %r13,%rbp
jne    13a0 <bondi_rel+0x250>
xor    %eax,%eax
add    $0xa8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
@M_PI = global i32 0
@.str = [40 x i8] c"M: %.8lg rs: %.8lg us: %.8lg as: %.8lg\0A\00"
@.str.1 = [32 x i8] c"rho0: %.8lg h0: %.8lg k: %.8lg\0A\00"
define i32 @bondi_rel(double %0, double %1, double %2, double %3, double* nocapture %4, double* nocapture %5, double* nocapture writeonly %6, double* nocapture writeonly %7, i32 %8) {
%10 = tail call double @bondi_rel_as(double %3, double %2)
%11 = fmul double %10, 3.000000e+00
%12 = tail call double @llvm.fmuladd.f64(double %11, double %10, double 1.000000e+00)
%13 = insertelement <2 x double> poison, double %3, i64 0
%14 = insertelement <2 x double> %13, double %10, i64 1
%15 = fmul <2 x double> %14, %14
%16 = extractelement <2 x double> %15, i64 1
%17 = fdiv double %16, %12
%18 = fmul double %17, 2.000000e+00
%19 = fdiv double %1, %18
%20 = tail call double @sqrt(double %17)
%21 = fneg double %20
%22 = fneg double %0
%23 = load i32, i32* @M_PI
%24 = shl nsw i32 %23, 2
%25 = sitofp i32 %24 to double
%26 = fmul double %19, %25
%27 = fmul double %19, %26
%28 = fmul double %27, %21
%29 = fdiv double %22, %28
%30 = fadd double %2, -1.000000e+00
%31 = fdiv double %16, %30
%32 = fsub double 1.000000e+00, %31
%33 = tail call double @pow(double %29, double %30)
%34 = fmul double %32, %33
%35 = fmul double %34, %2
%36 = tail call i32 (i8*, double, double, double, ...) @printf(i8* getelementptr inbounds ([40 x i8], [40 x i8]* @.str, i64 0, i64 0), double %1, double %19, double %21, double %10)
%37 = insertelement <2 x double> poison, double %30, i64 0
%38 = insertelement <2 x double> %37, double %35, i64 1
%39 = fdiv <2 x double> %15, %38
%40 = extractelement <2 x double> %39, i64 0
%41 = fsub double 1.000000e+00, %40
%42 = extractelement <2 x double> %39, i64 1
%43 = fmul double %42, %2
%44 = fmul double %41, %43
%45 = extractelement <2 x double> %15, i64 0
%46 = fdiv double %45, %44
%47 = fdiv double 1.000000e+00, %30
%48 = tail call double @pow(double %46, double %47)
%49 = fdiv double %2, %30
%50 = extractelement <2 x double> %39, i64 1
%51 = fmul double %49, %50
%52 = tail call double @pow(double %48, double %30)
%53 = tail call double @llvm.fmuladd.f64(double %51, double %52, double 1.000000e+00)
%54 = tail call i32 (i8*, double, double, double, ...) @printf(i8* getelementptr inbounds ([32 x i8], [32 x i8]* @.str.1, i64 0, i64 0), double %48, double %53, double %50)
%55 = icmp sgt i32 %8, 0
br i1 %55, label %56, label %79
56: ; preds = %9
%57 = zext i32 %8 to i64
br label %58
58: ; preds = %56, %58
%59 = phi i64 [ 0, %56 ], [ %77, %58 ]
%60 = getelementptr inbounds double, double* %4, i64 %59
%61 = load double, double* %60
%62 = tail call double @bondi_rel_solve(double %0, double %1, double %61, double %2, double %53, double %50, double %10)
%63 = load i32, i32* @M_PI
%64 = shl nsw i32 %63, 2
%65 = sitofp i32 %64 to double
%66 = load double, double* %60
%67 = fmul double %66, %65
%68 = fmul double %66, %67
%69 = fmul double %62, %68
%70 = fdiv double %22, %69
%71 = getelementptr inbounds double, double* %5, i64 %59
store double %70, double* %71
%72 = getelementptr inbounds double, double* %6, i64 %59
store double %62, double* %72
%73 = load double, double* %71
%74 = tail call double @pow(double %73, double %2)
%75 = fmul double %50, %74
%76 = getelementptr inbounds double, double* %7, i64 %59
store double %75, double* %76
%77 = add nuw nsw i64 %59, 1
%78 = icmp eq i64 %77, %57
br i1 %78, label %79, label %58
79: ; preds = %58, %9
ret i32 0
}