void shell_sysrem(n, arr) float arr[];
int n;
{
  int nn, m, j, i, lognb2;
  float t;
  lognb2 = (log((double)n) * 1.442695022 + 1.0e-5);
  m = n;
  for (nn = 1; nn <= lognb2; nn++) {
    m >>= 1;
    for (j = m + 1; j <= n; j++) {
      i = j - m;
      t = arr[j];
      while (i >= 1 && arr[i] > t) {
        arr[i + m] = arr[i];
        i -= m;
      }
      arr[i + m] = t;
    }
  }
}


x86 assembly (clang, -O0):
<shell_sysrem>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%rsp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
cvtsi2sdl -0x4(%rbp),%xmm0
call   1030 <log@plt>
movsd  0xecf(%rip),%xmm1
mulsd  %xmm1,%xmm0
movsd  0xecb(%rip),%xmm1
addsd  %xmm1,%xmm0
cvttsd2si %xmm0,%eax
mov    %eax,-0x24(%rbp)
mov    -0x4(%rbp),%eax
mov    %eax,-0x18(%rbp)
movl   $0x1,-0x14(%rbp)
mov    -0x14(%rbp),%eax
cmp    -0x24(%rbp),%eax
jg     1237 <shell_sysrem+0x127>
mov    -0x18(%rbp),%eax
sar    $0x1,%eax
mov    %eax,-0x18(%rbp)
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jg     1224 <shell_sysrem+0x114>
mov    -0x1c(%rbp),%eax
sub    -0x18(%rbp),%eax
mov    %eax,-0x20(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x28(%rbp)
xor    %eax,%eax
cmpl   $0x1,-0x20(%rbp)
mov    %al,-0x29(%rbp)
jl     11c0 <shell_sysrem+0xb0>
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
ucomiss -0x28(%rbp),%xmm0
seta   %al
mov    %al,-0x29(%rbp)
mov    -0x29(%rbp),%al
test   $0x1,%al
jne    11d0 <shell_sysrem+0xc0>
jmp    11ff <shell_sysrem+0xef>
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%ecx
add    -0x18(%rbp),%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x18(%rbp),%ecx
mov    -0x20(%rbp),%eax
sub    %ecx,%eax
mov    %eax,-0x20(%rbp)
jmp    119a <shell_sysrem+0x8a>
movss  -0x28(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%ecx
add    -0x18(%rbp),%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1173 <shell_sysrem+0x63>
jmp    1229 <shell_sysrem+0x119>
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    1155 <shell_sysrem+0x45>
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<shell_sysrem>:
push   %rbp
push   %r14
push   %rbx
mov    %rsi,%rbx
mov    %edi,%r14d
cvtsi2sd %edi,%xmm0
call   1030 <log@plt>
mulsd  0xed5(%rip),%xmm0
addsd  0xed5(%rip),%xmm0
cvttsd2si %xmm0,%r8d
test   %r8d,%r8d
jle    11b2 <shell_sysrem+0xa2>
movslq %r14d,%rcx
mov    $0x1,%r9d
mov    %r14d,%esi
jmp    115c <shell_sysrem+0x4c>
nopl   0x0(%rax,%rax,1)
lea    0x1(%r9),%eax
cmp    %r8d,%r9d
mov    %eax,%r9d
je     11b2 <shell_sysrem+0xa2>
sar    %esi
cmp    %r14d,%esi
jge    1150 <shell_sysrem+0x40>
movslq %esi,%rbp
jmp    117c <shell_sysrem+0x6c>
nopl   0x0(%rax,%rax,1)
cltq
movss  %xmm0,(%rbx,%rax,4)
cmp    %rcx,%rbp
je     1150 <shell_sysrem+0x40>
movss  0x4(%rbx,%rbp,4),%xmm0
add    $0x1,%rbp
mov    %ebp,%edi
sub    %esi,%edi
nopw   0x0(%rax,%rax,1)
lea    (%rsi,%rdi,1),%eax
mov    %edi,%edx
movss  (%rbx,%rdx,4),%xmm1
ucomiss %xmm0,%xmm1
jbe    1170 <shell_sysrem+0x60>
cltq
movss  %xmm1,(%rbx,%rax,4)
sub    %esi,%edi
test   %edi,%edi
jg     1190 <shell_sysrem+0x80>
add    %esi,%edi
mov    %edi,%eax
jmp    1170 <shell_sysrem+0x60>
pop    %rbx
pop    %r14
pop    %rbp
ret

LLVM IR (-O0):
define void @shell_sysrem(i32 %0, float* %1) {
%3 = alloca i32
%4 = alloca float*
%5 = alloca i32
%6 = alloca i32
%7 = alloca i32
%8 = alloca i32
%9 = alloca i32
%10 = alloca float
store i32 %0, i32* %3
store float* %1, float** %4
%11 = load i32, i32* %3
%12 = sitofp i32 %11 to double
%13 = call double @log(double %12)
%14 = call double @llvm.fmuladd.f64(double %13, double 0x3FF71547601978D6, double 1.000000e-05)
%15 = fptosi double %14 to i32
store i32 %15, i32* %9
%16 = load i32, i32* %3
store i32 %16, i32* %6
store i32 1, i32* %5
br label %17
17: ; preds = %79, %2
%18 = load i32, i32* %5
%19 = load i32, i32* %9
%20 = icmp sle i32 %18, %19
br i1 %20, label %21, label %82
21: ; preds = %17
%22 = load i32, i32* %6
%23 = ashr i32 %22, 1
store i32 %23, i32* %6
%24 = load i32, i32* %6
%25 = add nsw i32 %24, 1
store i32 %25, i32* %7
br label %26
26: ; preds = %75, %21
%27 = load i32, i32* %7
%28 = load i32, i32* %3
%29 = icmp sle i32 %27, %28
br i1 %29, label %30, label %78
30: ; preds = %26
%31 = load i32, i32* %7
%32 = load i32, i32* %6
%33 = sub nsw i32 %31, %32
store i32 %33, i32* %8
%34 = load float*, float** %4
%35 = load i32, i32* %7
%36 = sext i32 %35 to i64
%37 = getelementptr inbounds float, float* %34, i64 %36
%38 = load float, float* %37
store float %38, float* %10
br label %39
39: ; preds = %52, %30
%40 = load i32, i32* %8
%41 = icmp sge i32 %40, 1
br i1 %41, label %42, label %50
42: ; preds = %39
%43 = load float*, float** %4
%44 = load i32, i32* %8
%45 = sext i32 %44 to i64
%46 = getelementptr inbounds float, float* %43, i64 %45
%47 = load float, float* %46
%48 = load float, float* %10
%49 = fcmp ogt float %47, %48
br label %50
50: ; preds = %42, %39
%51 = phi i1 [ false, %39 ], [ %49, %42 ]
br i1 %51, label %52, label %67
52: ; preds = %50
%53 = load float*, float** %4
%54 = load i32, i32* %8
%55 = sext i32 %54 to i64
%56 = getelementptr inbounds float, float* %53, i64 %55
%57 = load float, float* %56
%58 = load float*, float** %4
%59 = load i32, i32* %8
%60 = load i32, i32* %6
%61 = add nsw i32 %59, %60
%62 = sext i32 %61 to i64
%63 = getelementptr inbounds float, float* %58, i64 %62
store float %57, float* %63
%64 = load i32, i32* %6
%65 = load i32, i32* %8
%66 = sub nsw i32 %65, %64
store i32 %66, i32* %8
br label %39
67: ; preds = %50
%68 = load float, float* %10
%69 = load float*, float** %4
%70 = load i32, i32* %8
%71 = load i32, i32* %6
%72 = add nsw i32 %70, %71
%73 = sext i32 %72 to i64
%74 = getelementptr inbounds float, float* %69, i64 %73
store float %68, float* %74
br label %75
75: ; preds = %67
%76 = load i32, i32* %7
%77 = add nsw i32 %76, 1
store i32 %77, i32* %7
br label %26
78: ; preds = %26
br label %79
79: ; preds = %78
%80 = load i32, i32* %5
%81 = add nsw i32 %80, 1
store i32 %81, i32* %5
br label %17
82: ; preds = %17
ret void
}

LLVM IR (-O2):
define void @shell_sysrem(i32 %0, float* nocapture %1) {
%3 = sitofp i32 %0 to double
%4 = tail call double @log(double %3)
%5 = tail call double @llvm.fmuladd.f64(double %4, double 0x3FF71547601978D6, double 1.000000e-05)
%6 = fptosi double %5 to i32
%7 = icmp slt i32 %6, 1
br i1 %7, label %45, label %8
8: ; preds = %2
%9 = sext i32 %0 to i64
br label %10
10: ; preds = %8, %42
%11 = phi i32 [ %43, %42 ], [ 1, %8 ]
%12 = phi i32 [ %13, %42 ], [ %0, %8 ]
%13 = ashr i32 %12, 1
%14 = icmp slt i32 %13, %0
br i1 %14, label %15, label %42
15: ; preds = %10
%16 = sext i32 %13 to i64
br label %17
17: ; preds = %37, %15
%18 = phi i64 [ %16, %15 ], [ %19, %37 ]
%19 = add nsw i64 %18, 1
%20 = getelementptr inbounds float, float* %1, i64 %19
%21 = load float, float* %20
%22 = trunc i64 %19 to i32
%23 = trunc i64 %19 to i32
%24 = sub i32 %23, %13
br label %25
25: ; preds = %17, %32
%26 = phi i32 [ %35, %32 ], [ %24, %17 ]
%27 = phi i32 [ %26, %32 ], [ %22, %17 ]
%28 = zext i32 %26 to i64
%29 = getelementptr inbounds float, float* %1, i64 %28
%30 = load float, float* %29
%31 = fcmp ogt float %30, %21
br i1 %31, label %32, label %37
32: ; preds = %25
%33 = sext i32 %27 to i64
%34 = getelementptr inbounds float, float* %1, i64 %33
store float %30, float* %34
%35 = sub nsw i32 %26, %13
%36 = icmp sgt i32 %35, 0
br i1 %36, label %25, label %37
37: ; preds = %25, %32
%38 = phi i32 [ %27, %25 ], [ %26, %32 ]
%39 = sext i32 %38 to i64
%40 = getelementptr inbounds float, float* %1, i64 %39
store float %21, float* %40
%41 = icmp eq i64 %19, %9
br i1 %41, label %42, label %17
42: ; preds = %37, %10
%43 = add nuw i32 %11, 1
%44 = icmp eq i32 %11, %6
br i1 %44, label %45, label %10
45: ; preds = %42, %2
ret void
}