int HPDinterval(double x[], int n, double HPD[2], double alpha) {
  int jL0 = (int)(n * alpha / 2), jU0 = (int)(n * (1 - alpha / 2)), jL, jU,
      jLb = jL0;
  double w0 = x[jU0] - x[jL0], w = w0;
  int debug = 0;
  HPD[0] = x[jL0];
  HPD[1] = x[jU0];
  if (n < 3)
    return (-1);
  for (jL = 0, jU = jL + (jU0 - jL0); jU < n; jL++, jU++) {
    if (x[jU] - x[jL] < w) {
      jLb = jL;
      w = x[jU] - x[jL];
    }
  }
  HPD[0] = x[jLb];
  HPD[1] = x[jLb + jU0 - jL0];
  return (0);
}


x86 assembly (clang, -O0):
<HPDinterval>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x10(%rbp)
mov    %esi,-0x14(%rbp)
mov    %rdx,-0x20(%rbp)
movsd  %xmm0,-0x28(%rbp)
cvtsi2sdl -0x14(%rbp),%xmm0
mulsd  -0x28(%rbp),%xmm0
movsd  0xee2(%rip),%xmm1
divsd  %xmm1,%xmm0
cvttsd2si %xmm0,%eax
mov    %eax,-0x2c(%rbp)
cvtsi2sdl -0x14(%rbp),%xmm0
movsd  -0x28(%rbp),%xmm2
movsd  0xec5(%rip),%xmm1
divsd  %xmm1,%xmm2
movsd  0xeb1(%rip),%xmm1
subsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
cvttsd2si %xmm0,%eax
mov    %eax,-0x30(%rbp)
mov    -0x2c(%rbp),%eax
mov    %eax,-0x3c(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
movl   $0x0,-0x54(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x10(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
movsd  %xmm0,0x8(%rax)
cmpl   $0x3,-0x14(%rbp)
jge    11d5 <HPDinterval+0xd5>
movl   $0xffffffff,-0x4(%rbp)
jmp    129a <HPDinterval+0x19a>
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
mov    -0x30(%rbp),%ecx
sub    -0x2c(%rbp),%ecx
add    %ecx,%eax
mov    %eax,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1260 <HPDinterval+0x160>
mov    -0x10(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm1
movsd  -0x50(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    1244 <HPDinterval+0x144>
mov    -0x34(%rbp),%eax
mov    %eax,-0x3c(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x50(%rbp)
jmp    1249 <HPDinterval+0x149>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    11ea <HPDinterval+0xea>
mov    -0x10(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x10(%rbp),%rax
mov    -0x3c(%rbp),%ecx
add    -0x30(%rbp),%ecx
sub    -0x2c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
movsd  %xmm0,0x8(%rax)
movl   $0x0,-0x4(%rbp)
mov    -0x4(%rbp),%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<HPDinterval>:
cvtsi2sd %esi,%xmm1
movapd %xmm1,%xmm2
mulsd  %xmm0,%xmm2
movsd  0xeec(%rip),%xmm3
mulsd  %xmm3,%xmm2
cvttsd2si %xmm2,%r9d
mulsd  %xmm3,%xmm0
movsd  0xedf(%rip),%xmm2
subsd  %xmm0,%xmm2
mulsd  %xmm1,%xmm2
cvttsd2si %xmm2,%ecx
movslq %ecx,%rax
movsd  (%rdi,%rax,8),%xmm0
movslq %r9d,%r10
movsd  (%rdi,%r10,8),%xmm1
movsd  %xmm1,(%rdx)
movsd  (%rdi,%rax,8),%xmm2
movsd  %xmm2,0x8(%rdx)
mov    $0xffffffff,%eax
cmp    $0x3,%esi
jl     120b <HPDinterval+0x10b>
mov    %ecx,%r8d
sub    %r9d,%r8d
cmp    %esi,%r8d
jge    11ef <HPDinterval+0xef>
subsd  %xmm1,%xmm0
movslq %r8d,%r11
add    %r9d,%esi
sub    %ecx,%esi
cmp    $0x1,%esi
jne    1186 <HPDinterval+0x86>
xor    %ecx,%ecx
jmp    11d3 <HPDinterval+0xd3>
mov    %esi,%r10d
and    $0xfffffffe,%r10d
xor    %ecx,%ecx
nop
movsd  (%rdi,%r11,8),%xmm1
subsd  (%rdi,%rcx,8),%xmm1
ucomisd %xmm1,%xmm0
minsd  %xmm0,%xmm1
movsd  0x8(%rdi,%r11,8),%xmm0
mov    %r9d,%eax
cmova  %ecx,%eax
subsd  0x8(%rdi,%rcx,8),%xmm0
ucomisd %xmm0,%xmm1
lea    0x1(%rcx),%r9d
cmovbe %eax,%r9d
minsd  %xmm1,%xmm0
add    $0x2,%rcx
add    $0x2,%r11
cmp    %r10,%rcx
jne    1190 <HPDinterval+0x90>
test   $0x1,%sil
je     11ec <HPDinterval+0xec>
movsd  (%rdi,%r11,8),%xmm1
subsd  (%rdi,%rcx,8),%xmm1
ucomisd %xmm1,%xmm0
cmova  %ecx,%r9d
movslq %r9d,%r10
movsd  (%rdi,%r10,8),%xmm0
movsd  %xmm0,(%rdx)
add    %r9d,%r8d
movslq %r8d,%rax
movsd  (%rdi,%rax,8),%xmm0
movsd  %xmm0,0x8(%rdx)
xor    %eax,%eax
ret

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

LLVM IR (-O2):
define i32 @HPDinterval(double* nocapture %0, i32 %1, double* nocapture writeonly %2, double %3) {
%5 = sitofp i32 %1 to double
%6 = fmul double %5, %3
%7 = fmul double %6, 5.000000e-01
%8 = fptosi double %7 to i32
%9 = fmul double %3, 5.000000e-01
%10 = fsub double 1.000000e+00, %9
%11 = fmul double %10, %5
%12 = fptosi double %11 to i32
%13 = sext i32 %12 to i64
%14 = getelementptr inbounds double, double* %0, i64 %13
%15 = load double, double* %14
%16 = sext i32 %8 to i64
%17 = getelementptr inbounds double, double* %0, i64 %16
%18 = load double, double* %17
store double %18, double* %2
%19 = load double, double* %14
%20 = getelementptr inbounds double, double* %2, i64 1
store double %19, double* %20
%21 = icmp slt i32 %1, 3
br i1 %21, label %93, label %22
22: ; preds = %4
%23 = sub i32 %12, %8
%24 = icmp slt i32 %23, %1
br i1 %24, label %25, label %84
25: ; preds = %22
%26 = fsub double %15, %18
%27 = sext i32 %23 to i64
%28 = add i32 %8, %1
%29 = sub i32 %28, %12
%30 = zext i32 %29 to i64
%31 = and i64 %30, 1
%32 = icmp eq i32 %29, 1
br i1 %32, label %65, label %33
33: ; preds = %25
%34 = and i64 %30, 4294967294
br label %35
35: ; preds = %35, %33
%36 = phi i64 [ %27, %33 ], [ %62, %35 ]
%37 = phi i64 [ 0, %33 ], [ %61, %35 ]
%38 = phi double [ %26, %33 ], [ %60, %35 ]
%39 = phi i32 [ %8, %33 ], [ %59, %35 ]
%40 = phi i64 [ 0, %33 ], [ %63, %35 ]
%41 = getelementptr inbounds double, double* %0, i64 %36
%42 = load double, double* %41
%43 = getelementptr inbounds double, double* %0, i64 %37
%44 = load double, double* %43
%45 = fsub double %42, %44
%46 = fcmp olt double %45, %38
%47 = trunc i64 %37 to i32
%48 = select i1 %46, i32 %47, i32 %39
%49 = select i1 %46, double %45, double %38
%50 = or i64 %37, 1
%51 = add nsw i64 %36, 1
%52 = getelementptr inbounds double, double* %0, i64 %51
%53 = load double, double* %52
%54 = getelementptr inbounds double, double* %0, i64 %50
%55 = load double, double* %54
%56 = fsub double %53, %55
%57 = fcmp olt double %56, %49
%58 = trunc i64 %50 to i32
%59 = select i1 %57, i32 %58, i32 %48
%60 = select i1 %57, double %56, double %49
%61 = add nuw nsw i64 %37, 2
%62 = add nsw i64 %36, 2
%63 = add i64 %40, 2
%64 = icmp eq i64 %63, %34
br i1 %64, label %65, label %35
65: ; preds = %35, %25
%66 = phi i32 [ undef, %25 ], [ %59, %35 ]
%67 = phi i64 [ %27, %25 ], [ %62, %35 ]
%68 = phi i64 [ 0, %25 ], [ %61, %35 ]
%69 = phi double [ %26, %25 ], [ %60, %35 ]
%70 = phi i32 [ %8, %25 ], [ %59, %35 ]
%71 = icmp eq i64 %31, 0
br i1 %71, label %81, label %72
72: ; preds = %65
%73 = getelementptr inbounds double, double* %0, i64 %67
%74 = load double, double* %73
%75 = getelementptr inbounds double, double* %0, i64 %68
%76 = load double, double* %75
%77 = fsub double %74, %76
%78 = fcmp olt double %77, %69
%79 = trunc i64 %68 to i32
%80 = select i1 %78, i32 %79, i32 %70
br label %81
81: ; preds = %65, %72
%82 = phi i32 [ %66, %65 ], [ %80, %72 ]
%83 = sext i32 %82 to i64
br label %84
84: ; preds = %81, %22
%85 = phi i64 [ %83, %81 ], [ %16, %22 ]
%86 = phi i32 [ %82, %81 ], [ %8, %22 ]
%87 = getelementptr inbounds double, double* %0, i64 %85
%88 = load double, double* %87
store double %88, double* %2
%89 = add i32 %23, %86
%90 = sext i32 %89 to i64
%91 = getelementptr inbounds double, double* %0, i64 %90
%92 = load double, double* %91
store double %92, double* %20
br label %93
93: ; preds = %4, %84
%94 = phi i32 [ 0, %84 ], [ -1, %4 ]
ret i32 %94
}