void silk_bwexpander_FLP(float *ar, const int d, const float chirp) {
  int i;
  float cfac = chirp;
  for (i = 0; i < d - 1; i++) {
    ar[i] *= cfac;
    cfac *= chirp;
  }
  ar[d - 1] *= cfac;
}


x86 assembly (clang, -O0):
<silk_bwexpander_FLP>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
movss  %xmm0,-0x10(%rbp)
movss  -0x10(%rbp),%xmm0
movss  %xmm0,-0x18(%rbp)
movl   $0x0,-0x14(%rbp)
mov    -0x14(%rbp),%eax
mov    -0xc(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    1166 <silk_bwexpander_FLP+0x66>
movss  -0x18(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x14(%rbp),%rcx
mulss  (%rax,%rcx,4),%xmm0
movss  %xmm0,(%rax,%rcx,4)
movss  -0x10(%rbp),%xmm0
mulss  -0x18(%rbp),%xmm0
movss  %xmm0,-0x18(%rbp)
mov    -0x14(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x14(%rbp)
jmp    1121 <silk_bwexpander_FLP+0x21>
movss  -0x18(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0xc(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mulss  (%rax,%rcx,4),%xmm0
movss  %xmm0,(%rax,%rcx,4)
pop    %rbp
ret

x86 assembly (clang, -O2):
<silk_bwexpander_FLP>:
lea    -0x1(%rsi),%eax
cmp    $0x2,%esi
jl     1120 <silk_bwexpander_FLP+0x20>
mov    %eax,%esi
lea    -0x1(%rsi),%rdx
mov    %esi,%ecx
and    $0x3,%ecx
cmp    $0x3,%rdx
jae    1128 <silk_bwexpander_FLP+0x28>
xor    %edx,%edx
movaps %xmm0,%xmm1
jmp    1187 <silk_bwexpander_FLP+0x87>
movaps %xmm0,%xmm1
jmp    11bb <silk_bwexpander_FLP+0xbb>
and    $0xfffffffc,%esi
xor    %edx,%edx
movaps %xmm0,%xmm1
movss  (%rdi,%rdx,4),%xmm2
mulss  %xmm1,%xmm2
movss  %xmm2,(%rdi,%rdx,4)
mulss  %xmm0,%xmm1
movss  0x4(%rdi,%rdx,4),%xmm2
mulss  %xmm1,%xmm2
movss  %xmm2,0x4(%rdi,%rdx,4)
mulss  %xmm0,%xmm1
movss  0x8(%rdi,%rdx,4),%xmm2
mulss  %xmm1,%xmm2
movss  %xmm2,0x8(%rdi,%rdx,4)
mulss  %xmm0,%xmm1
movss  0xc(%rdi,%rdx,4),%xmm2
mulss  %xmm1,%xmm2
movss  %xmm2,0xc(%rdi,%rdx,4)
mulss  %xmm0,%xmm1
add    $0x4,%rdx
cmp    %rdx,%rsi
jne    1130 <silk_bwexpander_FLP+0x30>
test   %rcx,%rcx
je     11bb <silk_bwexpander_FLP+0xbb>
lea    (%rdi,%rdx,4),%rdx
xor    %esi,%esi
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movss  (%rdx,%rsi,4),%xmm2
mulss  %xmm1,%xmm2
movss  %xmm2,(%rdx,%rsi,4)
mulss  %xmm0,%xmm1
add    $0x1,%rsi
cmp    %rsi,%rcx
jne    11a0 <silk_bwexpander_FLP+0xa0>
cltq
mulss  (%rdi,%rax,4),%xmm1
movss  %xmm1,(%rdi,%rax,4)
ret

LLVM IR (-O0):
define void @silk_bwexpander_FLP(float* %0, i32 %1, float %2) {
%4 = alloca float*
%5 = alloca i32
%6 = alloca float
%7 = alloca i32
%8 = alloca float
store float* %0, float** %4
store i32 %1, i32* %5
store float %2, float* %6
%9 = load float, float* %6
store float %9, float* %8
store i32 0, i32* %7
br label %10
10: ; preds = %26, %3
%11 = load i32, i32* %7
%12 = load i32, i32* %5
%13 = sub nsw i32 %12, 1
%14 = icmp slt i32 %11, %13
br i1 %14, label %15, label %29
15: ; preds = %10
%16 = load float, float* %8
%17 = load float*, float** %4
%18 = load i32, i32* %7
%19 = sext i32 %18 to i64
%20 = getelementptr inbounds float, float* %17, i64 %19
%21 = load float, float* %20
%22 = fmul float %21, %16
store float %22, float* %20
%23 = load float, float* %6
%24 = load float, float* %8
%25 = fmul float %24, %23
store float %25, float* %8
br label %26
26: ; preds = %15
%27 = load i32, i32* %7
%28 = add nsw i32 %27, 1
store i32 %28, i32* %7
br label %10
29: ; preds = %10
%30 = load float, float* %8
%31 = load float*, float** %4
%32 = load i32, i32* %5
%33 = sub nsw i32 %32, 1
%34 = sext i32 %33 to i64
%35 = getelementptr inbounds float, float* %31, i64 %34
%36 = load float, float* %35
%37 = fmul float %36, %30
store float %37, float* %35
ret void
}

LLVM IR (-O2):
define void @silk_bwexpander_FLP(float* nocapture %0, i32 %1, float %2) {
%4 = add i32 %1, -1
%5 = icmp sgt i32 %1, 1
br i1 %5, label %6, label %55
6: ; preds = %3
%7 = zext i32 %4 to i64
%8 = add nsw i64 %7, -1
%9 = and i64 %7, 3
%10 = icmp ult i64 %8, 3
br i1 %10, label %39, label %11
11: ; preds = %6
%12 = and i64 %7, 4294967292
br label %13
13: ; preds = %13, %11
%14 = phi i64 [ 0, %11 ], [ %36, %13 ]
%15 = phi float [ %2, %11 ], [ %35, %13 ]
%16 = phi i64 [ 0, %11 ], [ %37, %13 ]
%17 = getelementptr inbounds float, float* %0, i64 %14
%18 = load float, float* %17
%19 = fmul float %15, %18
store float %19, float* %17
%20 = fmul float %15, %2
%21 = or i64 %14, 1
%22 = getelementptr inbounds float, float* %0, i64 %21
%23 = load float, float* %22
%24 = fmul float %20, %23
store float %24, float* %22
%25 = fmul float %20, %2
%26 = or i64 %14, 2
%27 = getelementptr inbounds float, float* %0, i64 %26
%28 = load float, float* %27
%29 = fmul float %25, %28
store float %29, float* %27
%30 = fmul float %25, %2
%31 = or i64 %14, 3
%32 = getelementptr inbounds float, float* %0, i64 %31
%33 = load float, float* %32
%34 = fmul float %30, %33
store float %34, float* %32
%35 = fmul float %30, %2
%36 = add nuw nsw i64 %14, 4
%37 = add i64 %16, 4
%38 = icmp eq i64 %37, %12
br i1 %38, label %39, label %13
39: ; preds = %13, %6
%40 = phi float [ undef, %6 ], [ %35, %13 ]
%41 = phi i64 [ 0, %6 ], [ %36, %13 ]
%42 = phi float [ %2, %6 ], [ %35, %13 ]
%43 = icmp eq i64 %9, 0
br i1 %43, label %55, label %44
44: ; preds = %39, %44
%45 = phi i64 [ %52, %44 ], [ %41, %39 ]
%46 = phi float [ %51, %44 ], [ %42, %39 ]
%47 = phi i64 [ %53, %44 ], [ 0, %39 ]
%48 = getelementptr inbounds float, float* %0, i64 %45
%49 = load float, float* %48
%50 = fmul float %46, %49
store float %50, float* %48
%51 = fmul float %46, %2
%52 = add nuw nsw i64 %45, 1
%53 = add i64 %47, 1
%54 = icmp eq i64 %53, %9
br i1 %54, label %55, label %44
55: ; preds = %39, %44, %3
%56 = phi float [ %2, %3 ], [ %40, %39 ], [ %51, %44 ]
%57 = sext i32 %4 to i64
%58 = getelementptr inbounds float, float* %0, i64 %57
%59 = load float, float* %58
%60 = fmul float %56, %59
store float %60, float* %58
ret void
}