double JD(int y, int m, double d) {
  if (m == 1 || m == 2) {
    y--;
    m += 12;
  }
  int A = y / 100;
  int B = 2 - A + A / 4;
  double JD;
  JD = 36525 * (y + 4716) / 100 + 306 * (m + 1) / 10 + d + B - 1524.5;
  return JD;
}


x86 assembly (clang, -O0):
<JD>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
movsd  %xmm0,-0x10(%rbp)
cmpl   $0x1,-0x8(%rbp)
je     1123 <JD+0x23>
cmpl   $0x2,-0x8(%rbp)
jne    1135 <JD+0x35>
mov    -0x4(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x4(%rbp)
mov    -0x8(%rbp),%eax
add    $0xc,%eax
mov    %eax,-0x8(%rbp)
mov    -0x4(%rbp),%eax
mov    $0x64,%ecx
cltd
idiv   %ecx
mov    %eax,-0x14(%rbp)
mov    $0x2,%eax
sub    -0x14(%rbp),%eax
mov    %eax,-0x28(%rbp)
mov    -0x14(%rbp),%eax
mov    $0x4,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x28(%rbp),%eax
add    %ecx,%eax
mov    %eax,-0x18(%rbp)
mov    -0x4(%rbp),%eax
add    $0x126c,%eax
imul   $0x8ead,%eax,%eax
mov    $0x64,%ecx
cltd
idiv   %ecx
mov    %eax,-0x24(%rbp)
mov    -0x8(%rbp),%eax
add    $0x1,%eax
imul   $0x132,%eax,%eax
mov    $0xa,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x24(%rbp),%eax
add    %ecx,%eax
cvtsi2sd %eax,%xmm0
addsd  -0x10(%rbp),%xmm0
cvtsi2sdl -0x18(%rbp),%xmm1
addsd  %xmm1,%xmm0
movsd  0xe4f(%rip),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x20(%rbp)
movsd  -0x20(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<JD>:
lea    -0x1(%rsi),%ecx
lea    0xc(%rsi),%eax
cmp    $0x2,%ecx
cmovae %esi,%eax
sbb    $0x0,%edi
movslq %edi,%rdx
imul   $0xffffffffae147ae1,%rdx,%rcx
mov    %rcx,%rsi
shr    $0x3f,%rsi
sar    $0x25,%rcx
add    %esi,%ecx
imul   $0x51eb851f,%rdx,%rsi
mov    %rsi,%rdi
shr    $0x3f,%rdi
sar    $0x27,%rsi
add    %edi,%esi
add    %esi,%ecx
add    $0x2,%ecx
imul   $0x8ead,%edx,%edx
add    $0xa445afc,%edx
movslq %edx,%rdx
imul   $0x51eb851f,%rdx,%rdx
mov    %rdx,%rsi
shr    $0x3f,%rsi
sar    $0x25,%rdx
add    %esi,%edx
imul   $0x132,%eax,%eax
add    $0x132,%eax
cltq
imul   $0x66666667,%rax,%rax
mov    %rax,%rsi
shr    $0x3f,%rsi
sar    $0x22,%rax
add    %esi,%eax
add    %edx,%eax
cvtsi2sd %eax,%xmm1
cvtsi2sd %ecx,%xmm2
addsd  %xmm1,%xmm0
addsd  %xmm2,%xmm0
addsd  0xe63(%rip),%xmm0
ret

LLVM IR (-O0):
define double @JD(i32 %0, i32 %1, double %2) {
%4 = alloca i32
%5 = alloca i32
%6 = alloca double
%7 = alloca i32
%8 = alloca i32
%9 = alloca double
store i32 %0, i32* %4
store i32 %1, i32* %5
store double %2, double* %6
%10 = load i32, i32* %5
%11 = icmp eq i32 %10, 1
br i1 %11, label %15, label %12
12: ; preds = %3
%13 = load i32, i32* %5
%14 = icmp eq i32 %13, 2
br i1 %14, label %15, label %20
15: ; preds = %12, %3
%16 = load i32, i32* %4
%17 = add nsw i32 %16, -1
store i32 %17, i32* %4
%18 = load i32, i32* %5
%19 = add nsw i32 %18, 12
store i32 %19, i32* %5
br label %20
20: ; preds = %15, %12
%21 = load i32, i32* %4
%22 = sdiv i32 %21, 100
store i32 %22, i32* %7
%23 = load i32, i32* %7
%24 = sub nsw i32 2, %23
%25 = load i32, i32* %7
%26 = sdiv i32 %25, 4
%27 = add nsw i32 %24, %26
store i32 %27, i32* %8
%28 = load i32, i32* %4
%29 = add nsw i32 %28, 4716
%30 = mul nsw i32 36525, %29
%31 = sdiv i32 %30, 100
%32 = load i32, i32* %5
%33 = add nsw i32 %32, 1
%34 = mul nsw i32 306, %33
%35 = sdiv i32 %34, 10
%36 = add nsw i32 %31, %35
%37 = sitofp i32 %36 to double
%38 = load double, double* %6
%39 = fadd double %37, %38
%40 = load i32, i32* %8
%41 = sitofp i32 %40 to double
%42 = fadd double %39, %41
%43 = fsub double %42, 1.524500e+03
store double %43, double* %9
%44 = load double, double* %9
ret double %44
}

LLVM IR (-O2):
define double @JD(i32 %0, i32 %1, double %2) {
%4 = add i32 %1, -1
%5 = icmp ult i32 %4, 2
%6 = add nsw i32 %1, 12
%7 = select i1 %5, i32 %6, i32 %1
%8 = sext i1 %5 to i32
%9 = add nsw i32 %8, %0
%10 = sdiv i32 %9, -100
%11 = add nsw i32 %10, 2
%12 = sdiv i32 %9, 400
%13 = add nsw i32 %11, %12
%14 = mul i32 %9, 36525
%15 = add i32 %14, 172251900
%16 = sdiv i32 %15, 100
%17 = mul i32 %7, 306
%18 = add i32 %17, 306
%19 = sdiv i32 %18, 10
%20 = add nsw i32 %16, %19
%21 = sitofp i32 %20 to double
%22 = fadd double %21, %2
%23 = sitofp i32 %13 to double
%24 = fadd double %22, %23
%25 = fadd double %24, -1.524500e+03
ret double %25
}