void band_den_bin(int *n, int *nb, double *d, double *x, int *cnt) {
  int i, j, ii, jj, iij, nn = *n;
  double xmin, xmax, rang, dd;
  for (i = 0; i < *nb; i++)
    cnt[i] = 0;
  xmin = xmax = x[0];
  for (i = 1; i < nn; i++) {
    xmin = ((xmin) < (x[i]) ? (xmin) : (x[i]));
    xmax = ((xmax) > (x[i]) ? (xmax) : (x[i]));
  }
  rang = (xmax - xmin) * 1.01;
  *d = dd = rang / (*nb);
  for (i = 1; i < nn; i++) {
    ii = x[i] / dd;
    for (j = 0; j < i; j++) {
      jj = x[j] / dd;
      iij = ((ii - jj) > 0 ? (ii - jj) : -(ii - jj));
      cnt[iij]++;
    }
  }
}


x86 assembly (clang, -O0):
<band_den_bin>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    -0x8(%rbp),%rax
mov    (%rax),%eax
mov    %eax,-0x40(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
mov    -0x10(%rbp),%rcx
cmp    (%rcx),%eax
jge    1154 <band_den_bin+0x54>
mov    -0x28(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movl   $0x0,(%rax,%rcx,4)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1128 <band_den_bin+0x28>
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  %xmm0,-0x48(%rbp)
movl   $0x1,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x40(%rbp),%eax
jge    1211 <band_den_bin+0x111>
movsd  -0x48(%rbp),%xmm1
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
ucomisd %xmm1,%xmm0
jbe    11a4 <band_den_bin+0xa4>
movsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x68(%rbp)
jmp    11b6 <band_den_bin+0xb6>
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  -0x68(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x50(%rbp),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
ucomisd (%rax,%rcx,8),%xmm0
jbe    11e7 <band_den_bin+0xe7>
movsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x70(%rbp)
jmp    11f9 <band_den_bin+0xf9>
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x70(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    116d <band_den_bin+0x6d>
movsd  -0x50(%rbp),%xmm0
subsd  -0x48(%rbp),%xmm0
movsd  0xddd(%rip),%xmm1
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x58(%rbp),%xmm0
mov    -0x10(%rbp),%rax
cvtsi2sdl (%rax),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    -0x18(%rbp),%rax
movsd  %xmm0,(%rax)
movl   $0x1,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x40(%rbp),%eax
jge    1304 <band_den_bin+0x204>
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  -0x60(%rbp),%xmm0
cvttsd2si %xmm0,%eax
mov    %eax,-0x34(%rbp)
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x2c(%rbp),%eax
jge    12f1 <band_den_bin+0x1f1>
mov    -0x20(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  -0x60(%rbp),%xmm0
cvttsd2si %xmm0,%eax
mov    %eax,-0x38(%rbp)
mov    -0x34(%rbp),%eax
sub    -0x38(%rbp),%eax
cmp    $0x0,%eax
jle    12bf <band_den_bin+0x1bf>
mov    -0x34(%rbp),%eax
sub    -0x38(%rbp),%eax
mov    %eax,-0x74(%rbp)
jmp    12cc <band_den_bin+0x1cc>
mov    -0x34(%rbp),%ecx
sub    -0x38(%rbp),%ecx
xor    %eax,%eax
sub    %ecx,%eax
mov    %eax,-0x74(%rbp)
mov    -0x74(%rbp),%eax
mov    %eax,-0x3c(%rbp)
mov    -0x28(%rbp),%rax
movslq -0x3c(%rbp),%rcx
mov    (%rax,%rcx,4),%edx
add    $0x1,%edx
mov    %edx,(%rax,%rcx,4)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    127d <band_den_bin+0x17d>
jmp    12f6 <band_den_bin+0x1f6>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1251 <band_den_bin+0x151>
pop    %rbp
ret

x86 assembly (clang, -O2):
<band_den_bin>:
mov    (%rdi),%r9d
mov    (%rsi),%edi
test   %edi,%edi
jle    1124 <band_den_bin+0x24>
xor    %eax,%eax
nopl   0x0(%rax,%rax,1)
movl   $0x0,(%r8,%rax,4)
add    $0x1,%rax
movslq (%rsi),%rdi
cmp    %rdi,%rax
jl     1110 <band_den_bin+0x10>
movsd  (%rcx),%xmm0
cmp    $0x1,%r9d
jle    1156 <band_den_bin+0x56>
lea    -0x1(%r9),%rsi
lea    -0x2(%r9),%rax
mov    %esi,%r10d
and    $0x3,%r10d
cmp    $0x3,%rax
jae    116f <band_den_bin+0x6f>
mov    $0x1,%eax
movapd %xmm0,%xmm1
test   %r10,%r10
jne    11ca <band_den_bin+0xca>
jmp    11e6 <band_den_bin+0xe6>
subsd  %xmm0,%xmm0
mulsd  0xe9e(%rip),%xmm0
cvtsi2sd %edi,%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rdx)
ret
and    $0xfffffffffffffffc,%rsi
xor    %eax,%eax
movapd %xmm0,%xmm1
nopl   0x0(%rax)
movsd  0x8(%rcx,%rax,8),%xmm2
movsd  0x10(%rcx,%rax,8),%xmm3
minsd  %xmm2,%xmm1
maxsd  %xmm2,%xmm0
minsd  %xmm3,%xmm1
maxsd  %xmm3,%xmm0
movsd  0x18(%rcx,%rax,8),%xmm2
minsd  %xmm2,%xmm1
maxsd  %xmm2,%xmm0
movsd  0x20(%rcx,%rax,8),%xmm2
minsd  %xmm2,%xmm1
maxsd  %xmm2,%xmm0
add    $0x4,%rax
cmp    %rax,%rsi
jne    1180 <band_den_bin+0x80>
add    $0x1,%rax
test   %r10,%r10
je     11e6 <band_den_bin+0xe6>
lea    (%rcx,%rax,8),%rax
xor    %esi,%esi
movsd  (%rax,%rsi,8),%xmm2
minsd  %xmm2,%xmm1
maxsd  %xmm2,%xmm0
add    $0x1,%rsi
cmp    %rsi,%r10
jne    11d0 <band_den_bin+0xd0>
subsd  %xmm1,%xmm0
mulsd  0xe0e(%rip),%xmm0
xorps  %xmm1,%xmm1
cvtsi2sd %edi,%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rdx)
cmp    $0x2,%r9d
jl     124f <band_den_bin+0x14f>
mov    $0x1,%r10d
nopl   (%rax)
movsd  (%rcx,%r10,8),%xmm1
divsd  %xmm0,%xmm1
cvttsd2si %xmm1,%eax
xor    %esi,%esi
movsd  (%rcx,%rsi,8),%xmm1
divsd  %xmm0,%xmm1
cvttsd2si %xmm1,%edi
mov    %eax,%edx
sub    %edi,%edx
mov    %edx,%edi
neg    %edi
cmovs  %edx,%edi
addl   $0x1,(%r8,%rdi,4)
add    $0x1,%rsi
cmp    %rsi,%r10
jne    1220 <band_den_bin+0x120>
add    $0x1,%r10
cmp    %r9,%r10
jne    1210 <band_den_bin+0x110>
ret

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

LLVM IR (-O2):
define void @band_den_bin(i32* nocapture %0, i32* nocapture %1, double* nocapture writeonly %2, double* nocapture %3, i32* nocapture %4) {
%6 = load i32, i32* %0
%7 = load i32, i32* %1
%8 = icmp sgt i32 %7, 0
br i1 %8, label %9, label %16
9: ; preds = %5, %9
%10 = phi i64 [ %12, %9 ], [ 0, %5 ]
%11 = getelementptr inbounds i32, i32* %4, i64 %10
store i32 0, i32* %11
%12 = add nuw nsw i64 %10, 1
%13 = load i32, i32* %1
%14 = sext i32 %13 to i64
%15 = icmp slt i64 %12, %14
br i1 %15, label %9, label %16
16: ; preds = %9, %5
%17 = phi i32 [ %7, %5 ], [ %13, %9 ]
%18 = load double, double* %3
%19 = icmp sgt i32 %6, 1
br i1 %19, label %25, label %20
20: ; preds = %16
%21 = fsub double %18, %18
%22 = fmul double %21, 1.010000e+00
%23 = sitofp i32 %17 to double
%24 = fdiv double %22, %23
store double %24, double* %2
br label %121
25: ; preds = %16
%26 = zext i32 %6 to i64
%27 = add nsw i64 %26, -1
%28 = add nsw i64 %26, -2
%29 = and i64 %27, 3
%30 = icmp ult i64 %28, 3
br i1 %30, label %68, label %31
31: ; preds = %25
%32 = and i64 %27, -4
br label %33
33: ; preds = %33, %31
%34 = phi i64 [ 1, %31 ], [ %65, %33 ]
%35 = phi double [ %18, %31 ], [ %64, %33 ]
%36 = phi double [ %18, %31 ], [ %62, %33 ]
%37 = phi i64 [ 0, %31 ], [ %66, %33 ]
%38 = getelementptr inbounds double, double* %3, i64 %34
%39 = load double, double* %38
%40 = fcmp olt double %36, %39
%41 = select i1 %40, double %36, double %39
%42 = fcmp ogt double %35, %39
%43 = select i1 %42, double %35, double %39
%44 = add nuw nsw i64 %34, 1
%45 = getelementptr inbounds double, double* %3, i64 %44
%46 = load double, double* %45
%47 = fcmp olt double %41, %46
%48 = select i1 %47, double %41, double %46
%49 = fcmp ogt double %43, %46
%50 = select i1 %49, double %43, double %46
%51 = add nuw nsw i64 %34, 2
%52 = getelementptr inbounds double, double* %3, i64 %51
%53 = load double, double* %52
%54 = fcmp olt double %48, %53
%55 = select i1 %54, double %48, double %53
%56 = fcmp ogt double %50, %53
%57 = select i1 %56, double %50, double %53
%58 = add nuw nsw i64 %34, 3
%59 = getelementptr inbounds double, double* %3, i64 %58
%60 = load double, double* %59
%61 = fcmp olt double %55, %60
%62 = select i1 %61, double %55, double %60
%63 = fcmp ogt double %57, %60
%64 = select i1 %63, double %57, double %60
%65 = add nuw nsw i64 %34, 4
%66 = add i64 %37, 4
%67 = icmp eq i64 %66, %32
br i1 %67, label %68, label %33
68: ; preds = %33, %25
%69 = phi double [ undef, %25 ], [ %62, %33 ]
%70 = phi double [ undef, %25 ], [ %64, %33 ]
%71 = phi i64 [ 1, %25 ], [ %65, %33 ]
%72 = phi double [ %18, %25 ], [ %64, %33 ]
%73 = phi double [ %18, %25 ], [ %62, %33 ]
%74 = icmp eq i64 %29, 0
br i1 %74, label %89, label %75
75: ; preds = %68, %75
%76 = phi i64 [ %86, %75 ], [ %71, %68 ]
%77 = phi double [ %85, %75 ], [ %72, %68 ]
%78 = phi double [ %83, %75 ], [ %73, %68 ]
%79 = phi i64 [ %87, %75 ], [ 0, %68 ]
%80 = getelementptr inbounds double, double* %3, i64 %76
%81 = load double, double* %80
%82 = fcmp olt double %78, %81
%83 = select i1 %82, double %78, double %81
%84 = fcmp ogt double %77, %81
%85 = select i1 %84, double %77, double %81
%86 = add nuw nsw i64 %76, 1
%87 = add i64 %79, 1
%88 = icmp eq i64 %87, %29
br i1 %88, label %89, label %75
89: ; preds = %75, %68
%90 = phi double [ %69, %68 ], [ %83, %75 ]
%91 = phi double [ %70, %68 ], [ %85, %75 ]
%92 = fsub double %91, %90
%93 = fmul double %92, 1.010000e+00
%94 = sitofp i32 %17 to double
%95 = fdiv double %93, %94
store double %95, double* %2
br i1 %19, label %96, label %121
96: ; preds = %89
%97 = zext i32 %6 to i64
br label %98
98: ; preds = %96, %118
%99 = phi i64 [ 1, %96 ], [ %119, %118 ]
%100 = getelementptr inbounds double, double* %3, i64 %99
%101 = load double, double* %100
%102 = fdiv double %101, %95
%103 = fptosi double %102 to i32
br label %104
104: ; preds = %98, %104
%105 = phi i64 [ 0, %98 ], [ %116, %104 ]
%106 = getelementptr inbounds double, double* %3, i64 %105
%107 = load double, double* %106
%108 = fdiv double %107, %95
%109 = fptosi double %108 to i32
%110 = sub nsw i32 %103, %109
%111 = tail call i32 @llvm.abs.i32(i32 %110, i1 true)
%112 = zext i32 %111 to i64
%113 = getelementptr inbounds i32, i32* %4, i64 %112
%114 = load i32, i32* %113
%115 = add nsw i32 %114, 1
store i32 %115, i32* %113
%116 = add nuw nsw i64 %105, 1
%117 = icmp eq i64 %116, %99
br i1 %117, label %118, label %104
118: ; preds = %104
%119 = add nuw nsw i64 %99, 1
%120 = icmp eq i64 %119, %97
br i1 %120, label %121, label %98
121: ; preds = %118, %20, %89
ret void
}