void fSort(double *array, int *index, int nRef) {
  int nInt, ii, i, j, nf;
  double tf;
  nInt = nRef;
  do {
    nInt = nInt / 2;
    if ((nInt % 2) == 0) {
      nInt = nInt - 1;
    }
    for (ii = 0; ii < nRef - nInt; ii++) {
      i = ii;
      j = i + nInt;
      if (array[i] > array[j]) {
        tf = array[j];
        nf = index[j];
        do {
          array[j] = array[i];
          index[j] = index[i];
          j = i;
          i = i - nInt;
        } while ((i >= 0) && (array[i] > tf));
        array[j] = tf;
        index[j] = nf;
      }
    }
  } while (nInt > 1);
}


x86 assembly (clang, -O0):
<fSort>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    -0x14(%rbp),%eax
mov    %eax,-0x18(%rbp)
mov    -0x18(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
mov    %eax,-0x18(%rbp)
mov    -0x18(%rbp),%eax
mov    $0x2,%ecx
cltd
idiv   %ecx
cmp    $0x0,%edx
jne    1140 <fSort+0x40>
mov    -0x18(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x18(%rbp)
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
mov    -0x14(%rbp),%ecx
sub    -0x18(%rbp),%ecx
cmp    %ecx,%eax
jge    124b <fSort+0x14b>
mov    -0x1c(%rbp),%eax
mov    %eax,-0x20(%rbp)
mov    -0x20(%rbp),%eax
add    -0x18(%rbp),%eax
mov    %eax,-0x24(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
ucomisd (%rax,%rcx,8),%xmm0
jbe    1238 <fSort+0x138>
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x28(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
mov    (%rax,%rcx,4),%edx
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x20(%rbp),%eax
mov    %eax,-0x24(%rbp)
mov    -0x20(%rbp),%eax
sub    -0x18(%rbp),%eax
mov    %eax,-0x20(%rbp)
xor    %eax,%eax
cmpl   $0x0,-0x20(%rbp)
mov    %al,-0x31(%rbp)
jl     120d <fSort+0x10d>
mov    -0x8(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
ucomisd -0x30(%rbp),%xmm0
seta   %al
mov    %al,-0x31(%rbp)
mov    -0x31(%rbp),%al
test   $0x1,%al
jne    11a7 <fSort+0xa7>
movsd  -0x30(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%edx
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
jmp    123d <fSort+0x13d>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1147 <fSort+0x47>
jmp    1250 <fSort+0x150>
cmpl   $0x1,-0x18(%rbp)
jg     1115 <fSort+0x15>
pop    %rbp
ret

x86 assembly (clang, -O2):
<fSort>:
push   %rbp
push   %r15
push   %r14
push   %rbx
mov    %edx,%r15d
jmp    111a <fSort+0x1a>
nopl   0x0(%rax,%rax,1)
cmp    $0x1,%r15d
jle    11c3 <fSort+0xc3>
mov    %r15d,%ecx
shr    $0x1f,%ecx
add    %r15d,%ecx
sar    %ecx
mov    %ecx,%r10d
and    $0x1,%r10d
lea    (%rcx,%r10,1),%r15d
add    $0xffffffff,%r15d
mov    %edx,%ebx
sub    %r15d,%ebx
test   %ebx,%ebx
jle    1110 <fSort+0x10>
movslq %r15d,%r8
mov    %ebx,%r9d
add    %r10d,%ecx
mov    $0x1,%r10d
sub    %ecx,%r10d
xor    %r14d,%r14d
jmp    1172 <fSort+0x72>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movsd  %xmm0,(%rdi,%rcx,8)
mov    %r11d,(%rsi,%rcx,4)
add    $0x1,%r14
cmp    %r9,%r14
je     1110 <fSort+0x10>
lea    (%r14,%r8,1),%rcx
movsd  (%rdi,%r14,8),%xmm1
movsd  (%rdi,%rcx,8),%xmm0
ucomisd %xmm0,%xmm1
jbe    1169 <fSort+0x69>
mov    (%rsi,%rcx,4),%r11d
mov    %r14d,%ebx
xchg   %ax,%ax
lea    (%r15,%rbx,1),%ebp
movslq %ebx,%rcx
movsd  (%rdi,%rcx,8),%xmm1
movslq %ebp,%rbp
movsd  %xmm1,(%rdi,%rbp,8)
mov    (%rsi,%rcx,4),%eax
mov    %eax,(%rsi,%rbp,4)
mov    %r10d,%ebp
add    %ebx,%ebp
js     1160 <fSort+0x60>
sub    %r15d,%ebx
mov    %ebp,%eax
movsd  (%rdi,%rax,8),%xmm1
ucomisd %xmm0,%xmm1
ja     1190 <fSort+0x90>
jmp    1160 <fSort+0x60>
pop    %rbx
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @fSort(double* %0, i32* %1, i32 %2) {
%4 = alloca double*
%5 = alloca i32*
%6 = alloca i32
%7 = alloca i32
%8 = alloca i32
%9 = alloca i32
%10 = alloca i32
%11 = alloca i32
%12 = alloca double
store double* %0, double** %4
store i32* %1, i32** %5
store i32 %2, i32* %6
%13 = load i32, i32* %6
store i32 %13, i32* %7
br label %14
14: ; preds = %109, %3
%15 = load i32, i32* %7
%16 = sdiv i32 %15, 2
store i32 %16, i32* %7
%17 = load i32, i32* %7
%18 = srem i32 %17, 2
%19 = icmp eq i32 %18, 0
br i1 %19, label %20, label %23
20: ; preds = %14
%21 = load i32, i32* %7
%22 = sub nsw i32 %21, 1
store i32 %22, i32* %7
br label %23
23: ; preds = %20, %14
store i32 0, i32* %8
br label %24
24: ; preds = %105, %23
%25 = load i32, i32* %8
%26 = load i32, i32* %6
%27 = load i32, i32* %7
%28 = sub nsw i32 %26, %27
%29 = icmp slt i32 %25, %28
br i1 %29, label %30, label %108
30: ; preds = %24
%31 = load i32, i32* %8
store i32 %31, i32* %9
%32 = load i32, i32* %9
%33 = load i32, i32* %7
%34 = add nsw i32 %32, %33
store i32 %34, i32* %10
%35 = load double*, double** %4
%36 = load i32, i32* %9
%37 = sext i32 %36 to i64
%38 = getelementptr inbounds double, double* %35, i64 %37
%39 = load double, double* %38
%40 = load double*, double** %4
%41 = load i32, i32* %10
%42 = sext i32 %41 to i64
%43 = getelementptr inbounds double, double* %40, i64 %42
%44 = load double, double* %43
%45 = fcmp ogt double %39, %44
br i1 %45, label %46, label %104
46: ; preds = %30
%47 = load double*, double** %4
%48 = load i32, i32* %10
%49 = sext i32 %48 to i64
%50 = getelementptr inbounds double, double* %47, i64 %49
%51 = load double, double* %50
store double %51, double* %12
%52 = load i32*, i32** %5
%53 = load i32, i32* %10
%54 = sext i32 %53 to i64
%55 = getelementptr inbounds i32, i32* %52, i64 %54
%56 = load i32, i32* %55
store i32 %56, i32* %11
br label %57
57: ; preds = %91, %46
%58 = load double*, double** %4
%59 = load i32, i32* %9
%60 = sext i32 %59 to i64
%61 = getelementptr inbounds double, double* %58, i64 %60
%62 = load double, double* %61
%63 = load double*, double** %4
%64 = load i32, i32* %10
%65 = sext i32 %64 to i64
%66 = getelementptr inbounds double, double* %63, i64 %65
store double %62, double* %66
%67 = load i32*, i32** %5
%68 = load i32, i32* %9
%69 = sext i32 %68 to i64
%70 = getelementptr inbounds i32, i32* %67, i64 %69
%71 = load i32, i32* %70
%72 = load i32*, i32** %5
%73 = load i32, i32* %10
%74 = sext i32 %73 to i64
%75 = getelementptr inbounds i32, i32* %72, i64 %74
store i32 %71, i32* %75
%76 = load i32, i32* %9
store i32 %76, i32* %10
%77 = load i32, i32* %9
%78 = load i32, i32* %7
%79 = sub nsw i32 %77, %78
store i32 %79, i32* %9
br label %80
80: ; preds = %57
%81 = load i32, i32* %9
%82 = icmp sge i32 %81, 0
br i1 %82, label %83, label %91
83: ; preds = %80
%84 = load double*, double** %4
%85 = load i32, i32* %9
%86 = sext i32 %85 to i64
%87 = getelementptr inbounds double, double* %84, i64 %86
%88 = load double, double* %87
%89 = load double, double* %12
%90 = fcmp ogt double %88, %89
br label %91
91: ; preds = %83, %80
%92 = phi i1 [ false, %80 ], [ %90, %83 ]
br i1 %92, label %57, label %93
93: ; preds = %91
%94 = load double, double* %12
%95 = load double*, double** %4
%96 = load i32, i32* %10
%97 = sext i32 %96 to i64
%98 = getelementptr inbounds double, double* %95, i64 %97
store double %94, double* %98
%99 = load i32, i32* %11
%100 = load i32*, i32** %5
%101 = load i32, i32* %10
%102 = sext i32 %101 to i64
%103 = getelementptr inbounds i32, i32* %100, i64 %102
store i32 %99, i32* %103
br label %104
104: ; preds = %93, %30
br label %105
105: ; preds = %104
%106 = load i32, i32* %8
%107 = add nsw i32 %106, 1
store i32 %107, i32* %8
br label %24
108: ; preds = %24
br label %109
109: ; preds = %108
%110 = load i32, i32* %7
%111 = icmp sgt i32 %110, 1
br i1 %111, label %14, label %112
112: ; preds = %109
ret void
}

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