void shellsort(a, b, n) double *a;
short *b;
short n;
{
  short gap, i, j, itemp;
  double rtemp;
  gap = n / 2;
  while (gap > 0) {
    for (i = gap + 1; i <= n; i++) {
      j = i - gap;
      while (j > 0) {
        if (a[j - 1] > a[j + gap - 1]) {
          rtemp = a[j - 1];
          a[j - 1] = a[j + gap - 1];
          a[j + gap - 1] = rtemp;
          itemp = b[j - 1];
          b[j - 1] = b[j + gap - 1];
          b[j + gap - 1] = itemp;
        }
        j -= gap;
      }
    }
    gap /= 2;
  }
}


x86 assembly (clang, -O0):
<shellsort>:
push   %rbp
mov    %rsp,%rbp
mov    %dx,%ax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %ax,-0x12(%rbp)
movswl -0x12(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %ax,-0x14(%rbp)
movswl -0x14(%rbp),%eax
cmp    $0x0,%eax
jle    1294 <shellsort+0x194>
movswl -0x14(%rbp),%eax
add    $0x1,%eax
mov    %ax,-0x16(%rbp)
movswl -0x16(%rbp),%eax
movswl -0x12(%rbp),%ecx
cmp    %ecx,%eax
jg     127f <shellsort+0x17f>
movswl -0x16(%rbp),%eax
movswl -0x14(%rbp),%ecx
sub    %ecx,%eax
mov    %ax,-0x18(%rbp)
movswl -0x18(%rbp),%eax
cmp    $0x0,%eax
jle    1269 <shellsort+0x169>
mov    -0x8(%rbp),%rax
movswl -0x18(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
movswl -0x18(%rbp),%ecx
movswl -0x14(%rbp),%edx
add    %edx,%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
ucomisd (%rax,%rcx,8),%xmm0
jbe    1256 <shellsort+0x156>
mov    -0x8(%rbp),%rax
movswl -0x18(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0x8(%rbp),%rax
movswl -0x18(%rbp),%ecx
movswl -0x14(%rbp),%edx
add    %edx,%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
movswl -0x18(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x28(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movswl -0x18(%rbp),%ecx
movswl -0x14(%rbp),%edx
add    %edx,%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%rax
movswl -0x18(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,2),%ax
mov    %ax,-0x1a(%rbp)
mov    -0x10(%rbp),%rax
movswl -0x18(%rbp),%ecx
movswl -0x14(%rbp),%edx
add    %edx,%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,2),%dx
mov    -0x10(%rbp),%rax
movswl -0x18(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mov    %dx,(%rax,%rcx,2)
mov    -0x1a(%rbp),%dx
mov    -0x10(%rbp),%rax
movswl -0x18(%rbp),%ecx
movswl -0x14(%rbp),%esi
add    %esi,%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mov    %dx,(%rax,%rcx,2)
movswl -0x14(%rbp),%ecx
movswl -0x18(%rbp),%eax
sub    %ecx,%eax
mov    %ax,-0x18(%rbp)
jmp    1159 <shellsort+0x59>
jmp    126e <shellsort+0x16e>
mov    -0x16(%rbp),%ax
add    $0x1,%ax
mov    %ax,-0x16(%rbp)
jmp    113b <shellsort+0x3b>
movswl -0x14(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %ax,-0x14(%rbp)
jmp    1123 <shellsort+0x23>
pop    %rbp
ret

x86 assembly (clang, -O2):
<shellsort>:
push   %rbp
push   %rbx
mov    %edx,%eax
shl    $0x10,%eax
cmp    $0x10001,%eax
jl     1190 <shellsort+0x90>
movswl %dx,%r9d
mov    %r9d,%r11d
jmp    1126 <shellsort+0x26>
nopl   0x0(%rax,%rax,1)
cmp    $0x4,%r8d
jl     1190 <shellsort+0x90>
mov    %r11d,%r8d
shr    $0x1f,%r11d
add    %r8d,%r11d
sar    %r11d
lea    0x1(%r11),%r10d
jmp    1144 <shellsort+0x44>
nopl   0x0(%rax)
add    $0x1,%r10d
movswl %r10w,%eax
cmp    %eax,%r9d
jl     1120 <shellsort+0x20>
mov    %r10d,%edx
sub    %r11d,%edx
test   %dx,%dx
jle    1140 <shellsort+0x40>
movzwl %dx,%ecx
lea    -0x1(%rcx),%eax
movsd  (%rdi,%rax,8),%xmm0
add    %r11d,%ecx
add    $0xffffffff,%ecx
movsd  (%rdi,%rcx,8),%xmm1
ucomisd %xmm1,%xmm0
jbe    1150 <shellsort+0x50>
movsd  %xmm1,(%rdi,%rax,8)
movsd  %xmm0,(%rdi,%rcx,8)
movzwl (%rsi,%rax,2),%ebx
movzwl (%rsi,%rcx,2),%ebp
mov    %bp,(%rsi,%rax,2)
mov    %bx,(%rsi,%rcx,2)
jmp    1150 <shellsort+0x50>
pop    %rbx
pop    %rbp
ret

LLVM IR (-O0):
define void @shellsort(double* %0, i16* %1, i32 %2) {
%4 = alloca double*
%5 = alloca i16*
%6 = alloca i16
%7 = alloca i16
%8 = alloca i16
%9 = alloca i16
%10 = alloca i16
%11 = alloca double
%12 = trunc i32 %2 to i16
store double* %0, double** %4
store i16* %1, i16** %5
store i16 %12, i16* %6
%13 = load i16, i16* %6
%14 = sext i16 %13 to i32
%15 = sdiv i32 %14, 2
%16 = trunc i32 %15 to i16
store i16 %16, i16* %7
br label %17
17: ; preds = %140, %3
%18 = load i16, i16* %7
%19 = sext i16 %18 to i32
%20 = icmp sgt i32 %19, 0
br i1 %20, label %21, label %145
21: ; preds = %17
%22 = load i16, i16* %7
%23 = sext i16 %22 to i32
%24 = add nsw i32 %23, 1
%25 = trunc i32 %24 to i16
store i16 %25, i16* %8
br label %26
26: ; preds = %137, %21
%27 = load i16, i16* %8
%28 = sext i16 %27 to i32
%29 = load i16, i16* %6
%30 = sext i16 %29 to i32
%31 = icmp sle i32 %28, %30
br i1 %31, label %32, label %140
32: ; preds = %26
%33 = load i16, i16* %8
%34 = sext i16 %33 to i32
%35 = load i16, i16* %7
%36 = sext i16 %35 to i32
%37 = sub nsw i32 %34, %36
%38 = trunc i32 %37 to i16
store i16 %38, i16* %9
br label %39
39: ; preds = %129, %32
%40 = load i16, i16* %9
%41 = sext i16 %40 to i32
%42 = icmp sgt i32 %41, 0
br i1 %42, label %43, label %136
43: ; preds = %39
%44 = load double*, double** %4
%45 = load i16, i16* %9
%46 = sext i16 %45 to i32
%47 = sub nsw i32 %46, 1
%48 = sext i32 %47 to i64
%49 = getelementptr inbounds double, double* %44, i64 %48
%50 = load double, double* %49
%51 = load double*, double** %4
%52 = load i16, i16* %9
%53 = sext i16 %52 to i32
%54 = load i16, i16* %7
%55 = sext i16 %54 to i32
%56 = add nsw i32 %53, %55
%57 = sub nsw i32 %56, 1
%58 = sext i32 %57 to i64
%59 = getelementptr inbounds double, double* %51, i64 %58
%60 = load double, double* %59
%61 = fcmp ogt double %50, %60
br i1 %61, label %62, label %129
62: ; preds = %43
%63 = load double*, double** %4
%64 = load i16, i16* %9
%65 = sext i16 %64 to i32
%66 = sub nsw i32 %65, 1
%67 = sext i32 %66 to i64
%68 = getelementptr inbounds double, double* %63, i64 %67
%69 = load double, double* %68
store double %69, double* %11
%70 = load double*, double** %4
%71 = load i16, i16* %9
%72 = sext i16 %71 to i32
%73 = load i16, i16* %7
%74 = sext i16 %73 to i32
%75 = add nsw i32 %72, %74
%76 = sub nsw i32 %75, 1
%77 = sext i32 %76 to i64
%78 = getelementptr inbounds double, double* %70, i64 %77
%79 = load double, double* %78
%80 = load double*, double** %4
%81 = load i16, i16* %9
%82 = sext i16 %81 to i32
%83 = sub nsw i32 %82, 1
%84 = sext i32 %83 to i64
%85 = getelementptr inbounds double, double* %80, i64 %84
store double %79, double* %85
%86 = load double, double* %11
%87 = load double*, double** %4
%88 = load i16, i16* %9
%89 = sext i16 %88 to i32
%90 = load i16, i16* %7
%91 = sext i16 %90 to i32
%92 = add nsw i32 %89, %91
%93 = sub nsw i32 %92, 1
%94 = sext i32 %93 to i64
%95 = getelementptr inbounds double, double* %87, i64 %94
store double %86, double* %95
%96 = load i16*, i16** %5
%97 = load i16, i16* %9
%98 = sext i16 %97 to i32
%99 = sub nsw i32 %98, 1
%100 = sext i32 %99 to i64
%101 = getelementptr inbounds i16, i16* %96, i64 %100
%102 = load i16, i16* %101
store i16 %102, i16* %10
%103 = load i16*, i16** %5
%104 = load i16, i16* %9
%105 = sext i16 %104 to i32
%106 = load i16, i16* %7
%107 = sext i16 %106 to i32
%108 = add nsw i32 %105, %107
%109 = sub nsw i32 %108, 1
%110 = sext i32 %109 to i64
%111 = getelementptr inbounds i16, i16* %103, i64 %110
%112 = load i16, i16* %111
%113 = load i16*, i16** %5
%114 = load i16, i16* %9
%115 = sext i16 %114 to i32
%116 = sub nsw i32 %115, 1
%117 = sext i32 %116 to i64
%118 = getelementptr inbounds i16, i16* %113, i64 %117
store i16 %112, i16* %118
%119 = load i16, i16* %10
%120 = load i16*, i16** %5
%121 = load i16, i16* %9
%122 = sext i16 %121 to i32
%123 = load i16, i16* %7
%124 = sext i16 %123 to i32
%125 = add nsw i32 %122, %124
%126 = sub nsw i32 %125, 1
%127 = sext i32 %126 to i64
%128 = getelementptr inbounds i16, i16* %120, i64 %127
store i16 %119, i16* %128
br label %129
129: ; preds = %62, %43
%130 = load i16, i16* %7
%131 = sext i16 %130 to i32
%132 = load i16, i16* %9
%133 = sext i16 %132 to i32
%134 = sub nsw i32 %133, %131
%135 = trunc i32 %134 to i16
store i16 %135, i16* %9
br label %39
136: ; preds = %39
br label %137
137: ; preds = %136
%138 = load i16, i16* %8
%139 = add i16 %138, 1
store i16 %139, i16* %8
br label %26
140: ; preds = %26
%141 = load i16, i16* %7
%142 = sext i16 %141 to i32
%143 = sdiv i32 %142, 2
%144 = trunc i32 %143 to i16
store i16 %144, i16* %7
br label %17
145: ; preds = %17
ret void
}

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