#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
int /*<<< orphan*/ fqmedian(float *, int);
scalar_t__ pow(double, double);
float compute_fwhmrange(float *fwhm, int nfwhm, float maxvar, float minin,
                        float maxin, float *minout, float *maxout) {
  float *fwhmt, *fwhmt2, df, dfmin, fmin;
  int i, nw;
  fqmedian(fwhm, nfwhm);
  nw = nfwhm / 4;
  if (nw < 4)
    nw = 1;
  dfmin = 1e+30;
  fmin = 0.0;
  fwhmt = fwhm;
  fwhmt2 = fwhm + nw;
  for (i = nfwhm - nw; i--; fwhmt++, fwhmt2++) {
    if ((df = *fwhmt2 - *fwhmt) < dfmin) {
      dfmin = df;
      fmin = (*fwhmt2 + *fwhmt) / 2.0;
    }
  }
  if (nfwhm < 2)
    fmin = *fwhm;
  dfmin = (float)pow((double)maxvar + 1.0, 0.3333333);
  *minout = dfmin > 0.0 ? fmin / dfmin : 0.0;
  if (*minout < minin)
    *minout = minin;
  *maxout = fmin * dfmin * dfmin;
  if (*maxout > maxin)
    *maxout = maxin;
  return fmin;
}


x86 assembly (clang, -O0):
<compute_fwhmrange>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
movss  %xmm0,-0x10(%rbp)
movss  %xmm1,-0x14(%rbp)
movss  %xmm2,-0x18(%rbp)
mov    %rdx,-0x20(%rbp)
mov    %rcx,-0x28(%rbp)
mov    -0x8(%rbp),%rdi
mov    -0xc(%rbp),%esi
call   1040 <fqmedian@plt>
mov    -0xc(%rbp),%eax
mov    $0x4,%ecx
cltd
idiv   %ecx
mov    %eax,-0x4c(%rbp)
cmpl   $0x4,-0x4c(%rbp)
jge    1171 <compute_fwhmrange+0x51>
movl   $0x1,-0x4c(%rbp)
movss  0xe87(%rip),%xmm0
movss  %xmm0,-0x40(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x44(%rbp)
mov    -0x8(%rbp),%rax
mov    %rax,-0x30(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x4c(%rbp),%rcx
shl    $0x2,%rcx
add    %rcx,%rax
mov    %rax,-0x38(%rbp)
mov    -0xc(%rbp),%eax
sub    -0x4c(%rbp),%eax
mov    %eax,-0x48(%rbp)
mov    -0x48(%rbp),%eax
mov    %eax,%ecx
add    $0xffffffff,%ecx
mov    %ecx,-0x48(%rbp)
cmp    $0x0,%eax
je     1236 <compute_fwhmrange+0x116>
mov    -0x38(%rbp),%rax
movss  (%rax),%xmm1
mov    -0x30(%rbp),%rax
subss  (%rax),%xmm1
movss  %xmm1,-0x3c(%rbp)
movss  -0x40(%rbp),%xmm0
ucomiss %xmm1,%xmm0
jbe    1214 <compute_fwhmrange+0xf4>
movss  -0x3c(%rbp),%xmm0
movss  %xmm0,-0x40(%rbp)
mov    -0x38(%rbp),%rax
movss  (%rax),%xmm0
mov    -0x30(%rbp),%rax
addss  (%rax),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xe11(%rip),%xmm1
divsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x44(%rbp)
jmp    1219 <compute_fwhmrange+0xf9>
mov    -0x30(%rbp),%rax
add    $0x4,%rax
mov    %rax,-0x30(%rbp)
mov    -0x38(%rbp),%rax
add    $0x4,%rax
mov    %rax,-0x38(%rbp)
jmp    11aa <compute_fwhmrange+0x8a>
cmpl   $0x2,-0xc(%rbp)
jge    124d <compute_fwhmrange+0x12d>
mov    -0x8(%rbp),%rax
movss  (%rax),%xmm0
movss  %xmm0,-0x44(%rbp)
movss  -0x10(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xdb2(%rip),%xmm1
addsd  %xmm1,%xmm0
movsd  0xd9e(%rip),%xmm1
call   1030 <pow@plt>
cvtsi2ss %rax,%xmm0
movss  %xmm0,-0x40(%rbp)
movss  -0x40(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    12a7 <compute_fwhmrange+0x187>
movss  -0x44(%rbp),%xmm0
divss  -0x40(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  %xmm0,-0x58(%rbp)
jmp    12b4 <compute_fwhmrange+0x194>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x58(%rbp)
jmp    12b4 <compute_fwhmrange+0x194>
movsd  -0x58(%rbp),%xmm0
cvtsd2ss %xmm0,%xmm0
mov    -0x20(%rbp),%rax
movss  %xmm0,(%rax)
mov    -0x20(%rbp),%rax
movss  (%rax),%xmm1
movss  -0x14(%rbp),%xmm0
ucomiss %xmm1,%xmm0
jbe    12e8 <compute_fwhmrange+0x1c8>
movss  -0x14(%rbp),%xmm0
mov    -0x20(%rbp),%rax
movss  %xmm0,(%rax)
movss  -0x44(%rbp),%xmm0
mulss  -0x40(%rbp),%xmm0
mulss  -0x40(%rbp),%xmm0
mov    -0x28(%rbp),%rax
movss  %xmm0,(%rax)
mov    -0x28(%rbp),%rax
movss  (%rax),%xmm0
ucomiss -0x18(%rbp),%xmm0
jbe    131e <compute_fwhmrange+0x1fe>
movss  -0x18(%rbp),%xmm0
mov    -0x28(%rbp),%rax
movss  %xmm0,(%rax)
movss  -0x44(%rbp),%xmm0
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<compute_fwhmrange>:
push   %r15
push   %r14
push   %r12
push   %rbx
sub    $0x28,%rsp
mov    %rcx,%r14
mov    %rdx,%r15
movss  %xmm2,0xc(%rsp)
movss  %xmm1,0x8(%rsp)
movss  %xmm0,0x4(%rsp)
mov    %esi,%ebx
mov    %rdi,%r12
call   1040 <fqmedian@plt>
mov    $0x1,%edx
cmp    $0x10,%ebx
jl     1162 <compute_fwhmrange+0x42>
lea    0x3(%rbx),%edx
test   %ebx,%ebx
cmovns %ebx,%edx
sar    $0x2,%edx
mov    %ebx,%esi
sub    %edx,%esi
je     118e <compute_fwhmrange+0x6e>
movslq %edx,%rax
lea    (%r12,%rax,4),%rax
not    %edx
test   $0x1,%sil
jne    1196 <compute_fwhmrange+0x76>
movss  0xe81(%rip),%xmm0
xorps  %xmm7,%xmm7
mov    %r12,%rcx
add    %ebx,%edx
jne    11da <compute_fwhmrange+0xba>
jmp    1259 <compute_fwhmrange+0x139>
xorps  %xmm7,%xmm7
jmp    1259 <compute_fwhmrange+0x139>
add    $0xffffffff,%esi
movss  (%rax),%xmm1
movss  (%r12),%xmm0
movaps %xmm1,%xmm7
subss  %xmm0,%xmm7
addss  %xmm0,%xmm1
mulss  0xe4e(%rip),%xmm1
movss  0xe42(%rip),%xmm2
movaps %xmm7,%xmm0
minss  %xmm2,%xmm0
cmpltss %xmm2,%xmm7
andps  %xmm1,%xmm7
lea    0x4(%r12),%rcx
add    $0x4,%rax
add    %ebx,%edx
je     1259 <compute_fwhmrange+0x139>
mov    %esi,%edx
xor    %esi,%esi
movss  0xe1e(%rip),%xmm1
cs nopw 0x0(%rax,%rax,1)
movss  (%rax,%rsi,4),%xmm2
movss  0x4(%rax,%rsi,4),%xmm3
movss  (%rcx,%rsi,4),%xmm4
movss  0x4(%rcx,%rsi,4),%xmm5
movaps %xmm2,%xmm6
subss  %xmm4,%xmm6
addss  %xmm4,%xmm2
mulss  %xmm1,%xmm2
movaps %xmm6,%xmm4
minss  %xmm0,%xmm4
cmpltss %xmm0,%xmm6
andps  %xmm6,%xmm2
andnps %xmm7,%xmm6
orps   %xmm2,%xmm6
movaps %xmm3,%xmm7
subss  %xmm5,%xmm7
addss  %xmm5,%xmm3
mulss  %xmm1,%xmm3
movaps %xmm7,%xmm0
minss  %xmm4,%xmm0
cmpltss %xmm4,%xmm7
movaps %xmm7,%xmm2
andnps %xmm6,%xmm2
andps  %xmm3,%xmm7
orps   %xmm2,%xmm7
add    $0x2,%rsi
cmp    %esi,%edx
jne    11f0 <compute_fwhmrange+0xd0>
cmp    $0x1,%ebx
jg     1264 <compute_fwhmrange+0x144>
movss  (%r12),%xmm7
movaps %xmm7,0x10(%rsp)
movss  0x4(%rsp),%xmm0
cvtss2sd %xmm0,%xmm0
addsd  0xd8d(%rip),%xmm0
movsd  0xd8d(%rip),%xmm1
call   1030 <pow@plt>
xorps  %xmm3,%xmm3
cvtsi2ss %rax,%xmm3
test   %rax,%rax
jle    12a3 <compute_fwhmrange+0x183>
movaps 0x10(%rsp),%xmm0
movaps %xmm0,%xmm1
divss  %xmm3,%xmm1
jmp    12ab <compute_fwhmrange+0x18b>
xorps  %xmm1,%xmm1
movaps 0x10(%rsp),%xmm0
movss  0x8(%rsp),%xmm2
maxss  %xmm1,%xmm2
movss  %xmm2,(%r15)
movaps %xmm0,%xmm1
mulss  %xmm3,%xmm1
mulss  %xmm3,%xmm1
movss  0xc(%rsp),%xmm2
minss  %xmm1,%xmm2
movss  %xmm2,(%r14)
add    $0x28,%rsp
pop    %rbx
pop    %r12
pop    %r14
pop    %r15
ret

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

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