void sift_up(double **c, int j, int *heap, int *heap_inv, int i) {
  int p, t;
  while (i > 0) {
    p = (i - 1) / 2;
    if (c[heap[p]][j] <= c[heap[i]][j])
      return;
    heap_inv[heap[p]] = i;
    heap_inv[heap[i]] = p;
    t = heap[p];
    heap[p] = heap[i];
    heap[i] = t;
    i = p;
  }
}


x86 assembly (clang, -O0):
<sift_up>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8d,-0x24(%rbp)
cmpl   $0x0,-0x24(%rbp)
jle    11e4 <sift_up+0xe4>
mov    -0x24(%rbp),%eax
sub    $0x1,%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,-0x28(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x18(%rbp),%rcx
movslq -0x28(%rbp),%rdx
movslq (%rcx,%rdx,4),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0xc(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x8(%rbp),%rax
mov    -0x18(%rbp),%rcx
movslq -0x24(%rbp),%rdx
movslq (%rcx,%rdx,4),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0xc(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
ucomisd %xmm1,%xmm0
jb     117b <sift_up+0x7b>
jmp    11e4 <sift_up+0xe4>
mov    -0x24(%rbp),%edx
mov    -0x20(%rbp),%rax
mov    -0x18(%rbp),%rcx
movslq -0x28(%rbp),%rsi
movslq (%rcx,%rsi,4),%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x28(%rbp),%edx
mov    -0x20(%rbp),%rax
mov    -0x18(%rbp),%rcx
movslq -0x24(%rbp),%rsi
movslq (%rcx,%rsi,4),%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x18(%rbp),%rax
movslq -0x28(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x2c(%rbp)
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
mov    (%rax,%rcx,4),%edx
mov    -0x18(%rbp),%rax
movslq -0x28(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x2c(%rbp),%edx
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x28(%rbp),%eax
mov    %eax,-0x24(%rbp)
jmp    1117 <sift_up+0x17>
pop    %rbp
ret

x86 assembly (clang, -O2):
<sift_up>:
push   %r14
push   %rbx
test   %r8d,%r8d
jle    115c <sift_up+0x5c>
movslq %esi,%r9
nopl   0x0(%rax,%rax,1)
lea    -0x1(%r8),%r10d
mov    %r10d,%eax
shr    %eax
movslq (%rdx,%rax,4),%r11
mov    (%rdi,%r11,8),%r14
mov    %r8d,%esi
movslq (%rdx,%rsi,4),%rbx
mov    (%rdi,%rbx,8),%rbx
movsd  (%rbx,%r9,8),%xmm0
ucomisd (%r14,%r9,8),%xmm0
jae    115c <sift_up+0x5c>
mov    %r8d,(%rcx,%r11,4)
movslq (%rdx,%rsi,4),%rbx
mov    %eax,(%rcx,%rbx,4)
mov    (%rdx,%rax,4),%r8d
mov    (%rdx,%rsi,4),%ebx
mov    %ebx,(%rdx,%rax,4)
mov    %r8d,(%rdx,%rsi,4)
mov    %eax,%r8d
cmp    $0x1,%r10d
ja     1110 <sift_up+0x10>
pop    %rbx
pop    %r14
ret

LLVM IR (-O0):
define void @sift_up(double** %0, i32 %1, i32* %2, i32* %3, i32 %4) {
%6 = alloca double**
%7 = alloca i32
%8 = alloca i32*
%9 = alloca i32*
%10 = alloca i32
%11 = alloca i32
%12 = alloca i32
store double** %0, double*** %6
store i32 %1, i32* %7
store i32* %2, i32** %8
store i32* %3, i32** %9
store i32 %4, i32* %10
br label %13
13: ; preds = %48, %5
%14 = load i32, i32* %10
%15 = icmp sgt i32 %14, 0
br i1 %15, label %16, label %87
16: ; preds = %13
%17 = load i32, i32* %10
%18 = sub nsw i32 %17, 1
%19 = sdiv i32 %18, 2
store i32 %19, i32* %11
%20 = load double**, double*** %6
%21 = load i32*, i32** %8
%22 = load i32, i32* %11
%23 = sext i32 %22 to i64
%24 = getelementptr inbounds i32, i32* %21, i64 %23
%25 = load i32, i32* %24
%26 = sext i32 %25 to i64
%27 = getelementptr inbounds double*, double** %20, i64 %26
%28 = load double*, double** %27
%29 = load i32, i32* %7
%30 = sext i32 %29 to i64
%31 = getelementptr inbounds double, double* %28, i64 %30
%32 = load double, double* %31
%33 = load double**, double*** %6
%34 = load i32*, i32** %8
%35 = load i32, i32* %10
%36 = sext i32 %35 to i64
%37 = getelementptr inbounds i32, i32* %34, i64 %36
%38 = load i32, i32* %37
%39 = sext i32 %38 to i64
%40 = getelementptr inbounds double*, double** %33, i64 %39
%41 = load double*, double** %40
%42 = load i32, i32* %7
%43 = sext i32 %42 to i64
%44 = getelementptr inbounds double, double* %41, i64 %43
%45 = load double, double* %44
%46 = fcmp ole double %32, %45
br i1 %46, label %47, label %48
47: ; preds = %16
br label %87
48: ; preds = %16
%49 = load i32, i32* %10
%50 = load i32*, i32** %9
%51 = load i32*, i32** %8
%52 = load i32, i32* %11
%53 = sext i32 %52 to i64
%54 = getelementptr inbounds i32, i32* %51, i64 %53
%55 = load i32, i32* %54
%56 = sext i32 %55 to i64
%57 = getelementptr inbounds i32, i32* %50, i64 %56
store i32 %49, i32* %57
%58 = load i32, i32* %11
%59 = load i32*, i32** %9
%60 = load i32*, i32** %8
%61 = load i32, i32* %10
%62 = sext i32 %61 to i64
%63 = getelementptr inbounds i32, i32* %60, i64 %62
%64 = load i32, i32* %63
%65 = sext i32 %64 to i64
%66 = getelementptr inbounds i32, i32* %59, i64 %65
store i32 %58, i32* %66
%67 = load i32*, i32** %8
%68 = load i32, i32* %11
%69 = sext i32 %68 to i64
%70 = getelementptr inbounds i32, i32* %67, i64 %69
%71 = load i32, i32* %70
store i32 %71, i32* %12
%72 = load i32*, i32** %8
%73 = load i32, i32* %10
%74 = sext i32 %73 to i64
%75 = getelementptr inbounds i32, i32* %72, i64 %74
%76 = load i32, i32* %75
%77 = load i32*, i32** %8
%78 = load i32, i32* %11
%79 = sext i32 %78 to i64
%80 = getelementptr inbounds i32, i32* %77, i64 %79
store i32 %76, i32* %80
%81 = load i32, i32* %12
%82 = load i32*, i32** %8
%83 = load i32, i32* %10
%84 = sext i32 %83 to i64
%85 = getelementptr inbounds i32, i32* %82, i64 %84
store i32 %81, i32* %85
%86 = load i32, i32* %11
store i32 %86, i32* %10
br label %13
87: ; preds = %47, %13
ret void
}

LLVM IR (-O2):
define void @sift_up(double** nocapture %0, i32 %1, i32* nocapture %2, i32* nocapture writeonly %3, i32 %4) {
%6 = icmp sgt i32 %4, 0
br i1 %6, label %7, label %38
7: ; preds = %5
%8 = sext i32 %1 to i64
br label %9
9: ; preds = %7, %30
%10 = phi i32 [ %4, %7 ], [ %12, %30 ]
%11 = add nsw i32 %10, -1
%12 = lshr i32 %11, 1
%13 = zext i32 %12 to i64
%14 = getelementptr inbounds i32, i32* %2, i64 %13
%15 = load i32, i32* %14
%16 = sext i32 %15 to i64
%17 = getelementptr inbounds double*, double** %0, i64 %16
%18 = load double*, double** %17
%19 = getelementptr inbounds double, double* %18, i64 %8
%20 = load double, double* %19
%21 = zext i32 %10 to i64
%22 = getelementptr inbounds i32, i32* %2, i64 %21
%23 = load i32, i32* %22
%24 = sext i32 %23 to i64
%25 = getelementptr inbounds double*, double** %0, i64 %24
%26 = load double*, double** %25
%27 = getelementptr inbounds double, double* %26, i64 %8
%28 = load double, double* %27
%29 = fcmp ugt double %20, %28
br i1 %29, label %30, label %38
30: ; preds = %9
%31 = getelementptr inbounds i32, i32* %3, i64 %16
store i32 %10, i32* %31
%32 = load i32, i32* %22
%33 = sext i32 %32 to i64
%34 = getelementptr inbounds i32, i32* %3, i64 %33
store i32 %12, i32* %34
%35 = load i32, i32* %14
%36 = load i32, i32* %22
store i32 %36, i32* %14
store i32 %35, i32* %22
%37 = icmp ult i32 %11, 2
br i1 %37, label %38, label %9
38: ; preds = %9, %30, %5
ret void
}