int lpc_bwex(float *a, float *aw, float gamma, int p) {
  int i;
  float gk;
  for (i = 1, gk = gamma; i <= p; i++, gk *= gamma)
    aw[i] = a[i] * gk;
  return (0);
}


x86 assembly (clang, -O0):
<lpc_bwex>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
movss  %xmm0,-0x14(%rbp)
mov    %edx,-0x18(%rbp)
movl   $0x1,-0x1c(%rbp)
movss  -0x14(%rbp),%xmm0
movss  %xmm0,-0x20(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x18(%rbp),%eax
jg     116d <lpc_bwex+0x6d>
mov    -0x8(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mulss  -0x20(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
movss  -0x14(%rbp),%xmm0
mulss  -0x20(%rbp),%xmm0
movss  %xmm0,-0x20(%rbp)
jmp    1125 <lpc_bwex+0x25>
xor    %eax,%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<lpc_bwex>:
test   %edx,%edx
jle    11bb <lpc_bwex+0xbb>
mov    %edx,%edx
lea    -0x1(%rdx),%r8
mov    %edx,%eax
and    $0x3,%eax
mov    $0x1,%ecx
movaps %xmm0,%xmm1
cmp    $0x3,%r8
jb     118d <lpc_bwex+0x8d>
and    $0xfffffffc,%edx
xor    %ecx,%ecx
movaps %xmm0,%xmm1
nopl   0x0(%rax)
movss  0x4(%rdi,%rcx,4),%xmm2
mulss  %xmm1,%xmm2
movss  %xmm2,0x4(%rsi,%rcx,4)
mulss  %xmm0,%xmm1
movss  0x8(%rdi,%rcx,4),%xmm2
mulss  %xmm1,%xmm2
movss  %xmm2,0x8(%rsi,%rcx,4)
mulss  %xmm0,%xmm1
movss  0xc(%rdi,%rcx,4),%xmm2
mulss  %xmm1,%xmm2
movss  %xmm2,0xc(%rsi,%rcx,4)
mulss  %xmm0,%xmm1
movss  0x10(%rdi,%rcx,4),%xmm2
mulss  %xmm1,%xmm2
movss  %xmm2,0x10(%rsi,%rcx,4)
mulss  %xmm0,%xmm1
add    $0x4,%rcx
cmp    %rcx,%rdx
jne    1130 <lpc_bwex+0x30>
add    $0x1,%rcx
test   %rax,%rax
je     11bb <lpc_bwex+0xbb>
lea    (%rsi,%rcx,4),%rdx
lea    (%rdi,%rcx,4),%rcx
xor    %esi,%esi
nopl   0x0(%rax)
movss  (%rcx,%rsi,4),%xmm2
mulss  %xmm1,%xmm2
movss  %xmm2,(%rdx,%rsi,4)
mulss  %xmm0,%xmm1
add    $0x1,%rsi
cmp    %rsi,%rax
jne    11a0 <lpc_bwex+0xa0>
xor    %eax,%eax
ret

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

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