void jul2greg(double jd, int *m, int *d, int *y, int *h, int *mi, double *sec) {
  long j = jd;
  double tmp, frac = jd - j;
  frac += 0.0005;
  if (frac >= 0.5) {
    frac = frac - 0.5;
    j++;
  } else {
    frac = frac + 0.5;
  }
  j -= 1721119L;
  *y = (4L * j - 1L) / 146097L;
  j = 4L * j - 1L - 146097L * *y;
  *d = j / 4L;
  j = (4L * *d + 3L) / 1461L;
  *d = 4L * *d + 3L - 1461L * j;
  *d = (*d + 4L) / 4L;
  *m = (5L * *d - 3L) / 153L;
  *d = 5L * *d - 3 - 153L * *m;
  *d = (*d + 5L) / 5L;
  *y = 100L * *y + j;
  if (*m < 10)
    *m += 3;
  else {
    *m -= 9;
    *y += 1;
  }
  tmp = 3600.0 * (frac * 24.0);
  *h = (int)(tmp / 3600.0);
  tmp = tmp - *h * 3600.0;
  *mi = (int)(tmp / 60.0);
  *sec = tmp - *mi * 60.0;
}


x86 assembly (clang, -O0):
<jul2greg>:
push   %rbp
mov    %rsp,%rbp
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)
cvttsd2si -0x8(%rbp),%rax
mov    %rax,-0x40(%rbp)
movsd  -0x8(%rbp),%xmm0
cvtsi2sdq -0x40(%rbp),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  0xec1(%rip),%xmm0
addsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x50(%rbp),%xmm0
movsd  0xea2(%rip),%xmm1
ucomisd %xmm1,%xmm0
jb     118f <jul2greg+0x8f>
movsd  -0x50(%rbp),%xmm0
movsd  0xe8b(%rip),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x40(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x40(%rbp)
jmp    11a1 <jul2greg+0xa1>
movsd  0xe69(%rip),%xmm0
addsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x40(%rbp),%rax
sub    $0x1a431f,%rax
mov    %rax,-0x40(%rbp)
mov    -0x40(%rbp),%rax
shl    $0x2,%rax
sub    $0x1,%rax
mov    $0x23ab1,%ecx
cqto
idiv   %rcx
mov    %eax,%ecx
mov    -0x20(%rbp),%rax
mov    %ecx,(%rax)
mov    -0x40(%rbp),%rax
shl    $0x2,%rax
sub    $0x1,%rax
mov    -0x20(%rbp),%rcx
movslq (%rcx),%rcx
imul   $0x23ab1,%rcx,%rcx
sub    %rcx,%rax
mov    %rax,-0x40(%rbp)
mov    -0x40(%rbp),%rax
mov    $0x4,%ecx
cqto
idiv   %rcx
mov    %eax,%ecx
mov    -0x18(%rbp),%rax
mov    %ecx,(%rax)
mov    -0x18(%rbp),%rax
movslq (%rax),%rax
shl    $0x2,%rax
add    $0x3,%rax
mov    $0x5b5,%ecx
cqto
idiv   %rcx
mov    %rax,-0x40(%rbp)
mov    -0x18(%rbp),%rax
movslq (%rax),%rax
shl    $0x2,%rax
add    $0x3,%rax
imul   $0x5b5,-0x40(%rbp),%rcx
sub    %rcx,%rax
mov    %eax,%ecx
mov    -0x18(%rbp),%rax
mov    %ecx,(%rax)
mov    -0x18(%rbp),%rax
movslq (%rax),%rax
add    $0x4,%rax
mov    $0x4,%ecx
cqto
idiv   %rcx
mov    %eax,%ecx
mov    -0x18(%rbp),%rax
mov    %ecx,(%rax)
mov    -0x18(%rbp),%rax
movslq (%rax),%rax
imul   $0x5,%rax,%rax
sub    $0x3,%rax
mov    $0x99,%ecx
cqto
idiv   %rcx
mov    %eax,%ecx
mov    -0x10(%rbp),%rax
mov    %ecx,(%rax)
mov    -0x18(%rbp),%rax
movslq (%rax),%rax
imul   $0x5,%rax,%rax
sub    $0x3,%rax
mov    -0x10(%rbp),%rcx
movslq (%rcx),%rcx
imul   $0x99,%rcx,%rcx
sub    %rcx,%rax
mov    %eax,%ecx
mov    -0x18(%rbp),%rax
mov    %ecx,(%rax)
mov    -0x18(%rbp),%rax
movslq (%rax),%rax
add    $0x5,%rax
mov    $0x5,%ecx
cqto
idiv   %rcx
mov    %eax,%ecx
mov    -0x18(%rbp),%rax
mov    %ecx,(%rax)
mov    -0x20(%rbp),%rax
movslq (%rax),%rax
imul   $0x64,%rax,%rax
add    -0x40(%rbp),%rax
mov    %eax,%ecx
mov    -0x20(%rbp),%rax
mov    %ecx,(%rax)
mov    -0x10(%rbp),%rax
cmpl   $0xa,(%rax)
jge    12fa <jul2greg+0x1fa>
mov    -0x10(%rbp),%rax
mov    (%rax),%ecx
add    $0x3,%ecx
mov    %ecx,(%rax)
jmp    1310 <jul2greg+0x210>
mov    -0x10(%rbp),%rax
mov    (%rax),%ecx
sub    $0x9,%ecx
mov    %ecx,(%rax)
mov    -0x20(%rbp),%rax
mov    (%rax),%ecx
add    $0x1,%ecx
mov    %ecx,(%rax)
movsd  0xd08(%rip),%xmm1
mulsd  -0x50(%rbp),%xmm1
movsd  0xcf3(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x48(%rbp),%xmm0
movsd  0xcdd(%rip),%xmm1
divsd  %xmm1,%xmm0
cvttsd2si %xmm0,%ecx
mov    -0x28(%rbp),%rax
mov    %ecx,(%rax)
movsd  -0x48(%rbp),%xmm1
mov    -0x28(%rbp),%rax
cvtsi2sdl (%rax),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  0xca3(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x48(%rbp),%xmm0
movsd  0xc81(%rip),%xmm1
divsd  %xmm1,%xmm0
cvttsd2si %xmm0,%ecx
mov    -0x30(%rbp),%rax
mov    %ecx,(%rax)
movsd  -0x48(%rbp),%xmm1
mov    -0x30(%rbp),%rax
cvtsi2sdl (%rax),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  0xc47(%rip),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x38(%rbp),%rax
movsd  %xmm0,(%rax)
pop    %rbp
ret

x86 assembly (clang, -O2):
<jul2greg>:
push   %r14
push   %rbx
cvttsd2si %xmm0,%rbx
cvtsi2sd %rbx,%xmm1
subsd  %xmm1,%xmm0
addsd  0xee7(%rip),%xmm0
mov    %rdx,%r10
ucomisd 0xee4(%rip),%xmm0
sbb    $0xffffffffffffffff,%rbx
lea    -0x690c7d(,%rbx,4),%rax
ff
movabs $0x396b06bcc8f862ed,%rdx
06 6b 39
imul   %rdx
mov    %rdx,%rax
shr    $0x3f,%rax
sar    $0xf,%rdx
add    %rax,%rdx
movslq %edx,%rax
imul   $0xfffffffffffdc54f,%rax,%rax
lea    (%rax,%rbx,4),%r11
add    $0xffffffffff96f383,%r11
lea    (%rax,%rbx,4),%rax
add    $0xffffffffff96f386,%rax
test   %r11,%r11
cmovns %r11,%rax
mov    %edx,(%r10)
shr    $0x2,%rax
movslq %eax,%rbx
lea    0x3(,%rbx,4),%r11
movabs $0xb36d8397914268c9,%rdx
83 6d b3
mov    %r11,%rax
imul   %rdx
lea    (%rdx,%rbx,4),%r14
add    $0x3,%r14
mov    %r14,%rax
shr    $0x3f,%rax
shr    $0xa,%r14
add    %eax,%r14d
imul   $0xfffffa4b,%r14d,%eax
add    %eax,%r11d
movslq %r11d,%rbx
lea    0x4(%rbx),%rax
add    $0x7,%rbx
test   %rax,%rax
cmovns %rax,%rbx
sar    $0x2,%rbx
lea    (%rbx,%rbx,4),%rax
add    $0xfffffffffffffffd,%rax
movabs $0x358ae0358ae0358b,%rdx
e0 8a 35
imul   %rdx
mov    %ebx,(%rsi)
mov    %rdx,%rax
shr    $0x3f,%rax
shr    $0x5,%rdx
add    %eax,%edx
mov    %edx,(%rdi)
mov    (%rsi),%eax
lea    (%rax,%rax,4),%eax
imul   $0xffffff67,%edx,%edx
add    %edx,%eax
add    $0xfffffffd,%eax
cltq
add    $0x5,%rax
movabs $0x6666666666666667,%rdx
66 66 66
imul   %rdx
mov    %rdx,%rax
shr    $0x3f,%rax
shr    %rdx
add    %eax,%edx
mov    %edx,(%rsi)
imul   $0x64,(%r10),%eax
add    %r14d,%eax
mov    %eax,(%r10)
mov    (%rdi),%eax
cmp    $0x9,%eax
jg     1238 <jul2greg+0x138>
add    $0x3,%eax
mov    %eax,(%rdi)
jmp    1241 <jul2greg+0x141>
add    $0xfffffff7,%eax
mov    %eax,(%rdi)
addl   $0x1,(%r10)
xor    %eax,%eax
ucomisd 0xdbd(%rip),%xmm0
setb   %al
lea    0xdd3(%rip),%rdx
addsd  (%rdx,%rax,8),%xmm0
mulsd  0xdae(%rip),%xmm0
movsd  0xdae(%rip),%xmm1
mulsd  %xmm1,%xmm0
movapd %xmm0,%xmm2
divsd  %xmm1,%xmm2
cvttsd2si %xmm2,%eax
mov    %eax,(%rcx)
cvttpd2dq %xmm2,%xmm2
cvtdq2pd %xmm2,%xmm2
mulsd  %xmm1,%xmm2
subsd  %xmm2,%xmm0
movsd  0xd8c(%rip),%xmm1
movapd %xmm0,%xmm2
divsd  %xmm1,%xmm2
cvttsd2si %xmm2,%eax
cvttpd2dq %xmm2,%xmm2
cvtdq2pd %xmm2,%xmm2
mov    %eax,(%r8)
mulsd  %xmm1,%xmm2
subsd  %xmm2,%xmm0
movsd  %xmm0,(%r9)
pop    %rbx
pop    %r14
ret

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

LLVM IR (-O2):
define void @jul2greg(double %0, i32* nocapture %1, i32* nocapture %2, i32* nocapture %3, i32* nocapture writeonly %4, i32* nocapture writeonly %5, double* nocapture writeonly %6) {
%8 = fptosi double %0 to i64
%9 = sitofp i64 %8 to double
%10 = fsub double %0, %9
%11 = fadd double %10, 5.000000e-04
%12 = fcmp ult double %11, 5.000000e-01
%13 = xor i1 %12, true
%14 = zext i1 %13 to i64
%15 = add nsw i64 %14, %8
%16 = shl i64 %15, 2
%17 = add i64 %16, -6884477
%18 = sdiv i64 %17, 146097
%19 = trunc i64 %18 to i32
store i32 %19, i32* %3
%20 = shl i64 %18, 32
%21 = ashr exact i64 %20, 32
%22 = mul nsw i64 %21, -146097
%23 = add i64 %22, %17
%24 = sdiv i64 %23, 4
%25 = shl i64 %24, 32
%26 = ashr exact i64 %25, 30
%27 = or i64 %26, 3
%28 = sdiv i64 %27, 1461
%29 = mul nsw i64 %28, -1461
%30 = add nsw i64 %29, %27
%31 = shl i64 %30, 32
%32 = ashr exact i64 %31, 32
%33 = add nsw i64 %32, 4
%34 = sdiv i64 %33, 4
%35 = trunc i64 %34 to i32
store i32 %35, i32* %2
%36 = mul nsw i64 %34, 5
%37 = add nsw i64 %36, -3
%38 = sdiv i64 %37, 153
%39 = trunc i64 %38 to i32
store i32 %39, i32* %1
%40 = load i32, i32* %2
%41 = zext i32 %40 to i64
%42 = mul nuw nsw i64 %41, 5
%43 = mul nsw i64 %38, -153
%44 = add nsw i64 %43, 4294967293
%45 = add nuw nsw i64 %44, %42
%46 = shl i64 %45, 32
%47 = ashr exact i64 %46, 32
%48 = add nsw i64 %47, 5
%49 = sdiv i64 %48, 5
%50 = trunc i64 %49 to i32
store i32 %50, i32* %2
%51 = load i32, i32* %3
%52 = zext i32 %51 to i64
%53 = mul nuw nsw i64 %52, 100
%54 = add nsw i64 %53, %28
%55 = trunc i64 %54 to i32
store i32 %55, i32* %3
%56 = load i32, i32* %1
%57 = icmp slt i32 %56, 10
br i1 %57, label %58, label %60
58: ; preds = %7
%59 = add nsw i32 %56, 3
store i32 %59, i32* %1
br label %64
60: ; preds = %7
%61 = add nsw i32 %56, -9
store i32 %61, i32* %1
%62 = load i32, i32* %3
%63 = add nsw i32 %62, 1
store i32 %63, i32* %3
br label %64
64: ; preds = %60, %58
%65 = select i1 %12, double 5.000000e-01, double -5.000000e-01
%66 = fadd double %11, %65
%67 = fmul double %66, 2.400000e+01
%68 = fmul double %67, 3.600000e+03
%69 = fdiv double %68, 3.600000e+03
%70 = fptosi double %69 to i32
store i32 %70, i32* %4
%71 = sitofp i32 %70 to double
%72 = fneg double %71
%73 = tail call double @llvm.fmuladd.f64(double %72, double 3.600000e+03, double %68)
%74 = fdiv double %73, 6.000000e+01
%75 = fptosi double %74 to i32
store i32 %75, i32* %5
%76 = sitofp i32 %75 to double
%77 = fneg double %76
%78 = tail call double @llvm.fmuladd.f64(double %77, double 6.000000e+01, double %73)
store double %78, double* %6
ret void
}