long grib_julian_to_date(long jdate) {
  long x, y, d, m, e;
  long day, month, year;
  x = 4 * jdate - 6884477;
  y = (x / 146097) * 100;
  e = x % 146097;
  d = e / 4;
  x = 4 * d + 3;
  y = (x / 1461) + y;
  e = x % 1461;
  d = e / 4 + 1;
  x = 5 * d - 3;
  m = x / 153 + 1;
  e = x % 153;
  d = e / 5 + 1;
  if (m < 11)
    month = m + 2;
  else
    month = m - 10;
  day = d;
  year = y + m / 11;
  return year * 10000 + month * 100 + day;
}


x86 assembly (clang, -O0):
<grib_julian_to_date>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    -0x8(%rbp),%rax
shl    $0x2,%rax
sub    $0x690c7d,%rax
mov    %rax,-0x10(%rbp)
mov    -0x10(%rbp),%rax
mov    $0x23ab1,%ecx
cqto
idiv   %rcx
imul   $0x64,%rax,%rax
mov    %rax,-0x18(%rbp)
mov    -0x10(%rbp),%rax
mov    $0x23ab1,%ecx
cqto
idiv   %rcx
mov    %rdx,-0x30(%rbp)
mov    -0x30(%rbp),%rax
mov    $0x4,%ecx
cqto
idiv   %rcx
mov    %rax,-0x20(%rbp)
mov    -0x20(%rbp),%rax
shl    $0x2,%rax
add    $0x3,%rax
mov    %rax,-0x10(%rbp)
mov    -0x10(%rbp),%rax
mov    $0x5b5,%ecx
cqto
idiv   %rcx
add    -0x18(%rbp),%rax
mov    %rax,-0x18(%rbp)
mov    -0x10(%rbp),%rax
mov    $0x5b5,%ecx
cqto
idiv   %rcx
mov    %rdx,-0x30(%rbp)
mov    -0x30(%rbp),%rax
mov    $0x4,%ecx
cqto
idiv   %rcx
add    $0x1,%rax
mov    %rax,-0x20(%rbp)
imul   $0x5,-0x20(%rbp),%rax
sub    $0x3,%rax
mov    %rax,-0x10(%rbp)
mov    -0x10(%rbp),%rax
mov    $0x99,%ecx
cqto
idiv   %rcx
add    $0x1,%rax
mov    %rax,-0x28(%rbp)
mov    -0x10(%rbp),%rax
mov    $0x99,%ecx
cqto
idiv   %rcx
mov    %rdx,-0x30(%rbp)
mov    -0x30(%rbp),%rax
mov    $0x5,%ecx
cqto
idiv   %rcx
add    $0x1,%rax
mov    %rax,-0x20(%rbp)
cmpq   $0xb,-0x28(%rbp)
jge    1209 <grib_julian_to_date+0x109>
mov    -0x28(%rbp),%rax
add    $0x2,%rax
mov    %rax,-0x40(%rbp)
jmp    1215 <grib_julian_to_date+0x115>
mov    -0x28(%rbp),%rax
sub    $0xa,%rax
mov    %rax,-0x40(%rbp)
mov    -0x20(%rbp),%rax
mov    %rax,-0x38(%rbp)
mov    -0x18(%rbp),%rax
mov    %rax,-0x50(%rbp)
mov    -0x28(%rbp),%rax
mov    $0xb,%ecx
cqto
idiv   %rcx
mov    %rax,%rcx
mov    -0x50(%rbp),%rax
add    %rcx,%rax
mov    %rax,-0x48(%rbp)
imul   $0x2710,-0x48(%rbp),%rax
imul   $0x64,-0x40(%rbp),%rcx
add    %rcx,%rax
add    -0x38(%rbp),%rax
pop    %rbp
ret

x86 assembly (clang, -O2):
<grib_julian_to_date>:
lea    -0x690c7d(,%rdi,4),%rsi
ff
movabs $0x396b06bcc8f862ed,%rcx
06 6b 39
mov    %rsi,%rax
imul   %rcx
mov    %rdx,%r10
mov    %rdx,%rax
shr    $0x3f,%rax
sar    $0xf,%r10
add    %rax,%r10
imul   $0x23ab1,%r10d,%eax
sub    %eax,%esi
lea    0x3(%rsi),%eax
test   %esi,%esi
cmovns %esi,%eax
sar    $0x2,%eax
movslq %eax,%rdi
lea    0x3(,%rdi,4),%rsi
movabs $0xb36d8397914268c9,%rdx
83 6d b3
mov    %rsi,%rax
imul   %rdx
lea    (%rdx,%rdi,4),%rax
add    $0x3,%rax
mov    %rax,%rdx
shr    $0x3f,%rdx
sar    $0xa,%rax
add    %rdx,%rax
imul   $0x5b5,%eax,%edx
sub    %edx,%esi
lea    0x3(%rsi),%edx
test   %si,%si
cmovns %esi,%edx
movswl %dx,%edx
sar    $0x2,%edx
movslq %edx,%rdx
lea    (%rdx,%rdx,4),%r8
lea    0x2(%r8),%esi
movswl %si,%edi
imul   $0x358b,%edi,%edi
mov    %edi,%edx
shr    $0x1f,%edx
sar    $0x15,%edi
add    %edx,%edi
movslq %edi,%rdx
lea    -0x9(%rdx),%rdi
lea    0x3(%rdx),%r9
cmp    $0x5f8,%r8
cmovl  %r9,%rdi
imul   $0x99,%edx,%ecx
sub    %ecx,%esi
movswl %si,%ecx
imul   $0x6667,%ecx,%ecx
mov    %ecx,%esi
shr    $0x1f,%esi
sar    $0x11,%ecx
add    %esi,%ecx
movslq %ecx,%rcx
imul   $0x64,%rdi,%rsi
add    %rcx,%rsi
imul   $0x64,%r10,%rcx
add    %rax,%rcx
add    $0x1,%dl
movsbl %dl,%eax
imul   $0xffffffbb,%eax,%edx
shr    $0x8,%edx
add    %dl,%al
mov    %eax,%edx
shr    $0x7,%dl
sar    $0x3,%al
add    %dl,%al
movzbl %al,%eax
add    %rcx,%rax
imul   $0x2710,%rax,%rax
add    %rsi,%rax
add    $0x1,%rax
ret

LLVM IR (-O0):
define i64 @grib_julian_to_date(i64 %0) {
%2 = alloca i64
%3 = alloca i64
%4 = alloca i64
%5 = alloca i64
%6 = alloca i64
%7 = alloca i64
%8 = alloca i64
%9 = alloca i64
%10 = alloca i64
store i64 %0, i64* %2
%11 = load i64, i64* %2
%12 = mul nsw i64 4, %11
%13 = sub nsw i64 %12, 6884477
store i64 %13, i64* %3
%14 = load i64, i64* %3
%15 = sdiv i64 %14, 146097
%16 = mul nsw i64 %15, 100
store i64 %16, i64* %4
%17 = load i64, i64* %3
%18 = srem i64 %17, 146097
store i64 %18, i64* %7
%19 = load i64, i64* %7
%20 = sdiv i64 %19, 4
store i64 %20, i64* %5
%21 = load i64, i64* %5
%22 = mul nsw i64 4, %21
%23 = add nsw i64 %22, 3
store i64 %23, i64* %3
%24 = load i64, i64* %3
%25 = sdiv i64 %24, 1461
%26 = load i64, i64* %4
%27 = add nsw i64 %25, %26
store i64 %27, i64* %4
%28 = load i64, i64* %3
%29 = srem i64 %28, 1461
store i64 %29, i64* %7
%30 = load i64, i64* %7
%31 = sdiv i64 %30, 4
%32 = add nsw i64 %31, 1
store i64 %32, i64* %5
%33 = load i64, i64* %5
%34 = mul nsw i64 5, %33
%35 = sub nsw i64 %34, 3
store i64 %35, i64* %3
%36 = load i64, i64* %3
%37 = sdiv i64 %36, 153
%38 = add nsw i64 %37, 1
store i64 %38, i64* %6
%39 = load i64, i64* %3
%40 = srem i64 %39, 153
store i64 %40, i64* %7
%41 = load i64, i64* %7
%42 = sdiv i64 %41, 5
%43 = add nsw i64 %42, 1
store i64 %43, i64* %5
%44 = load i64, i64* %6
%45 = icmp slt i64 %44, 11
br i1 %45, label %46, label %49
46: ; preds = %1
%47 = load i64, i64* %6
%48 = add nsw i64 %47, 2
store i64 %48, i64* %9
br label %52
49: ; preds = %1
%50 = load i64, i64* %6
%51 = sub nsw i64 %50, 10
store i64 %51, i64* %9
br label %52
52: ; preds = %49, %46
%53 = load i64, i64* %5
store i64 %53, i64* %8
%54 = load i64, i64* %4
%55 = load i64, i64* %6
%56 = sdiv i64 %55, 11
%57 = add nsw i64 %54, %56
store i64 %57, i64* %10
%58 = load i64, i64* %10
%59 = mul nsw i64 %58, 10000
%60 = load i64, i64* %9
%61 = mul nsw i64 %60, 100
%62 = add nsw i64 %59, %61
%63 = load i64, i64* %8
%64 = add nsw i64 %62, %63
ret i64 %64
}

LLVM IR (-O2):
define i64 @grib_julian_to_date(i64 %0) {
%2 = shl nsw i64 %0, 2
%3 = add nsw i64 %2, -6884477
%4 = srem i64 %3, 146097
%5 = trunc i64 %4 to i32
%6 = sdiv i32 %5, 4
%7 = sext i32 %6 to i64
%8 = shl nsw i64 %7, 2
%9 = or i64 %8, 3
%10 = srem i64 %9, 1461
%11 = trunc i64 %10 to i16
%12 = sdiv i16 %11, 4
%13 = sext i16 %12 to i64
%14 = mul nsw i64 %13, 5
%15 = trunc i64 %14 to i16
%16 = add nsw i16 %15, 2
%17 = sdiv i16 %16, 153
%18 = sext i16 %17 to i64
%19 = icmp slt i64 %14, 1528
%20 = select i1 %19, i64 3, i64 -9
%21 = add nsw i64 %20, %18
%22 = srem i16 %16, 153
%23 = sdiv i16 %22, 5
%24 = add nsw i16 %23, 1
%25 = sext i16 %24 to i64
%26 = sdiv i64 %9, 1461
%27 = sdiv i64 %3, 146097
%28 = mul nsw i64 %27, 100
%29 = add nsw i64 %26, %28
%30 = trunc i16 %17 to i8
%31 = add nsw i8 %30, 1
%32 = sdiv i8 %31, 11
%33 = zext i8 %32 to i64
%34 = add nsw i64 %29, %33
%35 = mul nsw i64 %34, 10000
%36 = mul nsw i64 %21, 100
%37 = add nsw i64 %36, %25
%38 = add i64 %37, %35
ret i64 %38
}