void jd_to_gregorian(double jd, int *year, int *month, int *day, int *hour,
                     int *min, double *sec) {
  long ijd = (long)(jd + 0.5);
  double frac = jd + 0.5 - (double)ijd + 1.0e-10;
  long ka = ijd;
  long ialp;
  if (ijd >= 2299161L) {
    ialp = ((double)ijd - 1867216.25) / (36524.25);
    ka = ijd + 1L + ialp - (ialp >> 2);
  }
  long kb = ka + 1524L;
  long kc = ((double)kb - 122.1) / 365.25;
  long kd = (double)kc * 365.25;
  long ke = (double)(kb - kd) / 30.6001;
  *day = kb - kd - ((long)((double)ke * 30.6001));
  if (ke > 13L)
    *month = ke - 13L;
  else
    *month = ke - 1L;
  if ((*month == 2) && (*day > 28))
    *day = 29;
  if ((*month == 2) && (*day == 29) && (ke == 3L))
    *year = kc - 4716L;
  else if (*month > 2)
    *year = kc - 4716L;
  else
    *year = kc - 4715L;
  double dh = frac * 24.0;
  *hour = dh;
  double dm = (dh - (double)*hour) * 60.0;
  *min = dm;
  *sec = (dm - (double)*min) * 60.0;
}


x86 assembly (clang, -O0):
<jd_to_gregorian>:
push   %rbp
mov    %rsp,%rbp
push   %rax
movsd  %xmm0,-0x8(%rbp)
mov    %rdi,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %rdx,-0x20(%rbp)
mov    %rcx,-0x28(%rbp)
mov    %r8,-0x30(%rbp)
mov    %r9,-0x38(%rbp)
movsd  0xede(%rip),%xmm0
addsd  -0x8(%rbp),%xmm0
cvttsd2si %xmm0,%rax
mov    %rax,-0x40(%rbp)
movsd  0xec8(%rip),%xmm0
addsd  -0x8(%rbp),%xmm0
cvtsi2sdq -0x40(%rbp),%xmm1
subsd  %xmm1,%xmm0
movsd  0xea9(%rip),%xmm1
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x40(%rbp),%rax
mov    %rax,-0x50(%rbp)
cmpq   $0x231519,-0x40(%rbp)
jl     11b8 <jd_to_gregorian+0xb8>
cvtsi2sdq -0x40(%rbp),%xmm0
movsd  0xe94(%rip),%xmm1
subsd  %xmm1,%xmm0
movsd  0xe80(%rip),%xmm1
divsd  %xmm1,%xmm0
cvttsd2si %xmm0,%rax
mov    %rax,-0x58(%rbp)
mov    -0x40(%rbp),%rax
add    $0x1,%rax
add    -0x58(%rbp),%rax
mov    -0x58(%rbp),%rcx
sar    $0x2,%rcx
sub    %rcx,%rax
mov    %rax,-0x50(%rbp)
mov    -0x50(%rbp),%rax
add    $0x5f4,%rax
mov    %rax,-0x60(%rbp)
cvtsi2sdq -0x60(%rbp),%xmm0
movsd  0xe5c(%rip),%xmm1
subsd  %xmm1,%xmm0
movsd  0xe48(%rip),%xmm1
divsd  %xmm1,%xmm0
cvttsd2si %xmm0,%rax
mov    %rax,-0x68(%rbp)
cvtsi2sdq -0x68(%rbp),%xmm0
movsd  0xe2d(%rip),%xmm1
mulsd  %xmm1,%xmm0
cvttsd2si %xmm0,%rax
mov    %rax,-0x70(%rbp)
mov    -0x60(%rbp),%rax
sub    -0x70(%rbp),%rax
cvtsi2sd %rax,%xmm0
movsd  0xe03(%rip),%xmm1
divsd  %xmm1,%xmm0
cvttsd2si %xmm0,%rax
mov    %rax,-0x78(%rbp)
mov    -0x60(%rbp),%rax
sub    -0x70(%rbp),%rax
cvtsi2sdq -0x78(%rbp),%xmm0
movsd  0xde0(%rip),%xmm1
mulsd  %xmm1,%xmm0
cvttsd2si %xmm0,%rcx
sub    %rcx,%rax
mov    %eax,%ecx
mov    -0x20(%rbp),%rax
mov    %ecx,(%rax)
cmpq   $0xd,-0x78(%rbp)
jle    1274 <jd_to_gregorian+0x174>
mov    -0x78(%rbp),%rax
sub    $0xd,%rax
mov    %eax,%ecx
mov    -0x18(%rbp),%rax
mov    %ecx,(%rax)
jmp    1284 <jd_to_gregorian+0x184>
mov    -0x78(%rbp),%rax
sub    $0x1,%rax
mov    %eax,%ecx
mov    -0x18(%rbp),%rax
mov    %ecx,(%rax)
mov    -0x18(%rbp),%rax
cmpl   $0x2,(%rax)
jne    12a8 <jd_to_gregorian+0x1a8>
mov    -0x20(%rbp),%rax
cmpl   $0x1c,(%rax)
jle    12a8 <jd_to_gregorian+0x1a8>
mov    -0x20(%rbp),%rax
movl   $0x1d,(%rax)
mov    -0x18(%rbp),%rax
cmpl   $0x2,(%rax)
jne    12e4 <jd_to_gregorian+0x1e4>
mov    -0x20(%rbp),%rax
cmpl   $0x1d,(%rax)
jne    12e4 <jd_to_gregorian+0x1e4>
cmpq   $0x3,-0x78(%rbp)
jne    12e4 <jd_to_gregorian+0x1e4>
mov    -0x68(%rbp),%rax
sub    $0x126c,%rax
mov    %eax,%ecx
mov    -0x10(%rbp),%rax
mov    %ecx,(%rax)
jmp    131f <jd_to_gregorian+0x21f>
mov    -0x18(%rbp),%rax
cmpl   $0x2,(%rax)
jle    1308 <jd_to_gregorian+0x208>
mov    -0x68(%rbp),%rax
sub    $0x126c,%rax
mov    %eax,%ecx
mov    -0x10(%rbp),%rax
mov    %ecx,(%rax)
jmp    131a <jd_to_gregorian+0x21a>
mov    -0x68(%rbp),%rax
sub    $0x126b,%rax
mov    %eax,%ecx
mov    -0x10(%rbp),%rax
mov    %ecx,(%rax)
jmp    131f <jd_to_gregorian+0x21f>
movsd  0xd19(%rip),%xmm0
mulsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x80(%rbp)
cvttsd2si -0x80(%rbp),%ecx
mov    -0x28(%rbp),%rax
mov    %ecx,(%rax)
movsd  -0x80(%rbp),%xmm0
mov    -0x28(%rbp),%rax
cvtsi2sdl (%rax),%xmm1
subsd  %xmm1,%xmm0
movsd  0xce3(%rip),%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
cvttsd2si -0x88(%rbp),%ecx
ff
mov    -0x30(%rbp),%rax
mov    %ecx,(%rax)
movsd  -0x88(%rbp),%xmm0
ff
mov    -0x30(%rbp),%rax
cvtsi2sdl (%rax),%xmm1
subsd  %xmm1,%xmm0
movsd  0xcad(%rip),%xmm1
mulsd  %xmm1,%xmm0
mov    -0x38(%rbp),%rax
movsd  %xmm0,(%rax)
add    $0x8,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<jd_to_gregorian>:
push   %rbx
addsd  0xef7(%rip),%xmm0
cvttsd2si %xmm0,%rbx
cvtsi2sd %rbx,%xmm1
cmp    $0x231519,%rbx
jl     1143 <jd_to_gregorian+0x43>
movsd  0xee4(%rip),%xmm2
addsd  %xmm1,%xmm2
divsd  0xee0(%rip),%xmm2
cvttsd2si %xmm2,%r10
add    %r10,%rbx
sar    $0x2,%r10
sub    %r10,%rbx
add    $0x1,%rbx
add    $0x5f4,%rbx
xorps  %xmm2,%xmm2
cvtsi2sd %rbx,%xmm2
addsd  0xebe(%rip),%xmm2
movsd  0xebe(%rip),%xmm3
divsd  %xmm3,%xmm2
cvttsd2si %xmm2,%r10
xorps  %xmm2,%xmm2
cvtsi2sd %r10,%xmm2
mulsd  %xmm3,%xmm2
cvttsd2si %xmm2,%rax
sub    %rax,%rbx
xorps  %xmm2,%xmm2
cvtsi2sd %rbx,%xmm2
movsd  0xe99(%rip),%xmm3
divsd  %xmm3,%xmm2
cvttsd2si %xmm2,%r11
xorps  %xmm2,%xmm2
cvtsi2sd %r11,%xmm2
mulsd  %xmm3,%xmm2
cvttsd2si %xmm2,%rax
sub    %eax,%ebx
mov    %ebx,(%rdx)
lea    -0x1(%r11),%ebx
lea    -0xd(%r11),%eax
cmp    $0xe,%r11
cmovl  %ebx,%eax
mov    %eax,(%rsi)
cmp    $0x2,%eax
jne    11e2 <jd_to_gregorian+0xe2>
cmpl   $0x1d,(%rdx)
jl     11ec <jd_to_gregorian+0xec>
movl   $0x1d,(%rdx)
mov    (%rsi),%eax
cmp    $0x2,%eax
jne    11e2 <jd_to_gregorian+0xe2>
mov    $0xffffed94,%edx
cmp    $0x3,%r11
jne    11ec <jd_to_gregorian+0xec>
jmp    11f1 <jd_to_gregorian+0xf1>
mov    $0xffffed94,%edx
cmp    $0x2,%eax
jg     11f1 <jd_to_gregorian+0xf1>
mov    $0xffffed95,%edx
add    %edx,%r10d
mov    %r10d,(%rdi)
subsd  %xmm1,%xmm0
addsd  0xe2d(%rip),%xmm0
mulsd  0xe2d(%rip),%xmm0
cvttsd2si %xmm0,%eax
cvttpd2dq %xmm0,%xmm1
mov    %eax,(%rcx)
cvtdq2pd %xmm1,%xmm1
subsd  %xmm1,%xmm0
movsd  0xe1b(%rip),%xmm1
mulsd  %xmm1,%xmm0
cvttsd2si %xmm0,%eax
cvttpd2dq %xmm0,%xmm2
cvtdq2pd %xmm2,%xmm2
mov    %eax,(%r8)
subsd  %xmm2,%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,(%r9)
pop    %rbx
ret

LLVM IR (-O0):
define void @jd_to_gregorian(double %0, i32* %1, i32* %2, i32* %3, i32* %4, i32* %5, double* %6) {
%8 = alloca double
%9 = alloca i32*
%10 = alloca i32*
%11 = alloca i32*
%12 = alloca i32*
%13 = alloca i32*
%14 = alloca double*
%15 = alloca i64
%16 = alloca double
%17 = alloca i64
%18 = alloca i64
%19 = alloca i64
%20 = alloca i64
%21 = alloca i64
%22 = alloca i64
%23 = alloca double
%24 = alloca double
store double %0, double* %8
store i32* %1, i32** %9
store i32* %2, i32** %10
store i32* %3, i32** %11
store i32* %4, i32** %12
store i32* %5, i32** %13
store double* %6, double** %14
%25 = load double, double* %8
%26 = fadd double %25, 5.000000e-01
%27 = fptosi double %26 to i64
store i64 %27, i64* %15
%28 = load double, double* %8
%29 = fadd double %28, 5.000000e-01
%30 = load i64, i64* %15
%31 = sitofp i64 %30 to double
%32 = fsub double %29, %31
%33 = fadd double %32, 1.000000e-10
store double %33, double* %16
%34 = load i64, i64* %15
store i64 %34, i64* %17
%35 = load i64, i64* %15
%36 = icmp sge i64 %35, 2299161
br i1 %36, label %37, label %50
37: ; preds = %7
%38 = load i64, i64* %15
%39 = sitofp i64 %38 to double
%40 = fsub double %39, 0x413C7DD040000000
%41 = fdiv double %40, 0x40E1D58800000000
%42 = fptosi double %41 to i64
store i64 %42, i64* %18
%43 = load i64, i64* %15
%44 = add nsw i64 %43, 1
%45 = load i64, i64* %18
%46 = add nsw i64 %44, %45
%47 = load i64, i64* %18
%48 = ashr i64 %47, 2
%49 = sub nsw i64 %46, %48
store i64 %49, i64* %17
br label %50
50: ; preds = %37, %7
%51 = load i64, i64* %17
%52 = add nsw i64 %51, 1524
store i64 %52, i64* %19
%53 = load i64, i64* %19
%54 = sitofp i64 %53 to double
%55 = fsub double %54, 1.221000e+02
%56 = fdiv double %55, 3.652500e+02
%57 = fptosi double %56 to i64
store i64 %57, i64* %20
%58 = load i64, i64* %20
%59 = sitofp i64 %58 to double
%60 = fmul double %59, 3.652500e+02
%61 = fptosi double %60 to i64
store i64 %61, i64* %21
%62 = load i64, i64* %19
%63 = load i64, i64* %21
%64 = sub nsw i64 %62, %63
%65 = sitofp i64 %64 to double
%66 = fdiv double %65, 3.060010e+01
%67 = fptosi double %66 to i64
store i64 %67, i64* %22
%68 = load i64, i64* %19
%69 = load i64, i64* %21
%70 = sub nsw i64 %68, %69
%71 = load i64, i64* %22
%72 = sitofp i64 %71 to double
%73 = fmul double %72, 3.060010e+01
%74 = fptosi double %73 to i64
%75 = sub nsw i64 %70, %74
%76 = trunc i64 %75 to i32
%77 = load i32*, i32** %11
store i32 %76, i32* %77
%78 = load i64, i64* %22
%79 = icmp sgt i64 %78, 13
br i1 %79, label %80, label %85
80: ; preds = %50
%81 = load i64, i64* %22
%82 = sub nsw i64 %81, 13
%83 = trunc i64 %82 to i32
%84 = load i32*, i32** %10
store i32 %83, i32* %84
br label %90
85: ; preds = %50
%86 = load i64, i64* %22
%87 = sub nsw i64 %86, 1
%88 = trunc i64 %87 to i32
%89 = load i32*, i32** %10
store i32 %88, i32* %89
br label %90
90: ; preds = %85, %80
%91 = load i32*, i32** %10
%92 = load i32, i32* %91
%93 = icmp eq i32 %92, 2
br i1 %93, label %94, label %100
94: ; preds = %90
%95 = load i32*, i32** %11
%96 = load i32, i32* %95
%97 = icmp sgt i32 %96, 28
br i1 %97, label %98, label %100
98: ; preds = %94
%99 = load i32*, i32** %11
store i32 29, i32* %99
br label %100
100: ; preds = %98, %94, %90
%101 = load i32*, i32** %10
%102 = load i32, i32* %101
%103 = icmp eq i32 %102, 2
br i1 %103, label %104, label %116
104: ; preds = %100
%105 = load i32*, i32** %11
%106 = load i32, i32* %105
%107 = icmp eq i32 %106, 29
br i1 %107, label %108, label %116
108: ; preds = %104
%109 = load i64, i64* %22
%110 = icmp eq i64 %109, 3
br i1 %110, label %111, label %116
111: ; preds = %108
%112 = load i64, i64* %20
%113 = sub nsw i64 %112, 4716
%114 = trunc i64 %113 to i32
%115 = load i32*, i32** %9
store i32 %114, i32* %115
br label %131
116: ; preds = %108, %104, %100
%117 = load i32*, i32** %10
%118 = load i32, i32* %117
%119 = icmp sgt i32 %118, 2
br i1 %119, label %120, label %125
120: ; preds = %116
%121 = load i64, i64* %20
%122 = sub nsw i64 %121, 4716
%123 = trunc i64 %122 to i32
%124 = load i32*, i32** %9
store i32 %123, i32* %124
br label %130
125: ; preds = %116
%126 = load i64, i64* %20
%127 = sub nsw i64 %126, 4715
%128 = trunc i64 %127 to i32
%129 = load i32*, i32** %9
store i32 %128, i32* %129
br label %130
130: ; preds = %125, %120
br label %131
131: ; preds = %130, %111
%132 = load double, double* %16
%133 = fmul double %132, 2.400000e+01
store double %133, double* %23
%134 = load double, double* %23
%135 = fptosi double %134 to i32
%136 = load i32*, i32** %12
store i32 %135, i32* %136
%137 = load double, double* %23
%138 = load i32*, i32** %12
%139 = load i32, i32* %138
%140 = sitofp i32 %139 to double
%141 = fsub double %137, %140
%142 = fmul double %141, 6.000000e+01
store double %142, double* %24
%143 = load double, double* %24
%144 = fptosi double %143 to i32
%145 = load i32*, i32** %13
store i32 %144, i32* %145
%146 = load double, double* %24
%147 = load i32*, i32** %13
%148 = load i32, i32* %147
%149 = sitofp i32 %148 to double
%150 = fsub double %146, %149
%151 = fmul double %150, 6.000000e+01
%152 = load double*, double** %14
store double %151, double* %152
ret void
}

LLVM IR (-O2):
define void @jd_to_gregorian(double %0, i32* nocapture writeonly %1, i32* nocapture %2, i32* nocapture %3, i32* nocapture writeonly %4, i32* nocapture writeonly %5, double* nocapture writeonly %6) {
%8 = fadd double %0, 5.000000e-01
%9 = fptosi double %8 to i64
%10 = sitofp i64 %9 to double
%11 = icmp sgt i64 %9, 2299160
br i1 %11, label %12, label %20
12: ; preds = %7
%13 = fadd double %10, 0xC13C7DD040000000
%14 = fdiv double %13, 0x40E1D58800000000
%15 = fptosi double %14 to i64
%16 = add nuw nsw i64 %9, 1
%17 = add nsw i64 %16, %15
%18 = ashr i64 %15, 2
%19 = sub i64 %17, %18
br label %20
20: ; preds = %12, %7
%21 = phi i64 [ %19, %12 ], [ %9, %7 ]
%22 = add nsw i64 %21, 1524
%23 = sitofp i64 %22 to double
%24 = fadd double %23, -1.221000e+02
%25 = fdiv double %24, 3.652500e+02
%26 = fptosi double %25 to i64
%27 = sitofp i64 %26 to double
%28 = fmul double %27, 3.652500e+02
%29 = fptosi double %28 to i64
%30 = sub nsw i64 %22, %29
%31 = sitofp i64 %30 to double
%32 = fdiv double %31, 3.060010e+01
%33 = fptosi double %32 to i64
%34 = sitofp i64 %33 to double
%35 = fmul double %34, 3.060010e+01
%36 = fptosi double %35 to i64
%37 = sub nsw i64 %30, %36
%38 = trunc i64 %37 to i32
store i32 %38, i32* %3
%39 = icmp sgt i64 %33, 13
%40 = trunc i64 %33 to i32
%41 = select i1 %39, i32 -13, i32 -1
%42 = add i32 %41, %40
store i32 %42, i32* %2
%43 = icmp eq i32 %42, 2
br i1 %43, label %44, label %52
44: ; preds = %20
%45 = load i32, i32* %3
%46 = icmp sgt i32 %45, 28
br i1 %46, label %47, label %55
47: ; preds = %44
store i32 29, i32* %3
%48 = load i32, i32* %2
%49 = icmp eq i32 %48, 2
br i1 %49, label %50, label %52
50: ; preds = %47
%51 = icmp eq i64 %33, 3
br i1 %51, label %56, label %55
52: ; preds = %20, %47
%53 = phi i32 [ %48, %47 ], [ %42, %20 ]
%54 = icmp sgt i32 %53, 2
br i1 %54, label %56, label %55
55: ; preds = %44, %50, %52
br label %56
56: ; preds = %52, %50, %55
%57 = phi i32 [ -4715, %55 ], [ -4716, %50 ], [ -4716, %52 ]
%58 = trunc i64 %26 to i32
%59 = add i32 %57, %58
store i32 %59, i32* %1
%60 = fsub double %8, %10
%61 = fadd double %60, 1.000000e-10
%62 = fmul double %61, 2.400000e+01
%63 = fptosi double %62 to i32
store i32 %63, i32* %4
%64 = sitofp i32 %63 to double
%65 = fsub double %62, %64
%66 = fmul double %65, 6.000000e+01
%67 = fptosi double %66 to i32
store i32 %67, i32* %5
%68 = sitofp i32 %67 to double
%69 = fsub double %66, %68
%70 = fmul double %69, 6.000000e+01
store double %70, double* %6
ret void
}