#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 */
int /*<<< orphan*/ memset(float *, int /*<<< orphan*/, unsigned int);
void make_poly_from_roots(const double *roots, unsigned num_roots,
                          float *poly) {
  poly[0] = 1;
  poly[1] = -roots[0];
  memset(poly + 2, 0, (num_roots + 1 - 2) * sizeof(*poly));
  for (unsigned i = 1; i < num_roots; i++)
    for (unsigned j = num_roots; j > 0; j--)
      poly[j] -= poly[j - 1] * roots[i];
}


x86 assembly (clang, -O0):
<make_poly_from_roots>:
push   %rbp
mov    %rsp,%rbp
sub    $0x20,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %rdx,-0x18(%rbp)
mov    -0x18(%rbp),%rax
movss  0xed1(%rip),%xmm0
movss  %xmm0,(%rax)
mov    -0x8(%rbp),%rax
movsd  (%rax),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
cvtsd2ss %xmm0,%xmm0
mov    -0x18(%rbp),%rax
movss  %xmm0,0x4(%rax)
mov    -0x18(%rbp),%rdi
add    $0x8,%rdi
mov    -0xc(%rbp),%eax
add    $0x1,%eax
sub    $0x2,%eax
mov    %eax,%eax
shl    $0x2,%rax
mov    %eax,%edx
xor    %esi,%esi
call   1030 <memset@plt>
movl   $0x1,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jae    121c <make_poly_from_roots+0x10c>
mov    -0xc(%rbp),%eax
mov    %eax,-0x20(%rbp)
cmpl   $0x0,-0x20(%rbp)
jbe    1209 <make_poly_from_roots+0xf9>
mov    -0x18(%rbp),%rax
mov    -0x20(%rbp),%ecx
sub    $0x1,%ecx
mov    %ecx,%ecx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
mov    -0x8(%rbp),%rax
mov    -0x1c(%rbp),%ecx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x18(%rbp),%rax
mov    -0x20(%rbp),%ecx
movss  (%rax,%rcx,4),%xmm1
cvtss2sd %xmm1,%xmm1
movq   %xmm0,%rdx
movabs $0x8000000000000000,%rsi
xor    %rsi,%rdx
movq   %rdx,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x20(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x20(%rbp)
jmp    1198 <make_poly_from_roots+0x88>
jmp    120e <make_poly_from_roots+0xfe>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1186 <make_poly_from_roots+0x76>
add    $0x20,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<make_poly_from_roots>:
push   %r15
push   %r14
push   %rbx
mov    %rdx,%rbx
mov    %rdi,%r14
movl   $0x3f800000,(%rdx)
movsd  (%rdi),%xmm0
cvtsd2ss %xmm0,%xmm0
xorps  0xed0(%rip),%xmm0
mov    %esi,%r15d
movss  %xmm0,0x4(%rdx)
lea    0x8(%rdx),%rdi
lea    -0x4(,%r15,4),%edx
ff
xor    %esi,%esi
call   1030 <memset@plt>
cmp    $0x2,%r15d
jb     1222 <make_poly_from_roots+0x112>
mov    %r15d,%eax
lea    -0x1(%rax),%rcx
mov    $0x1,%edx
jmp    117d <make_poly_from_roots+0x6d>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
add    $0x1,%rdx
cmp    %rax,%rdx
je     1222 <make_poly_from_roots+0x112>
movsd  (%r14,%rdx,8),%xmm0
movss  (%rbx,%rax,4),%xmm1
test   $0x1,%al
jne    11a0 <make_poly_from_roots+0x90>
mov    %rax,%rsi
cmp    $0x1,%r15d
jne    11d0 <make_poly_from_roots+0xc0>
jmp    1170 <make_poly_from_roots+0x60>
nopw   0x0(%rax,%rax,1)
cvtss2sd %xmm1,%xmm1
movss  -0x4(%rbx,%rax,4),%xmm2
xorps  %xmm3,%xmm3
cvtss2sd %xmm2,%xmm3
mulsd  %xmm0,%xmm3
subsd  %xmm3,%xmm1
cvtsd2ss %xmm1,%xmm1
movss  %xmm1,(%rbx,%rax,4)
movaps %xmm2,%xmm1
mov    %rcx,%rsi
cmp    $0x1,%r15d
je     1170 <make_poly_from_roots+0x60>
xchg   %ax,%ax
cvtss2sd %xmm1,%xmm1
movss  -0x8(%rbx,%rsi,4),%xmm2
movss  -0x4(%rbx,%rsi,4),%xmm3
cvtss2sd %xmm3,%xmm3
movaps %xmm3,%xmm4
mulsd  %xmm0,%xmm4
subsd  %xmm4,%xmm1
cvtsd2ss %xmm1,%xmm1
movss  %xmm1,(%rbx,%rsi,4)
xorps  %xmm1,%xmm1
cvtss2sd %xmm2,%xmm1
mulsd  %xmm0,%xmm1
subsd  %xmm1,%xmm3
xorps  %xmm1,%xmm1
cvtsd2ss %xmm3,%xmm1
movss  %xmm1,-0x4(%rbx,%rsi,4)
movaps %xmm2,%xmm1
add    $0xfffffffffffffffe,%rsi
jne    11d0 <make_poly_from_roots+0xc0>
jmp    1170 <make_poly_from_roots+0x60>
pop    %rbx
pop    %r14
pop    %r15
ret

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

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