double julian_date(int y, int m, int d) {
  int A, B, C, D;
  double JD;
  if ((m == 1) || (m == 2)) {
    y -= 1;
    m += 12;
  }
  A = y / 100;
  B = 2 - A + (A / 4);
  C = (int)(365.25 * y);
  D = (int)(30.6001 * (m + 1));
  JD = B + C + D + d + 1720994.5;
  return JD;
}


x86 assembly (clang, -O0):
<julian_date>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %edx,-0xc(%rbp)
cmpl   $0x1,-0x8(%rbp)
je     1121 <julian_date+0x21>
cmpl   $0x2,-0x8(%rbp)
jne    1133 <julian_date+0x33>
mov    -0x4(%rbp),%eax
sub    $0x1,%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,-0x10(%rbp)
mov    $0x2,%eax
sub    -0x10(%rbp),%eax
mov    %eax,-0x2c(%rbp)
mov    -0x10(%rbp),%eax
mov    $0x4,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x2c(%rbp),%eax
add    %ecx,%eax
mov    %eax,-0x14(%rbp)
cvtsi2sdl -0x4(%rbp),%xmm1
movsd  0xea2(%rip),%xmm0
mulsd  %xmm1,%xmm0
cvttsd2si %xmm0,%eax
mov    %eax,-0x18(%rbp)
mov    -0x8(%rbp),%eax
add    $0x1,%eax
cvtsi2sd %eax,%xmm1
movsd  0xe7d(%rip),%xmm0
mulsd  %xmm1,%xmm0
cvttsd2si %xmm0,%eax
mov    %eax,-0x1c(%rbp)
mov    -0x14(%rbp),%eax
add    -0x18(%rbp),%eax
add    -0x1c(%rbp),%eax
add    -0xc(%rbp),%eax
cvtsi2sd %eax,%xmm0
movsd  0xe52(%rip),%xmm1
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x28(%rbp)
movsd  -0x28(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<julian_date>:
lea    -0x1(%rsi),%eax
lea    0xc(%rsi),%ecx
cmp    $0x2,%eax
cmovae %esi,%ecx
sbb    $0x0,%edi
movslq %edi,%r8
imul   $0xffffffffae147ae1,%r8,%rax
mov    %rax,%rdi
shr    $0x3f,%rdi
sar    $0x25,%rax
add    %edi,%eax
imul   $0x51eb851f,%r8,%rdi
mov    %rdi,%rsi
shr    $0x3f,%rsi
sar    $0x27,%rdi
cvtsi2sd %r8d,%xmm0
mulsd  0xebb(%rip),%xmm0
cvttsd2si %xmm0,%r8d
add    %esi,%edi
add    $0x1,%ecx
xorps  %xmm0,%xmm0
cvtsi2sd %ecx,%xmm0
mulsd  0xeaa(%rip),%xmm0
cvttsd2si %xmm0,%ecx
add    %edx,%eax
add    %edi,%eax
add    %r8d,%eax
add    %ecx,%eax
add    $0x2,%eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
addsd  0xe93(%rip),%xmm0
ret

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

LLVM IR (-O2):
define double @julian_date(i32 %0, i32 %1, i32 %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 = sdiv i32 %9, 400
%12 = sitofp i32 %9 to double
%13 = fmul double %12, 3.652500e+02
%14 = fptosi double %13 to i32
%15 = add nsw i32 %7, 1
%16 = sitofp i32 %15 to double
%17 = fmul double %16, 3.060010e+01
%18 = fptosi double %17 to i32
%19 = add i32 %2, 2
%20 = add i32 %19, %10
%21 = add i32 %20, %11
%22 = add i32 %21, %14
%23 = add i32 %22, %18
%24 = sitofp i32 %23 to double
%25 = fadd double %24, 0x413A42A280000000
ret double %25
}