#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 sqrt(double);
void bslipds(unsigned n, double t, double *Y, double *p, double *f) {
  double l1, l2, F1, Fx1, Fy1, Fz1, Fx2, Fy2, Fz2, F2;
  double xr1, zr1, xr2, zr2;
  xr1 = Y[0] - p[10];
  zr1 = Y[2] - p[12];
  l1 = sqrt(xr1 * xr1 + Y[1] * Y[1] + zr1 * zr1);
  if (p[16] == 1.) {
    F1 = -p[0] * (l1 - p[4]);
  } else {
    F1 = -p[1] * (l1 - p[5]);
  }
  if (F1 < 0.) {
    F1 = 0.;
  }
  xr2 = Y[0] - p[13];
  zr2 = Y[2] - p[15];
  l2 = sqrt(xr2 * xr2 + Y[1] * Y[1] + zr2 * zr2);
  if (p[16] == 1.) {
    F2 = -p[1] * (l2 - p[5]);
  } else {
    F2 = -p[0] * (l2 - p[4]);
  }
  if (F2 < 0.) {
    F2 = 0.;
  }
  Fx1 = xr1 / l1 * F1;
  Fy1 = Y[1] / l1 * F1;
  Fz1 = zr1 / l1 * F1;
  Fx2 = xr2 / l2 * F2;
  Fy2 = Y[1] / l2 * F2;
  Fz2 = zr2 / l2 * F2;
  f[0] = Y[3];
  f[1] = Y[4];
  f[2] = Y[5];
  f[3] = (Fx1 + Fx2) / p[8];
  f[4] = (Fy1 + Fy2) / p[8] + p[9];
  f[5] = (Fz1 + Fz2) / p[8];
}


x86 assembly (clang, -O0):
<bslipds>:
push   %rbp
mov    %rsp,%rbp
sub    $0xa0,%rsp
mov    %edi,-0x4(%rbp)
movsd  %xmm0,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %rdx,-0x20(%rbp)
mov    %rcx,-0x28(%rbp)
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x20(%rbp),%rax
subsd  0x50(%rax),%xmm0
movsd  %xmm0,-0x80(%rbp)
mov    -0x18(%rbp),%rax
movsd  0x10(%rax),%xmm0
mov    -0x20(%rbp),%rax
subsd  0x60(%rax),%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
movsd  -0x80(%rbp),%xmm1
movsd  -0x80(%rbp),%xmm2
mov    -0x18(%rbp),%rax
movsd  0x8(%rax),%xmm0
mov    -0x18(%rbp),%rax
mulsd  0x8(%rax),%xmm0
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
movsd  -0x88(%rbp),%xmm0
ff
movsd  -0x88(%rbp),%xmm2
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
call   1030 <sqrt@plt>
movsd  %xmm0,-0x30(%rbp)
mov    -0x20(%rbp),%rax
movsd  0x80(%rax),%xmm0
movsd  0xe47(%rip),%xmm1
ucomisd %xmm1,%xmm0
jne    1204 <bslipds+0xf4>
jp     1204 <bslipds+0xf4>
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  -0x30(%rbp),%xmm1
mov    -0x20(%rbp),%rax
subsd  0x20(%rax),%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
jmp    123b <bslipds+0x12b>
mov    -0x20(%rbp),%rax
movsd  0x8(%rax),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  -0x30(%rbp),%xmm1
mov    -0x20(%rbp),%rax
subsd  0x28(%rax),%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
xorps  %xmm0,%xmm0
ucomisd -0x40(%rbp),%xmm0
jbe    1251 <bslipds+0x141>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x18(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x20(%rbp),%rax
subsd  0x68(%rax),%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
mov    -0x18(%rbp),%rax
movsd  0x10(%rax),%xmm0
mov    -0x20(%rbp),%rax
subsd  0x78(%rax),%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
movsd  -0x90(%rbp),%xmm1
ff
movsd  -0x90(%rbp),%xmm2
ff
mov    -0x18(%rbp),%rax
movsd  0x8(%rax),%xmm0
mov    -0x18(%rbp),%rax
mulsd  0x8(%rax),%xmm0
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
movsd  -0x98(%rbp),%xmm0
ff
movsd  -0x98(%rbp),%xmm2
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
call   1030 <sqrt@plt>
movsd  %xmm0,-0x38(%rbp)
mov    -0x20(%rbp),%rax
movsd  0x80(%rax),%xmm0
movsd  0xd1c(%rip),%xmm1
ucomisd %xmm1,%xmm0
jne    1330 <bslipds+0x220>
jp     1330 <bslipds+0x220>
mov    -0x20(%rbp),%rax
movsd  0x8(%rax),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  -0x38(%rbp),%xmm1
mov    -0x20(%rbp),%rax
subsd  0x28(%rax),%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x78(%rbp)
jmp    1366 <bslipds+0x256>
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  -0x38(%rbp),%xmm1
mov    -0x20(%rbp),%rax
subsd  0x20(%rax),%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x78(%rbp)
xorps  %xmm0,%xmm0
ucomisd -0x78(%rbp),%xmm0
jbe    137c <bslipds+0x26c>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x78(%rbp)
movsd  -0x80(%rbp),%xmm0
divsd  -0x30(%rbp),%xmm0
mulsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x18(%rbp),%rax
movsd  0x8(%rax),%xmm0
divsd  -0x30(%rbp),%xmm0
mulsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x88(%rbp),%xmm0
ff
divsd  -0x30(%rbp),%xmm0
mulsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x90(%rbp),%xmm0
ff
divsd  -0x38(%rbp),%xmm0
mulsd  -0x78(%rbp),%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    -0x18(%rbp),%rax
movsd  0x8(%rax),%xmm0
divsd  -0x38(%rbp),%xmm0
mulsd  -0x78(%rbp),%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  -0x98(%rbp),%xmm0
ff
divsd  -0x38(%rbp),%xmm0
mulsd  -0x78(%rbp),%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    -0x18(%rbp),%rax
movsd  0x18(%rax),%xmm0
mov    -0x28(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x18(%rbp),%rax
movsd  0x20(%rax),%xmm0
mov    -0x28(%rbp),%rax
movsd  %xmm0,0x8(%rax)
mov    -0x18(%rbp),%rax
movsd  0x28(%rax),%xmm0
mov    -0x28(%rbp),%rax
movsd  %xmm0,0x10(%rax)
movsd  -0x48(%rbp),%xmm0
addsd  -0x60(%rbp),%xmm0
mov    -0x20(%rbp),%rax
divsd  0x40(%rax),%xmm0
mov    -0x28(%rbp),%rax
movsd  %xmm0,0x18(%rax)
movsd  -0x50(%rbp),%xmm0
addsd  -0x68(%rbp),%xmm0
mov    -0x20(%rbp),%rax
divsd  0x40(%rax),%xmm0
mov    -0x20(%rbp),%rax
addsd  0x48(%rax),%xmm0
mov    -0x28(%rbp),%rax
movsd  %xmm0,0x20(%rax)
movsd  -0x58(%rbp),%xmm0
addsd  -0x70(%rbp),%xmm0
mov    -0x20(%rbp),%rax
divsd  0x40(%rax),%xmm0
mov    -0x28(%rbp),%rax
movsd  %xmm0,0x28(%rax)
add    $0xa0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<bslipds>:
push   %r15
push   %r14
push   %rbx
sub    $0x30,%rsp
mov    %rcx,%r14
mov    %rdx,%rbx
mov    %rsi,%r15
movsd  (%rsi),%xmm1
subsd  0x50(%rdx),%xmm1
movsd  0x8(%rsi),%xmm0
movsd  0x10(%rsi),%xmm2
subsd  0x60(%rdx),%xmm2
mulsd  %xmm0,%xmm0
movsd  %xmm1,0x20(%rsp)
mulsd  %xmm1,%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm2,0x28(%rsp)
movapd %xmm2,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
xorpd  %xmm8,%xmm8
ucomisd %xmm8,%xmm0
jb     1171 <bslipds+0x61>
sqrtsd %xmm0,%xmm10
jmp    1180 <bslipds+0x70>
call   1030 <sqrt@plt>
xorpd  %xmm8,%xmm8
movapd %xmm0,%xmm10
movsd  0xe78(%rip),%xmm4
movsd  0x80(%rbx),%xmm0
cmpneqsd %xmm4,%xmm0
movq   %xmm0,%rax
and    $0x1,%eax
movsd  (%rbx,%rax,8),%xmm0
xorpd  0xe66(%rip),%xmm0
movapd %xmm10,%xmm1
subsd  0x20(%rbx,%rax,8),%xmm1
mulsd  %xmm0,%xmm1
movsd  (%r15),%xmm3
subsd  0x68(%rbx),%xmm3
movsd  0x8(%r15),%xmm0
movsd  0x10(%r15),%xmm7
subsd  0x78(%rbx),%xmm7
mulsd  %xmm0,%xmm0
movapd %xmm3,%xmm2
mulsd  %xmm3,%xmm2
addsd  %xmm0,%xmm2
movapd %xmm7,%xmm0
mulsd  %xmm7,%xmm0
addsd  %xmm2,%xmm0
ucomisd %xmm8,%xmm0
maxsd  %xmm1,%xmm8
jb     1202 <bslipds+0xf2>
sqrtsd %xmm0,%xmm0
jmp    1241 <bslipds+0x131>
movsd  %xmm8,0x18(%rsp)
movsd  %xmm7,0x10(%rsp)
movsd  %xmm3,0x8(%rsp)
movsd  %xmm10,(%rsp)
call   1030 <sqrt@plt>
movsd  0xdd8(%rip),%xmm4
movsd  (%rsp),%xmm10
movsd  0x8(%rsp),%xmm3
movsd  0x10(%rsp),%xmm7
movsd  0x18(%rsp),%xmm8
cmpeqsd 0x80(%rbx),%xmm4
movq   %xmm4,%rax
and    $0x1,%eax
movsd  (%rbx,%rax,8),%xmm1
xorpd  0xdb1(%rip),%xmm1
movapd %xmm0,%xmm2
subsd  0x20(%rbx,%rax,8),%xmm2
mulsd  %xmm1,%xmm2
xorpd  %xmm9,%xmm9
maxsd  %xmm2,%xmm9
movsd  0x20(%rsp),%xmm6
divsd  %xmm10,%xmm6
mulsd  %xmm8,%xmm6
movsd  0x8(%r15),%xmm1
movsd  0x18(%r15),%xmm2
movapd %xmm1,%xmm4
divsd  %xmm10,%xmm4
mulsd  %xmm8,%xmm4
movsd  0x28(%rsp),%xmm5
divsd  %xmm10,%xmm5
mulsd  %xmm8,%xmm5
divsd  %xmm0,%xmm3
mulsd  %xmm9,%xmm3
addsd  %xmm6,%xmm3
divsd  %xmm0,%xmm1
mulsd  %xmm9,%xmm1
addsd  %xmm4,%xmm1
divsd  %xmm0,%xmm7
mulsd  %xmm9,%xmm7
addsd  %xmm5,%xmm7
movsd  %xmm2,(%r14)
movsd  0x20(%r15),%xmm0
movsd  %xmm0,0x8(%r14)
movsd  0x28(%r15),%xmm0
movsd  %xmm0,0x10(%r14)
divsd  0x40(%rbx),%xmm3
movsd  %xmm3,0x18(%r14)
divsd  0x40(%rbx),%xmm1
addsd  0x48(%rbx),%xmm1
movsd  %xmm1,0x20(%r14)
divsd  0x40(%rbx),%xmm7
movsd  %xmm7,0x28(%r14)
add    $0x30,%rsp
pop    %rbx
pop    %r14
pop    %r15
ret

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

LLVM IR (-O2):
define void @bslipds(i32 %0, double %1, double* nocapture %2, double* nocapture %3, double* nocapture writeonly %4) {
%6 = load double, double* %2
%7 = getelementptr inbounds double, double* %3, i64 10
%8 = load double, double* %7
%9 = fsub double %6, %8
%10 = getelementptr inbounds double, double* %2, i64 2
%11 = load double, double* %10
%12 = getelementptr inbounds double, double* %3, i64 12
%13 = load double, double* %12
%14 = fsub double %11, %13
%15 = getelementptr inbounds double, double* %2, i64 1
%16 = load double, double* %15
%17 = fmul double %16, %16
%18 = tail call double @llvm.fmuladd.f64(double %9, double %9, double %17)
%19 = tail call double @llvm.fmuladd.f64(double %14, double %14, double %18)
%20 = tail call double @sqrt(double %19)
%21 = getelementptr inbounds double, double* %3, i64 16
%22 = load double, double* %21
%23 = fcmp oeq double %22, 1.000000e+00
%24 = xor i1 %23, true
%25 = zext i1 %24 to i64
%26 = getelementptr double, double* %3, i64 %25
%27 = select i1 %23, i64 4, i64 5
%28 = load double, double* %26
%29 = fneg double %28
%30 = getelementptr inbounds double, double* %3, i64 %27
%31 = load double, double* %30
%32 = fsub double %20, %31
%33 = fmul double %32, %29
%34 = fcmp olt double %33, 0.000000e+00
%35 = select i1 %34, double 0.000000e+00, double %33
%36 = load double, double* %2
%37 = getelementptr inbounds double, double* %3, i64 13
%38 = load double, double* %37
%39 = fsub double %36, %38
%40 = load double, double* %10
%41 = getelementptr inbounds double, double* %3, i64 15
%42 = load double, double* %41
%43 = fsub double %40, %42
%44 = load double, double* %15
%45 = fmul double %44, %44
%46 = tail call double @llvm.fmuladd.f64(double %39, double %39, double %45)
%47 = tail call double @llvm.fmuladd.f64(double %43, double %43, double %46)
%48 = tail call double @sqrt(double %47)
%49 = load double, double* %21
%50 = fcmp oeq double %49, 1.000000e+00
%51 = zext i1 %50 to i64
%52 = getelementptr double, double* %3, i64 %51
%53 = select i1 %50, i64 5, i64 4
%54 = load double, double* %52
%55 = fneg double %54
%56 = getelementptr inbounds double, double* %3, i64 %53
%57 = load double, double* %56
%58 = fsub double %48, %57
%59 = fmul double %58, %55
%60 = fcmp olt double %59, 0.000000e+00
%61 = select i1 %60, double 0.000000e+00, double %59
%62 = fdiv double %9, %20
%63 = fmul double %62, %35
%64 = load double, double* %15
%65 = fdiv double %64, %20
%66 = fmul double %35, %65
%67 = fdiv double %14, %20
%68 = fmul double %67, %35
%69 = fdiv double %39, %48
%70 = fmul double %69, %61
%71 = fdiv double %64, %48
%72 = fmul double %71, %61
%73 = fdiv double %43, %48
%74 = fmul double %73, %61
%75 = getelementptr inbounds double, double* %2, i64 3
%76 = load double, double* %75
store double %76, double* %4
%77 = getelementptr inbounds double, double* %2, i64 4
%78 = load double, double* %77
%79 = getelementptr inbounds double, double* %4, i64 1
store double %78, double* %79
%80 = getelementptr inbounds double, double* %2, i64 5
%81 = load double, double* %80
%82 = getelementptr inbounds double, double* %4, i64 2
store double %81, double* %82
%83 = fadd double %63, %70
%84 = getelementptr inbounds double, double* %3, i64 8
%85 = load double, double* %84
%86 = fdiv double %83, %85
%87 = getelementptr inbounds double, double* %4, i64 3
store double %86, double* %87
%88 = fadd double %66, %72
%89 = load double, double* %84
%90 = fdiv double %88, %89
%91 = getelementptr inbounds double, double* %3, i64 9
%92 = load double, double* %91
%93 = fadd double %90, %92
%94 = getelementptr inbounds double, double* %4, i64 4
store double %93, double* %94
%95 = fadd double %68, %74
%96 = load double, double* %84
%97 = fdiv double %95, %96
%98 = getelementptr inbounds double, double* %4, i64 5
store double %97, double* %98
ret void
}