double Julian(int year, int month, double day) {
  int a, b, c, e;
  if (month < 3) {
    year = year - 1;
    month = month + 12;
  }
  if (((year > 1582)) || ((year == 1582) && (month > 10)) ||
      ((year == 1582) && (month == 10) && (day > 15))) {
    a = year / 100;
    b = 2 - a + a / 4;
  }
  c = (long)(365.25 * year);
  e = (long)(30.6001 * (month + 1));
  return b + c + e + day + 1720994.5;
}


x86 assembly (clang, -O0):
<Julian>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
movsd  %xmm0,-0x10(%rbp)
cmpl   $0x3,-0x8(%rbp)
jge    112b <Julian+0x2b>
mov    -0x4(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x4(%rbp)
mov    -0x8(%rbp),%eax
add    $0xc,%eax
mov    %eax,-0x8(%rbp)
cmpl   $0x62e,-0x4(%rbp)
jg     117d <Julian+0x7d>
cmpl   $0x62e,-0x4(%rbp)
jne    114f <Julian+0x4f>
cmpl   $0xa,-0x8(%rbp)
jg     117d <Julian+0x7d>
cmpl   $0x62e,-0x4(%rbp)
jne    11ab <Julian+0xab>
cmpl   $0xa,-0x8(%rbp)
jne    11ab <Julian+0xab>
movsd  -0x10(%rbp),%xmm0
movsd  0xe8d(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    11ab <Julian+0xab>
mov    -0x4(%rbp),%eax
mov    $0x64,%ecx
cltd
idiv   %ecx
mov    %eax,-0x14(%rbp)
mov    $0x2,%eax
sub    -0x14(%rbp),%eax
mov    %eax,-0x24(%rbp)
mov    -0x14(%rbp),%eax
mov    $0x4,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x24(%rbp),%eax
add    %ecx,%eax
mov    %eax,-0x18(%rbp)
cvtsi2sdl -0x4(%rbp),%xmm1
movsd  0xe60(%rip),%xmm0
mulsd  %xmm1,%xmm0
cvttsd2si %xmm0,%rax
mov    %eax,-0x1c(%rbp)
mov    -0x8(%rbp),%eax
add    $0x1,%eax
cvtsi2sd %eax,%xmm1
movsd  0xe3a(%rip),%xmm0
mulsd  %xmm1,%xmm0
cvttsd2si %xmm0,%rax
mov    %eax,-0x20(%rbp)
mov    -0x18(%rbp),%eax
add    -0x1c(%rbp),%eax
add    -0x20(%rbp),%eax
cvtsi2sd %eax,%xmm0
addsd  -0x10(%rbp),%xmm0
movsd  0xe0c(%rip),%xmm1
addsd  %xmm1,%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<Julian>:
lea    0xc(%rsi),%eax
xor    %ecx,%ecx
cmp    $0x3,%esi
setl   %cl
cmovge %esi,%eax
sub    %ecx,%edi
cmp    $0x62e,%edi
jg     1136 <Julian+0x36>
jne    111f <Julian+0x1f>
cmp    $0xa,%eax
jg     1136 <Julian+0x36>
ucomisd 0xed9(%rip),%xmm0
jbe    114f <Julian+0x4f>
cmp    $0x62e,%edi
jne    114f <Julian+0x4f>
cmp    $0xa,%eax
jne    114f <Julian+0x4f>
mov    %edi,%ecx
imul   $0x51eb851f,%rcx,%rcx
mov    %rcx,%rdx
shr    $0x25,%rdx
shr    $0x27,%rcx
sub    %edx,%ecx
add    $0x2,%ecx
cvtsi2sd %edi,%xmm1
mulsd  0xead(%rip),%xmm1
cvttsd2si %xmm1,%rdx
add    $0x1,%eax
xorps  %xmm1,%xmm1
cvtsi2sd %eax,%xmm1
mulsd  0xe9e(%rip),%xmm1
cvttsd2si %xmm1,%rax
add    %edx,%eax
add    %ecx,%eax
xorps  %xmm1,%xmm1
cvtsi2sd %eax,%xmm1
addsd  %xmm0,%xmm1
addsd  0xe8a(%rip),%xmm1
movapd %xmm1,%xmm0
ret

LLVM IR (-O0):
define double @Julian(i32 %0, i32 %1, double %2) {
%4 = alloca i32
%5 = alloca i32
%6 = alloca double
%7 = alloca i32
%8 = alloca i32
%9 = alloca i32
%10 = alloca i32
store i32 %0, i32* %4
store i32 %1, i32* %5
store double %2, double* %6
%11 = load i32, i32* %5
%12 = icmp slt i32 %11, 3
br i1 %12, label %13, label %18
13: ; preds = %3
%14 = load i32, i32* %4
%15 = sub nsw i32 %14, 1
store i32 %15, i32* %4
%16 = load i32, i32* %5
%17 = add nsw i32 %16, 12
store i32 %17, i32* %5
br label %18
18: ; preds = %13, %3
%19 = load i32, i32* %4
%20 = icmp sgt i32 %19, 1582
br i1 %20, label %36, label %21
21: ; preds = %18
%22 = load i32, i32* %4
%23 = icmp eq i32 %22, 1582
br i1 %23, label %24, label %27
24: ; preds = %21
%25 = load i32, i32* %5
%26 = icmp sgt i32 %25, 10
br i1 %26, label %36, label %27
27: ; preds = %24, %21
%28 = load i32, i32* %4
%29 = icmp eq i32 %28, 1582
br i1 %29, label %30, label %44
30: ; preds = %27
%31 = load i32, i32* %5
%32 = icmp eq i32 %31, 10
br i1 %32, label %33, label %44
33: ; preds = %30
%34 = load double, double* %6
%35 = fcmp ogt double %34, 1.500000e+01
br i1 %35, label %36, label %44
36: ; preds = %33, %24, %18
%37 = load i32, i32* %4
%38 = sdiv i32 %37, 100
store i32 %38, i32* %7
%39 = load i32, i32* %7
%40 = sub nsw i32 2, %39
%41 = load i32, i32* %7
%42 = sdiv i32 %41, 4
%43 = add nsw i32 %40, %42
store i32 %43, i32* %8
br label %44
44: ; preds = %36, %33, %30, %27
%45 = load i32, i32* %4
%46 = sitofp i32 %45 to double
%47 = fmul double 3.652500e+02, %46
%48 = fptosi double %47 to i64
%49 = trunc i64 %48 to i32
store i32 %49, i32* %9
%50 = load i32, i32* %5
%51 = add nsw i32 %50, 1
%52 = sitofp i32 %51 to double
%53 = fmul double 3.060010e+01, %52
%54 = fptosi double %53 to i64
%55 = trunc i64 %54 to i32
store i32 %55, i32* %10
%56 = load i32, i32* %8
%57 = load i32, i32* %9
%58 = add nsw i32 %56, %57
%59 = load i32, i32* %10
%60 = add nsw i32 %58, %59
%61 = sitofp i32 %60 to double
%62 = load double, double* %6
%63 = fadd double %61, %62
%64 = fadd double %63, 0x413A42A280000000
ret double %64
}

LLVM IR (-O2):
define double @Julian(i32 %0, i32 %1, double %2) {
%4 = icmp slt i32 %1, 3
%5 = add nsw i32 %1, 12
%6 = select i1 %4, i32 %5, i32 %1
%7 = sext i1 %4 to i32
%8 = add nsw i32 %7, %0
%9 = icmp sgt i32 %8, 1582
br i1 %9, label %19, label %10
10: ; preds = %3
%11 = icmp eq i32 %8, 1582
%12 = icmp sgt i32 %6, 10
%13 = select i1 %11, i1 %12, i1 false
br i1 %13, label %19, label %14
14: ; preds = %10
%15 = icmp eq i32 %6, 10
%16 = select i1 %11, i1 %15, i1 false
%17 = fcmp ogt double %2, 1.500000e+01
%18 = and i1 %17, %16
br i1 %18, label %19, label %24
19: ; preds = %14, %10, %3
%20 = udiv i32 %8, 100
%21 = sub nsw i32 2, %20
%22 = udiv i32 %8, 400
%23 = add nsw i32 %21, %22
br label %24
24: ; preds = %19, %14
%25 = phi i32 [ %23, %19 ], [ undef, %14 ]
%26 = sitofp i32 %8 to double
%27 = fmul double %26, 3.652500e+02
%28 = fptosi double %27 to i64
%29 = trunc i64 %28 to i32
%30 = add nsw i32 %6, 1
%31 = sitofp i32 %30 to double
%32 = fmul double %31, 3.060010e+01
%33 = fptosi double %32 to i64
%34 = trunc i64 %33 to i32
%35 = add i32 %34, %29
%36 = add i32 %35, %25
%37 = sitofp i32 %36 to double
%38 = fadd double %37, %2
%39 = fadd double %38, 0x413A42A280000000
ret double %39
}