#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 */
float cos(float);
float sin(int);
void sinpoly(float x, float p[], int ndeg) {
  int j;
  float cosx;
  p[1] = 1.0;
  p[2] = sin(x / 229.0);
  if (ndeg > 2)
    p[3] = sin(2 * x / 229.0);
  if (ndeg > 3) {
    cosx = cos(x / 229.0);
    for (j = 3; j <= ndeg; j++) {
      p[j] = 2 * cosx * p[j - 1] - p[j - 2];
    }
  }
}


x86 assembly (clang, -O0):
<sinpoly>:
push   %rbp
mov    %rsp,%rbp
sub    $0x20,%rsp
movss  %xmm0,-0x4(%rbp)
mov    %rdi,-0x10(%rbp)
mov    %esi,-0x14(%rbp)
mov    -0x10(%rbp),%rax
movss  0xec8(%rip),%xmm0
movss  %xmm0,0x4(%rax)
movss  -0x4(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xeaa(%rip),%xmm1
divsd  %xmm1,%xmm0
cvttsd2si %xmm0,%edi
call   1030 <sin@plt>
mov    -0x10(%rbp),%rax
movss  %xmm0,0x8(%rax)
cmpl   $0x2,-0x14(%rbp)
jle    11a5 <sinpoly+0x85>
movss  0xe8e(%rip),%xmm0
mulss  -0x4(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xe71(%rip),%xmm1
divsd  %xmm1,%xmm0
cvttsd2si %xmm0,%edi
call   1030 <sin@plt>
mov    -0x10(%rbp),%rax
movss  %xmm0,0xc(%rax)
cmpl   $0x3,-0x14(%rbp)
jle    124b <sinpoly+0x12b>
movss  -0x4(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xe40(%rip),%xmm1
divsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
call   1040 <cos@plt>
movss  %xmm0,-0x1c(%rbp)
movl   $0x3,-0x18(%rbp)
mov    -0x18(%rbp),%eax
cmp    -0x14(%rbp),%eax
jg     1246 <sinpoly+0x126>
movss  0xe1f(%rip),%xmm0
mulss  -0x1c(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x18(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm2
mov    -0x10(%rbp),%rax
mov    -0x18(%rbp),%ecx
sub    $0x2,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm1
movd   %xmm1,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x18(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x18(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x18(%rbp)
jmp    11d9 <sinpoly+0xb9>
jmp    124b <sinpoly+0x12b>
add    $0x20,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<sinpoly>:
push   %rbp
push   %rbx
sub    $0x18,%rsp
mov    %esi,%ebp
mov    %rdi,%rbx
movl   $0x3f800000,0x4(%rdi)
movss  %xmm0,0xc(%rsp)
cvtss2sd %xmm0,%xmm0
divsd  0xebc(%rip),%xmm0
movsd  %xmm0,0x10(%rsp)
cvttsd2si %xmm0,%edi
call   1030 <sin@plt>
movss  %xmm0,0x8(%rbx)
cmp    $0x3,%ebp
jl     11ee <sinpoly+0xce>
movss  0xc(%rsp),%xmm0
addss  %xmm0,%xmm0
cvtss2sd %xmm0,%xmm0
divsd  0xe89(%rip),%xmm0
cvttsd2si %xmm0,%edi
call   1030 <sin@plt>
movss  %xmm0,0xc(%rbx)
cmp    $0x3,%ebp
je     11ee <sinpoly+0xce>
movsd  0x10(%rsp),%xmm0
cvtsd2ss %xmm0,%xmm0
call   1040 <cos@plt>
addss  %xmm0,%xmm0
movss  0x8(%rbx),%xmm1
mov    %ebp,%ecx
lea    -0x2(%rcx),%rdx
and    $0xfffffffffffffffe,%rdx
xor    %eax,%eax
xchg   %ax,%ax
mulss  %xmm0,%xmm1
subss  0x4(%rbx,%rax,4),%xmm1
movss  %xmm1,0xc(%rbx,%rax,4)
mulss  %xmm0,%xmm1
subss  0x8(%rbx,%rax,4),%xmm1
movss  %xmm1,0x10(%rbx,%rax,4)
add    $0x2,%rax
cmp    %rax,%rdx
jne    11b0 <sinpoly+0x90>
test   $0x1,%cl
je     11ee <sinpoly+0xce>
mulss  %xmm1,%xmm0
subss  0x4(%rbx,%rax,4),%xmm0
movss  %xmm0,0xc(%rbx,%rax,4)
add    $0x18,%rsp
pop    %rbx
pop    %rbp
ret

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

LLVM IR (-O2):
define void @sinpoly(float %0, float* nocapture %1, i32 %2) {
%4 = getelementptr inbounds float, float* %1, i64 1
store float 1.000000e+00, float* %4
%5 = fpext float %0 to double
%6 = fdiv double %5, 2.290000e+02
%7 = fptosi double %6 to i32
%8 = tail call float @sin(i32 %7)
%9 = getelementptr inbounds float, float* %1, i64 2
store float %8, float* %9
%10 = icmp sgt i32 %2, 2
br i1 %10, label %11, label %57
11: ; preds = %3
%12 = fmul float %0, 2.000000e+00
%13 = fpext float %12 to double
%14 = fdiv double %13, 2.290000e+02
%15 = fptosi double %14 to i32
%16 = tail call float @sin(i32 %15)
%17 = getelementptr inbounds float, float* %1, i64 3
store float %16, float* %17
%18 = icmp eq i32 %2, 3
br i1 %18, label %57, label %19
19: ; preds = %11
%20 = fptrunc double %6 to float
%21 = tail call float @cos(float %20)
%22 = fmul float %21, 2.000000e+00
%23 = load float, float* %9
%24 = zext i32 %2 to i64
%25 = and i64 %24, 1
%26 = add nsw i64 %24, -2
%27 = and i64 %26, -2
br label %28
28: ; preds = %28, %19
%29 = phi float [ %23, %19 ], [ %43, %28 ]
%30 = phi i64 [ 3, %19 ], [ %45, %28 ]
%31 = phi i64 [ 0, %19 ], [ %46, %28 ]
%32 = add nsw i64 %30, -2
%33 = getelementptr inbounds float, float* %1, i64 %32
%34 = load float, float* %33
%35 = fneg float %34
%36 = tail call float @llvm.fmuladd.f32(float %22, float %29, float %35)
%37 = getelementptr inbounds float, float* %1, i64 %30
store float %36, float* %37
%38 = add nuw nsw i64 %30, 1
%39 = add nsw i64 %30, -1
%40 = getelementptr inbounds float, float* %1, i64 %39
%41 = load float, float* %40
%42 = fneg float %41
%43 = tail call float @llvm.fmuladd.f32(float %22, float %36, float %42)
%44 = getelementptr inbounds float, float* %1, i64 %38
store float %43, float* %44
%45 = add nuw nsw i64 %30, 2
%46 = add i64 %31, 2
%47 = icmp eq i64 %46, %27
br i1 %47, label %48, label %28
48: ; preds = %28
%49 = icmp eq i64 %25, 0
br i1 %49, label %57, label %50
50: ; preds = %48
%51 = add nsw i64 %45, -2
%52 = getelementptr inbounds float, float* %1, i64 %51
%53 = load float, float* %52
%54 = fneg float %53
%55 = tail call float @llvm.fmuladd.f32(float %22, float %43, float %54)
%56 = getelementptr inbounds float, float* %1, i64 %45
store float %55, float* %56
br label %57
57: ; preds = %50, %48, %3, %11
ret void
}