#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 M_PI;
double sunmodel_mean_obliquity_of_ecliptic(double t) {
  double seconds, e0;
  seconds = 21.448 - t * (46.8150 + t * (0.00059 - t * 0.001813));
  e0 = (M_PI * (23.0 + (26.0 + seconds / 60.0) / 60.0) / 180.0);
  return e0;
}


x86 assembly (clang, -O0):
<sunmodel_mean_obliquity_of_ecliptic>:
push   %rbp
mov    %rsp,%rbp
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
movsd  -0x8(%rbp),%xmm1
movsd  -0x8(%rbp),%xmm2
movq   %xmm2,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm2
movsd  0xef9(%rip),%xmm3
mulsd  %xmm3,%xmm2
movsd  0xef5(%rip),%xmm3
addsd  %xmm3,%xmm2
mulsd  %xmm2,%xmm1
movsd  0xed5(%rip),%xmm2
addsd  %xmm2,%xmm1
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm1,%xmm0
movsd  0xea6(%rip),%xmm1
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x10(%rbp)
mov    0x2e5e(%rip),%rax
movsd  (%rax),%xmm0
movsd  -0x10(%rbp),%xmm1
movsd  0xe75(%rip),%xmm2
divsd  %xmm2,%xmm1
movsd  0xe71(%rip),%xmm2
addsd  %xmm1,%xmm2
movsd  0xe5d(%rip),%xmm1
divsd  %xmm1,%xmm2
movsd  0xe49(%rip),%xmm1
addsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
movsd  0xe31(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x18(%rbp)
movsd  -0x18(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<sunmodel_mean_obliquity_of_ecliptic>:
movsd  0xef8(%rip),%xmm1
mulsd  %xmm0,%xmm1
addsd  0xef4(%rip),%xmm1
mulsd  %xmm0,%xmm1
addsd  0xef0(%rip),%xmm1
mulsd  %xmm0,%xmm1
movsd  0xeec(%rip),%xmm0
subsd  %xmm1,%xmm0
movsd  0xee8(%rip),%xmm1
divsd  %xmm1,%xmm0
addsd  0xee4(%rip),%xmm0
mov    0x2e9d(%rip),%rax
divsd  %xmm1,%xmm0
addsd  0xed9(%rip),%xmm0
mulsd  (%rax),%xmm0
divsd  0xed5(%rip),%xmm0
ret

LLVM IR (-O0):
@M_PI = global double 0.000000e+00
define double @sunmodel_mean_obliquity_of_ecliptic(double %0) {
%2 = alloca double
%3 = alloca double
%4 = alloca double
store double %0, double* %2
%5 = load double, double* %2
%6 = load double, double* %2
%7 = load double, double* %2
%8 = fneg double %7
%9 = call double @llvm.fmuladd.f64(double %8, double 1.813000e-03, double 5.900000e-04)
%10 = call double @llvm.fmuladd.f64(double %6, double %9, double 4.681500e+01)
%11 = fneg double %5
%12 = call double @llvm.fmuladd.f64(double %11, double %10, double 2.144800e+01)
store double %12, double* %3
%13 = load double, double* @M_PI
%14 = load double, double* %3
%15 = fdiv double %14, 6.000000e+01
%16 = fadd double 2.600000e+01, %15
%17 = fdiv double %16, 6.000000e+01
%18 = fadd double 2.300000e+01, %17
%19 = fmul double %13, %18
%20 = fdiv double %19, 1.800000e+02
store double %20, double* %4
%21 = load double, double* %4
ret double %21
}

LLVM IR (-O2):
@M_PI = global double 0.000000e+00
define double @sunmodel_mean_obliquity_of_ecliptic(double %0) {
%2 = fneg double %0
%3 = tail call double @llvm.fmuladd.f64(double %2, double 1.813000e-03, double 5.900000e-04)
%4 = tail call double @llvm.fmuladd.f64(double %0, double %3, double 4.681500e+01)
%5 = tail call double @llvm.fmuladd.f64(double %2, double %4, double 2.144800e+01)
%6 = load double, double* @M_PI
%7 = fdiv double %5, 6.000000e+01
%8 = fadd double %7, 2.600000e+01
%9 = fdiv double %8, 6.000000e+01
%10 = fadd double %9, 2.300000e+01
%11 = fmul double %10, %6
%12 = fdiv double %11, 1.800000e+02
ret double %12
}