unsigned long gregdate_to_julday(year, month, day)
int year;
int month;
int day;
{
  long igreg = 15 + 31 * (10 + (12 * 1582));
  long iy;
  long ja;
  long jm;
  long jy;
  unsigned long julday;
  if (year == 0)
    year = 1;
  iy = year;
  if (year < 0)
    iy++;
  if (month > 2) {
    jy = iy;
    jm = month + 1;
  } else {
    jy = iy - 1;
    jm = month + 13;
  }
  julday = day + (int)(30.6001 * jm);
  if (jy >= 0) {
    julday += 365 * jy;
    julday += 0.25 * jy;
  } else {
    double xi = 365.25 * jy;
    if ((int)xi != xi)
      xi -= 1;
    julday += (int)xi;
  }
  julday += 1720995;
  if (day + (31 * (month + (12 * iy))) >= igreg) {
    ja = jy / 100;
    julday -= ja;
    julday += 2;
    julday += ja / 4;
  }
  return julday;
}


x86 assembly (clang, -O0):
<gregdate_to_julday>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %edx,-0xc(%rbp)
movq   $0x8fc1d,-0x18(%rbp)
cmpl   $0x0,-0x4(%rbp)
jne    1126 <gregdate_to_julday+0x26>
movl   $0x1,-0x4(%rbp)
movslq -0x4(%rbp),%rax
mov    %rax,-0x20(%rbp)
cmpl   $0x0,-0x4(%rbp)
jge    1144 <gregdate_to_julday+0x44>
mov    -0x20(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x20(%rbp)
cmpl   $0x2,-0x8(%rbp)
jle    1167 <gregdate_to_julday+0x67>
mov    -0x20(%rbp),%rax
mov    %rax,-0x38(%rbp)
mov    -0x8(%rbp),%eax
add    $0x1,%eax
cltq
mov    %rax,-0x30(%rbp)
jmp    117f <gregdate_to_julday+0x7f>
mov    -0x20(%rbp),%rax
sub    $0x1,%rax
mov    %rax,-0x38(%rbp)
mov    -0x8(%rbp),%eax
add    $0xd,%eax
cltq
mov    %rax,-0x30(%rbp)
mov    -0xc(%rbp),%eax
cvtsi2sdq -0x30(%rbp),%xmm1
movsd  0xe70(%rip),%xmm0
mulsd  %xmm1,%xmm0
cvttsd2si %xmm0,%ecx
add    %ecx,%eax
cltq
mov    %rax,-0x40(%rbp)
cmpq   $0x0,-0x38(%rbp)
jl     122d <gregdate_to_julday+0x12d>
mov    -0x38(%rbp),%rax
imul   $0x16d,%rax,%rcx
mov    -0x40(%rbp),%rax
add    %rcx,%rax
mov    %rax,-0x40(%rbp)
cvtsi2sdq -0x38(%rbp),%xmm0
movaps 0xe62(%rip),%xmm2
movq   -0x40(%rbp),%xmm1
punpckldq %xmm2,%xmm1
movapd 0xe61(%rip),%xmm2
subpd  %xmm2,%xmm1
movaps %xmm1,%xmm2
unpckhpd %xmm1,%xmm1
addsd  %xmm2,%xmm1
movsd  0xe22(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movaps %xmm0,%xmm1
cvttsd2si %xmm1,%rax
mov    %rax,%rdx
sar    $0x3f,%rdx
movsd  0xe0b(%rip),%xmm1
subsd  %xmm1,%xmm0
cvttsd2si %xmm0,%rcx
and    %rdx,%rcx
or     %rcx,%rax
mov    %rax,-0x40(%rbp)
jmp    1288 <gregdate_to_julday+0x188>
cvtsi2sdq -0x38(%rbp),%xmm1
movsd  0xdcd(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
cvttsd2si -0x48(%rbp),%eax
cvtsi2sd %eax,%xmm0
ucomisd -0x48(%rbp),%xmm0
jne    1263 <gregdate_to_julday+0x163>
jp     1263 <gregdate_to_julday+0x163>
jmp    1279 <gregdate_to_julday+0x179>
movsd  -0x48(%rbp),%xmm0
movsd  0xda0(%rip),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
cvttsd2si -0x48(%rbp),%eax
cltq
add    -0x40(%rbp),%rax
mov    %rax,-0x40(%rbp)
mov    -0x40(%rbp),%rax
add    $0x1a42a3,%rax
mov    %rax,-0x40(%rbp)
movslq -0xc(%rbp),%rax
movslq -0x8(%rbp),%rcx
imul   $0xc,-0x20(%rbp),%rdx
add    %rdx,%rcx
imul   $0x1f,%rcx,%rcx
add    %rcx,%rax
cmp    -0x18(%rbp),%rax
jl     12fa <gregdate_to_julday+0x1fa>
mov    -0x38(%rbp),%rax
mov    $0x64,%ecx
cqto
idiv   %rcx
mov    %rax,-0x28(%rbp)
mov    -0x28(%rbp),%rcx
mov    -0x40(%rbp),%rax
sub    %rcx,%rax
mov    %rax,-0x40(%rbp)
mov    -0x40(%rbp),%rax
add    $0x2,%rax
mov    %rax,-0x40(%rbp)
mov    -0x28(%rbp),%rax
mov    $0x4,%ecx
cqto
idiv   %rcx
add    -0x40(%rbp),%rax
mov    %rax,-0x40(%rbp)
mov    -0x40(%rbp),%rax
pop    %rbp
ret

x86 assembly (clang, -O2):
<gregdate_to_julday>:
test   %edi,%edi
mov    $0x1,%eax
cmovne %edi,%eax
movslq %eax,%rcx
mov    %ecx,%eax
shr    $0x1f,%eax
add    %rcx,%rax
lea    0xd(%rsi),%r8d
lea    0x1(%rsi),%edi
xor    %ecx,%ecx
cmp    $0x3,%esi
cmovl  %r8d,%edi
cvtsi2sd %edi,%xmm0
mulsd  0xecf(%rip),%xmm0
cvttsd2si %xmm0,%edi
setl   %cl
add    %edx,%edi
mov    %rax,%r9
sub    %rcx,%r9
movslq %edi,%r8
js     11a5 <gregdate_to_julday+0xa5>
imul   $0x16d,%r9,%rcx
add    %r8,%rcx
xorps  %xmm0,%xmm0
cvtsi2sd %r9,%xmm0
movq   %rcx,%xmm1
punpckldq 0xecc(%rip),%xmm1
subpd  0xed4(%rip),%xmm1
movapd %xmm1,%xmm2
unpckhpd %xmm1,%xmm2
mulsd  0xe9c(%rip),%xmm0
addsd  %xmm1,%xmm2
addsd  %xmm2,%xmm0
cvttsd2si %xmm0,%r8
mov    %r8,%rcx
sar    $0x3f,%rcx
subsd  0xe88(%rip),%xmm0
cvttsd2si %xmm0,%rdi
and    %rcx,%rdi
or     %r8,%rdi
jmp    11e4 <gregdate_to_julday+0xe4>
xorps  %xmm0,%xmm0
cvtsi2sd %r9,%xmm0
mulsd  0xe53(%rip),%xmm0
cvttpd2dq %xmm0,%xmm1
cvtdq2pd %xmm1,%xmm1
movsd  0xe4b(%rip),%xmm2
addsd  %xmm0,%xmm2
cmpneqsd %xmm0,%xmm1
andpd  %xmm1,%xmm2
andnpd %xmm0,%xmm1
orpd   %xmm2,%xmm1
cvttsd2si %xmm1,%ecx
movslq %ecx,%rdi
add    %r8,%rdi
add    $0x1a42a3,%rdi
movslq %edx,%rcx
movslq %esi,%rdx
lea    (%rax,%rax,2),%rax
lea    (%rdx,%rax,4),%rax
mov    %rax,%rdx
shl    $0x5,%rdx
sub    %rax,%rdx
add    %rcx,%rdx
cmp    $0x8fc1d,%rdx
jl     125e <gregdate_to_julday+0x15e>
movabs $0x5c28f5c28f5c28f5,%rcx
f5 28 5c
mov    %r9,%rax
imul   %rcx
mov    %rdx,%rsi
sub    %r9,%rsi
mov    %rsi,%rax
shr    $0x3f,%rax
sar    $0x6,%rsi
add    %rax,%rsi
movabs $0xa3d70a3d70a3d70b,%rcx
0a d7 a3
mov    %r9,%rax
imul   %rcx
add    %r9,%rdx
mov    %rdx,%rax
shr    $0x3f,%rax
sar    $0x8,%rdx
add    %rax,%rdx
add    %rsi,%rdx
add    %rdx,%rdi
add    $0x2,%rdi
mov    %rdi,%rax
ret

LLVM IR (-O0):
define i64 @gregdate_to_julday(i32 %0, i32 %1, i32 %2) {
%4 = alloca i32
%5 = alloca i32
%6 = alloca i32
%7 = alloca i64
%8 = alloca i64
%9 = alloca i64
%10 = alloca i64
%11 = alloca i64
%12 = alloca i64
%13 = alloca double
store i32 %0, i32* %4
store i32 %1, i32* %5
store i32 %2, i32* %6
store i64 588829, i64* %7
%14 = load i32, i32* %4
%15 = icmp eq i32 %14, 0
br i1 %15, label %16, label %17
16: ; preds = %3
store i32 1, i32* %4
br label %17
17: ; preds = %16, %3
%18 = load i32, i32* %4
%19 = sext i32 %18 to i64
store i64 %19, i64* %8
%20 = load i32, i32* %4
%21 = icmp slt i32 %20, 0
br i1 %21, label %22, label %25
22: ; preds = %17
%23 = load i64, i64* %8
%24 = add nsw i64 %23, 1
store i64 %24, i64* %8
br label %25
25: ; preds = %22, %17
%26 = load i32, i32* %5
%27 = icmp sgt i32 %26, 2
br i1 %27, label %28, label %33
28: ; preds = %25
%29 = load i64, i64* %8
store i64 %29, i64* %11
%30 = load i32, i32* %5
%31 = add nsw i32 %30, 1
%32 = sext i32 %31 to i64
store i64 %32, i64* %10
br label %39
33: ; preds = %25
%34 = load i64, i64* %8
%35 = sub nsw i64 %34, 1
store i64 %35, i64* %11
%36 = load i32, i32* %5
%37 = add nsw i32 %36, 13
%38 = sext i32 %37 to i64
store i64 %38, i64* %10
br label %39
39: ; preds = %33, %28
%40 = load i32, i32* %6
%41 = load i64, i64* %10
%42 = sitofp i64 %41 to double
%43 = fmul double 3.060010e+01, %42
%44 = fptosi double %43 to i32
%45 = add nsw i32 %40, %44
%46 = sext i32 %45 to i64
store i64 %46, i64* %12
%47 = load i64, i64* %11
%48 = icmp sge i64 %47, 0
br i1 %48, label %49, label %60
49: ; preds = %39
%50 = load i64, i64* %11
%51 = mul nsw i64 365, %50
%52 = load i64, i64* %12
%53 = add i64 %52, %51
store i64 %53, i64* %12
%54 = load i64, i64* %11
%55 = sitofp i64 %54 to double
%56 = load i64, i64* %12
%57 = uitofp i64 %56 to double
%58 = call double @llvm.fmuladd.f64(double 2.500000e-01, double %55, double %57)
%59 = fptoui double %58 to i64
store i64 %59, i64* %12
br label %78
60: ; preds = %39
%61 = load i64, i64* %11
%62 = sitofp i64 %61 to double
%63 = fmul double 3.652500e+02, %62
store double %63, double* %13
%64 = load double, double* %13
%65 = fptosi double %64 to i32
%66 = sitofp i32 %65 to double
%67 = load double, double* %13
%68 = fcmp une double %66, %67
br i1 %68, label %69, label %72
69: ; preds = %60
%70 = load double, double* %13
%71 = fsub double %70, 1.000000e+00
store double %71, double* %13
br label %72
72: ; preds = %69, %60
%73 = load double, double* %13
%74 = fptosi double %73 to i32
%75 = sext i32 %74 to i64
%76 = load i64, i64* %12
%77 = add i64 %76, %75
store i64 %77, i64* %12
br label %78
78: ; preds = %72, %49
%79 = load i64, i64* %12
%80 = add i64 %79, 1720995
store i64 %80, i64* %12
%81 = load i32, i32* %6
%82 = sext i32 %81 to i64
%83 = load i32, i32* %5
%84 = sext i32 %83 to i64
%85 = load i64, i64* %8
%86 = mul nsw i64 12, %85
%87 = add nsw i64 %84, %86
%88 = mul nsw i64 31, %87
%89 = add nsw i64 %82, %88
%90 = load i64, i64* %7
%91 = icmp sge i64 %89, %90
br i1 %91, label %92, label %104
92: ; preds = %78
%93 = load i64, i64* %11
%94 = sdiv i64 %93, 100
store i64 %94, i64* %9
%95 = load i64, i64* %9
%96 = load i64, i64* %12
%97 = sub i64 %96, %95
store i64 %97, i64* %12
%98 = load i64, i64* %12
%99 = add i64 %98, 2
store i64 %99, i64* %12
%100 = load i64, i64* %9
%101 = sdiv i64 %100, 4
%102 = load i64, i64* %12
%103 = add i64 %102, %101
store i64 %103, i64* %12
br label %104
104: ; preds = %92, %78
%105 = load i64, i64* %12
ret i64 %105
}

LLVM IR (-O2):
define i64 @gregdate_to_julday(i32 %0, i32 %1, i32 %2) {
%4 = icmp eq i32 %0, 0
%5 = select i1 %4, i32 1, i32 %0
%6 = sext i32 %5 to i64
%7 = lshr i32 %5, 31
%8 = zext i32 %7 to i64
%9 = add nsw i64 %8, %6
%10 = icmp sgt i32 %1, 2
%11 = select i1 %10, i32 1, i32 13
%12 = add nsw i32 %11, %1
%13 = xor i1 %10, true
%14 = sext i1 %13 to i64
%15 = add nsw i64 %9, %14
%16 = sitofp i32 %12 to double
%17 = fmul double %16, 3.060010e+01
%18 = fptosi double %17 to i32
%19 = add nsw i32 %18, %2
%20 = sext i32 %19 to i64
%21 = icmp sgt i64 %15, -1
br i1 %21, label %22, label %29
22: ; preds = %3
%23 = mul nuw nsw i64 %15, 365
%24 = add nsw i64 %23, %20
%25 = sitofp i64 %15 to double
%26 = uitofp i64 %24 to double
%27 = tail call double @llvm.fmuladd.f64(double %25, double 2.500000e-01, double %26)
%28 = fptoui double %27 to i64
br label %40
29: ; preds = %3
%30 = sitofp i64 %15 to double
%31 = fmul double %30, 3.652500e+02
%32 = fptosi double %31 to i32
%33 = sitofp i32 %32 to double
%34 = fcmp une double %31, %33
%35 = fadd double %31, -1.000000e+00
%36 = select i1 %34, double %35, double %31
%37 = fptosi double %36 to i32
%38 = sext i32 %37 to i64
%39 = add nsw i64 %38, %20
br label %40
40: ; preds = %29, %22
%41 = phi i64 [ %28, %22 ], [ %39, %29 ]
%42 = add i64 %41, 1720995
%43 = sext i32 %2 to i64
%44 = sext i32 %1 to i64
%45 = mul nsw i64 %9, 12
%46 = add nsw i64 %45, %44
%47 = mul nsw i64 %46, 31
%48 = add nsw i64 %47, %43
%49 = icmp sgt i64 %48, 588828
br i1 %49, label %50, label %56
50: ; preds = %40
%51 = sdiv i64 %15, -100
%52 = sdiv i64 %15, 400
%53 = add nsw i64 %51, 2
%54 = add nsw i64 %53, %52
%55 = add i64 %54, %42
br label %56
56: ; preds = %50, %40
%57 = phi i64 [ %55, %50 ], [ %42, %40 ]
ret i64 %57
}