double dmy2jd(double d, int m, int y) {
  int b;
  if (m <= 2) {
    y--;
    m += 12;
  }
  if (y < 1582 || (y == 1582 && (m < 10 || (m == 10 && d <= 4)))) {
    b = 0;
  } else {
    int a;
    a = y / 100;
    b = 2 - a + a / 4;
  }
  return (int)(365.25 * (y + 4716)) + (int)(30.6001 * (m + 1)) + d + b - 1524.5;
}


x86 assembly (clang, -O0):
<dmy2jd>:
push   %rbp
mov    %rsp,%rbp
movsd  %xmm0,-0x8(%rbp)
mov    %edi,-0xc(%rbp)
mov    %esi,-0x10(%rbp)
cmpl   $0x2,-0xc(%rbp)
jg     112b <dmy2jd+0x2b>
mov    -0x10(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x10(%rbp)
mov    -0xc(%rbp),%eax
add    $0xc,%eax
mov    %eax,-0xc(%rbp)
cmpl   $0x62e,-0x10(%rbp)
jl     116c <dmy2jd+0x6c>
cmpl   $0x62e,-0x10(%rbp)
jne    1178 <dmy2jd+0x78>
cmpl   $0xa,-0xc(%rbp)
jl     116c <dmy2jd+0x6c>
cmpl   $0xa,-0xc(%rbp)
jne    1178 <dmy2jd+0x78>
movsd  0xe9f(%rip),%xmm0
ucomisd -0x8(%rbp),%xmm0
jb     1178 <dmy2jd+0x78>
movl   $0x0,-0x14(%rbp)
jmp    11a6 <dmy2jd+0xa6>
mov    -0x10(%rbp),%eax
mov    $0x64,%ecx
cltd
idiv   %ecx
mov    %eax,-0x18(%rbp)
mov    $0x2,%eax
sub    -0x18(%rbp),%eax
mov    %eax,-0x1c(%rbp)
mov    -0x18(%rbp),%eax
mov    $0x4,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x1c(%rbp),%eax
add    %ecx,%eax
mov    %eax,-0x14(%rbp)
mov    -0x10(%rbp),%eax
add    $0x126c,%eax
cvtsi2sd %eax,%xmm1
movsd  0xe5e(%rip),%xmm0
mulsd  %xmm1,%xmm0
cvttsd2si %xmm0,%eax
mov    -0xc(%rbp),%ecx
add    $0x1,%ecx
cvtsi2sd %ecx,%xmm1
movsd  0xe3c(%rip),%xmm0
mulsd  %xmm1,%xmm0
cvttsd2si %xmm0,%ecx
add    %ecx,%eax
cvtsi2sd %eax,%xmm0
addsd  -0x8(%rbp),%xmm0
cvtsi2sdl -0x14(%rbp),%xmm1
addsd  %xmm1,%xmm0
movsd  0xe10(%rip),%xmm1
subsd  %xmm1,%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<dmy2jd>:
lea    0xc(%rdi),%eax
xor    %ecx,%ecx
cmp    $0x3,%edi
cmovge %edi,%eax
setl   %cl
sub    %ecx,%esi
xorpd  %xmm2,%xmm2
cmp    $0x62e,%esi
jl     1156 <dmy2jd+0x56>
jne    1136 <dmy2jd+0x36>
cmp    $0xa,%eax
jl     1156 <dmy2jd+0x56>
movsd  0xed5(%rip),%xmm1
ucomisd %xmm0,%xmm1
jb     1136 <dmy2jd+0x36>
cmp    $0xa,%eax
je     1156 <dmy2jd+0x56>
mov    %esi,%ecx
imul   $0x51eb851f,%rcx,%rcx
mov    %rcx,%rdx
shr    $0x25,%rdx
shr    $0x27,%rcx
sub    %edx,%ecx
add    $0x2,%ecx
xorps  %xmm2,%xmm2
cvtsi2sd %ecx,%xmm2
movd   %eax,%xmm1
movd   %esi,%xmm3
punpckldq %xmm1,%xmm3
paddd  0xea6(%rip),%xmm3
cvtdq2pd %xmm3,%xmm1
mulpd  0xeaa(%rip),%xmm1
cvttpd2dq %xmm1,%xmm1
pshufd $0x55,%xmm1,%xmm3
paddd  %xmm1,%xmm3
movd   %xmm3,%eax
xorps  %xmm1,%xmm1
cvtsi2sd %eax,%xmm1
addsd  %xmm0,%xmm1
addsd  %xmm2,%xmm1
addsd  0xe6a(%rip),%xmm1
movapd %xmm1,%xmm0
ret

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

LLVM IR (-O2):
define double @dmy2jd(double %0, i32 %1, i32 %2) {
%4 = icmp slt i32 %1, 3
%5 = add nsw i32 %1, 12
%6 = sext i1 %4 to i32
%7 = add nsw i32 %6, %2
%8 = select i1 %4, i32 %5, i32 %1
%9 = icmp slt i32 %7, 1582
br i1 %9, label %24, label %10
10: ; preds = %3
%11 = icmp eq i32 %7, 1582
br i1 %11, label %12, label %18
12: ; preds = %10
%13 = icmp slt i32 %8, 10
br i1 %13, label %24, label %14
14: ; preds = %12
%15 = icmp eq i32 %8, 10
%16 = fcmp ole double %0, 4.000000e+00
%17 = and i1 %16, %15
br i1 %17, label %24, label %18
18: ; preds = %14, %10
%19 = udiv i32 %7, 100
%20 = sub nsw i32 2, %19
%21 = udiv i32 %7, 400
%22 = add nsw i32 %20, %21
%23 = sitofp i32 %22 to double
br label %24
24: ; preds = %3, %12, %14, %18
%25 = phi double [ %23, %18 ], [ 0.000000e+00, %14 ], [ 0.000000e+00, %12 ], [ 0.000000e+00, %3 ]
%26 = insertelement <2 x i32> poison, i32 %7, i64 0
%27 = insertelement <2 x i32> %26, i32 %8, i64 1
%28 = add nsw <2 x i32> %27, <i32 4716, i32 1>
%29 = sitofp <2 x i32> %28 to <2 x double>
%30 = fmul <2 x double> %29, <double 3.652500e+02, double 3.060010e+01>
%31 = fptosi <2 x double> %30 to <2 x i32>
%32 = shufflevector <2 x i32> %31, <2 x i32> poison, <2 x i32> <i32 1, i32 undef>
%33 = add nsw <2 x i32> %32, %31
%34 = extractelement <2 x i32> %33, i64 0
%35 = sitofp i32 %34 to double
%36 = fadd double %35, %0
%37 = fadd double %36, %25
%38 = fadd double %37, -1.524500e+03
ret double %38
}