void lsort(long *x, int k) {
  int i, j, h;
  long iw;
  j = k / 3;
  h = 1;
  do
    h = 3 * h + 1;
  while (h < j);
  do {
    for (i = h; i < k; ++i) {
      iw = x[i];
      for (j = i; x[j - h] > iw;) {
        x[j] = x[j - h];
        if ((j -= h) < h)
          break;
      }
      x[j] = iw;
    }
    h /= 3;
  } while (h > 0);
}


x86 assembly (clang, -O0):
<lsort>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    -0xc(%rbp),%eax
mov    $0x3,%ecx
cltd
idiv   %ecx
mov    %eax,-0x14(%rbp)
movl   $0x1,-0x18(%rbp)
imul   $0x3,-0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
mov    -0x18(%rbp),%eax
cmp    -0x14(%rbp),%eax
jl     1120 <lsort+0x20>
jmp    113b <lsort+0x3b>
mov    -0x18(%rbp),%eax
mov    %eax,-0x10(%rbp)
mov    -0x10(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    11d7 <lsort+0xd7>
mov    -0x8(%rbp),%rax
movslq -0x10(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    %rax,-0x20(%rbp)
mov    -0x10(%rbp),%eax
mov    %eax,-0x14(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x14(%rbp),%ecx
sub    -0x18(%rbp),%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,8),%rax
cmp    -0x20(%rbp),%rax
jle    11b9 <lsort+0xb9>
mov    -0x8(%rbp),%rax
mov    -0x14(%rbp),%ecx
sub    -0x18(%rbp),%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,8),%rdx
mov    -0x8(%rbp),%rax
movslq -0x14(%rbp),%rcx
mov    %rdx,(%rax,%rcx,8)
mov    -0x18(%rbp),%ecx
mov    -0x14(%rbp),%eax
sub    %ecx,%eax
mov    %eax,-0x14(%rbp)
cmp    -0x18(%rbp),%eax
jge    11b4 <lsort+0xb4>
jmp    11b9 <lsort+0xb9>
jmp    1163 <lsort+0x63>
mov    -0x20(%rbp),%rdx
mov    -0x8(%rbp),%rax
movslq -0x14(%rbp),%rcx
mov    %rdx,(%rax,%rcx,8)
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    1141 <lsort+0x41>
mov    -0x18(%rbp),%eax
mov    $0x3,%ecx
cltd
idiv   %ecx
mov    %eax,-0x18(%rbp)
cmpl   $0x0,-0x18(%rbp)
jg     113b <lsort+0x3b>
pop    %rbp
ret

x86 assembly (clang, -O2):
<lsort>:
push   %r15
push   %r14
push   %rbx
movslq %esi,%r9
imul   $0x55555556,%r9,%rax
mov    %rax,%rcx
shr    $0x3f,%rcx
shr    $0x20,%rax
add    %ecx,%eax
mov    $0x1,%r14d
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
lea    (%r14,%r14,2),%r14d
add    $0x1,%r14d
cmp    %eax,%r14d
jl     1130 <lsort+0x30>
jmp    115d <lsort+0x5d>
nop
imul   $0x55555556,%r8,%rax
mov    %rax,%rcx
shr    $0x3f,%rcx
shr    $0x20,%rax
add    %ecx,%eax
cmp    $0x2,%r14d
mov    %eax,%r14d
jle    11ae <lsort+0xae>
movslq %r14d,%r8
cmp    %esi,%r14d
jge    1140 <lsort+0x40>
xor    %r10d,%r10d
mov    %r8,%r11
jmp    1184 <lsort+0x84>
nopl   (%rax)
mov    %rbx,%rdx
mov    %r15,(%rdi,%rdx,8)
add    $0x1,%r11
add    $0x1,%r10d
cmp    %r9,%r11
je     1140 <lsort+0x40>
mov    (%rdi,%r11,8),%r15
mov    %r10d,%eax
nopl   0x0(%rax,%rax,1)
lea    (%r14,%rax,1),%ebx
movslq %eax,%rdx
mov    (%rdi,%rdx,8),%rcx
movslq %ebx,%rbx
cmp    %r15,%rcx
jle    1170 <lsort+0x70>
mov    %rcx,(%rdi,%rbx,8)
sub    %r14d,%eax
jge    1190 <lsort+0x90>
jmp    1173 <lsort+0x73>
pop    %rbx
pop    %r14
pop    %r15
ret

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

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