void cal_mjd(mn, dy, yr, mjd) int mn, yr;
double dy;
double *mjd;
{
  static double last_mjd, last_dy;
  static int last_mn, last_yr;
  int b, d, m, y;
  long c;
  if (mn == last_mn && yr == last_yr && dy == last_dy) {
    *mjd = last_mjd;
    return;
  }
  m = mn;
  y = (yr < 0) ? yr + 1 : yr;
  if (mn < 3) {
    m += 12;
    y -= 1;
  }
  if (yr < 1582 || (yr == 1582 && (mn < 10 || (mn == 10 && dy < 15))))
    b = 0;
  else {
    int a;
    a = y / 100;
    b = 2 - a + a / 4;
  }
  if (y < 0)
    c = (long)((365.25 * y) - 0.75) - 694025L;
  else
    c = (long)(365.25 * y) - 694025L;
  d = 30.6001 * (m + 1);
  *mjd = b + c + d + dy - 0.5;
  last_mn = mn;
  last_dy = dy;
  last_yr = yr;
  last_mjd = *mjd;
}


x86 assembly (clang, -O0):
<cal_mjd>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
movsd  %xmm0,-0x10(%rbp)
mov    %esi,-0x14(%rbp)
mov    %rdx,-0x20(%rbp)
mov    -0x4(%rbp),%eax
cmp    0x2f1c(%rip),%eax
jne    115f <cal_mjd+0x5f>
mov    -0x14(%rbp),%eax
cmp    0x2f11(%rip),%eax
jne    115f <cal_mjd+0x5f>
movsd  -0x10(%rbp),%xmm0
ucomisd 0x2ef2(%rip),%xmm0
jne    115f <cal_mjd+0x5f>
jp     115f <cal_mjd+0x5f>
movsd  0x2ed6(%rip),%xmm0
mov    -0x20(%rbp),%rax
movsd  %xmm0,(%rax)
jmp    12f4 <cal_mjd+0x1f4>
mov    -0x4(%rbp),%eax
mov    %eax,-0x2c(%rbp)
cmpl   $0x0,-0x14(%rbp)
jge    117d <cal_mjd+0x7d>
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    1183 <cal_mjd+0x83>
mov    -0x14(%rbp),%eax
mov    %eax,-0x40(%rbp)
mov    -0x40(%rbp),%eax
mov    %eax,-0x30(%rbp)
cmpl   $0x3,-0x4(%rbp)
jge    11a5 <cal_mjd+0xa5>
mov    -0x2c(%rbp),%eax
add    $0xc,%eax
mov    %eax,-0x2c(%rbp)
mov    -0x30(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x30(%rbp)
cmpl   $0x62e,-0x14(%rbp)
jl     11e6 <cal_mjd+0xe6>
cmpl   $0x62e,-0x14(%rbp)
jne    11f2 <cal_mjd+0xf2>
cmpl   $0xa,-0x4(%rbp)
jl     11e6 <cal_mjd+0xe6>
cmpl   $0xa,-0x4(%rbp)
jne    11f2 <cal_mjd+0xf2>
movsd  0xe25(%rip),%xmm0
ucomisd -0x10(%rbp),%xmm0
jbe    11f2 <cal_mjd+0xf2>
movl   $0x0,-0x24(%rbp)
jmp    1220 <cal_mjd+0x120>
mov    -0x30(%rbp),%eax
mov    $0x64,%ecx
cltd
idiv   %ecx
mov    %eax,-0x3c(%rbp)
mov    $0x2,%eax
sub    -0x3c(%rbp),%eax
mov    %eax,-0x44(%rbp)
mov    -0x3c(%rbp),%eax
mov    $0x4,%ecx
cltd
idiv   %ecx
mov    %eax,%ecx
mov    -0x44(%rbp),%eax
add    %ecx,%eax
mov    %eax,-0x24(%rbp)
cmpl   $0x0,-0x30(%rbp)
jge    125b <cal_mjd+0x15b>
cvtsi2sdl -0x30(%rbp),%xmm0
movsd  0xdd1(%rip),%xmm1
mulsd  %xmm1,%xmm0
movsd  0xdcd(%rip),%xmm1
addsd  %xmm1,%xmm0
cvttsd2si %xmm0,%rax
sub    $0xa9709,%rax
mov    %rax,-0x38(%rbp)
jmp    127b <cal_mjd+0x17b>
cvtsi2sdl -0x30(%rbp),%xmm1
movsd  0xda0(%rip),%xmm0
mulsd  %xmm1,%xmm0
cvttsd2si %xmm0,%rax
sub    $0xa9709,%rax
mov    %rax,-0x38(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
cvtsi2sd %eax,%xmm1
movsd  0xd93(%rip),%xmm0
mulsd  %xmm1,%xmm0
cvttsd2si %xmm0,%eax
mov    %eax,-0x28(%rbp)
movslq -0x24(%rbp),%rax
add    -0x38(%rbp),%rax
movslq -0x28(%rbp),%rcx
add    %rcx,%rax
cvtsi2sd %rax,%xmm0
addsd  -0x10(%rbp),%xmm0
movsd  0xd5f(%rip),%xmm1
subsd  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x4(%rbp),%eax
mov    %eax,0x2d6a(%rip)
movsd  -0x10(%rbp),%xmm0
movsd  %xmm0,0x2d55(%rip)
mov    -0x14(%rbp),%eax
mov    %eax,0x2d58(%rip)
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,0x2d34(%rip)
pop    %rbp
ret

x86 assembly (clang, -O2):
<cal_mjd>:
cmp    %edi,0x2f32(%rip)
jne    112d <cal_mjd+0x2d>
cmp    %esi,0x2f2e(%rip)
jne    112d <cal_mjd+0x2d>
movsd  0x2f18(%rip),%xmm1
ucomisd %xmm0,%xmm1
jne    112d <cal_mjd+0x2d>
jp     112d <cal_mjd+0x2d>
movsd  0x2f00(%rip),%xmm0
movsd  %xmm0,(%rdx)
ret
mov    %esi,%eax
shr    $0x1f,%eax
lea    0xc(%rdi),%r9d
xor    %ecx,%ecx
cmp    $0x3,%edi
cmovge %edi,%r9d
setl   %cl
mov    %esi,%r10d
sub    %ecx,%r10d
add    %eax,%r10d
cmp    $0x62e,%esi
jl     116a <cal_mjd+0x6a>
jne    116f <cal_mjd+0x6f>
cmp    $0xa,%edi
jl     116a <cal_mjd+0x6a>
jne    116f <cal_mjd+0x6f>
movsd  0xe9c(%rip),%xmm1
ucomisd %xmm0,%xmm1
jbe    116f <cal_mjd+0x6f>
xor    %r8d,%r8d
jmp    11a3 <cal_mjd+0xa3>
movslq %r10d,%r8
imul   $0xffffffffae147ae1,%r8,%rcx
mov    %rcx,%rax
shr    $0x3f,%rax
sar    $0x25,%rcx
add    %eax,%ecx
imul   $0x51eb851f,%r8,%rax
mov    %rax,%r8
shr    $0x3f,%r8
sar    $0x27,%rax
add    %r8d,%eax
add    %ecx,%eax
add    $0x2,%eax
movslq %eax,%r8
xorps  %xmm1,%xmm1
cvtsi2sd %r10d,%xmm1
mulsd  0xe55(%rip),%xmm1
test   %r10d,%r10d
jns    11c0 <cal_mjd+0xc0>
addsd  0xe50(%rip),%xmm1
cvttsd2si %xmm1,%rax
add    $0x1,%r9d
xorps  %xmm1,%xmm1
cvtsi2sd %r9d,%xmm1
mulsd  0xe3f(%rip),%xmm1
cvttsd2si %xmm1,%ecx
movslq %ecx,%rcx
add    %r8,%rcx
add    %rcx,%rax
add    $0xfffffffffff568f7,%rax
xorps  %xmm1,%xmm1
cvtsi2sd %rax,%xmm1
addsd  %xmm0,%xmm1
addsd  0xe20(%rip),%xmm1
movsd  %xmm1,(%rdx)
mov    %edi,0x2e2e(%rip)
movsd  %xmm0,0x2e1e(%rip)
mov    %esi,0x2e24(%rip)
movsd  %xmm1,0x2e08(%rip)
ret

LLVM IR (-O0):
@cal_mjd.last_mjd = internal global double 0.000000e+00
@cal_mjd.last_dy = internal global double 0.000000e+00
@cal_mjd.last_mn = internal global i32 0
@cal_mjd.last_yr = internal global i32 0
define void @cal_mjd(i32 %0, double %1, i32 %2, double* %3) {
%5 = alloca i32
%6 = alloca double
%7 = alloca i32
%8 = alloca double*
%9 = alloca i32
%10 = alloca i32
%11 = alloca i32
%12 = alloca i32
%13 = alloca i64
%14 = alloca i32
store i32 %0, i32* %5
store double %1, double* %6
store i32 %2, i32* %7
store double* %3, double** %8
%15 = load i32, i32* %5
%16 = load i32, i32* @cal_mjd.last_mn
%17 = icmp eq i32 %15, %16
br i1 %17, label %18, label %29
18: ; preds = %4
%19 = load i32, i32* %7
%20 = load i32, i32* @cal_mjd.last_yr
%21 = icmp eq i32 %19, %20
br i1 %21, label %22, label %29
22: ; preds = %18
%23 = load double, double* %6
%24 = load double, double* @cal_mjd.last_dy
%25 = fcmp oeq double %23, %24
br i1 %25, label %26, label %29
26: ; preds = %22
%27 = load double, double* @cal_mjd.last_mjd
%28 = load double*, double** %8
store double %27, double* %28
br label %109
29: ; preds = %22, %18, %4
%30 = load i32, i32* %5
store i32 %30, i32* %11
%31 = load i32, i32* %7
%32 = icmp slt i32 %31, 0
br i1 %32, label %33, label %36
33: ; preds = %29
%34 = load i32, i32* %7
%35 = add nsw i32 %34, 1
br label %38
36: ; preds = %29
%37 = load i32, i32* %7
br label %38
38: ; preds = %36, %33
%39 = phi i32 [ %35, %33 ], [ %37, %36 ]
store i32 %39, i32* %12
%40 = load i32, i32* %5
%41 = icmp slt i32 %40, 3
br i1 %41, label %42, label %47
42: ; preds = %38
%43 = load i32, i32* %11
%44 = add nsw i32 %43, 12
store i32 %44, i32* %11
%45 = load i32, i32* %12
%46 = sub nsw i32 %45, 1
store i32 %46, i32* %12
br label %47
47: ; preds = %42, %38
%48 = load i32, i32* %7
%49 = icmp slt i32 %48, 1582
br i1 %49, label %62, label %50
50: ; preds = %47
%51 = load i32, i32* %7
%52 = icmp eq i32 %51, 1582
br i1 %52, label %53, label %63
53: ; preds = %50
%54 = load i32, i32* %5
%55 = icmp slt i32 %54, 10
br i1 %55, label %62, label %56
56: ; preds = %53
%57 = load i32, i32* %5
%58 = icmp eq i32 %57, 10
br i1 %58, label %59, label %63
59: ; preds = %56
%60 = load double, double* %6
%61 = fcmp olt double %60, 1.500000e+01
br i1 %61, label %62, label %63
62: ; preds = %59, %53, %47
store i32 0, i32* %9
br label %71
63: ; preds = %59, %56, %50
%64 = load i32, i32* %12
%65 = sdiv i32 %64, 100
store i32 %65, i32* %14
%66 = load i32, i32* %14
%67 = sub nsw i32 2, %66
%68 = load i32, i32* %14
%69 = sdiv i32 %68, 4
%70 = add nsw i32 %67, %69
store i32 %70, i32* %9
br label %71
71: ; preds = %63, %62
%72 = load i32, i32* %12
%73 = icmp slt i32 %72, 0
br i1 %73, label %74, label %80
74: ; preds = %71
%75 = load i32, i32* %12
%76 = sitofp i32 %75 to double
%77 = call double @llvm.fmuladd.f64(double 3.652500e+02, double %76, double -7.500000e-01)
%78 = fptosi double %77 to i64
%79 = sub nsw i64 %78, 694025
store i64 %79, i64* %13
br label %86
80: ; preds = %71
%81 = load i32, i32* %12
%82 = sitofp i32 %81 to double
%83 = fmul double 3.652500e+02, %82
%84 = fptosi double %83 to i64
%85 = sub nsw i64 %84, 694025
store i64 %85, i64* %13
br label %86
86: ; preds = %80, %74
%87 = load i32, i32* %11
%88 = add nsw i32 %87, 1
%89 = sitofp i32 %88 to double
%90 = fmul double 3.060010e+01, %89
%91 = fptosi double %90 to i32
store i32 %91, i32* %10
%92 = load i32, i32* %9
%93 = sext i32 %92 to i64
%94 = load i64, i64* %13
%95 = add nsw i64 %93, %94
%96 = load i32, i32* %10
%97 = sext i32 %96 to i64
%98 = add nsw i64 %95, %97
%99 = sitofp i64 %98 to double
%100 = load double, double* %6
%101 = fadd double %99, %100
%102 = fsub double %101, 5.000000e-01
%103 = load double*, double** %8
store double %102, double* %103
%104 = load i32, i32* %5
store i32 %104, i32* @cal_mjd.last_mn
%105 = load double, double* %6
store double %105, double* @cal_mjd.last_dy
%106 = load i32, i32* %7
store i32 %106, i32* @cal_mjd.last_yr
%107 = load double*, double** %8
%108 = load double, double* %107
store double %108, double* @cal_mjd.last_mjd
br label %109
109: ; preds = %86, %26
ret void
}

LLVM IR (-O2):
@cal_mjd.last_mjd = internal unnamed_addr global double 0.000000e+00
@cal_mjd.last_dy = internal unnamed_addr global double 0.000000e+00
@cal_mjd.last_mn = internal unnamed_addr global i32 0
@cal_mjd.last_yr = internal unnamed_addr global i32 0
define void @cal_mjd(i32 %0, double %1, i32 %2, double* nocapture %3) {
%5 = load i32, i32* @cal_mjd.last_mn
%6 = icmp eq i32 %5, %0
%7 = load i32, i32* @cal_mjd.last_yr
%8 = icmp eq i32 %7, %2
%9 = select i1 %6, i1 %8, i1 false
%10 = load double, double* @cal_mjd.last_dy
%11 = fcmp oeq double %10, %1
%12 = select i1 %9, i1 %11, i1 false
br i1 %12, label %13, label %15
13: ; preds = %4
%14 = load double, double* @cal_mjd.last_mjd
store double %14, double* %3
br label %57
15: ; preds = %4
%16 = lshr i32 %2, 31
%17 = icmp slt i32 %0, 3
%18 = add nsw i32 %0, 12
%19 = select i1 %17, i32 %18, i32 %0
%20 = sext i1 %17 to i32
%21 = add i32 %20, %2
%22 = add i32 %21, %16
%23 = icmp slt i32 %2, 1582
br i1 %23, label %38, label %24
24: ; preds = %15
%25 = icmp eq i32 %2, 1582
br i1 %25, label %26, label %32
26: ; preds = %24
%27 = icmp slt i32 %0, 10
br i1 %27, label %38, label %28
28: ; preds = %26
%29 = icmp eq i32 %0, 10
%30 = fcmp olt double %1, 1.500000e+01
%31 = and i1 %29, %30
br i1 %31, label %38, label %32
32: ; preds = %28, %24
%33 = sdiv i32 %22, -100
%34 = add nsw i32 %33, 2
%35 = sdiv i32 %22, 400
%36 = add nsw i32 %34, %35
%37 = sext i32 %36 to i64
br label %38
38: ; preds = %15, %26, %28, %32
%39 = phi i64 [ %37, %32 ], [ 0, %28 ], [ 0, %26 ], [ 0, %15 ]
%40 = icmp slt i32 %22, 0
%41 = sitofp i32 %22 to double
%42 = tail call double @llvm.fmuladd.f64(double %41, double 3.652500e+02, double -7.500000e-01)
%43 = fmul double %41, 3.652500e+02
%44 = select i1 %40, double %42, double %43
%45 = fptosi double %44 to i64
%46 = add nsw i32 %19, 1
%47 = sitofp i32 %46 to double
%48 = fmul double %47, 3.060010e+01
%49 = fptosi double %48 to i32
%50 = sext i32 %49 to i64
%51 = add nsw i64 %50, -694025
%52 = add nsw i64 %51, %39
%53 = add i64 %52, %45
%54 = sitofp i64 %53 to double
%55 = fadd double %54, %1
%56 = fadd double %55, -5.000000e-01
store double %56, double* %3
store i32 %0, i32* @cal_mjd.last_mn
store double %1, double* @cal_mjd.last_dy
store i32 %2, i32* @cal_mjd.last_yr
store double %56, double* @cal_mjd.last_mjd
br label %57
57: ; preds = %38, %13
ret void
}