#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 */
scalar_t__ sqrt(float);
float fa_get_subband_sqrtpower(float *X, int kmin, int kmax) {
  int k;
  float tmp;
  float Px;
  Px = 0.0;
  for (k = kmin; k <= kmax; k++) {
    tmp = ((X[k]) > 0 ? (X[k]) : (-(X[k])));
    Px += (sqrt(tmp));
  }
  return Px;
}


x86 assembly (clang, -O0):
<fa_get_subband_sqrtpower>:
push   %rbp
mov    %rsp,%rbp
sub    $0x20,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x1c(%rbp)
mov    -0xc(%rbp),%eax
mov    %eax,-0x14(%rbp)
mov    -0x14(%rbp),%eax
cmp    -0x10(%rbp),%eax
jg     11bc <fa_get_subband_sqrtpower+0xac>
mov    -0x8(%rbp),%rax
movslq -0x14(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
xorps  %xmm1,%xmm1
ucomiss %xmm1,%xmm0
jbe    116c <fa_get_subband_sqrtpower+0x5c>
mov    -0x8(%rbp),%rax
movslq -0x14(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x20(%rbp)
jmp    118b <fa_get_subband_sqrtpower+0x7b>
mov    -0x8(%rbp),%rax
movslq -0x14(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
movss  %xmm0,-0x20(%rbp)
movss  -0x20(%rbp),%xmm0
movss  %xmm0,-0x18(%rbp)
movss  -0x18(%rbp),%xmm0
call   1030 <sqrt@plt>
cvtsi2ss %rax,%xmm0
addss  -0x1c(%rbp),%xmm0
movss  %xmm0,-0x1c(%rbp)
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    1130 <fa_get_subband_sqrtpower+0x20>
movss  -0x1c(%rbp),%xmm0
add    $0x20,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<fa_get_subband_sqrtpower>:
push   %rbp
push   %r15
push   %r14
push   %r12
push   %rbx
sub    $0x10,%rsp
xorps  %xmm0,%xmm0
cmp    %esi,%edx
jl     120f <fa_get_subband_sqrtpower+0xff>
mov    %edx,%r14d
mov    %esi,%r15d
mov    %rdi,%rbx
movslq %esi,%rbp
lea    0x1(%r14),%r12d
mov    %r12d,%eax
sub    %esi,%eax
test   $0x1,%al
jne    114d <fa_get_subband_sqrtpower+0x3d>
xorps  %xmm0,%xmm0
cmp    %r15d,%r14d
jne    1187 <fa_get_subband_sqrtpower+0x77>
jmp    120f <fa_get_subband_sqrtpower+0xff>
movss  (%rbx,%rbp,4),%xmm1
movaps 0xea7(%rip),%xmm2
xorps  %xmm1,%xmm2
xorps  %xmm0,%xmm0
cmpltss %xmm1,%xmm0
andps  %xmm0,%xmm1
andnps %xmm2,%xmm0
orps   %xmm1,%xmm0
call   1030 <sqrt@plt>
xorps  %xmm0,%xmm0
cvtsi2ss %rax,%xmm0
add    $0x1,%rbp
cmp    %r15d,%r14d
je     120f <fa_get_subband_sqrtpower+0xff>
mov    %r12d,%r14d
nopw   0x0(%rax,%rax,1)
movss  %xmm0,0xc(%rsp)
movss  (%rbx,%rbp,4),%xmm1
movaps %xmm1,%xmm2
movaps 0xe5b(%rip),%xmm3
xorps  %xmm3,%xmm2
xorps  %xmm0,%xmm0
cmpltss %xmm1,%xmm0
andps  %xmm0,%xmm1
andnps %xmm2,%xmm0
orps   %xmm1,%xmm0
call   1030 <sqrt@plt>
xorps  %xmm0,%xmm0
cvtsi2ss %rax,%xmm0
addss  0xc(%rsp),%xmm0
movss  %xmm0,0xc(%rsp)
movss  0x4(%rbx,%rbp,4),%xmm1
movaps %xmm1,%xmm2
xorps  0xe1e(%rip),%xmm2
xorps  %xmm0,%xmm0
cmpltss %xmm1,%xmm0
andps  %xmm0,%xmm1
andnps %xmm2,%xmm0
orps   %xmm1,%xmm0
call   1030 <sqrt@plt>
xorps  %xmm0,%xmm0
cvtsi2ss %rax,%xmm0
addss  0xc(%rsp),%xmm0
add    $0x2,%rbp
cmp    %ebp,%r14d
jne    1190 <fa_get_subband_sqrtpower+0x80>
add    $0x10,%rsp
pop    %rbx
pop    %r12
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define float @fa_get_subband_sqrtpower(float* %0, i32 %1, i32 %2) {
%4 = alloca float*
%5 = alloca i32
%6 = alloca i32
%7 = alloca i32
%8 = alloca float
%9 = alloca float
store float* %0, float** %4
store i32 %1, i32* %5
store i32 %2, i32* %6
store float 0.000000e+00, float* %9
%10 = load i32, i32* %5
store i32 %10, i32* %7
br label %11
11: ; preds = %42, %3
%12 = load i32, i32* %7
%13 = load i32, i32* %6
%14 = icmp sle i32 %12, %13
br i1 %14, label %15, label %45
15: ; preds = %11
%16 = load float*, float** %4
%17 = load i32, i32* %7
%18 = sext i32 %17 to i64
%19 = getelementptr inbounds float, float* %16, i64 %18
%20 = load float, float* %19
%21 = fcmp ogt float %20, 0.000000e+00
br i1 %21, label %22, label %28
22: ; preds = %15
%23 = load float*, float** %4
%24 = load i32, i32* %7
%25 = sext i32 %24 to i64
%26 = getelementptr inbounds float, float* %23, i64 %25
%27 = load float, float* %26
br label %35
28: ; preds = %15
%29 = load float*, float** %4
%30 = load i32, i32* %7
%31 = sext i32 %30 to i64
%32 = getelementptr inbounds float, float* %29, i64 %31
%33 = load float, float* %32
%34 = fneg float %33
br label %35
35: ; preds = %28, %22
%36 = phi float [ %27, %22 ], [ %34, %28 ]
store float %36, float* %8
%37 = load float, float* %8
%38 = call i64 @sqrt(float %37)
%39 = sitofp i64 %38 to float
%40 = load float, float* %9
%41 = fadd float %40, %39
store float %41, float* %9
br label %42
42: ; preds = %35
%43 = load i32, i32* %7
%44 = add nsw i32 %43, 1
store i32 %44, i32* %7
br label %11
45: ; preds = %11
%46 = load float, float* %9
ret float %46
}

LLVM IR (-O2):
define float @fa_get_subband_sqrtpower(float* nocapture %0, i32 %1, i32 %2) {
%4 = icmp slt i32 %2, %1
br i1 %4, label %48, label %5
5: ; preds = %3
%6 = sext i32 %1 to i64
%7 = add i32 %2, 1
%8 = sub i32 %7, %1
%9 = and i32 %8, 1
%10 = icmp eq i32 %9, 0
br i1 %10, label %20, label %11
11: ; preds = %5
%12 = getelementptr inbounds float, float* %0, i64 %6
%13 = load float, float* %12
%14 = fcmp ogt float %13, 0.000000e+00
%15 = fneg float %13
%16 = select i1 %14, float %13, float %15
%17 = tail call i64 @sqrt(float %16)
%18 = sitofp i64 %17 to float
%19 = add nsw i64 %6, 1
br label %20
20: ; preds = %11, %5
%21 = phi float [ undef, %5 ], [ %18, %11 ]
%22 = phi i64 [ %6, %5 ], [ %19, %11 ]
%23 = phi float [ 0.000000e+00, %5 ], [ %18, %11 ]
%24 = icmp eq i32 %2, %1
br i1 %24, label %48, label %25
25: ; preds = %20, %25
%26 = phi i64 [ %45, %25 ], [ %22, %20 ]
%27 = phi float [ %44, %25 ], [ %23, %20 ]
%28 = getelementptr inbounds float, float* %0, i64 %26
%29 = load float, float* %28
%30 = fcmp ogt float %29, 0.000000e+00
%31 = fneg float %29
%32 = select i1 %30, float %29, float %31
%33 = tail call i64 @sqrt(float %32)
%34 = sitofp i64 %33 to float
%35 = fadd float %27, %34
%36 = add nsw i64 %26, 1
%37 = getelementptr inbounds float, float* %0, i64 %36
%38 = load float, float* %37
%39 = fcmp ogt float %38, 0.000000e+00
%40 = fneg float %38
%41 = select i1 %39, float %38, float %40
%42 = tail call i64 @sqrt(float %41)
%43 = sitofp i64 %42 to float
%44 = fadd float %35, %43
%45 = add nsw i64 %26, 2
%46 = trunc i64 %45 to i32
%47 = icmp eq i32 %7, %46
br i1 %47, label %48, label %25
48: ; preds = %20, %25, %3
%49 = phi float [ 0.000000e+00, %3 ], [ %21, %20 ], [ %44, %25 ]
ret float %49
}