#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 acos(int);
double cos(double);
double fabs(double);
double sgn(double);
double sin(double);
double EccentricAnomalyFast(double dMeanA, double dEcc, double tol,
                            int maxiter) {
  double dEccA;
  double di_1, di_2, di_3 = 1.0, fi, fi_1, fi_2, fi_3;
  double lo = -2 * acos(-1.);
  double up = 2 * acos(-1.);
  double next;
  int iter;
  if (dEcc == 0.)
    return dMeanA;
  dEccA = dMeanA + sgn(sin(dMeanA)) * 0.85 * dEcc;
  for (iter = 1; iter <= maxiter; iter++) {
    fi = dEccA - dEcc * sin(dEccA) - dMeanA;
    if (fi > 0)
      up = dEccA;
    else
      lo = dEccA;
    fi_1 = 1.0 - dEcc * cos(dEccA);
    fi_2 = dEcc * sin(dEccA);
    fi_3 = dEcc * cos(dEccA);
    di_1 = -fi / fi_1;
    di_2 = -fi / (fi_1 + 0.5 * di_1 * fi_2);
    di_3 = -fi / (fi_1 + 0.5 * di_2 * fi_2 + 1. / 6. * di_2 * di_2 * fi_3);
    next = dEccA + di_3;
    if (fabs(dEccA - next) < tol)
      break;
    if ((next > lo) && (next < up))
      dEccA = next;
    else
      dEccA = (lo + up) / 2.;
    if ((fabs(dEccA - lo) < tol) || (fabs(dEccA - up) < tol))
      break;
  }
  if (iter >= maxiter)
    return -1.;
  else
    return dEccA;
}


x86 assembly (clang, -O0):
<EccentricAnomalyFast>:
push   %rbp
mov    %rsp,%rbp
sub    $0xc0,%rsp
movsd  %xmm0,-0x10(%rbp)
movsd  %xmm1,-0x18(%rbp)
movsd  %xmm2,-0x20(%rbp)
mov    %edi,-0x24(%rbp)
movsd  0xe9b(%rip),%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    $0xffffffff,%edi
call   1030 <acos@plt>
imul   $0xfffffffe,%eax,%eax
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    $0xffffffff,%edi
call   1030 <acos@plt>
shl    $0x1,%eax
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x78(%rbp)
movsd  -0x18(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    11c0 <EccentricAnomalyFast+0x80>
jp     11c0 <EccentricAnomalyFast+0x80>
movsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    152e <EccentricAnomalyFast+0x3ee>
movsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
movsd  -0x10(%rbp),%xmm0
call   1040 <sin@plt>
call   1050 <sgn@plt>
movsd  -0x90(%rbp),%xmm1
ff
movsd  0xe1c(%rip),%xmm2
mulsd  %xmm2,%xmm0
movsd  -0x18(%rbp),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
movl   $0x1,-0x84(%rbp)
mov    -0x84(%rbp),%eax
cmp    -0x24(%rbp),%eax
jg     1503 <EccentricAnomalyFast+0x3c3>
movsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0xa0(%rbp)
ff
movsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
movsd  -0x30(%rbp),%xmm0
call   1040 <sin@plt>
movsd  -0xa0(%rbp),%xmm1
ff
movaps %xmm0,%xmm2
movsd  -0x98(%rbp),%xmm0
ff
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
subsd  -0x10(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x50(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    129c <EccentricAnomalyFast+0x15c>
movsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x78(%rbp)
jmp    12a6 <EccentricAnomalyFast+0x166>
movsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0xb8(%rbp)
ff
movsd  -0x30(%rbp),%xmm0
call   1060 <cos@plt>
movaps %xmm0,%xmm1
movsd  -0xb8(%rbp),%xmm0
ff
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm1,%xmm0
movsd  0xd15(%rip),%xmm1
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0xb0(%rbp)
movsd  -0x30(%rbp),%xmm0
call   1040 <sin@plt>
movaps %xmm0,%xmm1
movsd  -0xb0(%rbp),%xmm0
ff
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  -0x18(%rbp),%xmm0
movsd  %xmm0,-0xa8(%rbp)
ff
movsd  -0x30(%rbp),%xmm0
call   1060 <cos@plt>
movaps %xmm0,%xmm1
movsd  -0xa8(%rbp),%xmm0
ff
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  -0x50(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
divsd  -0x58(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x50(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  -0x58(%rbp),%xmm2
movsd  0xc7f(%rip),%xmm1
mulsd  -0x38(%rbp),%xmm1
movsd  -0x60(%rbp),%xmm3
mulsd  %xmm3,%xmm1
addsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x50(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  -0x58(%rbp),%xmm1
movsd  0xc3b(%rip),%xmm2
mulsd  -0x40(%rbp),%xmm2
movsd  -0x60(%rbp),%xmm3
mulsd  %xmm3,%xmm2
addsd  %xmm1,%xmm2
movsd  0xc19(%rip),%xmm1
mulsd  -0x40(%rbp),%xmm1
mulsd  -0x40(%rbp),%xmm1
movsd  -0x68(%rbp),%xmm3
mulsd  %xmm3,%xmm1
addsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x30(%rbp),%xmm0
addsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x80(%rbp)
movsd  -0x30(%rbp),%xmm1
subsd  -0x80(%rbp),%xmm1
movaps 0xbf9(%rip),%xmm0
pand   %xmm0,%xmm1
movsd  -0x20(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    144f <EccentricAnomalyFast+0x30f>
jmp    1503 <EccentricAnomalyFast+0x3c3>
movsd  -0x80(%rbp),%xmm0
ucomisd -0x70(%rbp),%xmm0
jbe    1482 <EccentricAnomalyFast+0x342>
movsd  -0x80(%rbp),%xmm1
movsd  -0x78(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    1482 <EccentricAnomalyFast+0x342>
movsd  -0x80(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
jmp    149d <EccentricAnomalyFast+0x35d>
movsd  -0x70(%rbp),%xmm0
addsd  -0x78(%rbp),%xmm0
movsd  0xb8c(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x30(%rbp),%xmm1
subsd  -0x70(%rbp),%xmm1
movaps 0xb82(%rip),%xmm0
pand   %xmm0,%xmm1
movsd  -0x20(%rbp),%xmm0
ucomisd %xmm1,%xmm0
ja     14e5 <EccentricAnomalyFast+0x3a5>
movsd  -0x30(%rbp),%xmm1
subsd  -0x78(%rbp),%xmm1
movaps 0xb5e(%rip),%xmm0
pand   %xmm0,%xmm1
movsd  -0x20(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    14ea <EccentricAnomalyFast+0x3aa>
jmp    1503 <EccentricAnomalyFast+0x3c3>
jmp    14ef <EccentricAnomalyFast+0x3af>
mov    -0x84(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x84(%rbp)
jmp    120c <EccentricAnomalyFast+0xcc>
mov    -0x84(%rbp),%eax
cmp    -0x24(%rbp),%eax
jl     1524 <EccentricAnomalyFast+0x3e4>
movsd  0xb0e(%rip),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    152e <EccentricAnomalyFast+0x3ee>
movsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
add    $0xc0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<EccentricAnomalyFast>:
push   %rbp
push   %r15
push   %r14
push   %rbx
sub    $0x88,%rsp
mov    %edi,%r14d
movsd  %xmm2,0x38(%rsp)
movaps %xmm1,0x10(%rsp)
movsd  %xmm0,(%rsp)
mov    $0xffffffff,%edi
call   1030 <acos@plt>
mov    %eax,%ebp
mov    $0xffffffff,%edi
call   1030 <acos@plt>
movapd 0x10(%rsp),%xmm1
xorpd  %xmm0,%xmm0
ucomisd %xmm0,%xmm1
jne    1192 <EccentricAnomalyFast+0x52>
jp     1192 <EccentricAnomalyFast+0x52>
movsd  (%rsp),%xmm0
jmp    13f4 <EccentricAnomalyFast+0x2b4>
mov    %eax,%ebx
movsd  (%rsp),%xmm0
call   1040 <sin@plt>
call   1050 <sgn@plt>
movapd 0x10(%rsp),%xmm1
movapd %xmm0,%xmm5
mulsd  0xe4b(%rip),%xmm5
mulsd  %xmm1,%xmm5
addsd  (%rsp),%xmm5
test   %r14d,%r14d
jle    13de <EccentricAnomalyFast+0x29e>
add    %ebp,%ebp
neg    %ebp
cvtsi2sd %ebp,%xmm6
add    %ebx,%ebx
cvtsi2sd %ebx,%xmm4
movapd 0xe53(%rip),%xmm0
xorpd  %xmm1,%xmm0
movapd %xmm0,0x50(%rsp)
lea    0x1(%r14),%r15d
mov    %r14d,%ebp
neg    %ebp
mov    $0x1,%ebx
cs nopw 0x0(%rax,%rax,1)
nop
movsd  %xmm4,0x8(%rsp)
movapd %xmm6,0x60(%rsp)
movapd %xmm5,0x20(%rsp)
movapd %xmm5,%xmm0
call   1040 <sin@plt>
movapd 0x20(%rsp),%xmm1
mulsd  0x50(%rsp),%xmm0
addsd  %xmm1,%xmm0
subsd  (%rsp),%xmm0
movapd %xmm0,0x70(%rsp)
ucomisd 0xdea(%rip),%xmm0
movapd %xmm1,%xmm0
ja     124a <EccentricAnomalyFast+0x10a>
movsd  0x8(%rsp),%xmm0
movsd  %xmm0,0x40(%rsp)
movapd %xmm1,%xmm0
call   1060 <cos@plt>
mulsd  0x50(%rsp),%xmm0
addsd  0xda1(%rip),%xmm0
movsd  %xmm0,0x8(%rsp)
movapd 0x20(%rsp),%xmm0
call   1040 <sin@plt>
mulsd  0x10(%rsp),%xmm0
movsd  %xmm0,0x48(%rsp)
movapd 0x20(%rsp),%xmm0
call   1060 <cos@plt>
movapd 0x20(%rsp),%xmm5
mulsd  0x10(%rsp),%xmm0
movapd 0x70(%rsp),%xmm4
movapd %xmm4,%xmm1
xorpd  0xd83(%rip),%xmm1
movapd %xmm4,%xmm2
movsd  0x8(%rsp),%xmm6
divsd  %xmm6,%xmm2
movsd  0xd4d(%rip),%xmm7
mulsd  %xmm7,%xmm2
movsd  0x48(%rsp),%xmm8
mulsd  %xmm8,%xmm2
movapd %xmm6,%xmm3
subsd  %xmm2,%xmm3
divsd  %xmm3,%xmm1
movapd %xmm1,%xmm2
mulsd  0xd2d(%rip),%xmm2
mulsd  %xmm1,%xmm2
mulsd  %xmm7,%xmm1
mulsd  %xmm8,%xmm1
addsd  %xmm6,%xmm1
mulsd  %xmm0,%xmm2
addsd  %xmm1,%xmm2
movapd %xmm4,%xmm1
divsd  %xmm2,%xmm1
movapd %xmm5,%xmm0
subsd  %xmm1,%xmm0
movapd %xmm5,%xmm1
subsd  %xmm0,%xmm1
movapd 0xd1b(%rip),%xmm8
andpd  %xmm8,%xmm1
movsd  0x38(%rsp),%xmm9
ucomisd %xmm1,%xmm9
ja     13e3 <EccentricAnomalyFast+0x2a3>
movapd %xmm4,%xmm2
xorpd  %xmm1,%xmm1
cmpltsd %xmm4,%xmm1
movapd 0x60(%rsp),%xmm6
andpd  %xmm1,%xmm6
andnpd %xmm5,%xmm1
orpd   %xmm1,%xmm6
movapd %xmm6,%xmm1
movsd  0x40(%rsp),%xmm4
addsd  %xmm4,%xmm1
mulsd  %xmm7,%xmm1
movapd %xmm6,%xmm2
cmpltsd %xmm0,%xmm2
movapd %xmm0,%xmm5
cmpltsd %xmm4,%xmm5
movapd %xmm5,%xmm3
andnpd %xmm1,%xmm3
andpd  %xmm0,%xmm5
orpd   %xmm3,%xmm5
andpd  %xmm2,%xmm5
andnpd %xmm1,%xmm2
orpd   %xmm2,%xmm5
movapd %xmm5,%xmm0
subsd  %xmm6,%xmm0
andpd  %xmm8,%xmm0
ucomisd %xmm0,%xmm9
ja     13e3 <EccentricAnomalyFast+0x2a3>
movapd %xmm5,%xmm0
subsd  %xmm4,%xmm0
andpd  %xmm8,%xmm0
ucomisd %xmm0,%xmm9
ja     13e3 <EccentricAnomalyFast+0x2a3>
lea    (%rbx,%rbp,1),%eax
add    $0x1,%eax
mov    %ebx,%ecx
add    $0x1,%ecx
mov    %ecx,%ebx
cmp    $0x1,%eax
jne    1200 <EccentricAnomalyFast+0xc0>
mov    %r15d,%ebx
jmp    13e3 <EccentricAnomalyFast+0x2a3>
mov    $0x1,%ebx
cmp    %r14d,%ebx
jl     13f0 <EccentricAnomalyFast+0x2b0>
movsd  0xc30(%rip),%xmm5
movapd %xmm5,%xmm0
add    $0x88,%rsp
pop    %rbx
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define double @EccentricAnomalyFast(double %0, double %1, double %2, i32 %3) {
%5 = alloca double
%6 = alloca double
%7 = alloca double
%8 = alloca double
%9 = alloca i32
%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 i32
store double %0, double* %6
store double %1, double* %7
store double %2, double* %8
store i32 %3, i32* %9
store double 1.000000e+00, double* %13
%22 = call i32 @acos(i32 -1)
%23 = mul nsw i32 -2, %22
%24 = sitofp i32 %23 to double
store double %24, double* %18
%25 = call i32 @acos(i32 -1)
%26 = mul nsw i32 2, %25
%27 = sitofp i32 %26 to double
store double %27, double* %19
%28 = load double, double* %7
%29 = fcmp oeq double %28, 0.000000e+00
br i1 %29, label %30, label %32
30: ; preds = %4
%31 = load double, double* %6
store double %31, double* %5
br label %150
32: ; preds = %4
%33 = load double, double* %6
%34 = load double, double* %6
%35 = call double @sin(double %34)
%36 = call double @sgn(double %35)
%37 = fmul double %36, 8.500000e-01
%38 = load double, double* %7
%39 = call double @llvm.fmuladd.f64(double %37, double %38, double %33)
store double %39, double* %10
store i32 1, i32* %21
br label %40
40: ; preds = %140, %32
%41 = load i32, i32* %21
%42 = load i32, i32* %9
%43 = icmp sle i32 %41, %42
br i1 %43, label %44, label %143
44: ; preds = %40
%45 = load double, double* %10
%46 = load double, double* %7
%47 = load double, double* %10
%48 = call double @sin(double %47)
%49 = fneg double %46
%50 = call double @llvm.fmuladd.f64(double %49, double %48, double %45)
%51 = load double, double* %6
%52 = fsub double %50, %51
store double %52, double* %14
%53 = load double, double* %14
%54 = fcmp ogt double %53, 0.000000e+00
br i1 %54, label %55, label %57
55: ; preds = %44
%56 = load double, double* %10
store double %56, double* %19
br label %59
57: ; preds = %44
%58 = load double, double* %10
store double %58, double* %18
br label %59
59: ; preds = %57, %55
%60 = load double, double* %7
%61 = load double, double* %10
%62 = call double @cos(double %61)
%63 = fneg double %60
%64 = call double @llvm.fmuladd.f64(double %63, double %62, double 1.000000e+00)
store double %64, double* %15
%65 = load double, double* %7
%66 = load double, double* %10
%67 = call double @sin(double %66)
%68 = fmul double %65, %67
store double %68, double* %16
%69 = load double, double* %7
%70 = load double, double* %10
%71 = call double @cos(double %70)
%72 = fmul double %69, %71
store double %72, double* %17
%73 = load double, double* %14
%74 = fneg double %73
%75 = load double, double* %15
%76 = fdiv double %74, %75
store double %76, double* %11
%77 = load double, double* %14
%78 = fneg double %77
%79 = load double, double* %15
%80 = load double, double* %11
%81 = fmul double 5.000000e-01, %80
%82 = load double, double* %16
%83 = call double @llvm.fmuladd.f64(double %81, double %82, double %79)
%84 = fdiv double %78, %83
store double %84, double* %12
%85 = load double, double* %14
%86 = fneg double %85
%87 = load double, double* %15
%88 = load double, double* %12
%89 = fmul double 5.000000e-01, %88
%90 = load double, double* %16
%91 = call double @llvm.fmuladd.f64(double %89, double %90, double %87)
%92 = load double, double* %12
%93 = fmul double 0x3FC5555555555555, %92
%94 = load double, double* %12
%95 = fmul double %93, %94
%96 = load double, double* %17
%97 = call double @llvm.fmuladd.f64(double %95, double %96, double %91)
%98 = fdiv double %86, %97
store double %98, double* %13
%99 = load double, double* %10
%100 = load double, double* %13
%101 = fadd double %99, %100
store double %101, double* %20
%102 = load double, double* %10
%103 = load double, double* %20
%104 = fsub double %102, %103
%105 = call double @llvm.fabs.f64(double %104)
%106 = load double, double* %8
%107 = fcmp olt double %105, %106
br i1 %107, label %108, label %109
108: ; preds = %59
br label %143
109: ; preds = %59
%110 = load double, double* %20
%111 = load double, double* %18
%112 = fcmp ogt double %110, %111
br i1 %112, label %113, label %119
113: ; preds = %109
%114 = load double, double* %20
%115 = load double, double* %19
%116 = fcmp olt double %114, %115
br i1 %116, label %117, label %119
117: ; preds = %113
%118 = load double, double* %20
store double %118, double* %10
br label %124
119: ; preds = %113, %109
%120 = load double, double* %18
%121 = load double, double* %19
%122 = fadd double %120, %121
%123 = fdiv double %122, 2.000000e+00
store double %123, double* %10
br label %124
124: ; preds = %119, %117
%125 = load double, double* %10
%126 = load double, double* %18
%127 = fsub double %125, %126
%128 = call double @llvm.fabs.f64(double %127)
%129 = load double, double* %8
%130 = fcmp olt double %128, %129
br i1 %130, label %138, label %131
131: ; preds = %124
%132 = load double, double* %10
%133 = load double, double* %19
%134 = fsub double %132, %133
%135 = call double @llvm.fabs.f64(double %134)
%136 = load double, double* %8
%137 = fcmp olt double %135, %136
br i1 %137, label %138, label %139
138: ; preds = %131, %124
br label %143
139: ; preds = %131
br label %140
140: ; preds = %139
%141 = load i32, i32* %21
%142 = add nsw i32 %141, 1
store i32 %142, i32* %21
br label %40
143: ; preds = %138, %108, %40
%144 = load i32, i32* %21
%145 = load i32, i32* %9
%146 = icmp sge i32 %144, %145
br i1 %146, label %147, label %148
147: ; preds = %143
store double -1.000000e+00, double* %5
br label %150
148: ; preds = %143
%149 = load double, double* %10
store double %149, double* %5
br label %150
150: ; preds = %148, %147, %30
%151 = load double, double* %5
ret double %151
}

LLVM IR (-O2):
define double @EccentricAnomalyFast(double %0, double %1, double %2, i32 %3) {
%5 = tail call i32 @acos(i32 -1)
%6 = tail call i32 @acos(i32 -1)
%7 = fcmp oeq double %1, 0.000000e+00
br i1 %7, label %75, label %8
8: ; preds = %4
%9 = tail call double @sin(double %0)
%10 = tail call double @sgn(double %9)
%11 = fmul double %10, 8.500000e-01
%12 = tail call double @llvm.fmuladd.f64(double %11, double %1, double %0)
%13 = icmp slt i32 %3, 1
br i1 %13, label %70, label %14
14: ; preds = %8
%15 = mul nsw i32 %5, -2
%16 = sitofp i32 %15 to double
%17 = shl nsw i32 %6, 1
%18 = sitofp i32 %17 to double
%19 = fneg double %1
%20 = add nuw i32 %3, 1
br label %21
21: ; preds = %14, %67
%22 = phi i32 [ 1, %14 ], [ %68, %67 ]
%23 = phi double [ %18, %14 ], [ %31, %67 ]
%24 = phi double [ %16, %14 ], [ %30, %67 ]
%25 = phi double [ %12, %14 ], [ %59, %67 ]
%26 = tail call double @sin(double %25)
%27 = tail call double @llvm.fmuladd.f64(double %19, double %26, double %25)
%28 = fsub double %27, %0
%29 = fcmp ogt double %28, 0.000000e+00
%30 = select i1 %29, double %24, double %25
%31 = select i1 %29, double %25, double %23
%32 = tail call double @cos(double %25)
%33 = tail call double @llvm.fmuladd.f64(double %19, double %32, double 1.000000e+00)
%34 = tail call double @sin(double %25)
%35 = fmul double %34, %1
%36 = tail call double @cos(double %25)
%37 = fmul double %36, %1
%38 = fneg double %28
%39 = fdiv double %38, %33
%40 = fmul double %39, 5.000000e-01
%41 = tail call double @llvm.fmuladd.f64(double %40, double %35, double %33)
%42 = fdiv double %38, %41
%43 = fmul double %42, 5.000000e-01
%44 = tail call double @llvm.fmuladd.f64(double %43, double %35, double %33)
%45 = fmul double %42, 0x3FC5555555555555
%46 = fmul double %42, %45
%47 = tail call double @llvm.fmuladd.f64(double %46, double %37, double %44)
%48 = fdiv double %28, %47
%49 = fsub double %25, %48
%50 = fsub double %25, %49
%51 = tail call double @llvm.fabs.f64(double %50)
%52 = fcmp olt double %51, %2
br i1 %52, label %70, label %53
53: ; preds = %21
%54 = fcmp ogt double %49, %30
%55 = fcmp olt double %49, %31
%56 = select i1 %54, i1 %55, i1 false
%57 = fadd double %30, %31
%58 = fmul double %57, 5.000000e-01
%59 = select i1 %56, double %49, double %58
%60 = fsub double %59, %30
%61 = tail call double @llvm.fabs.f64(double %60)
%62 = fcmp olt double %61, %2
br i1 %62, label %70, label %63
63: ; preds = %53
%64 = fsub double %59, %31
%65 = tail call double @llvm.fabs.f64(double %64)
%66 = fcmp olt double %65, %2
br i1 %66, label %70, label %67
67: ; preds = %63
%68 = add nuw i32 %22, 1
%69 = icmp eq i32 %22, %3
br i1 %69, label %70, label %21
70: ; preds = %67, %21, %63, %53, %8
%71 = phi i32 [ 1, %8 ], [ %22, %53 ], [ %22, %63 ], [ %22, %21 ], [ %20, %67 ]
%72 = phi double [ %12, %8 ], [ %59, %53 ], [ %59, %63 ], [ %25, %21 ], [ %59, %67 ]
%73 = icmp slt i32 %71, %3
%74 = select i1 %73, double %72, double -1.000000e+00
br label %75
75: ; preds = %70, %4
%76 = phi double [ %0, %4 ], [ %74, %70 ]
ret double %76
}