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

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
}