float lpc_aejw(float *a, float w, int p) {
  int i;
  float c_re, c_im;
  float cs, sn, tmp;
  if (p == 0)
    return (1.);
  cs = (float)cos((double)w);
  sn = -(float)sin((double)w);
  c_re = cs * a[p];
  c_im = sn * a[p];
  for (i = p - 1; i > 0; i--) {
    c_re += a[i];
    c_im = cs * (tmp = c_im) + sn * c_re;
    c_re = cs * c_re - sn * tmp;
  }
  c_re += 1.0;
  return (((c_re) * (c_re)) + ((c_im) * (c_im)));
}


x86 assembly (clang, -O0):
<lpc_aejw>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%rsp
mov    %rdi,-0x10(%rbp)
movss  %xmm0,-0x14(%rbp)
mov    %esi,-0x18(%rbp)
cmpl   $0x0,-0x18(%rbp)
jne    1150 <lpc_aejw+0x30>
movss  0xec2(%rip),%xmm0
movss  %xmm0,-0x4(%rbp)
jmp    1284 <lpc_aejw+0x164>
movss  -0x14(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
call   1040 <cos@plt>
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x28(%rbp)
movss  -0x14(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
call   1030 <sin@plt>
cvtsd2ss %xmm0,%xmm0
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
movss  %xmm0,-0x2c(%rbp)
movss  -0x28(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x18(%rbp),%rcx
mulss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x20(%rbp)
movss  -0x2c(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x18(%rbp),%rcx
mulss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x24(%rbp)
mov    -0x18(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x1c(%rbp)
cmpl   $0x0,-0x1c(%rbp)
jle    1245 <lpc_aejw+0x125>
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
addss  -0x20(%rbp),%xmm0
movss  %xmm0,-0x20(%rbp)
movss  -0x28(%rbp),%xmm0
movss  -0x24(%rbp),%xmm2
movss  %xmm2,-0x30(%rbp)
movss  -0x2c(%rbp),%xmm1
mulss  -0x20(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x24(%rbp)
movss  -0x28(%rbp),%xmm0
movss  -0x20(%rbp),%xmm2
movss  -0x2c(%rbp),%xmm1
mulss  -0x30(%rbp),%xmm1
movd   %xmm1,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x20(%rbp)
mov    -0x1c(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x1c(%rbp)
jmp    11c2 <lpc_aejw+0xa2>
movss  -0x20(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xdaa(%rip),%xmm1
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x20(%rbp)
movss  -0x20(%rbp),%xmm0
movss  -0x20(%rbp),%xmm2
movss  -0x24(%rbp),%xmm1
mulss  -0x24(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x4(%rbp)
movss  -0x4(%rbp),%xmm0
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<lpc_aejw>:
push   %rbp
push   %rbx
sub    $0x18,%rsp
test   %esi,%esi
je     11c5 <lpc_aejw+0xa5>
mov    %esi,%ebp
mov    %rdi,%rbx
cvtss2sd %xmm0,%xmm0
movsd  %xmm0,0x10(%rsp)
call   1040 <cos@plt>
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,0xc(%rsp)
movsd  0x10(%rsp),%xmm0
call   1030 <sin@plt>
movss  0xc(%rsp),%xmm6
cvtsd2ss %xmm0,%xmm1
movaps 0xe98(%rip),%xmm2
xorps  %xmm1,%xmm2
movslq %ebp,%rax
movss  (%rbx,%rax,4),%xmm3
movaps %xmm3,%xmm0
mulss  %xmm6,%xmm0
mulss  %xmm2,%xmm3
cmp    $0x2,%ebp
jl     11cf <lpc_aejw+0xaf>
mov    %ebp,%eax
add    $0x1,%rax
nopl   0x0(%rax)
lea    -0x2(%rax),%ecx
addss  (%rbx,%rcx,4),%xmm0
movaps %xmm0,%xmm5
mulss  %xmm2,%xmm5
movaps %xmm6,%xmm4
mulss  %xmm3,%xmm4
addss  %xmm5,%xmm4
mulss  %xmm1,%xmm3
mulss  %xmm6,%xmm0
addss  %xmm3,%xmm0
add    $0xffffffffffffffff,%rax
movaps %xmm4,%xmm3
cmp    $0x2,%rax
ja     1190 <lpc_aejw+0x70>
jmp    11d2 <lpc_aejw+0xb2>
movss  0xe43(%rip),%xmm0
jmp    11e6 <lpc_aejw+0xc6>
movaps %xmm3,%xmm4
addss  0xe36(%rip),%xmm0
mulss  %xmm4,%xmm4
mulss  %xmm0,%xmm0
addss  %xmm4,%xmm0
add    $0x18,%rsp
pop    %rbx
pop    %rbp
ret

LLVM IR (-O0):
define float @lpc_aejw(float* %0, float %1, i32 %2) {
%4 = alloca float
%5 = alloca float*
%6 = alloca float
%7 = alloca i32
%8 = alloca i32
%9 = alloca float
%10 = alloca float
%11 = alloca float
%12 = alloca float
%13 = alloca float
store float* %0, float** %5
store float %1, float* %6
store i32 %2, i32* %7
%14 = load i32, i32* %7
%15 = icmp eq i32 %14, 0
br i1 %15, label %16, label %17
16: ; preds = %3
store float 1.000000e+00, float* %4
br label %81
17: ; preds = %3
%18 = load float, float* %6
%19 = fpext float %18 to double
%20 = call double @cos(double %19)
%21 = fptrunc double %20 to float
store float %21, float* %11
%22 = load float, float* %6
%23 = fpext float %22 to double
%24 = call double @sin(double %23)
%25 = fptrunc double %24 to float
%26 = fneg float %25
store float %26, float* %12
%27 = load float, float* %11
%28 = load float*, float** %5
%29 = load i32, i32* %7
%30 = sext i32 %29 to i64
%31 = getelementptr inbounds float, float* %28, i64 %30
%32 = load float, float* %31
%33 = fmul float %27, %32
store float %33, float* %9
%34 = load float, float* %12
%35 = load float*, float** %5
%36 = load i32, i32* %7
%37 = sext i32 %36 to i64
%38 = getelementptr inbounds float, float* %35, i64 %37
%39 = load float, float* %38
%40 = fmul float %34, %39
store float %40, float* %10
%41 = load i32, i32* %7
%42 = sub nsw i32 %41, 1
store i32 %42, i32* %8
br label %43
43: ; preds = %67, %17
%44 = load i32, i32* %8
%45 = icmp sgt i32 %44, 0
br i1 %45, label %46, label %70
46: ; preds = %43
%47 = load float*, float** %5
%48 = load i32, i32* %8
%49 = sext i32 %48 to i64
%50 = getelementptr inbounds float, float* %47, i64 %49
%51 = load float, float* %50
%52 = load float, float* %9
%53 = fadd float %52, %51
store float %53, float* %9
%54 = load float, float* %11
%55 = load float, float* %10
store float %55, float* %13
%56 = load float, float* %12
%57 = load float, float* %9
%58 = fmul float %56, %57
%59 = call float @llvm.fmuladd.f32(float %54, float %55, float %58)
store float %59, float* %10
%60 = load float, float* %11
%61 = load float, float* %9
%62 = load float, float* %12
%63 = load float, float* %13
%64 = fmul float %62, %63
%65 = fneg float %64
%66 = call float @llvm.fmuladd.f32(float %60, float %61, float %65)
store float %66, float* %9
br label %67
67: ; preds = %46
%68 = load i32, i32* %8
%69 = add nsw i32 %68, -1
store i32 %69, i32* %8
br label %43
70: ; preds = %43
%71 = load float, float* %9
%72 = fpext float %71 to double
%73 = fadd double %72, 1.000000e+00
%74 = fptrunc double %73 to float
store float %74, float* %9
%75 = load float, float* %9
%76 = load float, float* %9
%77 = load float, float* %10
%78 = load float, float* %10
%79 = fmul float %77, %78
%80 = call float @llvm.fmuladd.f32(float %75, float %76, float %79)
store float %80, float* %4
br label %81
81: ; preds = %70, %16
%82 = load float, float* %4
ret float %82
}

LLVM IR (-O2):
define float @lpc_aejw(float* nocapture %0, float %1, i32 %2) {
%4 = icmp eq i32 %2, 0
br i1 %4, label %40, label %5
5: ; preds = %3
%6 = fpext float %1 to double
%7 = tail call double @cos(double %6)
%8 = fptrunc double %7 to float
%9 = tail call double @sin(double %6)
%10 = fptrunc double %9 to float
%11 = fneg float %10
%12 = sext i32 %2 to i64
%13 = getelementptr inbounds float, float* %0, i64 %12
%14 = load float, float* %13
%15 = fmul float %14, %8
%16 = fmul float %14, %11
%17 = icmp sgt i32 %2, 1
br i1 %17, label %18, label %34
18: ; preds = %5
%19 = zext i32 %2 to i64
br label %20
20: ; preds = %18, %20
%21 = phi i64 [ %19, %18 ], [ %24, %20 ]
%22 = phi float [ %15, %18 ], [ %32, %20 ]
%23 = phi float [ %16, %18 ], [ %30, %20 ]
%24 = add nsw i64 %21, -1
%25 = and i64 %24, 4294967295
%26 = getelementptr inbounds float, float* %0, i64 %25
%27 = load float, float* %26
%28 = fadd float %22, %27
%29 = fmul float %28, %11
%30 = tail call float @llvm.fmuladd.f32(float %8, float %23, float %29)
%31 = fmul float %23, %10
%32 = tail call float @llvm.fmuladd.f32(float %8, float %28, float %31)
%33 = icmp ugt i64 %21, 2
br i1 %33, label %20, label %34
34: ; preds = %20, %5
%35 = phi float [ %16, %5 ], [ %30, %20 ]
%36 = phi float [ %15, %5 ], [ %32, %20 ]
%37 = fadd float %36, 1.000000e+00
%38 = fmul float %35, %35
%39 = tail call float @llvm.fmuladd.f32(float %37, float %37, float %38)
br label %40
40: ; preds = %3, %34
%41 = phi float [ %39, %34 ], [ 1.000000e+00, %3 ]
ret float %41
}