float fa_get_subband_power(float *X, int kmin, int kmax) {
  int k;
  float Px;
  Px = 0.0;
  for (k = kmin; k <= kmax; k++)
    Px += X[k] * X[k];
  return Px;
}


x86 assembly (clang, -O0):
<fa_get_subband_power>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x18(%rbp)
mov    -0xc(%rbp),%eax
mov    %eax,-0x14(%rbp)
mov    -0x14(%rbp),%eax
cmp    -0x10(%rbp),%eax
jg     1162 <fa_get_subband_power+0x62>
mov    -0x8(%rbp),%rax
movslq -0x14(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x14(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
movss  -0x18(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x18(%rbp)
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    111c <fa_get_subband_power+0x1c>
movss  -0x18(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<fa_get_subband_power>:
xorps  %xmm0,%xmm0
mov    %edx,%r8d
sub    %esi,%r8d
jl     117f <fa_get_subband_power+0x7f>
movslq %esi,%rax
add    $0x1,%edx
mov    %edx,%ecx
sub    %esi,%ecx
xorps  %xmm0,%xmm0
and    $0x3,%ecx
je     1136 <fa_get_subband_power+0x36>
nopl   (%rax)
movss  (%rdi,%rax,4),%xmm1
mulss  %xmm1,%xmm1
addss  %xmm1,%xmm0
add    $0x1,%rax
add    $0xffffffff,%ecx
jne    1120 <fa_get_subband_power+0x20>
cmp    $0x3,%r8d
jb     117f <fa_get_subband_power+0x7f>
mov    %edx,%ecx
xchg   %ax,%ax
movss  (%rdi,%rax,4),%xmm1
movss  0x4(%rdi,%rax,4),%xmm2
mulss  %xmm1,%xmm1
addss  %xmm0,%xmm1
mulss  %xmm2,%xmm2
addss  %xmm1,%xmm2
movss  0x8(%rdi,%rax,4),%xmm1
mulss  %xmm1,%xmm1
addss  %xmm2,%xmm1
movss  0xc(%rdi,%rax,4),%xmm0
mulss  %xmm0,%xmm0
addss  %xmm1,%xmm0
add    $0x4,%rax
cmp    %eax,%ecx
jne    1140 <fa_get_subband_power+0x40>
ret

LLVM IR (-O0):
define float @fa_get_subband_power(float* %0, i32 %1, i32 %2) {
%4 = alloca float*
%5 = alloca i32
%6 = alloca i32
%7 = alloca i32
%8 = alloca float
store float* %0, float** %4
store i32 %1, i32* %5
store i32 %2, i32* %6
store float 0.000000e+00, float* %8
%9 = load i32, i32* %5
store i32 %9, i32* %7
br label %10
10: ; preds = %27, %3
%11 = load i32, i32* %7
%12 = load i32, i32* %6
%13 = icmp sle i32 %11, %12
br i1 %13, label %14, label %30
14: ; preds = %10
%15 = load float*, float** %4
%16 = load i32, i32* %7
%17 = sext i32 %16 to i64
%18 = getelementptr inbounds float, float* %15, i64 %17
%19 = load float, float* %18
%20 = load float*, float** %4
%21 = load i32, i32* %7
%22 = sext i32 %21 to i64
%23 = getelementptr inbounds float, float* %20, i64 %22
%24 = load float, float* %23
%25 = load float, float* %8
%26 = call float @llvm.fmuladd.f32(float %19, float %24, float %25)
store float %26, float* %8
br label %27
27: ; preds = %14
%28 = load i32, i32* %7
%29 = add nsw i32 %28, 1
store i32 %29, i32* %7
br label %10
30: ; preds = %10
%31 = load float, float* %8
ret float %31
}

LLVM IR (-O2):
define float @fa_get_subband_power(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 = sub i32 %2, %1
%10 = and i32 %8, 3
%11 = icmp eq i32 %10, 0
br i1 %11, label %22, label %12
12: ; preds = %5, %12
%13 = phi i64 [ %19, %12 ], [ %6, %5 ]
%14 = phi float [ %18, %12 ], [ 0.000000e+00, %5 ]
%15 = phi i32 [ %20, %12 ], [ 0, %5 ]
%16 = getelementptr inbounds float, float* %0, i64 %13
%17 = load float, float* %16
%18 = tail call float @llvm.fmuladd.f32(float %17, float %17, float %14)
%19 = add nsw i64 %13, 1
%20 = add i32 %15, 1
%21 = icmp eq i32 %20, %10
br i1 %21, label %22, label %12
22: ; preds = %12, %5
%23 = phi float [ undef, %5 ], [ %18, %12 ]
%24 = phi i64 [ %6, %5 ], [ %19, %12 ]
%25 = phi float [ 0.000000e+00, %5 ], [ %18, %12 ]
%26 = icmp ult i32 %9, 3
br i1 %26, label %48, label %27
27: ; preds = %22, %27
%28 = phi i64 [ %45, %27 ], [ %24, %22 ]
%29 = phi float [ %44, %27 ], [ %25, %22 ]
%30 = getelementptr inbounds float, float* %0, i64 %28
%31 = load float, float* %30
%32 = tail call float @llvm.fmuladd.f32(float %31, float %31, float %29)
%33 = add nsw i64 %28, 1
%34 = getelementptr inbounds float, float* %0, i64 %33
%35 = load float, float* %34
%36 = tail call float @llvm.fmuladd.f32(float %35, float %35, float %32)
%37 = add nsw i64 %28, 2
%38 = getelementptr inbounds float, float* %0, i64 %37
%39 = load float, float* %38
%40 = tail call float @llvm.fmuladd.f32(float %39, float %39, float %36)
%41 = add nsw i64 %28, 3
%42 = getelementptr inbounds float, float* %0, i64 %41
%43 = load float, float* %42
%44 = tail call float @llvm.fmuladd.f32(float %43, float %43, float %40)
%45 = add nsw i64 %28, 4
%46 = trunc i64 %45 to i32
%47 = icmp eq i32 %7, %46
br i1 %47, label %48, label %27
48: ; preds = %22, %27, %3
%49 = phi float [ 0.000000e+00, %3 ], [ %23, %22 ], [ %44, %27 ]
ret float %49
}