#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*/ memcpy(double *, double *, int);
void poly_coefs(double *x, double *beta, int k, double *phi) {
  int j;
  int ell;
  memcpy(phi, beta, (k + 1) * sizeof(double));
  for (j = 1; j <= k; j++) {
    for (ell = k; ell >= j; ell--) {
      phi[ell] = (phi[ell] - phi[ell - 1]) / (x[ell] - x[ell - j]);
    }
  }
}


x86 assembly (clang, -O0):
<poly_coefs>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %rcx,-0x20(%rbp)
mov    -0x20(%rbp),%rdi
mov    -0x10(%rbp),%rsi
mov    -0x14(%rbp),%eax
add    $0x1,%eax
cltq
shl    $0x3,%rax
mov    %eax,%edx
call   1030 <memcpy@plt>
movl   $0x1,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x14(%rbp),%eax
jg     11d7 <poly_coefs+0xc7>
mov    -0x14(%rbp),%eax
mov    %eax,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x24(%rbp),%eax
jl     11c4 <poly_coefs+0xb4>
mov    -0x20(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x28(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x8(%rbp),%rax
mov    -0x28(%rbp),%ecx
sub    -0x24(%rbp),%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm1
divsd  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
movslq -0x28(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x28(%rbp)
jmp    115b <poly_coefs+0x4b>
jmp    11c9 <poly_coefs+0xb9>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1149 <poly_coefs+0x39>
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<poly_coefs>:
push   %r15
push   %r14
push   %rbx
mov    %rcx,%rbx
mov    %edx,%r14d
mov    %rdi,%r15
lea    0x8(,%r14,8),%edx
mov    %rcx,%rdi
call   1030 <memcpy@plt>
test   %r14d,%r14d
jle    119b <poly_coefs+0x8b>
mov    %r14d,%eax
add    $0x1,%r14d
lea    0x1(%rax),%rcx
mov    $0x1,%edx
mov    %r15,%rsi
cs nopw 0x0(%rax,%rax,1)
movsd  (%rbx,%rax,8),%xmm0
mov    %rcx,%rdi
nopl   0x0(%rax,%rax,1)
movsd  -0x10(%rbx,%rdi,8),%xmm1
subsd  %xmm1,%xmm0
movsd  -0x8(%r15,%rdi,8),%xmm2
subsd  -0x10(%rsi,%rdi,8),%xmm2
divsd  %xmm2,%xmm0
movsd  %xmm0,-0x8(%rbx,%rdi,8)
add    $0xffffffffffffffff,%rdi
movapd %xmm1,%xmm0
cmp    %rdx,%rdi
jg     1160 <poly_coefs+0x50>
add    $0x1,%rdx
add    $0xfffffffffffffff8,%rsi
cmp    %r14,%rdx
jne    1150 <poly_coefs+0x40>
pop    %rbx
pop    %r14
pop    %r15
ret

LLVM IR (-O0):
define void @poly_coefs(double* %0, double* %1, i32 %2, double* %3) {
%5 = alloca double*
%6 = alloca double*
%7 = alloca i32
%8 = alloca double*
%9 = alloca i32
%10 = alloca i32
store double* %0, double** %5
store double* %1, double** %6
store i32 %2, i32* %7
store double* %3, double** %8
%11 = load double*, double** %8
%12 = load double*, double** %6
%13 = load i32, i32* %7
%14 = add nsw i32 %13, 1
%15 = sext i32 %14 to i64
%16 = mul i64 %15, 8
%17 = trunc i64 %16 to i32
%18 = call i32 @memcpy(double* %11, double* %12, i32 %17)
store i32 1, i32* %9
br label %19
19: ; preds = %64, %4
%20 = load i32, i32* %9
%21 = load i32, i32* %7
%22 = icmp sle i32 %20, %21
br i1 %22, label %23, label %67
23: ; preds = %19
%24 = load i32, i32* %7
store i32 %24, i32* %10
br label %25
25: ; preds = %60, %23
%26 = load i32, i32* %10
%27 = load i32, i32* %9
%28 = icmp sge i32 %26, %27
br i1 %28, label %29, label %63
29: ; preds = %25
%30 = load double*, double** %8
%31 = load i32, i32* %10
%32 = sext i32 %31 to i64
%33 = getelementptr inbounds double, double* %30, i64 %32
%34 = load double, double* %33
%35 = load double*, double** %8
%36 = load i32, i32* %10
%37 = sub nsw i32 %36, 1
%38 = sext i32 %37 to i64
%39 = getelementptr inbounds double, double* %35, i64 %38
%40 = load double, double* %39
%41 = fsub double %34, %40
%42 = load double*, double** %5
%43 = load i32, i32* %10
%44 = sext i32 %43 to i64
%45 = getelementptr inbounds double, double* %42, i64 %44
%46 = load double, double* %45
%47 = load double*, double** %5
%48 = load i32, i32* %10
%49 = load i32, i32* %9
%50 = sub nsw i32 %48, %49
%51 = sext i32 %50 to i64
%52 = getelementptr inbounds double, double* %47, i64 %51
%53 = load double, double* %52
%54 = fsub double %46, %53
%55 = fdiv double %41, %54
%56 = load double*, double** %8
%57 = load i32, i32* %10
%58 = sext i32 %57 to i64
%59 = getelementptr inbounds double, double* %56, i64 %58
store double %55, double* %59
br label %60
60: ; preds = %29
%61 = load i32, i32* %10
%62 = add nsw i32 %61, -1
store i32 %62, i32* %10
br label %25
63: ; preds = %25
br label %64
64: ; preds = %63
%65 = load i32, i32* %9
%66 = add nsw i32 %65, 1
store i32 %66, i32* %9
br label %19
67: ; preds = %19
ret void
}

LLVM IR (-O2):
define void @poly_coefs(double* nocapture %0, double* %1, i32 %2, double* %3) {
%5 = shl i32 %2, 3
%6 = add i32 %5, 8
%7 = tail call i32 @memcpy(double* %3, double* %1, i32 %6)
%8 = icmp slt i32 %2, 1
br i1 %8, label %36, label %9
9: ; preds = %4
%10 = zext i32 %2 to i64
%11 = add nuw i32 %2, 1
%12 = zext i32 %11 to i64
%13 = getelementptr inbounds double, double* %3, i64 %10
br label %14
14: ; preds = %33, %9
%15 = phi i64 [ 1, %9 ], [ %34, %33 ]
%16 = load double, double* %13
br label %17
17: ; preds = %14, %17
%18 = phi double [ %16, %14 ], [ %23, %17 ]
%19 = phi i64 [ %10, %14 ], [ %21, %17 ]
%20 = getelementptr inbounds double, double* %3, i64 %19
%21 = add nsw i64 %19, -1
%22 = getelementptr inbounds double, double* %3, i64 %21
%23 = load double, double* %22
%24 = fsub double %18, %23
%25 = getelementptr inbounds double, double* %0, i64 %19
%26 = load double, double* %25
%27 = sub nsw i64 %19, %15
%28 = getelementptr inbounds double, double* %0, i64 %27
%29 = load double, double* %28
%30 = fsub double %26, %29
%31 = fdiv double %24, %30
store double %31, double* %20
%32 = icmp sgt i64 %19, %15
br i1 %32, label %17, label %33
33: ; preds = %17
%34 = add nuw nsw i64 %15, 1
%35 = icmp eq i64 %34, %12
br i1 %35, label %36, label %14
36: ; preds = %33, %4
ret void
}