#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 */
typedef int /*<<< orphan*/ z5;
typedef int /*<<< orphan*/ z4;
typedef int /*<<< orphan*/ z3;
typedef int /*<<< orphan*/ z2;
typedef int /*<<< orphan*/ x1;
typedef int /*<<< orphan*/ slopes1;
typedef int /*<<< orphan*/ pz5;
typedef int /*<<< orphan*/ pz4;
typedef int /*<<< orphan*/ pz3;
typedef int /*<<< orphan*/ pz2;
typedef int /*<<< orphan*/ fx5;
typedef int /*<<< orphan*/ fx4;
typedef int /*<<< orphan*/ fx3;
typedef int /*<<< orphan*/ fx2;
typedef int /*<<< orphan*/ fx1;

/* Variables and functions */
int /*<<< orphan*/ DividedDifferences(double *, double *, int, double *);
int /*<<< orphan*/ SimpleNewtonPolynomialInterpolation(double *, double *, int,
                                                       double *, double *, int);
int /*<<< orphan*/ printvdnl(double *, int);
int main() {
  double x1[] = {1.0, 2.0, 3.0, 4.0};
  double fx1[] = {6.0, 9.0, 2.0, 5.0};
  double slopes1[(sizeof(x1) / sizeof(x1[0]))];
  double x2[] = {-4.0, -3.5, -2.5, -1.5, 0.0, 1.5, 2.5, 3.5, 4.0};
  double fx2[] = {-3.0, 1.5, 2.5, 1.5, 0.5, 1.5, 2.5, 1.5, -3.0};
  double z2[] = {-3.0, -2.0, -1.0};
  double pz2[(sizeof(z2) / sizeof(z2[0]))];
  double x3[] = {0.0, 10.0};
  double fx3[] = {0.0, 10.0};
  double z3[] = {0.0, 1.0, 2.0, 3.0, 4.0, 100.0};
  double pz3[(sizeof(z3) / sizeof(z3[0]))];
  double x4[] = {-2.0, -1.0, -0.5, 0.0, 1.0, 2.5, 3.0};
  double fx4[] = {4.0, 1.0, 0.25, 0.0, 1.0, 6.25, 9.0};
  double z4[] = {-1.5, -0.9, -0.8, 0.5, 2.0, 10.0};
  double pz4[(sizeof(z4) / sizeof(z4[0]))];
  double x5[] = {-3.0, -2.0, 2.0, 3.0};
  double fx5[] = {-5.0, -1.1, 1.9, 4.8};
  double z5[] = {-2.5, 0.0, 2.5};
  double pz5[(sizeof(z5) / sizeof(z5[0]))];
  DividedDifferences(x1, fx1, (sizeof(fx1) / sizeof(fx1[0])), slopes1);
  printvdnl(slopes1, (sizeof(slopes1) / sizeof(slopes1[0])));
  SimpleNewtonPolynomialInterpolation(x2, fx2, (sizeof(fx2) / sizeof(fx2[0])),
                                      z2, pz2, (sizeof(pz2) / sizeof(pz2[0])));
  printvdnl(pz2, (sizeof(pz2) / sizeof(pz2[0])));
  SimpleNewtonPolynomialInterpolation(x3, fx3, (sizeof(fx3) / sizeof(fx3[0])),
                                      z3, pz3, (sizeof(pz3) / sizeof(pz3[0])));
  printvdnl(pz3, (sizeof(pz3) / sizeof(pz3[0])));
  SimpleNewtonPolynomialInterpolation(x4, fx4, (sizeof(fx4) / sizeof(fx4[0])),
                                      z4, pz4, (sizeof(pz4) / sizeof(pz4[0])));
  printvdnl(pz4, (sizeof(pz4) / sizeof(pz4[0])));
  SimpleNewtonPolynomialInterpolation(x5, fx5, (sizeof(fx5) / sizeof(fx5[0])),
                                      z5, pz5, (sizeof(pz5) / sizeof(pz5[0])));
  printvdnl(pz5, (sizeof(pz5) / sizeof(pz5[0])));
  return 0;
}


x86 assembly (clang, -O0):
<main>:
push   %rbp
mov    %rsp,%rbp
sub    $0x330,%rsp
movl   $0x0,-0x4(%rbp)
mov    0xf27(%rip),%rax
mov    %rax,-0x30(%rbp)
mov    0xf24(%rip),%rax
mov    %rax,-0x28(%rbp)
mov    0xf21(%rip),%rax
mov    %rax,-0x20(%rbp)
mov    0xf1e(%rip),%rax
mov    %rax,-0x18(%rbp)
mov    0xf1b(%rip),%rax
mov    %rax,-0x50(%rbp)
mov    0xf18(%rip),%rax
mov    %rax,-0x48(%rbp)
mov    0xf15(%rip),%rax
mov    %rax,-0x40(%rbp)
mov    0xf12(%rip),%rax
mov    %rax,-0x38(%rbp)
lea    -0xc0(%rbp),%rdi
lea    0xf48(%rip),%rsi
mov    $0x48,%edx
call   1060 <memcpy@plt>
lea    -0x110(%rbp),%rdi
lea    0xf80(%rip),%rsi
mov    $0x48,%edx
call   1060 <memcpy@plt>
mov    0xfbf(%rip),%rax
mov    %rax,-0x130(%rbp)
mov    0xfb9(%rip),%rax
mov    %rax,-0x128(%rbp)
mov    0xfb3(%rip),%rax
mov    %rax,-0x120(%rbp)
mov    0xfd5(%rip),%rax
mov    %rax,-0x160(%rbp)
mov    0xfcf(%rip),%rax
mov    %rax,-0x158(%rbp)
mov    0xfb9(%rip),%rax
mov    %rax,-0x170(%rbp)
mov    0xfb3(%rip),%rax
mov    %rax,-0x168(%rbp)
lea    -0x1a0(%rbp),%rdi
xor    %esi,%esi
mov    $0x30,%edx
call   1050 <memset@plt>
movsd  0xdd9(%rip),%xmm0
movsd  %xmm0,-0x198(%rbp)
ff
movsd  0xd91(%rip),%xmm0
movsd  %xmm0,-0x190(%rbp)
ff
movsd  0xdd1(%rip),%xmm0
movsd  %xmm0,-0x188(%rbp)
ff
movsd  0xdb9(%rip),%xmm0
movsd  %xmm0,-0x180(%rbp)
ff
movsd  0xdd9(%rip),%xmm0
movsd  %xmm0,-0x178(%rbp)
ff
lea    -0x210(%rbp),%rdi
xor    %esi,%esi
mov    $0x38,%edx
call   1050 <memset@plt>
movsd  0xdae(%rip),%xmm0
movsd  %xmm0,-0x210(%rbp)
ff
movsd  0xd96(%rip),%xmm0
movsd  %xmm0,-0x208(%rbp)
ff
movsd  0xd7e(%rip),%xmm0
movsd  %xmm0,-0x200(%rbp)
ff
movsd  0xd46(%rip),%xmm0
movsd  %xmm0,-0x1f0(%rbp)
ff
movsd  0xd56(%rip),%xmm0
movsd  %xmm0,-0x1e8(%rbp)
ff
movsd  0xd3e(%rip),%xmm0
movsd  %xmm0,-0x1e0(%rbp)
ff
lea    -0x250(%rbp),%rdi
xor    %esi,%esi
mov    $0x38,%edx
call   1050 <memset@plt>
movsd  0xd13(%rip),%xmm0
movsd  %xmm0,-0x250(%rbp)
ff
movsd  0xcf3(%rip),%xmm0
movsd  %xmm0,-0x248(%rbp)
ff
movsd  0xceb(%rip),%xmm0
movsd  %xmm0,-0x240(%rbp)
ff
movsd  0xcd3(%rip),%xmm0
movsd  %xmm0,-0x230(%rbp)
ff
movsd  0xcbb(%rip),%xmm0
movsd  %xmm0,-0x228(%rbp)
ff
movsd  0xca3(%rip),%xmm0
movsd  %xmm0,-0x220(%rbp)
ff
lea    -0x280(%rbp),%rdi
xor    %esi,%esi
mov    $0x30,%edx
call   1050 <memset@plt>
movsd  0xc78(%rip),%xmm0
movsd  %xmm0,-0x280(%rbp)
ff
movsd  0xc60(%rip),%xmm0
movsd  %xmm0,-0x278(%rbp)
ff
movsd  0xc48(%rip),%xmm0
movsd  %xmm0,-0x270(%rbp)
ff
movsd  0xc30(%rip),%xmm0
movsd  %xmm0,-0x268(%rbp)
ff
movsd  0xc18(%rip),%xmm0
movsd  %xmm0,-0x260(%rbp)
ff
movsd  0xc00(%rip),%xmm0
movsd  %xmm0,-0x258(%rbp)
ff
mov    0xcc1(%rip),%rax
mov    %rax,-0x2d0(%rbp)
mov    0xcbb(%rip),%rax
mov    %rax,-0x2c8(%rbp)
mov    0xcb5(%rip),%rax
mov    %rax,-0x2c0(%rbp)
mov    0xcaf(%rip),%rax
mov    %rax,-0x2b8(%rbp)
mov    0xca9(%rip),%rax
mov    %rax,-0x2f0(%rbp)
mov    0xca3(%rip),%rax
mov    %rax,-0x2e8(%rbp)
mov    0xc9d(%rip),%rax
mov    %rax,-0x2e0(%rbp)
mov    0xc97(%rip),%rax
mov    %rax,-0x2d8(%rbp)
mov    0xd51(%rip),%rax
mov    %rax,-0x310(%rbp)
mov    0xd4b(%rip),%rax
mov    %rax,-0x308(%rbp)
mov    0xd45(%rip),%rax
mov    %rax,-0x300(%rbp)
lea    -0x30(%rbp),%rdi
lea    -0x50(%rbp),%rsi
lea    -0x70(%rbp),%rcx
mov    $0x4,%edx
call   1040 <DividedDifferences@plt>
lea    -0x70(%rbp),%rdi
mov    $0x4,%esi
call   1070 <printvdnl@plt>
lea    -0xc0(%rbp),%rdi
lea    -0x110(%rbp),%rsi
lea    -0x130(%rbp),%rcx
lea    -0x150(%rbp),%r8
mov    $0x9,%edx
mov    $0x3,%r9d
call   1030 <SimpleNewtonPolynomialInterpolation@plt>
lea    -0x150(%rbp),%rdi
mov    $0x3,%esi
call   1070 <printvdnl@plt>
lea    -0x160(%rbp),%rdi
lea    -0x170(%rbp),%rsi
lea    -0x1a0(%rbp),%rcx
lea    -0x1d0(%rbp),%r8
mov    $0x2,%edx
mov    $0x6,%r9d
call   1030 <SimpleNewtonPolynomialInterpolation@plt>
lea    -0x1d0(%rbp),%rdi
mov    $0x6,%esi
call   1070 <printvdnl@plt>
lea    -0x210(%rbp),%rdi
lea    -0x250(%rbp),%rsi
lea    -0x280(%rbp),%rcx
lea    -0x2b0(%rbp),%r8
mov    $0x7,%edx
mov    $0x6,%r9d
call   1030 <SimpleNewtonPolynomialInterpolation@plt>
lea    -0x2b0(%rbp),%rdi
mov    $0x6,%esi
call   1070 <printvdnl@plt>
lea    -0x2d0(%rbp),%rdi
lea    -0x2f0(%rbp),%rsi
lea    -0x310(%rbp),%rcx
lea    -0x330(%rbp),%r8
mov    $0x4,%edx
mov    $0x3,%r9d
call   1030 <SimpleNewtonPolynomialInterpolation@plt>
lea    -0x330(%rbp),%rdi
mov    $0x3,%esi
call   1070 <printvdnl@plt>
xor    %eax,%eax
add    $0x330,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<main>:
push   %rbx
sub    $0x320,%rsp
movaps 0xf71(%rip),%xmm0
movaps %xmm0,0x180(%rsp)
movaps 0xf52(%rip),%xmm0
movaps %xmm0,0x170(%rsp)
movaps 0xf73(%rip),%xmm0
movaps %xmm0,0x160(%rsp)
movaps 0xf54(%rip),%xmm0
movaps %xmm0,0x150(%rsp)
mov    0xfe5(%rip),%rax
mov    %rax,0x250(%rsp)
movaps 0xfc6(%rip),%xmm0
movaps %xmm0,0x240(%rsp)
movaps 0xfa7(%rip),%xmm0
movaps %xmm0,0x230(%rsp)
movaps 0xf88(%rip),%xmm0
movaps %xmm0,0x220(%rsp)
movaps 0xf69(%rip),%xmm0
movaps %xmm0,0x210(%rsp)
mov    0xfea(%rip),%rax
mov    %rax,0x200(%rsp)
movaps 0xfcb(%rip),%xmm0
movaps %xmm0,0x1f0(%rsp)
movaps 0xfac(%rip),%xmm0
movaps %xmm0,0x1e0(%rsp)
movaps 0xf8d(%rip),%xmm0
movaps %xmm0,0x1d0(%rsp)
movaps 0xf6e(%rip),%xmm0
movaps %xmm0,0x1c0(%rsp)
mov    0xfbf(%rip),%rax
mov    %rax,0x50(%rsp)
movaps 0xfa3(%rip),%xmm0
movaps %xmm0,0x40(%rsp)
movaps 0xe67(%rip),%xmm0
movaps %xmm0,0x30(%rsp)
movaps %xmm0,0x20(%rsp)
xorps  %xmm0,%xmm0
movaps %xmm0,0xe0(%rsp)
movaps 0xdbb(%rip),%xmm1
movups %xmm1,0xe8(%rsp)
movaps 0xdbc(%rip),%xmm1
movups %xmm1,0xf8(%rsp)
movabs $0x4059000000000000,%rax
mov    %rax,0x108(%rsp)
movaps %xmm0,0xb0(%rsp)
movaps 0xda3(%rip),%xmm1
movaps %xmm1,0xa0(%rsp)
movabs $0xbfe0000000000000,%rax
mov    %rax,0xb0(%rsp)
movaps 0xd92(%rip),%xmm1
movaps %xmm1,0xc0(%rsp)
movabs $0x4008000000000000,%rax
mov    %rax,0xd0(%rsp)
movaps %xmm0,0x70(%rsp)
movaps 0xd7c(%rip),%xmm0
movaps %xmm0,0x60(%rsp)
movabs $0x3fd0000000000000,%rax
mov    %rax,0x70(%rsp)
movaps 0xd71(%rip),%xmm0
movaps %xmm0,0x80(%rsp)
movabs $0x4022000000000000,%rax
mov    %rax,0x90(%rsp)
movaps 0xd60(%rip),%xmm0
movaps %xmm0,0x190(%rsp)
movaps 0xd61(%rip),%xmm0
movaps %xmm0,0x1a0(%rsp)
movaps 0xd62(%rip),%xmm0
movaps %xmm0,0x1b0(%rsp)
movaps 0xdc3(%rip),%xmm0
movaps %xmm0,0x140(%rsp)
movaps 0xda4(%rip),%xmm0
movaps %xmm0,0x130(%rsp)
movaps 0xdc5(%rip),%xmm0
movaps %xmm0,0x120(%rsp)
movaps 0xda6(%rip),%xmm0
movaps %xmm0,0x110(%rsp)
mov    0xe87(%rip),%rax
mov    %rax,0x10(%rsp)
movaps 0xe6b(%rip),%xmm0
movaps %xmm0,(%rsp)
lea    0x170(%rsp),%rdi
lea    0x150(%rsp),%rsi
lea    0x2a0(%rsp),%rbx
mov    $0x4,%edx
mov    %rbx,%rcx
call   1040 <DividedDifferences@plt>
mov    %rbx,%rdi
mov    $0x4,%esi
call   1050 <printvdnl@plt>
lea    0x210(%rsp),%rdi
lea    0x1c0(%rsp),%rsi
lea    0x40(%rsp),%rcx
lea    0x280(%rsp),%rbx
mov    $0x9,%edx
mov    %rbx,%r8
mov    $0x3,%r9d
call   1030 <SimpleNewtonPolynomialInterpolation@plt>
mov    %rbx,%rdi
mov    $0x3,%esi
call   1050 <printvdnl@plt>
lea    0x30(%rsp),%rdi
lea    0x20(%rsp),%rsi
lea    0xe0(%rsp),%rcx
lea    0x2f0(%rsp),%rbx
mov    $0x2,%edx
mov    %rbx,%r8
mov    $0x6,%r9d
call   1030 <SimpleNewtonPolynomialInterpolation@plt>
mov    %rbx,%rdi
mov    $0x6,%esi
call   1050 <printvdnl@plt>
lea    0xa0(%rsp),%rdi
lea    0x60(%rsp),%rsi
lea    0x190(%rsp),%rcx
lea    0x2c0(%rsp),%rbx
mov    $0x7,%edx
mov    %rbx,%r8
mov    $0x6,%r9d
call   1030 <SimpleNewtonPolynomialInterpolation@plt>
mov    %rbx,%rdi
mov    $0x6,%esi
call   1050 <printvdnl@plt>
lea    0x130(%rsp),%rdi
lea    0x110(%rsp),%rsi
mov    %rsp,%rcx
lea    0x260(%rsp),%rbx
mov    $0x4,%edx
mov    %rbx,%r8
mov    $0x3,%r9d
call   1030 <SimpleNewtonPolynomialInterpolation@plt>
mov    %rbx,%rdi
mov    $0x3,%esi
call   1050 <printvdnl@plt>
xor    %eax,%eax
add    $0x320,%rsp
pop    %rbx
ret

LLVM IR (-O0):
@__const.main.x1 = [4 x double] [double 1.000000e+00, double 2.000000e+00, double 3.000000e+00, double 4.000000e+00]
@__const.main.fx1 = [4 x double] [double 6.000000e+00, double 9.000000e+00, double 2.000000e+00, double 5.000000e+00]
@__const.main.x2 = [9 x double] [double -4.000000e+00, double -3.500000e+00, double -2.500000e+00, double -1.500000e+00, double 0.000000e+00, double 1.500000e+00, double 2.500000e+00, double 3.500000e+00, double 4.000000e+00]
@__const.main.fx2 = [9 x double] [double -3.000000e+00, double 1.500000e+00, double 2.500000e+00, double 1.500000e+00, double 5.000000e-01, double 1.500000e+00, double 2.500000e+00, double 1.500000e+00, double -3.000000e+00]
@__const.main.z2 = [3 x double] [double -3.000000e+00, double -2.000000e+00, double -1.000000e+00]
@__const.main.x3 = [2 x double] [double 0.000000e+00, double 1.000000e+01]
@__const.main.fx3 = [2 x double] [double 0.000000e+00, double 1.000000e+01]
@__const.main.x5 = [4 x double] [double -3.000000e+00, double -2.000000e+00, double 2.000000e+00, double 3.000000e+00]
@__const.main.fx5 = [4 x double] [double -5.000000e+00, double -1.100000e+00, double 1.900000e+00, double 4.800000e+00]
@__const.main.z5 = [3 x double] [double -2.500000e+00, double 0.000000e+00, double 2.500000e+00]
define i32 @main() {
%1 = alloca i32
%2 = alloca [4 x double]
%3 = alloca [4 x double]
%4 = alloca [4 x double]
%5 = alloca [9 x double]
%6 = alloca [9 x double]
%7 = alloca [3 x double]
%8 = alloca [3 x double]
%9 = alloca [2 x double]
%10 = alloca [2 x double]
%11 = alloca [6 x double]
%12 = alloca [6 x double]
%13 = alloca [7 x double]
%14 = alloca [7 x double]
%15 = alloca [6 x double]
%16 = alloca [6 x double]
%17 = alloca [4 x double]
%18 = alloca [4 x double]
%19 = alloca [3 x double]
%20 = alloca [3 x double]
store i32 0, i32* %1
%21 = bitcast [4 x double]* %2 to i8*
%22 = bitcast [4 x double]* %3 to i8*
%23 = bitcast [9 x double]* %5 to i8*
%24 = bitcast [9 x double]* %6 to i8*
%25 = bitcast [3 x double]* %7 to i8*
%26 = bitcast [2 x double]* %9 to i8*
%27 = bitcast [2 x double]* %10 to i8*
%28 = bitcast [6 x double]* %11 to i8*
%29 = bitcast i8* %28 to [6 x double]*
%30 = getelementptr inbounds [6 x double], [6 x double]* %29, i32 0, i32 1
store double 1.000000e+00, double* %30
%31 = getelementptr inbounds [6 x double], [6 x double]* %29, i32 0, i32 2
store double 2.000000e+00, double* %31
%32 = getelementptr inbounds [6 x double], [6 x double]* %29, i32 0, i32 3
store double 3.000000e+00, double* %32
%33 = getelementptr inbounds [6 x double], [6 x double]* %29, i32 0, i32 4
store double 4.000000e+00, double* %33
%34 = getelementptr inbounds [6 x double], [6 x double]* %29, i32 0, i32 5
store double 1.000000e+02, double* %34
%35 = bitcast [7 x double]* %13 to i8*
%36 = bitcast i8* %35 to [7 x double]*
%37 = getelementptr inbounds [7 x double], [7 x double]* %36, i32 0, i32 0
store double -2.000000e+00, double* %37
%38 = getelementptr inbounds [7 x double], [7 x double]* %36, i32 0, i32 1
store double -1.000000e+00, double* %38
%39 = getelementptr inbounds [7 x double], [7 x double]* %36, i32 0, i32 2
store double -5.000000e-01, double* %39
%40 = getelementptr inbounds [7 x double], [7 x double]* %36, i32 0, i32 4
store double 1.000000e+00, double* %40
%41 = getelementptr inbounds [7 x double], [7 x double]* %36, i32 0, i32 5
store double 2.500000e+00, double* %41
%42 = getelementptr inbounds [7 x double], [7 x double]* %36, i32 0, i32 6
store double 3.000000e+00, double* %42
%43 = bitcast [7 x double]* %14 to i8*
%44 = bitcast i8* %43 to [7 x double]*
%45 = getelementptr inbounds [7 x double], [7 x double]* %44, i32 0, i32 0
store double 4.000000e+00, double* %45
%46 = getelementptr inbounds [7 x double], [7 x double]* %44, i32 0, i32 1
store double 1.000000e+00, double* %46
%47 = getelementptr inbounds [7 x double], [7 x double]* %44, i32 0, i32 2
store double 2.500000e-01, double* %47
%48 = getelementptr inbounds [7 x double], [7 x double]* %44, i32 0, i32 4
store double 1.000000e+00, double* %48
%49 = getelementptr inbounds [7 x double], [7 x double]* %44, i32 0, i32 5
store double 6.250000e+00, double* %49
%50 = getelementptr inbounds [7 x double], [7 x double]* %44, i32 0, i32 6
store double 9.000000e+00, double* %50
%51 = bitcast [6 x double]* %15 to i8*
%52 = bitcast i8* %51 to [6 x double]*
%53 = getelementptr inbounds [6 x double], [6 x double]* %52, i32 0, i32 0
store double -1.500000e+00, double* %53
%54 = getelementptr inbounds [6 x double], [6 x double]* %52, i32 0, i32 1
store double -9.000000e-01, double* %54
%55 = getelementptr inbounds [6 x double], [6 x double]* %52, i32 0, i32 2
store double -8.000000e-01, double* %55
%56 = getelementptr inbounds [6 x double], [6 x double]* %52, i32 0, i32 3
store double 5.000000e-01, double* %56
%57 = getelementptr inbounds [6 x double], [6 x double]* %52, i32 0, i32 4
store double 2.000000e+00, double* %57
%58 = getelementptr inbounds [6 x double], [6 x double]* %52, i32 0, i32 5
store double 1.000000e+01, double* %58
%59 = bitcast [4 x double]* %17 to i8*
%60 = bitcast [4 x double]* %18 to i8*
%61 = bitcast [3 x double]* %19 to i8*
%62 = getelementptr inbounds [4 x double], [4 x double]* %2, i64 0, i64 0
%63 = getelementptr inbounds [4 x double], [4 x double]* %3, i64 0, i64 0
%64 = getelementptr inbounds [4 x double], [4 x double]* %4, i64 0, i64 0
%65 = call i32 @DividedDifferences(double* %62, double* %63, i32 4, double* %64)
%66 = getelementptr inbounds [4 x double], [4 x double]* %4, i64 0, i64 0
%67 = call i32 @printvdnl(double* %66, i32 4)
%68 = getelementptr inbounds [9 x double], [9 x double]* %5, i64 0, i64 0
%69 = getelementptr inbounds [9 x double], [9 x double]* %6, i64 0, i64 0
%70 = getelementptr inbounds [3 x double], [3 x double]* %7, i64 0, i64 0
%71 = getelementptr inbounds [3 x double], [3 x double]* %8, i64 0, i64 0
%72 = call i32 @SimpleNewtonPolynomialInterpolation(double* %68, double* %69, i32 9, double* %70, double* %71, i32 3)
%73 = getelementptr inbounds [3 x double], [3 x double]* %8, i64 0, i64 0
%74 = call i32 @printvdnl(double* %73, i32 3)
%75 = getelementptr inbounds [2 x double], [2 x double]* %9, i64 0, i64 0
%76 = getelementptr inbounds [2 x double], [2 x double]* %10, i64 0, i64 0
%77 = getelementptr inbounds [6 x double], [6 x double]* %11, i64 0, i64 0
%78 = getelementptr inbounds [6 x double], [6 x double]* %12, i64 0, i64 0
%79 = call i32 @SimpleNewtonPolynomialInterpolation(double* %75, double* %76, i32 2, double* %77, double* %78, i32 6)
%80 = getelementptr inbounds [6 x double], [6 x double]* %12, i64 0, i64 0
%81 = call i32 @printvdnl(double* %80, i32 6)
%82 = getelementptr inbounds [7 x double], [7 x double]* %13, i64 0, i64 0
%83 = getelementptr inbounds [7 x double], [7 x double]* %14, i64 0, i64 0
%84 = getelementptr inbounds [6 x double], [6 x double]* %15, i64 0, i64 0
%85 = getelementptr inbounds [6 x double], [6 x double]* %16, i64 0, i64 0
%86 = call i32 @SimpleNewtonPolynomialInterpolation(double* %82, double* %83, i32 7, double* %84, double* %85, i32 6)
%87 = getelementptr inbounds [6 x double], [6 x double]* %16, i64 0, i64 0
%88 = call i32 @printvdnl(double* %87, i32 6)
%89 = getelementptr inbounds [4 x double], [4 x double]* %17, i64 0, i64 0
%90 = getelementptr inbounds [4 x double], [4 x double]* %18, i64 0, i64 0
%91 = getelementptr inbounds [3 x double], [3 x double]* %19, i64 0, i64 0
%92 = getelementptr inbounds [3 x double], [3 x double]* %20, i64 0, i64 0
%93 = call i32 @SimpleNewtonPolynomialInterpolation(double* %89, double* %90, i32 4, double* %91, double* %92, i32 3)
%94 = getelementptr inbounds [3 x double], [3 x double]* %20, i64 0, i64 0
%95 = call i32 @printvdnl(double* %94, i32 3)
ret i32 0
}

LLVM IR (-O2):
@__const.main.x1 = [4 x double] [double 1.000000e+00, double 2.000000e+00, double 3.000000e+00, double 4.000000e+00]
@__const.main.fx1 = [4 x double] [double 6.000000e+00, double 9.000000e+00, double 2.000000e+00, double 5.000000e+00]
@__const.main.x2 = [9 x double] [double -4.000000e+00, double -3.500000e+00, double -2.500000e+00, double -1.500000e+00, double 0.000000e+00, double 1.500000e+00, double 2.500000e+00, double 3.500000e+00, double 4.000000e+00]
@__const.main.fx2 = [9 x double] [double -3.000000e+00, double 1.500000e+00, double 2.500000e+00, double 1.500000e+00, double 5.000000e-01, double 1.500000e+00, double 2.500000e+00, double 1.500000e+00, double -3.000000e+00]
@__const.main.z2 = [3 x double] [double -3.000000e+00, double -2.000000e+00, double -1.000000e+00]
@__const.main.fx3 = [2 x double] [double 0.000000e+00, double 1.000000e+01]
@__const.main.x5 = [4 x double] [double -3.000000e+00, double -2.000000e+00, double 2.000000e+00, double 3.000000e+00]
@__const.main.fx5 = [4 x double] [double -5.000000e+00, double -1.100000e+00, double 1.900000e+00, double 4.800000e+00]
@__const.main.z5 = [3 x double] [double -2.500000e+00, double 0.000000e+00, double 2.500000e+00]
define i32 @main() {
%1 = alloca [4 x double]
%2 = alloca [4 x double]
%3 = alloca [4 x double]
%4 = alloca [9 x double]
%5 = alloca [9 x double]
%6 = alloca [3 x double]
%7 = alloca [3 x double]
%8 = alloca [2 x double]
%9 = alloca [2 x double]
%10 = alloca [6 x double]
%11 = alloca [6 x double]
%12 = alloca [7 x double]
%13 = alloca [7 x double]
%14 = alloca [6 x double]
%15 = alloca [6 x double]
%16 = alloca [4 x double]
%17 = alloca [4 x double]
%18 = alloca [3 x double]
%19 = alloca [3 x double]
%20 = bitcast [4 x double]* %1 to i8*
%21 = bitcast [4 x double]* %2 to i8*
%22 = bitcast [4 x double]* %3 to i8*
%23 = bitcast [9 x double]* %4 to i8*
%24 = bitcast [9 x double]* %5 to i8*
%25 = bitcast [3 x double]* %6 to i8*
%26 = bitcast [3 x double]* %7 to i8*
%27 = bitcast [2 x double]* %8 to i8*
%28 = bitcast [2 x double]* %9 to i8*
%29 = bitcast [6 x double]* %10 to i8*
%30 = getelementptr inbounds [6 x double], [6 x double]* %10, i64 0, i64 1
%31 = bitcast double* %30 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 2.000000e+00>, <2 x double>* %31
%32 = getelementptr inbounds [6 x double], [6 x double]* %10, i64 0, i64 3
%33 = bitcast double* %32 to <2 x double>*
store <2 x double> <double 3.000000e+00, double 4.000000e+00>, <2 x double>* %33
%34 = getelementptr inbounds [6 x double], [6 x double]* %10, i64 0, i64 5
store double 1.000000e+02, double* %34
%35 = bitcast [6 x double]* %11 to i8*
%36 = bitcast [7 x double]* %12 to i8*
%37 = getelementptr inbounds [7 x double], [7 x double]* %12, i64 0, i64 2
%38 = bitcast double* %37 to i8*
%39 = getelementptr inbounds [7 x double], [7 x double]* %12, i64 0, i64 0
%40 = bitcast [7 x double]* %12 to <2 x double>*
store <2 x double> <double -2.000000e+00, double -1.000000e+00>, <2 x double>* %40
%41 = getelementptr inbounds [7 x double], [7 x double]* %12, i64 0, i64 2
store double -5.000000e-01, double* %41
%42 = getelementptr inbounds [7 x double], [7 x double]* %12, i64 0, i64 4
%43 = bitcast double* %42 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 2.500000e+00>, <2 x double>* %43
%44 = getelementptr inbounds [7 x double], [7 x double]* %12, i64 0, i64 6
store double 3.000000e+00, double* %44
%45 = bitcast [7 x double]* %13 to i8*
%46 = getelementptr inbounds [7 x double], [7 x double]* %13, i64 0, i64 2
%47 = bitcast double* %46 to i8*
%48 = getelementptr inbounds [7 x double], [7 x double]* %13, i64 0, i64 0
%49 = bitcast [7 x double]* %13 to <2 x double>*
store <2 x double> <double 4.000000e+00, double 1.000000e+00>, <2 x double>* %49
%50 = getelementptr inbounds [7 x double], [7 x double]* %13, i64 0, i64 2
store double 2.500000e-01, double* %50
%51 = getelementptr inbounds [7 x double], [7 x double]* %13, i64 0, i64 4
%52 = bitcast double* %51 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 6.250000e+00>, <2 x double>* %52
%53 = getelementptr inbounds [7 x double], [7 x double]* %13, i64 0, i64 6
store double 9.000000e+00, double* %53
%54 = bitcast [6 x double]* %14 to i8*
%55 = getelementptr inbounds [6 x double], [6 x double]* %14, i64 0, i64 0
%56 = bitcast [6 x double]* %14 to <2 x double>*
store <2 x double> <double -1.500000e+00, double -9.000000e-01>, <2 x double>* %56
%57 = getelementptr inbounds [6 x double], [6 x double]* %14, i64 0, i64 2
%58 = bitcast double* %57 to <2 x double>*
store <2 x double> <double -8.000000e-01, double 5.000000e-01>, <2 x double>* %58
%59 = getelementptr inbounds [6 x double], [6 x double]* %14, i64 0, i64 4
%60 = bitcast double* %59 to <2 x double>*
store <2 x double> <double 2.000000e+00, double 1.000000e+01>, <2 x double>* %60
%61 = bitcast [6 x double]* %15 to i8*
%62 = bitcast [4 x double]* %16 to i8*
%63 = bitcast [4 x double]* %17 to i8*
%64 = bitcast [3 x double]* %18 to i8*
%65 = bitcast [3 x double]* %19 to i8*
%66 = getelementptr inbounds [4 x double], [4 x double]* %1, i64 0, i64 0
%67 = getelementptr inbounds [4 x double], [4 x double]* %2, i64 0, i64 0
%68 = getelementptr inbounds [4 x double], [4 x double]* %3, i64 0, i64 0
%69 = call i32 @DividedDifferences(double* %66, double* %67, i32 4, double* %68)
%70 = call i32 @printvdnl(double* %68, i32 4)
%71 = getelementptr inbounds [9 x double], [9 x double]* %4, i64 0, i64 0
%72 = getelementptr inbounds [9 x double], [9 x double]* %5, i64 0, i64 0
%73 = getelementptr inbounds [3 x double], [3 x double]* %6, i64 0, i64 0
%74 = getelementptr inbounds [3 x double], [3 x double]* %7, i64 0, i64 0
%75 = call i32 @SimpleNewtonPolynomialInterpolation(double* %71, double* %72, i32 9, double* %73, double* %74, i32 3)
%76 = call i32 @printvdnl(double* %74, i32 3)
%77 = getelementptr inbounds [2 x double], [2 x double]* %8, i64 0, i64 0
%78 = getelementptr inbounds [2 x double], [2 x double]* %9, i64 0, i64 0
%79 = getelementptr inbounds [6 x double], [6 x double]* %10, i64 0, i64 0
%80 = getelementptr inbounds [6 x double], [6 x double]* %11, i64 0, i64 0
%81 = call i32 @SimpleNewtonPolynomialInterpolation(double* %77, double* %78, i32 2, double* %79, double* %80, i32 6)
%82 = call i32 @printvdnl(double* %80, i32 6)
%83 = getelementptr inbounds [6 x double], [6 x double]* %15, i64 0, i64 0
%84 = call i32 @SimpleNewtonPolynomialInterpolation(double* %39, double* %48, i32 7, double* %55, double* %83, i32 6)
%85 = call i32 @printvdnl(double* %83, i32 6)
%86 = getelementptr inbounds [4 x double], [4 x double]* %16, i64 0, i64 0
%87 = getelementptr inbounds [4 x double], [4 x double]* %17, i64 0, i64 0
%88 = getelementptr inbounds [3 x double], [3 x double]* %18, i64 0, i64 0
%89 = getelementptr inbounds [3 x double], [3 x double]* %19, i64 0, i64 0
%90 = call i32 @SimpleNewtonPolynomialInterpolation(double* %86, double* %87, i32 4, double* %88, double* %89, i32 3)
%91 = call i32 @printvdnl(double* %89, i32 3)
ret i32 0
}