#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 JulianDayFromGregorian(int, int, float);
int cos(int);
double pow(double, int);
int /*<<< orphan*/ printf(char *, double, double, double, double, double);
int sin(int);
double tan(double);
double EquationOfTime(int year, int month, float day) {
  double T, L0, M, y, epsilon, epsilon0, e;
  T = (JulianDayFromGregorian(year, month, day) - 2451545.0) / 36525;
  L0 = 280.46645 + 36000.76983 * T + 0.0003032 * T * T;
  M = 357.52772 + 35999.05034 * T - 0.0001603 * T * T - (T * T * T) / 300000;
  const double CTE1 = 23 + 26.0 / 60 + 21.448 / 3600;
  const double CTE2 = 4680.93 / 3600;
  double U = T / 100;
  epsilon0 = CTE1 - CTE2 * U - 1.55 * pow(U, 2) + 1999.25 * pow(U, 3) -
             51.38 * pow(U, 4) - 249.67 * pow(U, 5) - 39.05 * pow(U, 6) +
             7.12 * pow(U, 7) + 27.87 * pow(U, 8) + 5.79 * pow(U, 9) +
             2.45 * pow(U, 10);
  while (L0 > 0)
    L0 -= 360;
  while (L0 < 0)
    L0 += 360;
  while (M > 0)
    M -= 360;
  while (M < 0)
    M += 360;
  y = tan(epsilon / 2) * tan(epsilon / 2);
  printf("JD\t= %.2f\nT\t= %.15f\nL0\t= %.15f\nM\t= %.15f\neps0\t= %.15f\n",
         JulianDayFromGregorian(year, month, day), T, L0, M, epsilon0);
  return y * sin(2 * L0) - 2 * e * sin(M) + 4 * e * y * sin(M) * cos(2 * L0) -
         (1 / 2) * y * y * sin(4 * L0) - (5 / 4) * e * e * sin(2 * M);
}


x86 assembly (clang, -O0):
<EquationOfTime>:
push   %rbp
mov    %rsp,%rbp
sub    $0x100,%rsp
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
movss  %xmm0,-0xc(%rbp)
mov    -0x4(%rbp),%edi
mov    -0x8(%rbp),%esi
movss  -0xc(%rbp),%xmm0
call   1030 <JulianDayFromGregorian@plt>
movsd  0xf1a(%rip),%xmm1
subsd  %xmm1,%xmm0
movsd  0xf06(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x18(%rbp)
movsd  -0x18(%rbp),%xmm1
movsd  0xee0(%rip),%xmm0
mulsd  %xmm0,%xmm1
movsd  0xedc(%rip),%xmm0
addsd  %xmm0,%xmm1
movsd  0xec0(%rip),%xmm0
mulsd  -0x18(%rbp),%xmm0
movsd  -0x18(%rbp),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x20(%rbp)
movsd  -0x18(%rbp),%xmm1
movsd  0xe8c(%rip),%xmm0
mulsd  %xmm0,%xmm1
movsd  0xe88(%rip),%xmm0
addsd  %xmm0,%xmm1
movsd  0xe6c(%rip),%xmm0
mulsd  -0x18(%rbp),%xmm0
movsd  -0x18(%rbp),%xmm2
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  -0x18(%rbp),%xmm1
mulsd  -0x18(%rbp),%xmm1
mulsd  -0x18(%rbp),%xmm1
movsd  0xe24(%rip),%xmm2
divsd  %xmm2,%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x28(%rbp)
movsd  0xdf7(%rip),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  0xdfa(%rip),%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x18(%rbp),%xmm0
movsd  0xde0(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  -0x60(%rbp),%xmm0
movsd  0xdba(%rip),%xmm1
mulsd  %xmm1,%xmm0
movsd  0xdb6(%rip),%xmm1
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xa8(%rbp)
ff
movsd  -0x60(%rbp),%xmm0
mov    $0x2,%edi
call   1080 <pow@plt>
movsd  -0xa8(%rbp),%xmm1
ff
movsd  0xd7b(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xa0(%rbp)
ff
movsd  -0x60(%rbp),%xmm0
mov    $0x3,%edi
call   1080 <pow@plt>
movsd  -0xa0(%rbp),%xmm1
ff
movsd  0xd44(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
movsd  -0x60(%rbp),%xmm0
mov    $0x4,%edi
call   1080 <pow@plt>
movsd  -0x98(%rbp),%xmm1
ff
movsd  0xd0d(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
movsd  -0x60(%rbp),%xmm0
mov    $0x5,%edi
call   1080 <pow@plt>
movsd  -0x90(%rbp),%xmm1
ff
movsd  0xcd6(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
movsd  -0x60(%rbp),%xmm0
mov    $0x6,%edi
call   1080 <pow@plt>
movsd  -0x88(%rbp),%xmm1
ff
movsd  0xc9f(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x80(%rbp)
movsd  -0x60(%rbp),%xmm0
mov    $0x7,%edi
call   1080 <pow@plt>
movsd  -0x80(%rbp),%xmm1
movsd  0xc6e(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x78(%rbp)
movsd  -0x60(%rbp),%xmm0
mov    $0x8,%edi
call   1080 <pow@plt>
movsd  -0x78(%rbp),%xmm1
movsd  0xc3d(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x60(%rbp),%xmm0
mov    $0x9,%edi
call   1080 <pow@plt>
movsd  -0x70(%rbp),%xmm1
movsd  0xc0c(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  -0x60(%rbp),%xmm0
mov    $0xa,%edi
call   1080 <pow@plt>
movsd  -0x68(%rbp),%xmm1
movsd  0xbdb(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x20(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    145f <EquationOfTime+0x2ff>
movsd  -0x20(%rbp),%xmm0
movsd  0xc77(%rip),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x20(%rbp)
jmp    1432 <EquationOfTime+0x2d2>
jmp    1464 <EquationOfTime+0x304>
xorps  %xmm0,%xmm0
ucomisd -0x20(%rbp),%xmm0
jbe    1489 <EquationOfTime+0x329>
movsd  0xc4e(%rip),%xmm0
addsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x20(%rbp)
jmp    1464 <EquationOfTime+0x304>
jmp    148e <EquationOfTime+0x32e>
movsd  -0x28(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    14bb <EquationOfTime+0x35b>
movsd  -0x28(%rbp),%xmm0
movsd  0xc1b(%rip),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x28(%rbp)
jmp    148e <EquationOfTime+0x32e>
jmp    14c0 <EquationOfTime+0x360>
xorps  %xmm0,%xmm0
ucomisd -0x28(%rbp),%xmm0
jbe    14e5 <EquationOfTime+0x385>
movsd  0xbf2(%rip),%xmm0
addsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0x28(%rbp)
jmp    14c0 <EquationOfTime+0x360>
movsd  -0x38(%rbp),%xmm0
movsd  0xbbe(%rip),%xmm1
divsd  %xmm1,%xmm0
call   1060 <tan@plt>
movsd  %xmm0,-0x100(%rbp)
ff
movsd  -0x38(%rbp),%xmm0
movsd  0xba0(%rip),%xmm1
divsd  %xmm1,%xmm0
call   1060 <tan@plt>
movaps %xmm0,%xmm1
movsd  -0x100(%rbp),%xmm0
ff
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x4(%rbp),%edi
mov    -0x8(%rbp),%esi
movss  -0xc(%rbp),%xmm0
call   1030 <JulianDayFromGregorian@plt>
movsd  -0x18(%rbp),%xmm1
movsd  -0x20(%rbp),%xmm2
movsd  -0x28(%rbp),%xmm3
movsd  -0x40(%rbp),%xmm4
lea    0xb78(%rip),%rdi
call   1040 <printf@plt>
movsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0xf0(%rbp)
ff
movsd  0xb3e(%rip),%xmm0
mulsd  -0x20(%rbp),%xmm0
cvttsd2si %xmm0,%edi
call   1050 <sin@plt>
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0xe8(%rbp)
ff
movsd  0xb1c(%rip),%xmm0
mulsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0xf8(%rbp)
ff
cvttsd2si -0x28(%rbp),%edi
call   1050 <sin@plt>
movsd  -0xf8(%rbp),%xmm1
ff
movsd  -0xf0(%rbp),%xmm0
ff
movsd  -0xe8(%rbp),%xmm2
ff
cvtsi2sd %eax,%xmm3
mulsd  %xmm3,%xmm1
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xd0(%rbp)
ff
movsd  0xac6(%rip),%xmm0
mulsd  -0x48(%rbp),%xmm0
mulsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0xe0(%rbp)
ff
cvttsd2si -0x28(%rbp),%edi
call   1050 <sin@plt>
movsd  -0xe0(%rbp),%xmm0
ff
cvtsi2sd %eax,%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0xd8(%rbp)
ff
movsd  0xa7a(%rip),%xmm0
mulsd  -0x20(%rbp),%xmm0
cvttsd2si %xmm0,%edi
call   1070 <cos@plt>
movsd  -0xd8(%rbp),%xmm0
ff
movsd  -0xd0(%rbp),%xmm1
ff
cvtsi2sd %eax,%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xc0(%rbp)
ff
xorps  %xmm0,%xmm0
mulsd  -0x30(%rbp),%xmm0
mulsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0xc8(%rbp)
ff
movsd  0xa3b(%rip),%xmm0
mulsd  -0x20(%rbp),%xmm0
cvttsd2si %xmm0,%edi
call   1050 <sin@plt>
movsd  -0xc8(%rbp),%xmm0
ff
movsd  -0xc0(%rbp),%xmm1
ff
cvtsi2sd %eax,%xmm2
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xb0(%rbp)
ff
movsd  0x9e2(%rip),%xmm0
mulsd  -0x48(%rbp),%xmm0
mulsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0xb8(%rbp)
ff
movsd  0x9c0(%rip),%xmm0
mulsd  -0x28(%rbp),%xmm0
cvttsd2si %xmm0,%edi
call   1050 <sin@plt>
movsd  -0xb8(%rbp),%xmm0
ff
movsd  -0xb0(%rbp),%xmm1
ff
cvtsi2sd %eax,%xmm2
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
add    $0x100,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<EquationOfTime>:
push   %rbp
push   %rbx
sub    $0x78,%rsp
mov    %esi,%ebx
mov    %edi,%ebp
movss  %xmm0,0x24(%rsp)
call   1030 <JulianDayFromGregorian@plt>
addsd  0xe93(%rip),%xmm0
divsd  0xe93(%rip),%xmm0
movsd  0xe93(%rip),%xmm1
mulsd  %xmm0,%xmm1
addsd  0xe8f(%rip),%xmm1
movsd  0xe8f(%rip),%xmm2
mulsd  %xmm0,%xmm2
mulsd  %xmm0,%xmm2
addsd  %xmm1,%xmm2
movsd  %xmm2,0x18(%rsp)
movsd  0xe7d(%rip),%xmm2
mulsd  %xmm0,%xmm2
addsd  0xe79(%rip),%xmm2
movsd  0xe79(%rip),%xmm1
mulsd  %xmm0,%xmm1
mulsd  %xmm0,%xmm1
addsd  %xmm2,%xmm1
movsd  %xmm1,0x10(%rsp)
movapd %xmm0,%xmm1
mulsd  %xmm0,%xmm1
mulsd  %xmm0,%xmm1
divsd  0xe5b(%rip),%xmm1
movsd  %xmm1,0x28(%rsp)
movsd  %xmm0,0x70(%rsp)
divsd  0xe4f(%rip),%xmm0
movsd  %xmm0,0x8(%rsp)
mov    $0x2,%edi
call   1070 <pow@plt>
movsd  %xmm0,0x68(%rsp)
movsd  0x8(%rsp),%xmm0
mov    $0x3,%edi
call   1070 <pow@plt>
movsd  %xmm0,0x60(%rsp)
movsd  0x8(%rsp),%xmm0
mov    $0x4,%edi
call   1070 <pow@plt>
movsd  %xmm0,0x58(%rsp)
movsd  0x8(%rsp),%xmm0
mov    $0x5,%edi
call   1070 <pow@plt>
movsd  %xmm0,0x50(%rsp)
movsd  0x8(%rsp),%xmm0
mov    $0x6,%edi
call   1070 <pow@plt>
movsd  %xmm0,0x48(%rsp)
movsd  0x8(%rsp),%xmm0
mov    $0x7,%edi
call   1070 <pow@plt>
movsd  %xmm0,0x40(%rsp)
movsd  0x8(%rsp),%xmm0
mov    $0x8,%edi
call   1070 <pow@plt>
movsd  %xmm0,0x38(%rsp)
movsd  0x8(%rsp),%xmm0
mov    $0x9,%edi
call   1070 <pow@plt>
movsd  %xmm0,0x30(%rsp)
movsd  0x8(%rsp),%xmm0
mov    $0xa,%edi
call   1070 <pow@plt>
movsd  0x18(%rsp),%xmm2
xorpd  %xmm4,%xmm4
ucomisd %xmm4,%xmm2
jbe    12ea <EquationOfTime+0x19a>
movsd  0xd7f(%rip),%xmm1
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
addsd  %xmm1,%xmm2
ucomisd %xmm4,%xmm2
ja     12e0 <EquationOfTime+0x190>
movsd  0x10(%rsp),%xmm3
addsd  0x28(%rsp),%xmm3
ucomisd %xmm2,%xmm4
jbe    131a <EquationOfTime+0x1ca>
movsd  0xd54(%rip),%xmm4
xorpd  %xmm1,%xmm1
nopl   0x0(%rax,%rax,1)
addsd  %xmm4,%xmm2
ucomisd %xmm2,%xmm1
ja     1310 <EquationOfTime+0x1c0>
xorpd  %xmm4,%xmm4
ucomisd %xmm4,%xmm3
jbe    133a <EquationOfTime+0x1ea>
movsd  0xd24(%rip),%xmm1
nopl   0x0(%rax)
addsd  %xmm1,%xmm3
ucomisd %xmm4,%xmm3
ja     1330 <EquationOfTime+0x1e0>
ucomisd %xmm3,%xmm4
movsd  %xmm2,0x18(%rsp)
jbe    136a <EquationOfTime+0x21a>
movsd  0xd0a(%rip),%xmm2
xorpd  %xmm1,%xmm1
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
addsd  %xmm2,%xmm3
ucomisd %xmm3,%xmm1
ja     1360 <EquationOfTime+0x210>
movsd  %xmm3,0x10(%rsp)
movsd  0x8(%rsp),%xmm2
mulsd  0xce2(%rip),%xmm2
addsd  0xce2(%rip),%xmm2
movsd  0x68(%rsp),%xmm1
mulsd  0xcdc(%rip),%xmm1
addsd  %xmm2,%xmm1
movsd  0x60(%rsp),%xmm3
mulsd  0xcd2(%rip),%xmm3
movsd  0x58(%rsp),%xmm2
mulsd  0xccc(%rip),%xmm2
addsd  %xmm1,%xmm3
addsd  %xmm3,%xmm2
movsd  0x50(%rsp),%xmm3
mulsd  0xcbe(%rip),%xmm3
addsd  %xmm2,%xmm3
movsd  0x48(%rsp),%xmm1
mulsd  0xcb4(%rip),%xmm1
addsd  %xmm3,%xmm1
movsd  0x40(%rsp),%xmm3
mulsd  0xcaa(%rip),%xmm3
movsd  0x38(%rsp),%xmm2
mulsd  0xca4(%rip),%xmm2
addsd  %xmm1,%xmm3
addsd  %xmm3,%xmm2
movsd  0x30(%rsp),%xmm1
mulsd  0xc96(%rip),%xmm1
addsd  %xmm2,%xmm1
mulsd  0xc92(%rip),%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,0x8(%rsp)
mov    %ebp,%edi
mov    %ebx,%esi
movss  0x24(%rsp),%xmm0
call   1030 <JulianDayFromGregorian@plt>
lea    0xc8a(%rip),%rdi
movsd  0x70(%rsp),%xmm1
movsd  0x18(%rsp),%xmm2
movsd  0x10(%rsp),%xmm3
movsd  0x8(%rsp),%xmm4
call   1040 <printf@plt>
movsd  0x18(%rsp),%xmm0
addsd  %xmm0,%xmm0
cvttsd2si %xmm0,%ebx
mov    %ebx,%edi
call   1050 <sin@plt>
cvttsd2si 0x10(%rsp),%ebp
mov    %ebp,%edi
call   1050 <sin@plt>
mov    %ebp,%edi
call   1050 <sin@plt>
mov    %ebx,%edi
call   1060 <cos@plt>
movsd  0x18(%rsp),%xmm0
mulsd  0xc1f(%rip),%xmm0
cvttsd2si %xmm0,%edi
call   1050 <sin@plt>
movsd  0x10(%rsp),%xmm0
addsd  %xmm0,%xmm0
cvttsd2si %xmm0,%edi
call   1050 <sin@plt>
movsd  0xc03(%rip),%xmm0
add    $0x78,%rsp
pop    %rbx
pop    %rbp
ret

LLVM IR (-O0):
@.str = [55 x i8] c"JD\09= %.2f\0AT\09= %.15f\0AL0\09= %.15f\0AM\09= %.15f\0Aeps0\09= %.15f\0A\00"
define double @EquationOfTime(i32 %0, i32 %1, float %2) {
%4 = alloca i32
%5 = alloca i32
%6 = alloca float
%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
store i32 %0, i32* %4
store i32 %1, i32* %5
store float %2, float* %6
%17 = load i32, i32* %4
%18 = load i32, i32* %5
%19 = load float, float* %6
%20 = call double @JulianDayFromGregorian(i32 %17, i32 %18, float %19)
%21 = fsub double %20, 0x4142B42C80000000
%22 = fdiv double %21, 3.652500e+04
store double %22, double* %7
%23 = load double, double* %7
%24 = call double @llvm.fmuladd.f64(double 0x40E19418A272862F, double %23, double 0x4071877694467382)
%25 = load double, double* %7
%26 = fmul double 3.032000e-04, %25
%27 = load double, double* %7
%28 = call double @llvm.fmuladd.f64(double %26, double %27, double %24)
store double %28, double* %8
%29 = load double, double* %7
%30 = call double @llvm.fmuladd.f64(double 0x40E193E19C62A1B6, double %29, double 0x407658718A86D71F)
%31 = load double, double* %7
%32 = fmul double 1.603000e-04, %31
%33 = load double, double* %7
%34 = fneg double %32
%35 = call double @llvm.fmuladd.f64(double %34, double %33, double %30)
%36 = load double, double* %7
%37 = load double, double* %7
%38 = fmul double %36, %37
%39 = load double, double* %7
%40 = fmul double %38, %39
%41 = fdiv double %40, 3.000000e+05
%42 = fsub double %35, %41
store double %42, double* %9
store double 0x4037707561DBA54E, double* %14
store double 0x3FF4CDDBAEA04A47, double* %15
%43 = load double, double* %7
%44 = fdiv double %43, 1.000000e+02
store double %44, double* %16
%45 = load double, double* %16
%46 = call double @llvm.fmuladd.f64(double 0xBFF4CDDBAEA04A47, double %45, double 0x4037707561DBA54E)
%47 = load double, double* %16
%48 = call double @pow(double %47, i32 2)
%49 = call double @llvm.fmuladd.f64(double -1.550000e+00, double %48, double %46)
%50 = load double, double* %16
%51 = call double @pow(double %50, i32 3)
%52 = call double @llvm.fmuladd.f64(double 1.999250e+03, double %51, double %49)
%53 = load double, double* %16
%54 = call double @pow(double %53, i32 4)
%55 = call double @llvm.fmuladd.f64(double -5.138000e+01, double %54, double %52)
%56 = load double, double* %16
%57 = call double @pow(double %56, i32 5)
%58 = call double @llvm.fmuladd.f64(double -2.496700e+02, double %57, double %55)
%59 = load double, double* %16
%60 = call double @pow(double %59, i32 6)
%61 = call double @llvm.fmuladd.f64(double -3.905000e+01, double %60, double %58)
%62 = load double, double* %16
%63 = call double @pow(double %62, i32 7)
%64 = call double @llvm.fmuladd.f64(double 7.120000e+00, double %63, double %61)
%65 = load double, double* %16
%66 = call double @pow(double %65, i32 8)
%67 = call double @llvm.fmuladd.f64(double 2.787000e+01, double %66, double %64)
%68 = load double, double* %16
%69 = call double @pow(double %68, i32 9)
%70 = call double @llvm.fmuladd.f64(double 5.790000e+00, double %69, double %67)
%71 = load double, double* %16
%72 = call double @pow(double %71, i32 10)
%73 = call double @llvm.fmuladd.f64(double 2.450000e+00, double %72, double %70)
store double %73, double* %12
br label %74
74: ; preds = %77, %3
%75 = load double, double* %8
%76 = fcmp ogt double %75, 0.000000e+00
br i1 %76, label %77, label %80
77: ; preds = %74
%78 = load double, double* %8
%79 = fsub double %78, 3.600000e+02
store double %79, double* %8
br label %74
80: ; preds = %74
br label %81
81: ; preds = %84, %80
%82 = load double, double* %8
%83 = fcmp olt double %82, 0.000000e+00
br i1 %83, label %84, label %87
84: ; preds = %81
%85 = load double, double* %8
%86 = fadd double %85, 3.600000e+02
store double %86, double* %8
br label %81
87: ; preds = %81
br label %88
88: ; preds = %91, %87
%89 = load double, double* %9
%90 = fcmp ogt double %89, 0.000000e+00
br i1 %90, label %91, label %94
91: ; preds = %88
%92 = load double, double* %9
%93 = fsub double %92, 3.600000e+02
store double %93, double* %9
br label %88
94: ; preds = %88
br label %95
95: ; preds = %98, %94
%96 = load double, double* %9
%97 = fcmp olt double %96, 0.000000e+00
br i1 %97, label %98, label %101
98: ; preds = %95
%99 = load double, double* %9
%100 = fadd double %99, 3.600000e+02
store double %100, double* %9
br label %95
101: ; preds = %95
%102 = load double, double* %11
%103 = fdiv double %102, 2.000000e+00
%104 = call double @tan(double %103)
%105 = load double, double* %11
%106 = fdiv double %105, 2.000000e+00
%107 = call double @tan(double %106)
%108 = fmul double %104, %107
store double %108, double* %10
%109 = load i32, i32* %4
%110 = load i32, i32* %5
%111 = load float, float* %6
%112 = call double @JulianDayFromGregorian(i32 %109, i32 %110, float %111)
%113 = load double, double* %7
%114 = load double, double* %8
%115 = load double, double* %9
%116 = load double, double* %12
%117 = call i32 @printf(i8* getelementptr inbounds ([55 x i8], [55 x i8]* @.str, i64 0, i64 0), double %112, double %113, double %114, double %115, double %116)
%118 = load double, double* %10
%119 = load double, double* %8
%120 = fmul double 2.000000e+00, %119
%121 = fptosi double %120 to i32
%122 = call i32 @sin(i32 %121)
%123 = sitofp i32 %122 to double
%124 = load double, double* %13
%125 = fmul double 2.000000e+00, %124
%126 = load double, double* %9
%127 = fptosi double %126 to i32
%128 = call i32 @sin(i32 %127)
%129 = sitofp i32 %128 to double
%130 = fmul double %125, %129
%131 = fneg double %130
%132 = call double @llvm.fmuladd.f64(double %118, double %123, double %131)
%133 = load double, double* %13
%134 = fmul double 4.000000e+00, %133
%135 = load double, double* %10
%136 = fmul double %134, %135
%137 = load double, double* %9
%138 = fptosi double %137 to i32
%139 = call i32 @sin(i32 %138)
%140 = sitofp i32 %139 to double
%141 = fmul double %136, %140
%142 = load double, double* %8
%143 = fmul double 2.000000e+00, %142
%144 = fptosi double %143 to i32
%145 = call i32 @cos(i32 %144)
%146 = sitofp i32 %145 to double
%147 = call double @llvm.fmuladd.f64(double %141, double %146, double %132)
%148 = load double, double* %10
%149 = fmul double 0.000000e+00, %148
%150 = load double, double* %10
%151 = fmul double %149, %150
%152 = load double, double* %8
%153 = fmul double 4.000000e+00, %152
%154 = fptosi double %153 to i32
%155 = call i32 @sin(i32 %154)
%156 = sitofp i32 %155 to double
%157 = fneg double %151
%158 = call double @llvm.fmuladd.f64(double %157, double %156, double %147)
%159 = load double, double* %13
%160 = fmul double 1.000000e+00, %159
%161 = load double, double* %13
%162 = fmul double %160, %161
%163 = load double, double* %9
%164 = fmul double 2.000000e+00, %163
%165 = fptosi double %164 to i32
%166 = call i32 @sin(i32 %165)
%167 = sitofp i32 %166 to double
%168 = fneg double %162
%169 = call double @llvm.fmuladd.f64(double %168, double %167, double %158)
ret double %169
}

LLVM IR (-O2):
@.str = [55 x i8] c"JD\09= %.2f\0AT\09= %.15f\0AL0\09= %.15f\0AM\09= %.15f\0Aeps0\09= %.15f\0A\00"
define double @EquationOfTime(i32 %0, i32 %1, float %2) {
%4 = tail call double @JulianDayFromGregorian(i32 %0, i32 %1, float %2)
%5 = fadd double %4, 0xC142B42C80000000
%6 = fdiv double %5, 3.652500e+04
%7 = tail call double @llvm.fmuladd.f64(double %6, double 0x40E19418A272862F, double 0x4071877694467382)
%8 = fmul double %6, 3.032000e-04
%9 = tail call double @llvm.fmuladd.f64(double %8, double %6, double %7)
%10 = tail call double @llvm.fmuladd.f64(double %6, double 0x40E193E19C62A1B6, double 0x407658718A86D71F)
%11 = fmul double %6, -1.603000e-04
%12 = tail call double @llvm.fmuladd.f64(double %11, double %6, double %10)
%13 = fmul double %6, %6
%14 = fmul double %6, %13
%15 = fdiv double %14, 3.000000e+05
%16 = fsub double %12, %15
%17 = fdiv double %6, 1.000000e+02
%18 = tail call double @pow(double %17, i32 2)
%19 = tail call double @pow(double %17, i32 3)
%20 = tail call double @pow(double %17, i32 4)
%21 = tail call double @pow(double %17, i32 5)
%22 = tail call double @pow(double %17, i32 6)
%23 = tail call double @pow(double %17, i32 7)
%24 = tail call double @pow(double %17, i32 8)
%25 = tail call double @pow(double %17, i32 9)
%26 = tail call double @pow(double %17, i32 10)
%27 = fcmp ogt double %9, 0.000000e+00
br i1 %27, label %31, label %28
28: ; preds = %31, %3
%29 = phi double [ %9, %3 ], [ %33, %31 ]
%30 = fcmp olt double %29, 0.000000e+00
br i1 %30, label %38, label %35
31: ; preds = %3, %31
%32 = phi double [ %33, %31 ], [ %9, %3 ]
%33 = fadd double %32, -3.600000e+02
%34 = fcmp ogt double %33, 0.000000e+00
br i1 %34, label %31, label %28
35: ; preds = %38, %28
%36 = phi double [ %29, %28 ], [ %40, %38 ]
%37 = fcmp ogt double %16, 0.000000e+00
br i1 %37, label %45, label %42
38: ; preds = %28, %38
%39 = phi double [ %40, %38 ], [ %29, %28 ]
%40 = fadd double %39, 3.600000e+02
%41 = fcmp olt double %40, 0.000000e+00
br i1 %41, label %38, label %35
42: ; preds = %45, %35
%43 = phi double [ %16, %35 ], [ %47, %45 ]
%44 = fcmp olt double %43, 0.000000e+00
br i1 %44, label %49, label %53
45: ; preds = %35, %45
%46 = phi double [ %47, %45 ], [ %16, %35 ]
%47 = fadd double %46, -3.600000e+02
%48 = fcmp ogt double %47, 0.000000e+00
br i1 %48, label %45, label %42
49: ; preds = %42, %49
%50 = phi double [ %51, %49 ], [ %43, %42 ]
%51 = fadd double %50, 3.600000e+02
%52 = fcmp olt double %51, 0.000000e+00
br i1 %52, label %49, label %53
53: ; preds = %49, %42
%54 = phi double [ %43, %42 ], [ %51, %49 ]
%55 = tail call double @llvm.fmuladd.f64(double %17, double 0xBFF4CDDBAEA04A47, double 0x4037707561DBA54E)
%56 = tail call double @llvm.fmuladd.f64(double %18, double -1.550000e+00, double %55)
%57 = tail call double @llvm.fmuladd.f64(double %19, double 1.999250e+03, double %56)
%58 = tail call double @llvm.fmuladd.f64(double %20, double -5.138000e+01, double %57)
%59 = tail call double @llvm.fmuladd.f64(double %21, double -2.496700e+02, double %58)
%60 = tail call double @llvm.fmuladd.f64(double %22, double -3.905000e+01, double %59)
%61 = tail call double @llvm.fmuladd.f64(double %23, double 7.120000e+00, double %60)
%62 = tail call double @llvm.fmuladd.f64(double %24, double 2.787000e+01, double %61)
%63 = tail call double @llvm.fmuladd.f64(double %25, double 5.790000e+00, double %62)
%64 = tail call double @llvm.fmuladd.f64(double %26, double 2.450000e+00, double %63)
%65 = tail call double @JulianDayFromGregorian(i32 %0, i32 %1, float %2)
%66 = tail call i32 @printf(i8* getelementptr inbounds ([55 x i8], [55 x i8]* @.str, i64 0, i64 0), double %65, double %6, double %36, double %54, double %64)
%67 = fmul double %36, 2.000000e+00
%68 = fptosi double %67 to i32
%69 = tail call i32 @sin(i32 %68)
%70 = fptosi double %54 to i32
%71 = tail call i32 @sin(i32 %70)
%72 = tail call i32 @sin(i32 %70)
%73 = tail call i32 @cos(i32 %68)
%74 = fmul double %36, 4.000000e+00
%75 = fptosi double %74 to i32
%76 = tail call i32 @sin(i32 %75)
%77 = fmul double %54, 2.000000e+00
%78 = fptosi double %77 to i32
%79 = tail call i32 @sin(i32 %78)
ret double 0x7FF8000000000000
}