double julday(month, day, year)
int month;
int day;
int year;
{
  int days[12] = {0, 31, 59, 90, 120, 151, 181, 212, 243, 273, 304, 334};
  int leap_year =
      (((year % 4) == 0) && (((year % 100) != 0) || ((year % 400) == 0)));
  double day_in_year = (days[month - 1] + day + (month > 2 ? leap_year : 0));
  return ((double)year + (day_in_year / (365.0 + leap_year)));
}


x86 assembly (clang, -O0):
<julday>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %edx,-0xc(%rbp)
lea    -0x40(%rbp),%rdi
lea    0xee4(%rip),%rsi
mov    $0x30,%edx
call   1030 <memcpy@plt>
mov    -0xc(%rbp),%eax
mov    $0x4,%ecx
cltd
idiv   %ecx
xor    %eax,%eax
cmp    $0x0,%edx
mov    %al,-0x51(%rbp)
jne    1182 <julday+0x72>
mov    -0xc(%rbp),%eax
mov    $0x64,%ecx
cltd
idiv   %ecx
mov    $0x1,%al
cmp    $0x0,%edx
mov    %al,-0x52(%rbp)
jne    117c <julday+0x6c>
mov    -0xc(%rbp),%eax
mov    $0x190,%ecx
cltd
idiv   %ecx
cmp    $0x0,%edx
sete   %al
mov    %al,-0x52(%rbp)
mov    -0x52(%rbp),%al
mov    %al,-0x51(%rbp)
mov    -0x51(%rbp),%al
and    $0x1,%al
movzbl %al,%eax
mov    %eax,-0x44(%rbp)
mov    -0x4(%rbp),%eax
sub    $0x1,%eax
cltq
mov    -0x40(%rbp,%rax,4),%eax
add    -0x8(%rbp),%eax
mov    %eax,-0x58(%rbp)
cmpl   $0x2,-0x4(%rbp)
jle    11b4 <julday+0xa4>
mov    -0x44(%rbp),%eax
mov    %eax,-0x5c(%rbp)
jmp    11be <julday+0xae>
xor    %eax,%eax
mov    %eax,-0x5c(%rbp)
jmp    11be <julday+0xae>
mov    -0x58(%rbp),%eax
mov    -0x5c(%rbp),%ecx
add    %ecx,%eax
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x50(%rbp)
cvtsi2sdl -0xc(%rbp),%xmm0
movsd  -0x50(%rbp),%xmm1
cvtsi2sdl -0x44(%rbp),%xmm3
movsd  0xe1a(%rip),%xmm2
addsd  %xmm3,%xmm2
divsd  %xmm2,%xmm1
addsd  %xmm1,%xmm0
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<julday>:
xor    %r10d,%r10d
mov    $0x0,%ecx
test   $0x3,%dl
jne    1140 <julday+0x40>
imul   $0xc28f5c29,%edx,%r8d
add    $0x51eb850,%r8d
mov    %r8d,%r9d
ror    $0x2,%r9d
mov    $0x1,%ecx
cmp    $0x28f5c28,%r9d
ja     1140 <julday+0x40>
ror    $0x4,%r8d
xor    %ecx,%ecx
cmp    $0xa3d70b,%r8d
setb   %cl
movslq %edi,%rdi
lea    0xeb6(%rip),%rax
add    -0x4(%rax,%rdi,4),%esi
cmp    $0x3,%edi
cmovge %ecx,%r10d
add    %esi,%r10d
cvtsi2sd %r10d,%xmm0
cvtsi2sd %edx,%xmm1
add    $0x16d,%ecx
cvtsi2sd %ecx,%xmm2
divsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
ret

LLVM IR (-O0):
@__const.julday.days = [12 x i32] [i32 0, i32 31, i32 59, i32 90, i32 120, i32 151, i32 181, i32 212, i32 243, i32 273, i32 304, i32 334]
define double @julday(i32 %0, i32 %1, i32 %2) {
%4 = alloca i32
%5 = alloca i32
%6 = alloca i32
%7 = alloca [12 x i32]
%8 = alloca i32
%9 = alloca double
store i32 %0, i32* %4
store i32 %1, i32* %5
store i32 %2, i32* %6
%10 = bitcast [12 x i32]* %7 to i8*
%11 = load i32, i32* %6
%12 = srem i32 %11, 4
%13 = icmp eq i32 %12, 0
br i1 %13, label %14, label %24
14: ; preds = %3
%15 = load i32, i32* %6
%16 = srem i32 %15, 100
%17 = icmp ne i32 %16, 0
br i1 %17, label %22, label %18
18: ; preds = %14
%19 = load i32, i32* %6
%20 = srem i32 %19, 400
%21 = icmp eq i32 %20, 0
br label %22
22: ; preds = %18, %14
%23 = phi i1 [ true, %14 ], [ %21, %18 ]
br label %24
24: ; preds = %22, %3
%25 = phi i1 [ false, %3 ], [ %23, %22 ]
%26 = zext i1 %25 to i32
store i32 %26, i32* %8
%27 = load i32, i32* %4
%28 = sub nsw i32 %27, 1
%29 = sext i32 %28 to i64
%30 = getelementptr inbounds [12 x i32], [12 x i32]* %7, i64 0, i64 %29
%31 = load i32, i32* %30
%32 = load i32, i32* %5
%33 = add nsw i32 %31, %32
%34 = load i32, i32* %4
%35 = icmp sgt i32 %34, 2
br i1 %35, label %36, label %38
36: ; preds = %24
%37 = load i32, i32* %8
br label %39
38: ; preds = %24
br label %39
39: ; preds = %38, %36
%40 = phi i32 [ %37, %36 ], [ 0, %38 ]
%41 = add nsw i32 %33, %40
%42 = sitofp i32 %41 to double
store double %42, double* %9
%43 = load i32, i32* %6
%44 = sitofp i32 %43 to double
%45 = load double, double* %9
%46 = load i32, i32* %8
%47 = sitofp i32 %46 to double
%48 = fadd double 3.650000e+02, %47
%49 = fdiv double %45, %48
%50 = fadd double %44, %49
ret double %50
}

LLVM IR (-O2):
@__const.julday.days = [12 x i32] [i32 0, i32 31, i32 59, i32 90, i32 120, i32 151, i32 181, i32 212, i32 243, i32 273, i32 304, i32 334]
define double @julday(i32 %0, i32 %1, i32 %2) {
%4 = and i32 %2, 3
%5 = icmp eq i32 %4, 0
br i1 %5, label %6, label %13
6: ; preds = %3
%7 = srem i32 %2, 100
%8 = icmp eq i32 %7, 0
br i1 %8, label %9, label %13
9: ; preds = %6
%10 = srem i32 %2, 400
%11 = icmp eq i32 %10, 0
%12 = zext i1 %11 to i32
br label %13
13: ; preds = %6, %9, %3
%14 = phi i32 [ 0, %3 ], [ 1, %6 ], [ %12, %9 ]
%15 = add nsw i32 %0, -1
%16 = sext i32 %15 to i64
%17 = getelementptr inbounds [12 x i32], [12 x i32]* @__const.julday.days, i64 0, i64 %16
%18 = load i32, i32* %17
%19 = add nsw i32 %18, %1
%20 = icmp sgt i32 %0, 2
%21 = select i1 %20, i32 %14, i32 0
%22 = add nsw i32 %19, %21
%23 = sitofp i32 %22 to double
%24 = sitofp i32 %2 to double
%25 = add nuw nsw i32 %14, 365
%26 = sitofp i32 %25 to double
%27 = fdiv double %23, %26
%28 = fadd double %27, %24
ret double %28
}