// File: matthewarcus[P]mmps/equations.cpp
// Function: equationoftime(double)
double equationoftime(double date)
{
  double p = perihelion;
  double m = date - periheliontime;
  double s = 
    -591.7 * sin(2 * (m+p))
    -459.6 * sin(m) +
     +19.8 * sin(m + 2 * p)
     -19.8 * sin(3 * m + 2 * p)
     -12.8 * sin(4 * (m + p))
      -4.8 * sin(2 * m)
      +0.9 * sin(3 * m + 4 * p)
      -0.9 * sin(5 * m + 4 * p)
      -0.5 * sin(4 * m + 2 * p)
      -0.4 * sin(6 * (m + p));
    
  // fprintf(stderr, "EOT: %f\n", s / 60.0);
  return s; // seconds
}

// Assembly:
pushq %rbp
movq %rsp, %rbp
subq $0x70, %rsp
movsd %xmm0, -0x8(%rbp)
movsd 0xd993(%rip), %xmm0     # 0x1d2f8
movsd %xmm0, -0x10(%rbp)
movsd -0x8(%rbp), %xmm0
movsd 0x5e69(%rip), %xmm1     # 0x157e0
subsd %xmm1, %xmm0
movsd %xmm0, -0x18(%rbp)
movsd -0x18(%rbp), %xmm1
addsd -0x10(%rbp), %xmm1
movsd 0x5a8e(%rip), %xmm0     # 0x15420
mulsd %xmm1, %xmm0
callq 0x5080
movsd %xmm0, -0x68(%rbp)
movsd -0x18(%rbp), %xmm0
callq 0x5080
movaps %xmm0, %xmm2
movsd -0x68(%rbp), %xmm0
movsd 0x5e1e(%rip), %xmm1     # 0x157d8
mulsd %xmm2, %xmm1
movq %xmm1, %rax
movabsq $-0x8000000000000000, %rcx # imm = 0x8000000000000000
xorq %rcx, %rax
movq %rax, %xmm1
movsd 0x5df3(%rip), %xmm2     # 0x157d0
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd %xmm0, -0x60(%rbp)
movsd -0x18(%rbp), %xmm1
movsd -0x10(%rbp), %xmm0
addsd %xmm0, %xmm0
addsd %xmm1, %xmm0
callq 0x5080
movsd -0x60(%rbp), %xmm1
movsd 0x5dba(%rip), %xmm2     # 0x157c8
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd %xmm0, -0x58(%rbp)
movsd -0x18(%rbp), %xmm0
movsd 0x59f8(%rip), %xmm1     # 0x15420
mulsd -0x10(%rbp), %xmm1
movsd 0x5d43(%rip), %xmm2     # 0x15778
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
callq 0x5080
movsd -0x58(%rbp), %xmm1
movsd 0x5d71(%rip), %xmm2     # 0x157c0
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd %xmm0, -0x50(%rbp)
movsd -0x18(%rbp), %xmm1
addsd -0x10(%rbp), %xmm1
movsd 0x59c2(%rip), %xmm0     # 0x15430
mulsd %xmm1, %xmm0
callq 0x5080
movsd -0x50(%rbp), %xmm1
movsd 0x5d34(%rip), %xmm2     # 0x157b8
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd %xmm0, -0x48(%rbp)
movsd 0x5987(%rip), %xmm0     # 0x15420
mulsd -0x18(%rbp), %xmm0
callq 0x5080
movsd -0x48(%rbp), %xmm1
movsd 0x5d00(%rip), %xmm2     # 0x157b0
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd %xmm0, -0x40(%rbp)
movsd -0x18(%rbp), %xmm0
movsd 0x5966(%rip), %xmm1     # 0x15430
mulsd -0x10(%rbp), %xmm1
movsd 0x5ca1(%rip), %xmm2     # 0x15778
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
callq 0x5080
movsd -0x40(%rbp), %xmm1
movsd 0x5cb7(%rip), %xmm2     # 0x157a8
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd %xmm0, -0x38(%rbp)
movsd -0x18(%rbp), %xmm0
movsd 0x5925(%rip), %xmm1     # 0x15430
mulsd -0x10(%rbp), %xmm1
movsd 0x5c68(%rip), %xmm2     # 0x15780
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
callq 0x5080
movsd -0x38(%rbp), %xmm1
movsd 0x5c6e(%rip), %xmm2     # 0x157a0
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd %xmm0, -0x30(%rbp)
movsd -0x18(%rbp), %xmm0
movsd 0x58d4(%rip), %xmm1     # 0x15420
mulsd -0x10(%rbp), %xmm1
movsd 0x58d7(%rip), %xmm2     # 0x15430
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
callq 0x5080
movsd -0x30(%rbp), %xmm1
movsd 0x58dd(%rip), %xmm2     # 0x15450
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd %xmm0, -0x28(%rbp)
movsd -0x18(%rbp), %xmm1
addsd -0x10(%rbp), %xmm1
movsd 0x58ce(%rip), %xmm0     # 0x15460
mulsd %xmm1, %xmm0
callq 0x5080
movsd -0x28(%rbp), %xmm1
movsd 0x5bf0(%rip), %xmm2     # 0x15798
mulsd %xmm2, %xmm0
addsd %xmm1, %xmm0
movsd %xmm0, -0x20(%rbp)
movsd -0x20(%rbp), %xmm0
addq $0x70, %rsp
popq %rbp
retq