unsigned long YMDtoJulian(unsigned int year, unsigned int month,
                          unsigned int day) {
  short a, b = 0;
  short work_year = year;
  short work_month = month;
  short work_day = day;
  if (work_year < 0) {
    work_year++;
  }
  if (work_month <= 2) {
    work_year--;
    work_month += (short)12;
  }
  if (work_year * 10000. + work_month * 100. + work_day >= 15821015.) {
    a = (short)(work_year / 100.);
    b = (short)(2 - a + a / 4);
  }
  return (long)(365.25 * work_year) + (long)(30.6001 * (work_month + 1)) +
         work_day + 1720994L + b;
}


x86 assembly (clang, -O0):
<YMDtoJulian>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %edx,-0xc(%rbp)
movw   $0x0,-0x10(%rbp)
mov    -0x4(%rbp),%eax
mov    %ax,-0x12(%rbp)
mov    -0x8(%rbp),%eax
mov    %ax,-0x14(%rbp)
mov    -0xc(%rbp),%eax
mov    %ax,-0x16(%rbp)
movswl -0x12(%rbp),%eax
cmp    $0x0,%eax
jge    1141 <YMDtoJulian+0x41>
mov    -0x12(%rbp),%ax
add    $0x1,%ax
mov    %ax,-0x12(%rbp)
movswl -0x14(%rbp),%eax
cmp    $0x2,%eax
jg     1165 <YMDtoJulian+0x65>
mov    -0x12(%rbp),%ax
add    $0xffff,%ax
mov    %ax,-0x12(%rbp)
movswl -0x14(%rbp),%eax
add    $0xc,%eax
mov    %ax,-0x14(%rbp)
movswl -0x12(%rbp),%eax
cvtsi2sd %eax,%xmm0
movswl -0x14(%rbp),%eax
cvtsi2sd %eax,%xmm1
movsd  0xe93(%rip),%xmm2
mulsd  %xmm2,%xmm1
movsd  0xe7f(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movswl -0x16(%rbp),%eax
cvtsi2sd %eax,%xmm1
addsd  %xmm1,%xmm0
movsd  0xe5b(%rip),%xmm1
ucomisd %xmm1,%xmm0
jb     11f0 <YMDtoJulian+0xf0>
movswl -0x12(%rbp),%eax
cvtsi2sd %eax,%xmm0
movsd  0xe51(%rip),%xmm1
divsd  %xmm1,%xmm0
cvttsd2si %xmm0,%eax
mov    %ax,-0xe(%rbp)
movswl -0xe(%rbp),%ecx
mov    $0x2,%eax
sub    %ecx,%eax
mov    %eax,-0x1c(%rbp)
movswl -0xe(%rbp),%eax
mov    $0x4,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x1c(%rbp),%eax
add    %ecx,%eax
mov    %ax,-0x10(%rbp)
movswl -0x12(%rbp),%eax
cvtsi2sd %eax,%xmm1
movsd  0xe20(%rip),%xmm0
mulsd  %xmm1,%xmm0
cvttsd2si %xmm0,%rax
movswl -0x14(%rbp),%ecx
add    $0x1,%ecx
cvtsi2sd %ecx,%xmm1
movsd  0xdfc(%rip),%xmm0
mulsd  %xmm1,%xmm0
cvttsd2si %xmm0,%rcx
add    %rcx,%rax
movswq -0x16(%rbp),%rcx
add    %rcx,%rax
add    $0x1a42a2,%rax
movswq -0x10(%rbp),%rcx
add    %rcx,%rax
pop    %rbp
ret

x86 assembly (clang, -O2):
<YMDtoJulian>:
movzwl %di,%ecx
shr    $0xf,%ecx
add    %edi,%ecx
mov    %esi,%edi
shl    $0x10,%edi
lea    0xc0000(%rdi),%r8d
xor    %eax,%eax
cmp    $0x3,%si
setl   %al
cmovge %edi,%r8d
sub    %eax,%ecx
movswl %cx,%eax
cvtsi2sd %eax,%xmm0
sar    $0x10,%r8d
cvtsi2sd %r8d,%xmm1
mulsd  0xec5(%rip),%xmm1
movsd  0xec5(%rip),%xmm2
mulsd  %xmm0,%xmm2
addsd  %xmm1,%xmm2
movswl %dx,%eax
xorps  %xmm1,%xmm1
cvtsi2sd %eax,%xmm1
addsd  %xmm2,%xmm1
ucomisd 0xeaf(%rip),%xmm1
jae    1167 <YMDtoJulian+0x67>
xor    %ecx,%ecx
jmp    118f <YMDtoJulian+0x8f>
movapd %xmm0,%xmm1
divsd  0xe8d(%rip),%xmm1
cvttsd2si %xmm1,%eax
lea    0x3(%rax),%ecx
test   %ax,%ax
cmovns %eax,%ecx
movswl %cx,%ecx
sar    $0x2,%ecx
sub    %eax,%ecx
add    $0x2,%ecx
movswq %cx,%rcx
mulsd  0xe81(%rip),%xmm0
cvttsd2si %xmm0,%rax
add    $0x1,%r8d
xorps  %xmm0,%xmm0
cvtsi2sd %r8d,%xmm0
mulsd  0xe70(%rip),%xmm0
cvttsd2si %xmm0,%rsi
movswq %dx,%rdx
add    %rax,%rdx
add    %rsi,%rdx
lea    (%rcx,%rdx,1),%rax
add    $0x1a42a2,%rax
ret

LLVM IR (-O0):
define i64 @YMDtoJulian(i32 %0, i32 %1, i32 %2) {
%4 = alloca i32
%5 = alloca i32
%6 = alloca i32
%7 = alloca i16
%8 = alloca i16
%9 = alloca i16
%10 = alloca i16
%11 = alloca i16
store i32 %0, i32* %4
store i32 %1, i32* %5
store i32 %2, i32* %6
store i16 0, i16* %8
%12 = load i32, i32* %4
%13 = trunc i32 %12 to i16
store i16 %13, i16* %9
%14 = load i32, i32* %5
%15 = trunc i32 %14 to i16
store i16 %15, i16* %10
%16 = load i32, i32* %6
%17 = trunc i32 %16 to i16
store i16 %17, i16* %11
%18 = load i16, i16* %9
%19 = sext i16 %18 to i32
%20 = icmp slt i32 %19, 0
br i1 %20, label %21, label %24
21: ; preds = %3
%22 = load i16, i16* %9
%23 = add i16 %22, 1
store i16 %23, i16* %9
br label %24
24: ; preds = %21, %3
%25 = load i16, i16* %10
%26 = sext i16 %25 to i32
%27 = icmp sle i32 %26, 2
br i1 %27, label %28, label %35
28: ; preds = %24
%29 = load i16, i16* %9
%30 = add i16 %29, -1
store i16 %30, i16* %9
%31 = load i16, i16* %10
%32 = sext i16 %31 to i32
%33 = add nsw i32 %32, 12
%34 = trunc i32 %33 to i16
store i16 %34, i16* %10
br label %35
35: ; preds = %28, %24
%36 = load i16, i16* %9
%37 = sext i16 %36 to i32
%38 = sitofp i32 %37 to double
%39 = load i16, i16* %10
%40 = sext i16 %39 to i32
%41 = sitofp i32 %40 to double
%42 = fmul double %41, 1.000000e+02
%43 = call double @llvm.fmuladd.f64(double %38, double 1.000000e+04, double %42)
%44 = load i16, i16* %11
%45 = sext i16 %44 to i32
%46 = sitofp i32 %45 to double
%47 = fadd double %43, %46
%48 = fcmp oge double %47, 0x416E2D1AE0000000
br i1 %48, label %49, label %63
49: ; preds = %35
%50 = load i16, i16* %9
%51 = sext i16 %50 to i32
%52 = sitofp i32 %51 to double
%53 = fdiv double %52, 1.000000e+02
%54 = fptosi double %53 to i16
store i16 %54, i16* %7
%55 = load i16, i16* %7
%56 = sext i16 %55 to i32
%57 = sub nsw i32 2, %56
%58 = load i16, i16* %7
%59 = sext i16 %58 to i32
%60 = sdiv i32 %59, 4
%61 = add nsw i32 %57, %60
%62 = trunc i32 %61 to i16
store i16 %62, i16* %8
br label %63
63: ; preds = %49, %35
%64 = load i16, i16* %9
%65 = sext i16 %64 to i32
%66 = sitofp i32 %65 to double
%67 = fmul double 3.652500e+02, %66
%68 = fptosi double %67 to i64
%69 = load i16, i16* %10
%70 = sext i16 %69 to i32
%71 = add nsw i32 %70, 1
%72 = sitofp i32 %71 to double
%73 = fmul double 3.060010e+01, %72
%74 = fptosi double %73 to i64
%75 = add nsw i64 %68, %74
%76 = load i16, i16* %11
%77 = sext i16 %76 to i64
%78 = add nsw i64 %75, %77
%79 = add nsw i64 %78, 1720994
%80 = load i16, i16* %8
%81 = sext i16 %80 to i64
%82 = add nsw i64 %79, %81
ret i64 %82
}

LLVM IR (-O2):
define i64 @YMDtoJulian(i32 %0, i32 %1, i32 %2) {
%4 = trunc i32 %0 to i16
%5 = lshr i16 %4, 15
%6 = add i16 %5, %4
%7 = trunc i32 %1 to i16
%8 = icmp slt i16 %7, 3
%9 = sext i1 %8 to i16
%10 = add i16 %6, %9
%11 = sitofp i16 %10 to double
%12 = shl i32 %1, 16
%13 = add i32 %12, 786432
%14 = select i1 %8, i32 %13, i32 %12
%15 = ashr exact i32 %14, 16
%16 = sitofp i32 %15 to double
%17 = fmul double %16, 1.000000e+02
%18 = tail call double @llvm.fmuladd.f64(double %11, double 1.000000e+04, double %17)
%19 = shl i32 %2, 16
%20 = ashr exact i32 %19, 16
%21 = sitofp i32 %20 to double
%22 = fadd double %18, %21
%23 = fcmp ult double %22, 0x416E2D1AE0000000
br i1 %23, label %34, label %24
24: ; preds = %3
%25 = fdiv double %11, 1.000000e+02
%26 = fptosi double %25 to i16
%27 = zext i16 %26 to i64
%28 = sub nsw i64 2, %27
%29 = sdiv i16 %26, 4
%30 = zext i16 %29 to i64
%31 = add nsw i64 %28, %30
%32 = shl i64 %31, 48
%33 = ashr exact i64 %32, 48
br label %34
34: ; preds = %24, %3
%35 = phi i64 [ %33, %24 ], [ 0, %3 ]
%36 = zext i32 %2 to i64
%37 = fmul double %11, 3.652500e+02
%38 = fptosi double %37 to i64
%39 = add nsw i32 %15, 1
%40 = sitofp i32 %39 to double
%41 = fmul double %40, 3.060010e+01
%42 = fptosi double %41 to i64
%43 = shl i64 %36, 48
%44 = ashr exact i64 %43, 48
%45 = add nsw i64 %44, 1720994
%46 = add i64 %45, %38
%47 = add i64 %46, %42
%48 = add i64 %47, %35
ret i64 %48
}