double gregorian_calendar_to_jd(int y, int m, int d, int h, int mi, int s) {
  y += 8000;
  if (m < 3) {
    y--;
    m += 12;
  }
  return (y * 365) + (y / 4) - (y / 100) + (y / 400) - 1200820 +
         (m * 153 + 3) / 5 - 92 + d - 1 + h / 24.0 + mi / 1440.0 + s / 86400.0 -
         0.5;
}


x86 assembly (clang, -O0):
<gregorian_calendar_to_jd>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %edx,-0xc(%rbp)
mov    %ecx,-0x10(%rbp)
mov    %r8d,-0x14(%rbp)
mov    %r9d,-0x18(%rbp)
mov    -0x4(%rbp),%eax
add    $0x1f40,%eax
mov    %eax,-0x4(%rbp)
cmpl   $0x3,-0x8(%rbp)
jge    113f <gregorian_calendar_to_jd+0x3f>
mov    -0x4(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x4(%rbp)
mov    -0x8(%rbp),%eax
add    $0xc,%eax
mov    %eax,-0x8(%rbp)
imul   $0x16d,-0x4(%rbp),%eax
mov    %eax,-0x28(%rbp)
mov    -0x4(%rbp),%eax
mov    $0x4,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x28(%rbp),%eax
add    %ecx,%eax
mov    %eax,-0x24(%rbp)
mov    -0x4(%rbp),%eax
mov    $0x64,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x24(%rbp),%eax
sub    %ecx,%eax
mov    %eax,-0x20(%rbp)
mov    -0x4(%rbp),%eax
mov    $0x190,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x20(%rbp),%eax
add    %ecx,%eax
sub    $0x1252b4,%eax
mov    %eax,-0x1c(%rbp)
imul   $0x99,-0x8(%rbp),%eax
add    $0x3,%eax
mov    $0x5,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x1c(%rbp),%eax
add    %ecx,%eax
sub    $0x5c,%eax
add    -0xc(%rbp),%eax
sub    $0x1,%eax
cvtsi2sd %eax,%xmm0
cvtsi2sdl -0x10(%rbp),%xmm1
movsd  0xe58(%rip),%xmm2
divsd  %xmm2,%xmm1
addsd  %xmm1,%xmm0
cvtsi2sdl -0x14(%rbp),%xmm1
movsd  0xe3b(%rip),%xmm2
divsd  %xmm2,%xmm1
addsd  %xmm1,%xmm0
cvtsi2sdl -0x18(%rbp),%xmm1
movsd  0xe1e(%rip),%xmm2
divsd  %xmm2,%xmm1
addsd  %xmm1,%xmm0
movsd  0xe06(%rip),%xmm1
subsd  %xmm1,%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<gregorian_calendar_to_jd>:
push   %r14
push   %rbx
lea    0xc(%rsi),%r11d
xor    %eax,%eax
cmp    $0x3,%esi
setl   %al
cmovge %esi,%r11d
sub    %eax,%edi
lea    0x1f40(%rdi),%eax
imul   $0x16d,%eax,%r10d
add    $0x1f43,%edi
test   %eax,%eax
cmovns %eax,%edi
sar    $0x2,%edi
cltq
imul   $0xffffffffae147ae1,%rax,%r14
mov    %r14,%rbx
shr    $0x3f,%rbx
sar    $0x25,%r14
add    %ebx,%r14d
imul   $0x51eb851f,%rax,%rax
mov    %rax,%rbx
shr    $0x3f,%rbx
sar    $0x27,%rax
add    %ebx,%eax
imul   $0x99,%r11d,%ebx
add    $0x3,%ebx
movslq %ebx,%rbx
imul   $0x66666667,%rbx,%rbx
mov    %rbx,%rsi
shr    $0x3f,%rsi
sar    $0x21,%rbx
add    %esi,%ebx
add    %edx,%edi
add    %r10d,%edi
add    %r14d,%edi
add    %eax,%edi
lea    (%rdi,%rbx,1),%eax
add    $0xffedacef,%eax
cvtsi2sd %eax,%xmm0
movd   %r8d,%xmm1
movd   %r9d,%xmm2
punpckldq %xmm1,%xmm2
xorps  %xmm1,%xmm1
cvtsi2sd %ecx,%xmm1
divsd  0xe51(%rip),%xmm1
cvtdq2pd %xmm2,%xmm2
addsd  %xmm0,%xmm1
divpd  0xe51(%rip),%xmm2
movapd %xmm2,%xmm0
unpckhpd %xmm2,%xmm0
addsd  %xmm1,%xmm0
addsd  %xmm2,%xmm0
addsd  0xe31(%rip),%xmm0
pop    %rbx
pop    %r14
ret

LLVM IR (-O0):
define double @gregorian_calendar_to_jd(i32 %0, i32 %1, i32 %2, i32 %3, i32 %4, i32 %5) {
%7 = alloca i32
%8 = alloca i32
%9 = alloca i32
%10 = alloca i32
%11 = alloca i32
%12 = alloca i32
store i32 %0, i32* %7
store i32 %1, i32* %8
store i32 %2, i32* %9
store i32 %3, i32* %10
store i32 %4, i32* %11
store i32 %5, i32* %12
%13 = load i32, i32* %7
%14 = add nsw i32 %13, 8000
store i32 %14, i32* %7
%15 = load i32, i32* %8
%16 = icmp slt i32 %15, 3
br i1 %16, label %17, label %22
17: ; preds = %6
%18 = load i32, i32* %7
%19 = add nsw i32 %18, -1
store i32 %19, i32* %7
%20 = load i32, i32* %8
%21 = add nsw i32 %20, 12
store i32 %21, i32* %8
br label %22
22: ; preds = %17, %6
%23 = load i32, i32* %7
%24 = mul nsw i32 %23, 365
%25 = load i32, i32* %7
%26 = sdiv i32 %25, 4
%27 = add nsw i32 %24, %26
%28 = load i32, i32* %7
%29 = sdiv i32 %28, 100
%30 = sub nsw i32 %27, %29
%31 = load i32, i32* %7
%32 = sdiv i32 %31, 400
%33 = add nsw i32 %30, %32
%34 = sub nsw i32 %33, 1200820
%35 = load i32, i32* %8
%36 = mul nsw i32 %35, 153
%37 = add nsw i32 %36, 3
%38 = sdiv i32 %37, 5
%39 = add nsw i32 %34, %38
%40 = sub nsw i32 %39, 92
%41 = load i32, i32* %9
%42 = add nsw i32 %40, %41
%43 = sub nsw i32 %42, 1
%44 = sitofp i32 %43 to double
%45 = load i32, i32* %10
%46 = sitofp i32 %45 to double
%47 = fdiv double %46, 2.400000e+01
%48 = fadd double %44, %47
%49 = load i32, i32* %11
%50 = sitofp i32 %49 to double
%51 = fdiv double %50, 1.440000e+03
%52 = fadd double %48, %51
%53 = load i32, i32* %12
%54 = sitofp i32 %53 to double
%55 = fdiv double %54, 8.640000e+04
%56 = fadd double %52, %55
%57 = fsub double %56, 5.000000e-01
ret double %57
}

LLVM IR (-O2):
define double @gregorian_calendar_to_jd(i32 %0, i32 %1, i32 %2, i32 %3, i32 %4, i32 %5) {
%7 = icmp slt i32 %1, 3
%8 = add nsw i32 %1, 12
%9 = select i1 %7, i32 %8, i32 %1
%10 = select i1 %7, i32 7999, i32 8000
%11 = add nsw i32 %10, %0
%12 = mul nsw i32 %11, 365
%13 = sdiv i32 %11, 4
%14 = sdiv i32 %11, -100
%15 = sdiv i32 %11, 400
%16 = mul nsw i32 %9, 153
%17 = add nsw i32 %16, 3
%18 = sdiv i32 %17, 5
%19 = add i32 %2, -1200913
%20 = add i32 %19, %13
%21 = add i32 %20, %12
%22 = add i32 %21, %14
%23 = add i32 %22, %15
%24 = add i32 %23, %18
%25 = sitofp i32 %24 to double
%26 = sitofp i32 %3 to double
%27 = fdiv double %26, 2.400000e+01
%28 = fadd double %27, %25
%29 = insertelement <2 x i32> poison, i32 %5, i64 0
%30 = insertelement <2 x i32> %29, i32 %4, i64 1
%31 = sitofp <2 x i32> %30 to <2 x double>
%32 = fdiv <2 x double> %31, <double 8.640000e+04, double 1.440000e+03>
%33 = extractelement <2 x double> %32, i64 1
%34 = fadd double %33, %28
%35 = extractelement <2 x double> %32, i64 0
%36 = fadd double %35, %34
%37 = fadd double %36, -5.000000e-01
ret double %37
}