#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 JDNow();
double Map360(double);
double cos(int);
double NObliquity(void) {
  double deps, lsun, lmoon, omega;
  double jdnow, dt, t;
  dt = 35.0;
  dt = dt + 32.184;
  dt = dt / (36525 * 24. * 60. * 60.);
  jdnow = JDNow();
  t = (jdnow - 2451545.0) / 36525.0;
  t = t + dt;
  omega = 125.04452 - 1934.136261 * t + 0.0020708 * t * t + t * t * t / 450000.;
  lsun = Map360(280.4665 + 3600.7698 * t);
  lmoon = Map360(218.3165 + 481267.8813 * t);
  omega = omega * 3.14159265358 / 180.;
  lsun = lsun * 3.14159265358 / 180.;
  lmoon = lmoon * 3.14159265358 / 180.;
  deps = 9.20 * cos(omega) + 0.57 * cos(2. * lsun) + 0.1 * cos(2. * lmoon) -
         0.09 * cos(2. * omega);
  deps = deps / 3600.;
  return (deps);
}


x86 assembly (clang, -O0):
<NObliquity>:
push   %rbp
mov    %rsp,%rbp
sub    $0x50,%rsp
movsd  0xf60(%rip),%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  0xf4b(%rip),%xmm0
addsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x30(%rbp),%xmm0
movsd  0xf2c(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    $0x0,%al
call   1040 <JDNow@plt>
movsd  %xmm0,-0x28(%rbp)
movsd  -0x28(%rbp),%xmm0
movsd  0xf02(%rip),%xmm1
subsd  %xmm1,%xmm0
movsd  0xeee(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x38(%rbp),%xmm0
addsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x38(%rbp),%xmm1
movsd  0xeb9(%rip),%xmm0
mulsd  %xmm0,%xmm1
movsd  0xeb5(%rip),%xmm0
addsd  %xmm0,%xmm1
movsd  0xe99(%rip),%xmm0
mulsd  -0x38(%rbp),%xmm0
movsd  -0x38(%rbp),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  -0x38(%rbp),%xmm1
mulsd  -0x38(%rbp),%xmm1
mulsd  -0x38(%rbp),%xmm1
movsd  0xe68(%rip),%xmm2
divsd  %xmm2,%xmm1
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x20(%rbp)
movsd  -0x38(%rbp),%xmm0
movsd  0xe3e(%rip),%xmm1
mulsd  %xmm1,%xmm0
movsd  0xe3a(%rip),%xmm1
addsd  %xmm1,%xmm0
call   1030 <Map360@plt>
movsd  %xmm0,-0x10(%rbp)
movsd  -0x38(%rbp),%xmm0
movsd  0xe07(%rip),%xmm1
mulsd  %xmm1,%xmm0
movsd  0xe03(%rip),%xmm1
addsd  %xmm1,%xmm0
call   1030 <Map360@plt>
movsd  %xmm0,-0x18(%rbp)
movsd  0xddd(%rip),%xmm0
mulsd  -0x20(%rbp),%xmm0
movsd  0xdc8(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x20(%rbp)
movsd  0xdbf(%rip),%xmm0
mulsd  -0x10(%rbp),%xmm0
movsd  0xdaa(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x10(%rbp)
movsd  0xda1(%rip),%xmm0
mulsd  -0x18(%rbp),%xmm0
movsd  0xd8c(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x18(%rbp)
cvttsd2si -0x20(%rbp),%edi
call   1050 <cos@plt>
movsd  %xmm0,-0x50(%rbp)
movsd  0xd4c(%rip),%xmm0
mulsd  -0x10(%rbp),%xmm0
cvttsd2si %xmm0,%edi
call   1050 <cos@plt>
movaps %xmm0,%xmm2
movsd  -0x50(%rbp),%xmm0
movsd  0xd46(%rip),%xmm1
mulsd  %xmm2,%xmm1
movsd  0xd32(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  0xd0d(%rip),%xmm0
mulsd  -0x18(%rbp),%xmm0
cvttsd2si %xmm0,%edi
call   1050 <cos@plt>
movsd  -0x48(%rbp),%xmm1
movsd  0xcfa(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  0xcdd(%rip),%xmm0
mulsd  -0x20(%rbp),%xmm0
cvttsd2si %xmm0,%edi
call   1050 <cos@plt>
movsd  -0x40(%rbp),%xmm1
movsd  0xcba(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
movsd  0xc98(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
add    $0x50,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<NObliquity>:
sub    $0x28,%rsp
xor    %eax,%eax
call   1040 <JDNow@plt>
addsd  0xebd(%rip),%xmm0
divsd  0xebd(%rip),%xmm0
addsd  0xebd(%rip),%xmm0
movsd  0xebd(%rip),%xmm2
mulsd  %xmm0,%xmm2
addsd  0xeb9(%rip),%xmm2
movsd  0xeb9(%rip),%xmm1
mulsd  %xmm0,%xmm1
mulsd  %xmm0,%xmm1
addsd  %xmm2,%xmm1
movapd %xmm0,%xmm2
mulsd  %xmm0,%xmm2
mulsd  %xmm0,%xmm2
movapd %xmm0,%xmm3
movsd  %xmm0,0x10(%rsp)
divsd  0xe97(%rip),%xmm2
addsd  %xmm1,%xmm2
movsd  %xmm2,0x18(%rsp)
movsd  0xe8d(%rip),%xmm0
mulsd  %xmm3,%xmm0
addsd  0xe89(%rip),%xmm0
call   1030 <Map360@plt>
movsd  %xmm0,0x8(%rsp)
movsd  0x10(%rsp),%xmm0
mulsd  0xe78(%rip),%xmm0
addsd  0xe78(%rip),%xmm0
call   1030 <Map360@plt>
movsd  0xe73(%rip),%xmm3
movsd  0x18(%rsp),%xmm2
mulsd  %xmm3,%xmm2
movsd  0xe69(%rip),%xmm1
divsd  %xmm1,%xmm2
movsd  %xmm2,0x18(%rsp)
movsd  0x8(%rsp),%xmm4
mulsd  %xmm3,%xmm4
divsd  %xmm1,%xmm4
movsd  %xmm4,0x8(%rsp)
mulsd  %xmm3,%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,0x10(%rsp)
cvttsd2si %xmm2,%edi
call   1050 <cos@plt>
movsd  %xmm0,0x20(%rsp)
movsd  0x8(%rsp),%xmm0
addsd  %xmm0,%xmm0
cvttsd2si %xmm0,%edi
call   1050 <cos@plt>
mulsd  0xe1b(%rip),%xmm0
movsd  0x20(%rsp),%xmm1
mulsd  0xe15(%rip),%xmm1
addsd  %xmm1,%xmm0
movsd  %xmm0,0x8(%rsp)
movsd  0x10(%rsp),%xmm0
addsd  %xmm0,%xmm0
cvttsd2si %xmm0,%edi
call   1050 <cos@plt>
mulsd  0xdf8(%rip),%xmm0
addsd  0x8(%rsp),%xmm0
movsd  %xmm0,0x8(%rsp)
movsd  0x18(%rsp),%xmm0
addsd  %xmm0,%xmm0
cvttsd2si %xmm0,%edi
call   1050 <cos@plt>
mulsd  0xdd9(%rip),%xmm0
addsd  0x8(%rsp),%xmm0
divsd  0xdd3(%rip),%xmm0
add    $0x28,%rsp
ret

LLVM IR (-O0):
define double @NObliquity() {
%1 = alloca double
%2 = alloca double
%3 = alloca double
%4 = alloca double
%5 = alloca double
%6 = alloca double
%7 = alloca double
store double 3.500000e+01, double* %6
%8 = load double, double* %6
%9 = fadd double %8, 3.218400e+01
store double %9, double* %6
%10 = load double, double* %6
%11 = fdiv double %10, 3.155760e+09
store double %11, double* %6
%12 = call double (...) @JDNow()
store double %12, double* %5
%13 = load double, double* %5
%14 = fsub double %13, 0x4142B42C80000000
%15 = fdiv double %14, 3.652500e+04
store double %15, double* %7
%16 = load double, double* %7
%17 = load double, double* %6
%18 = fadd double %16, %17
store double %18, double* %7
%19 = load double, double* %7
%20 = call double @llvm.fmuladd.f64(double 0xC09E388B8800EAE2, double %19, double 0x405F42D96A6A0126)
%21 = load double, double* %7
%22 = fmul double 2.070800e-03, %21
%23 = load double, double* %7
%24 = call double @llvm.fmuladd.f64(double %22, double %23, double %20)
%25 = load double, double* %7
%26 = load double, double* %7
%27 = fmul double %25, %26
%28 = load double, double* %7
%29 = fmul double %27, %28
%30 = fdiv double %29, 4.500000e+05
%31 = fadd double %24, %30
store double %31, double* %4
%32 = load double, double* %7
%33 = call double @llvm.fmuladd.f64(double 0x40AC218A2339C0EC, double %32, double 0x40718776C8B43958)
%34 = call double @Map360(double %33)
store double %34, double* %2
%35 = load double, double* %7
%36 = call double @llvm.fmuladd.f64(double 0x411D5FCF867381D8, double %35, double 0x406B4A20C49BA5E3)
%37 = call double @Map360(double %36)
store double %37, double* %3
%38 = load double, double* %4
%39 = fmul double %38, 0x400921FB5443D6F4
%40 = fdiv double %39, 1.800000e+02
store double %40, double* %4
%41 = load double, double* %2
%42 = fmul double %41, 0x400921FB5443D6F4
%43 = fdiv double %42, 1.800000e+02
store double %43, double* %2
%44 = load double, double* %3
%45 = fmul double %44, 0x400921FB5443D6F4
%46 = fdiv double %45, 1.800000e+02
store double %46, double* %3
%47 = load double, double* %4
%48 = fptosi double %47 to i32
%49 = call double @cos(i32 %48)
%50 = load double, double* %2
%51 = fmul double 2.000000e+00, %50
%52 = fptosi double %51 to i32
%53 = call double @cos(i32 %52)
%54 = fmul double 5.700000e-01, %53
%55 = call double @llvm.fmuladd.f64(double 0x4022666666666666, double %49, double %54)
%56 = load double, double* %3
%57 = fmul double 2.000000e+00, %56
%58 = fptosi double %57 to i32
%59 = call double @cos(i32 %58)
%60 = call double @llvm.fmuladd.f64(double 1.000000e-01, double %59, double %55)
%61 = load double, double* %4
%62 = fmul double 2.000000e+00, %61
%63 = fptosi double %62 to i32
%64 = call double @cos(i32 %63)
%65 = call double @llvm.fmuladd.f64(double 0xBFB70A3D70A3D70A, double %64, double %60)
store double %65, double* %1
%66 = load double, double* %1
%67 = fdiv double %66, 3.600000e+03
store double %67, double* %1
%68 = load double, double* %1
ret double %68
}

LLVM IR (-O2):
define double @NObliquity() {
%1 = tail call double (...) @JDNow()
%2 = fadd double %1, 0xC142B42C80000000
%3 = fdiv double %2, 3.652500e+04
%4 = fadd double %3, 0x3E56DBF70F0D61B1
%5 = tail call double @llvm.fmuladd.f64(double %4, double 0xC09E388B8800EAE2, double 0x405F42D96A6A0126)
%6 = fmul double %4, 2.070800e-03
%7 = tail call double @llvm.fmuladd.f64(double %6, double %4, double %5)
%8 = fmul double %4, %4
%9 = fmul double %4, %8
%10 = fdiv double %9, 4.500000e+05
%11 = fadd double %7, %10
%12 = tail call double @llvm.fmuladd.f64(double %4, double 0x40AC218A2339C0EC, double 0x40718776C8B43958)
%13 = tail call double @Map360(double %12)
%14 = tail call double @llvm.fmuladd.f64(double %4, double 0x411D5FCF867381D8, double 0x406B4A20C49BA5E3)
%15 = tail call double @Map360(double %14)
%16 = fmul double %11, 0x400921FB5443D6F4
%17 = fdiv double %16, 1.800000e+02
%18 = fmul double %13, 0x400921FB5443D6F4
%19 = fdiv double %18, 1.800000e+02
%20 = fmul double %15, 0x400921FB5443D6F4
%21 = fdiv double %20, 1.800000e+02
%22 = fptosi double %17 to i32
%23 = tail call double @cos(i32 %22)
%24 = fmul double %19, 2.000000e+00
%25 = fptosi double %24 to i32
%26 = tail call double @cos(i32 %25)
%27 = fmul double %26, 5.700000e-01
%28 = tail call double @llvm.fmuladd.f64(double %23, double 0x4022666666666666, double %27)
%29 = fmul double %21, 2.000000e+00
%30 = fptosi double %29 to i32
%31 = tail call double @cos(i32 %30)
%32 = tail call double @llvm.fmuladd.f64(double %31, double 1.000000e-01, double %28)
%33 = fmul double %17, 2.000000e+00
%34 = fptosi double %33 to i32
%35 = tail call double @cos(i32 %34)
%36 = tail call double @llvm.fmuladd.f64(double %35, double 0xBFB70A3D70A3D70A, double %32)
%37 = fdiv double %36, 3.600000e+03
ret double %37
}