double jd(ny, nm, nd, UT)
int ny, nm, nd;
double UT;
{
  double A, B, C, D, JD, day;
  day = nd + UT / 24.0;
  if ((nm == 1) || (nm == 2)) {
    ny = ny - 1;
    nm = nm + 12;
  }
  if (((double)ny + nm / 12.0 + day / 365.25) >=
      (1582.0 + 10.0 / 12.0 + 15.0 / 365.25)) {
    A = ((int)(ny / 100.0));
    B = 2.0 - A + (int)(A / 4.0);
  } else {
    B = 0.0;
  }
  if (ny < 0.0) {
    C = (int)((365.25 * (double)ny) - 0.75);
  } else {
    C = (int)(365.25 * (double)ny);
  }
  D = (int)(30.6001 * (double)(nm + 1));
  JD = B + C + D + day + 1720994.5;
  return (JD);
}


x86 assembly (clang, -O0):
<jd>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %edx,-0xc(%rbp)
movsd  %xmm0,-0x18(%rbp)
cvtsi2sdl -0xc(%rbp),%xmm0
movsd  -0x18(%rbp),%xmm1
movsd  0xedc(%rip),%xmm2
divsd  %xmm2,%xmm1
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
cmpl   $0x1,-0x8(%rbp)
je     1145 <jd+0x45>
cmpl   $0x2,-0x8(%rbp)
jne    1157 <jd+0x57>
mov    -0x4(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x4(%rbp)
mov    -0x8(%rbp),%eax
add    $0xc,%eax
mov    %eax,-0x8(%rbp)
cvtsi2sdl -0x4(%rbp),%xmm0
cvtsi2sdl -0x8(%rbp),%xmm1
movsd  0xeaf(%rip),%xmm2
divsd  %xmm2,%xmm1
addsd  %xmm1,%xmm0
movsd  -0x48(%rbp),%xmm1
movsd  0xe92(%rip),%xmm2
divsd  %xmm2,%xmm1
addsd  %xmm1,%xmm0
movsd  0xe7a(%rip),%xmm1
ucomisd %xmm1,%xmm0
jb     11ea <jd+0xea>
cvtsi2sdl -0x4(%rbp),%xmm0
movsd  0xe8b(%rip),%xmm1
divsd  %xmm1,%xmm0
cvttsd2si %xmm0,%eax
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x20(%rbp)
movsd  0xe6a(%rip),%xmm0
subsd  -0x20(%rbp),%xmm0
movsd  -0x20(%rbp),%xmm1
movsd  0xe50(%rip),%xmm2
divsd  %xmm2,%xmm1
cvttsd2si %xmm1,%eax
cvtsi2sd %eax,%xmm1
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x28(%rbp)
jmp    11f2 <jd+0xf2>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x28(%rbp)
cvtsi2sdl -0x4(%rbp),%xmm1
xorps  %xmm0,%xmm0
ucomisd %xmm1,%xmm0
jbe    1233 <jd+0x133>
cvtsi2sdl -0x4(%rbp),%xmm0
movsd  0xdff(%rip),%xmm1
mulsd  %xmm1,%xmm0
movsd  0xe1b(%rip),%xmm1
addsd  %xmm1,%xmm0
cvttsd2si %xmm0,%eax
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x30(%rbp)
jmp    1251 <jd+0x151>
cvtsi2sdl -0x4(%rbp),%xmm1
movsd  0xdd0(%rip),%xmm0
mulsd  %xmm1,%xmm0
cvttsd2si %xmm0,%eax
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x8(%rbp),%eax
add    $0x1,%eax
cvtsi2sd %eax,%xmm1
movsd  0xde5(%rip),%xmm0
mulsd  %xmm1,%xmm0
cvttsd2si %xmm0,%eax
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x28(%rbp),%xmm0
addsd  -0x30(%rbp),%xmm0
addsd  -0x38(%rbp),%xmm0
addsd  -0x48(%rbp),%xmm0
movsd  0xdb0(%rip),%xmm1
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x40(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<jd>:
cvtsi2sd %edx,%xmm2
divsd  0xef4(%rip),%xmm0
lea    -0x1(%rsi),%ecx
lea    0xc(%rsi),%eax
cmp    $0x2,%ecx
cmovae %esi,%eax
sbb    $0x0,%edi
cvtsi2sd %edi,%xmm1
cvtsi2sd %eax,%xmm3
addsd  %xmm2,%xmm0
movapd %xmm0,%xmm2
unpcklpd %xmm3,%xmm2
divpd  0xf19(%rip),%xmm2
movapd %xmm2,%xmm3
unpckhpd %xmm2,%xmm3
addsd  %xmm1,%xmm3
addsd  %xmm2,%xmm3
xorpd  %xmm2,%xmm2
ucomisd 0xeb5(%rip),%xmm3
jb     1189 <jd+0x89>
movapd %xmm1,%xmm2
divsd  0xeaf(%rip),%xmm2
cvttpd2dq %xmm2,%xmm2
cvtdq2pd %xmm2,%xmm2
movsd  0xea7(%rip),%xmm3
subsd  %xmm2,%xmm3
mulsd  0xea3(%rip),%xmm2
cvttpd2dq %xmm2,%xmm2
cvtdq2pd %xmm2,%xmm2
addsd  %xmm3,%xmm2
mulsd  0xe97(%rip),%xmm1
test   %edi,%edi
jns    119d <jd+0x9d>
addsd  0xe93(%rip),%xmm1
cvttpd2dq %xmm1,%xmm1
cvtdq2pd %xmm1,%xmm1
add    $0x1,%eax
xorps  %xmm3,%xmm3
cvtsi2sd %eax,%xmm3
mulsd  0xe81(%rip),%xmm3
cvttpd2dq %xmm3,%xmm3
cvtdq2pd %xmm3,%xmm3
addsd  %xmm1,%xmm2
addsd  %xmm3,%xmm2
addsd  %xmm2,%xmm0
addsd  0xe6d(%rip),%xmm0
ret

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

LLVM IR (-O2):
define double @jd(i32 %0, i32 %1, i32 %2, double %3) {
%5 = sitofp i32 %2 to double
%6 = fdiv double %3, 2.400000e+01
%7 = fadd double %6, %5
%8 = add i32 %1, -1
%9 = icmp ult i32 %8, 2
%10 = add nsw i32 %1, 12
%11 = select i1 %9, i32 %10, i32 %1
%12 = sext i1 %9 to i32
%13 = add nsw i32 %12, %0
%14 = sitofp i32 %13 to double
%15 = sitofp i32 %11 to double
%16 = insertelement <2 x double> poison, double %7, i64 0
%17 = insertelement <2 x double> %16, double %15, i64 1
%18 = fdiv <2 x double> %17, <double 3.652500e+02, double 1.200000e+01>
%19 = extractelement <2 x double> %18, i64 1
%20 = fadd double %19, %14
%21 = extractelement <2 x double> %18, i64 0
%22 = fadd double %21, %20
%23 = fcmp ult double %22, 0x4098BB7F63002CDB
br i1 %23, label %33, label %24
24: ; preds = %4
%25 = fdiv double %14, 1.000000e+02
%26 = fptosi double %25 to i32
%27 = sitofp i32 %26 to double
%28 = fsub double 2.000000e+00, %27
%29 = fmul double %27, 2.500000e-01
%30 = fptosi double %29 to i32
%31 = sitofp i32 %30 to double
%32 = fadd double %28, %31
br label %33
33: ; preds = %4, %24
%34 = phi double [ %32, %24 ], [ 0.000000e+00, %4 ]
%35 = icmp slt i32 %13, 0
%36 = tail call double @llvm.fmuladd.f64(double %14, double 3.652500e+02, double -7.500000e-01)
%37 = fmul double %14, 3.652500e+02
%38 = select i1 %35, double %36, double %37
%39 = fptosi double %38 to i32
%40 = sitofp i32 %39 to double
%41 = add nsw i32 %11, 1
%42 = sitofp i32 %41 to double
%43 = fmul double %42, 3.060010e+01
%44 = fptosi double %43 to i32
%45 = sitofp i32 %44 to double
%46 = fadd double %34, %40
%47 = fadd double %46, %45
%48 = fadd double %7, %47
%49 = fadd double %48, 0x413A42A280000000
ret double %49
}