#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 get_weight(int, int, int, int, float *);
int /*<<< orphan*/ printf(char *, int /*<<< orphan*/, int /*<<< orphan*/,
                          int /*<<< orphan*/, int);
void fornberg_weights(int M, int p, int *alpha, float *delta) {
  int N = p;
  float c1, c2, c3;
  int min;
  int i, j, n, v, m;
  delta[(N + 1) * (N + 1) * 0 + (N + 1) * 0 + 0] = 1;
  c1 = 1;
  printf("delta %d %d %d=%f\n", 0, 0, 0, get_weight(N, 0, 0, 0, delta));
  for (n = 1; n <= N; n++) {
    c2 = 1;
    for (v = 0; v <= n - 1; v++) {
      c3 = alpha[n] - alpha[v];
      c2 = c2 * c3;
      if (n <= M) {
        delta[(N + 1) * (N + 1) * n + (N + 1) * (n - 1) + v] = 0;
      }
      if (n < M)
        min = n;
      else
        min = M;
      for (m = 0; m <= min; m++) {
        if (m == 0)
          delta[(N + 1) * (N + 1) * m + (N + 1) * n + v] =
              (alpha[n] * get_weight(N, m, n - 1, v, delta)) / c3;
        else
          delta[(N + 1) * (N + 1) * m + (N + 1) * n + v] =
              (alpha[n] * get_weight(N, m, n - 1, v, delta) -
               m * get_weight(N, m - 1, n - 1, v, delta)) /
              c3;
      }
    }
    for (m = 0; m <= min; m++) {
      if (m == 0)
        delta[(N + 1) * (N + 1) * m + (N + 1) * n + n] =
            (c1 / c2) *
            (-1 * alpha[n - 1] * get_weight(N, m, n - 1, n - 1, delta));
      else
        delta[(N + 1) * (N + 1) * m + (N + 1) * n + n] =
            (c1 / c2) * (m * get_weight(N, m - 1, n - 1, n - 1, delta) -
                         alpha[n - 1] * get_weight(N, m, n - 1, n - 1, delta));
    }
    c1 = c2;
  }
}


x86 assembly (clang, -O0):
<fornberg_weights>:
push   %rbp
mov    %rsp,%rbp
sub    $0x70,%rsp
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %rdx,-0x10(%rbp)
mov    %rcx,-0x18(%rbp)
mov    -0x8(%rbp),%eax
mov    %eax,-0x1c(%rbp)
mov    -0x18(%rbp),%rax
mov    -0x1c(%rbp),%ecx
add    $0x1,%ecx
mov    -0x1c(%rbp),%edx
add    $0x1,%edx
imul   %edx,%ecx
imul   $0x0,%ecx,%ecx
mov    -0x1c(%rbp),%edx
add    $0x1,%edx
imul   $0x0,%edx,%edx
add    %edx,%ecx
add    $0x0,%ecx
movslq %ecx,%rcx
movss  0xe95(%rip),%xmm0
movss  %xmm0,(%rax,%rcx,4)
movss  0xe88(%rip),%xmm0
movss  %xmm0,-0x20(%rbp)
mov    -0x1c(%rbp),%edi
mov    -0x18(%rbp),%r8
xor    %ecx,%ecx
mov    %ecx,%esi
mov    %ecx,%edx
call   1030 <get_weight@plt>
mov    %eax,%r8d
lea    0xe6b(%rip),%rdi
xor    %ecx,%ecx
mov    %ecx,%esi
mov    %ecx,%edx
call   1040 <printf@plt>
movl   $0x1,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jg     1531 <fornberg_weights+0x411>
movss  0xe41(%rip),%xmm0
movss  %xmm0,-0x24(%rbp)
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
mov    -0x38(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jg     13ac <fornberg_weights+0x28c>
mov    -0x10(%rbp),%rax
movslq -0x38(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    -0x10(%rbp),%rcx
movslq -0x3c(%rbp),%rdx
sub    (%rcx,%rdx,4),%eax
cvtsi2ss %eax,%xmm0
movss  %xmm0,-0x28(%rbp)
movss  -0x24(%rbp),%xmm0
mulss  -0x28(%rbp),%xmm0
movss  %xmm0,-0x24(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x4(%rbp),%eax
jg     124c <fornberg_weights+0x12c>
mov    -0x18(%rbp),%rax
mov    -0x1c(%rbp),%ecx
add    $0x1,%ecx
mov    -0x1c(%rbp),%edx
add    $0x1,%edx
imul   %edx,%ecx
imul   -0x38(%rbp),%ecx
mov    -0x1c(%rbp),%edx
add    $0x1,%edx
mov    -0x38(%rbp),%esi
sub    $0x1,%esi
imul   %esi,%edx
add    %edx,%ecx
add    -0x3c(%rbp),%ecx
movslq %ecx,%rcx
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x38(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    1263 <fornberg_weights+0x143>
mov    -0x38(%rbp),%eax
mov    %eax,-0x2c(%rbp)
jmp    1269 <fornberg_weights+0x149>
mov    -0x4(%rbp),%eax
mov    %eax,-0x2c(%rbp)
movl   $0x0,-0x40(%rbp)
mov    -0x40(%rbp),%eax
cmp    -0x2c(%rbp),%eax
jg     1399 <fornberg_weights+0x279>
cmpl   $0x0,-0x40(%rbp)
jne    12f0 <fornberg_weights+0x1d0>
mov    -0x10(%rbp),%rax
movslq -0x38(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x44(%rbp)
mov    -0x1c(%rbp),%edi
mov    -0x40(%rbp),%esi
mov    -0x38(%rbp),%edx
sub    $0x1,%edx
mov    -0x3c(%rbp),%ecx
mov    -0x18(%rbp),%r8
call   1030 <get_weight@plt>
mov    %eax,%ecx
mov    -0x44(%rbp),%eax
imul   %ecx,%eax
cvtsi2ss %eax,%xmm0
divss  -0x28(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x1c(%rbp),%ecx
add    $0x1,%ecx
mov    -0x1c(%rbp),%edx
add    $0x1,%edx
imul   %edx,%ecx
imul   -0x40(%rbp),%ecx
mov    -0x1c(%rbp),%edx
add    $0x1,%edx
imul   -0x38(%rbp),%edx
add    %edx,%ecx
add    -0x3c(%rbp),%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
jmp    1386 <fornberg_weights+0x266>
mov    -0x10(%rbp),%rax
movslq -0x38(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x50(%rbp)
mov    -0x1c(%rbp),%edi
mov    -0x40(%rbp),%esi
mov    -0x38(%rbp),%edx
sub    $0x1,%edx
mov    -0x3c(%rbp),%ecx
mov    -0x18(%rbp),%r8
call   1030 <get_weight@plt>
mov    %eax,%ecx
mov    -0x50(%rbp),%eax
imul   %ecx,%eax
mov    %eax,-0x48(%rbp)
mov    -0x40(%rbp),%eax
mov    %eax,-0x4c(%rbp)
mov    -0x1c(%rbp),%edi
mov    -0x40(%rbp),%esi
sub    $0x1,%esi
mov    -0x38(%rbp),%edx
sub    $0x1,%edx
mov    -0x3c(%rbp),%ecx
mov    -0x18(%rbp),%r8
call   1030 <get_weight@plt>
mov    -0x4c(%rbp),%ecx
mov    %eax,%edx
mov    -0x48(%rbp),%eax
imul   %edx,%ecx
sub    %ecx,%eax
cvtsi2ss %eax,%xmm0
divss  -0x28(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x1c(%rbp),%ecx
add    $0x1,%ecx
mov    -0x1c(%rbp),%edx
add    $0x1,%edx
imul   %edx,%ecx
imul   -0x40(%rbp),%ecx
mov    -0x1c(%rbp),%edx
add    $0x1,%edx
imul   -0x38(%rbp),%edx
add    %edx,%ecx
add    -0x3c(%rbp),%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
jmp    138b <fornberg_weights+0x26b>
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    1270 <fornberg_weights+0x150>
jmp    139e <fornberg_weights+0x27e>
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    11cb <fornberg_weights+0xab>
movl   $0x0,-0x40(%rbp)
mov    -0x40(%rbp),%eax
cmp    -0x2c(%rbp),%eax
jg     1519 <fornberg_weights+0x3f9>
cmpl   $0x0,-0x40(%rbp)
jne    1452 <fornberg_weights+0x332>
movss  -0x20(%rbp),%xmm0
divss  -0x24(%rbp),%xmm0
movss  %xmm0,-0x58(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x38(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
imul   $0xffffffff,(%rax,%rcx,4),%eax
mov    %eax,-0x54(%rbp)
mov    -0x1c(%rbp),%edi
mov    -0x40(%rbp),%esi
mov    -0x38(%rbp),%edx
sub    $0x1,%edx
mov    -0x38(%rbp),%ecx
sub    $0x1,%ecx
mov    -0x18(%rbp),%r8
call   1030 <get_weight@plt>
movss  -0x58(%rbp),%xmm0
mov    %eax,%ecx
mov    -0x54(%rbp),%eax
imul   %ecx,%eax
cvtsi2ss %eax,%xmm1
mulss  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
mov    -0x1c(%rbp),%ecx
add    $0x1,%ecx
mov    -0x1c(%rbp),%edx
add    $0x1,%edx
imul   %edx,%ecx
imul   -0x40(%rbp),%ecx
mov    -0x1c(%rbp),%edx
add    $0x1,%edx
imul   -0x38(%rbp),%edx
add    %edx,%ecx
add    -0x38(%rbp),%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
jmp    1506 <fornberg_weights+0x3e6>
movss  -0x20(%rbp),%xmm0
divss  -0x24(%rbp),%xmm0
movss  %xmm0,-0x60(%rbp)
mov    -0x40(%rbp),%eax
mov    %eax,-0x68(%rbp)
mov    -0x1c(%rbp),%edi
mov    -0x40(%rbp),%esi
sub    $0x1,%esi
mov    -0x38(%rbp),%edx
sub    $0x1,%edx
mov    -0x38(%rbp),%ecx
sub    $0x1,%ecx
mov    -0x18(%rbp),%r8
call   1030 <get_weight@plt>
mov    %eax,%ecx
mov    -0x68(%rbp),%eax
imul   %ecx,%eax
mov    %eax,-0x5c(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x38(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x64(%rbp)
mov    -0x1c(%rbp),%edi
mov    -0x40(%rbp),%esi
mov    -0x38(%rbp),%edx
sub    $0x1,%edx
mov    -0x38(%rbp),%ecx
sub    $0x1,%ecx
mov    -0x18(%rbp),%r8
call   1030 <get_weight@plt>
mov    -0x64(%rbp),%ecx
movss  -0x60(%rbp),%xmm0
mov    %eax,%edx
mov    -0x5c(%rbp),%eax
imul   %edx,%ecx
sub    %ecx,%eax
cvtsi2ss %eax,%xmm1
mulss  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
mov    -0x1c(%rbp),%ecx
add    $0x1,%ecx
mov    -0x1c(%rbp),%edx
add    $0x1,%edx
imul   %edx,%ecx
imul   -0x40(%rbp),%ecx
mov    -0x1c(%rbp),%edx
add    $0x1,%edx
imul   -0x38(%rbp),%edx
add    %edx,%ecx
add    -0x38(%rbp),%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
jmp    150b <fornberg_weights+0x3eb>
mov    -0x40(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x40(%rbp)
jmp    13b3 <fornberg_weights+0x293>
movss  -0x24(%rbp),%xmm0
movss  %xmm0,-0x20(%rbp)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    11ab <fornberg_weights+0x8b>
add    $0x70,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<fornberg_weights>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0xa8,%rsp
mov    %rcx,%r15
mov    %rdx,0x30(%rsp)
mov    %esi,%r12d
mov    %edi,%ebx
movl   $0x3f800000,(%rcx)
mov    %esi,%edi
xor    %esi,%esi
xor    %edx,%edx
xor    %ecx,%ecx
mov    %r15,%r8
call   1030 <get_weight@plt>
lea    0xea9(%rip),%rdi
xor    %esi,%esi
xor    %edx,%edx
xor    %ecx,%ecx
mov    %eax,%r8d
call   1040 <printf@plt>
test   %r12d,%r12d
jle    1535 <fornberg_weights+0x415>
lea    0x1(%r12),%eax
mov    %eax,%ecx
imul   %ecx,%ecx
mov    %eax,0x28(%rsp)
cltq
movslq %ebx,%rdx
mov    %rdx,0x70(%rsp)
mov    %rax,0x90(%rsp)
mov    %eax,%eax
mov    %rax,0x88(%rsp)
lea    0x1(%r12,%rcx,1),%eax
mov    %rax,0x48(%rsp)
lea    0x2(%r12,%rcx,1),%eax
mov    %eax,0x24(%rsp)
lea    0x2(%r12),%eax
mov    %eax,0x5c(%rsp)
movss  0xe40(%rip),%xmm0
mov    $0x1,%r14d
mov    %r15,0x40(%rsp)
mov    %rcx,0x50(%rsp)
mov    %r12,0x80(%rsp)
mov    0x80(%rsp),%r13
jmp    121d <fornberg_weights+0xfd>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
add    $0x1,%r14
mov    0x48(%rsp),%rax
add    0x28(%rsp),%eax
mov    %rax,0x48(%rsp)
mov    0x5c(%rsp),%eax
add    %eax,0x24(%rsp)
cmp    0x88(%rsp),%r14
mov    0x50(%rsp),%rcx
je     1535 <fornberg_weights+0x415>
mov    0x70(%rsp),%rax
cmp    %r14d,%eax
mov    %r14d,%edx
cmovl  %eax,%edx
movss  %xmm0,0x8(%rsp)
mov    %rdx,0x68(%rsp)
lea    0x1(%rdx),%esi
add    $0xffffffffffffffff,%rsi
lea    -0x1(%r14),%rdx
cmp    %rax,%r14
mov    %eax,%ebx
cmovl  %r14d,%ebx
mov    %r14,%rax
imul   %rcx,%rax
mov    %rdx,%rdi
mov    0x90(%rsp),%rcx
imul   %rcx,%rdi
add    %rax,%rdi
mov    %rdi,0x98(%rsp)
mov    %r14,%rax
imul   %rcx,%rax
mov    %rax,0xa0(%rsp)
xor    %edi,%edi
movss  0xd7c(%rip),%xmm0
mov    %r14,0x18(%rsp)
mov    %rsi,0x78(%rsp)
mov    %rdx,0x10(%rsp)
mov    %ebx,0x60(%rsp)
jmp    12ad <fornberg_weights+0x18d>
nopl   0x0(%rax)
add    $0x1,%rdi
cmp    %r14,%rdi
je     11f0 <fornberg_weights+0xd0>
mov    0x30(%rsp),%rcx
mov    (%rcx,%r14,4),%ebp
mov    %ebp,%eax
sub    (%rcx,%rdi,4),%eax
cmp    0x70(%rsp),%r14
jg     12d5 <fornberg_weights+0x1b5>
mov    0x98(%rsp),%rcx
add    %rdi,%rcx
movl   $0x0,(%r15,%rcx,4)
xorps  %xmm1,%xmm1
cvtsi2ss %eax,%xmm1
mulss  %xmm1,%xmm0
test   %ebx,%ebx
js     12a0 <fornberg_weights+0x180>
movss  %xmm0,0xc(%rsp)
mov    0xa0(%rsp),%rax
lea    (%rdi,%rax,1),%r14
mov    %rdi,%rcx
mov    %r12d,%edi
xor    %esi,%esi
mov    %rcx,0x38(%rsp)
mov    %r15,%r8
movss  %xmm1,0x2c(%rsp)
call   1030 <get_weight@plt>
imul   %ebp,%eax
xorps  %xmm0,%xmm0
cvtsi2ss %eax,%xmm0
divss  0x2c(%rsp),%xmm0
movss  %xmm0,(%r15,%r14,4)
cmpl   $0x0,0x68(%rsp)
je     1400 <fornberg_weights+0x2e0>
mov    0x48(%rsp),%rax
mov    0x38(%rsp),%rcx
lea    (%rax,%rcx,1),%r12d
xor    %ebp,%ebp
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
mov    0x30(%rsp),%rax
mov    0x18(%rsp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,0x64(%rsp)
lea    0x1(%rbp),%r14
mov    %r13d,%edi
mov    %r14d,%esi
mov    0x10(%rsp),%rdx
mov    0x38(%rsp),%rbx
mov    %ebx,%ecx
mov    %r15,%r8
call   1030 <get_weight@plt>
mov    %eax,%r15d
imul   0x64(%rsp),%r15d
mov    %r13d,%edi
mov    %ebp,%esi
mov    0x10(%rsp),%rdx
mov    %ebx,%ecx
mov    0x40(%rsp),%r8
call   1030 <get_weight@plt>
imul   %r14d,%eax
sub    %eax,%r15d
xorps  %xmm0,%xmm0
cvtsi2ss %r15d,%xmm0
mov    0x40(%rsp),%r15
divss  0x2c(%rsp),%xmm0
movslq %r12d,%r12
movss  %xmm0,(%r15,%r12,4)
add    0x50(%rsp),%r12d
mov    %r14,%rbp
cmp    %r14,0x78(%rsp)
jne    1350 <fornberg_weights+0x230>
mov    0x38(%rsp),%rdi
add    $0x1,%rdi
mov    0x18(%rsp),%r14
cmp    %r14,%rdi
mov    0x80(%rsp),%r12
movss  0xc(%rsp),%xmm0
mov    0x10(%rsp),%rdx
mov    0x60(%rsp),%ebx
jne    12ad <fornberg_weights+0x18d>
jmp    1430 <fornberg_weights+0x310>
mov    0x38(%rsp),%rdi
add    $0x1,%rdi
mov    0x18(%rsp),%r14
cmp    %r14,%rdi
movss  0xc(%rsp),%xmm0
mov    0x10(%rsp),%rdx
jne    12ad <fornberg_weights+0x18d>
jmp    1438 <fornberg_weights+0x318>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
test   %ebx,%ebx
js     11f0 <fornberg_weights+0xd0>
movss  0x8(%rsp),%xmm1
divss  %xmm0,%xmm1
movss  %xmm1,0x8(%rsp)
mov    0x28(%rsp),%eax
imul   %r14d,%eax
movslq %eax,%rbx
add    %r14,%rbx
mov    0x30(%rsp),%rax
mov    -0x4(%rax,%r14,4),%ebp
mov    %r12d,%edi
xor    %esi,%esi
mov    %edx,%ecx
mov    %r15,%r8
call   1030 <get_weight@plt>
mov    0x10(%rsp),%rcx
imul   %ebp,%eax
neg    %eax
xorps  %xmm0,%xmm0
cvtsi2ss %eax,%xmm0
mulss  0x8(%rsp),%xmm0
movss  %xmm0,(%r15,%rbx,4)
movss  0xc(%rsp),%xmm0
cmpl   $0x0,0x68(%rsp)
je     11f0 <fornberg_weights+0xd0>
xor    %ebx,%ebx
mov    0x24(%rsp),%eax
mov    %eax,%ebp
cs nopw 0x0(%rax,%rax,1)
nop
mov    %r12d,%edi
mov    %ebx,%esi
mov    %ecx,%edx
mov    %r15,%r8
call   1030 <get_weight@plt>
mov    %eax,%r14d
add    $0x1,%rbx
imul   %ebx,%r14d
mov    0x18(%rsp),%rax
mov    0x30(%rsp),%rcx
mov    -0x4(%rcx,%rax,4),%r15d
mov    %r12d,%edi
mov    %ebx,%esi
mov    0x10(%rsp),%rdx
mov    0x10(%rsp),%rcx
mov    0x40(%rsp),%r8
call   1030 <get_weight@plt>
mov    0x10(%rsp),%rcx
imul   %r15d,%eax
mov    0x40(%rsp),%r15
sub    %eax,%r14d
xorps  %xmm0,%xmm0
cvtsi2ss %r14d,%xmm0
mov    0x18(%rsp),%r14
mulss  0x8(%rsp),%xmm0
movslq %ebp,%rbp
movss  %xmm0,(%r15,%rbp,4)
movss  0xc(%rsp),%xmm0
add    0x50(%rsp),%ebp
cmp    %rbx,0x78(%rsp)
jne    14b0 <fornberg_weights+0x390>
jmp    11f0 <fornberg_weights+0xd0>
add    $0xa8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
@.str = [19 x i8] c"delta %d %d %d=%f\0A\00"
define void @fornberg_weights(i32 %0, i32 %1, i32* %2, float* %3) {
%5 = alloca i32
%6 = alloca i32
%7 = alloca i32*
%8 = alloca float*
%9 = alloca i32
%10 = alloca float
%11 = alloca float
%12 = alloca float
%13 = alloca i32
%14 = alloca i32
%15 = alloca i32
%16 = alloca i32
%17 = alloca i32
%18 = alloca i32
store i32 %0, i32* %5
store i32 %1, i32* %6
store i32* %2, i32** %7
store float* %3, float** %8
%19 = load i32, i32* %6
store i32 %19, i32* %9
%20 = load float*, float** %8
%21 = load i32, i32* %9
%22 = add nsw i32 %21, 1
%23 = load i32, i32* %9
%24 = add nsw i32 %23, 1
%25 = mul nsw i32 %22, %24
%26 = mul nsw i32 %25, 0
%27 = load i32, i32* %9
%28 = add nsw i32 %27, 1
%29 = mul nsw i32 %28, 0
%30 = add nsw i32 %26, %29
%31 = add nsw i32 %30, 0
%32 = sext i32 %31 to i64
%33 = getelementptr inbounds float, float* %20, i64 %32
store float 1.000000e+00, float* %33
store float 1.000000e+00, float* %10
%34 = load i32, i32* %9
%35 = load float*, float** %8
%36 = call i32 @get_weight(i32 %34, i32 0, i32 0, i32 0, float* %35)
%37 = call i32 @printf(i8* getelementptr inbounds ([19 x i8], [19 x i8]* @.str, i64 0, i64 0), i32 0, i32 0, i32 0, i32 %36)
store i32 1, i32* %16
br label %38
38: ; preds = %292, %4
%39 = load i32, i32* %16
%40 = load i32, i32* %9
%41 = icmp sle i32 %39, %40
br i1 %41, label %42, label %295
42: ; preds = %38
store float 1.000000e+00, float* %11
store i32 0, i32* %17
br label %43
43: ; preds = %186, %42
%44 = load i32, i32* %17
%45 = load i32, i32* %16
%46 = sub nsw i32 %45, 1
%47 = icmp sle i32 %44, %46
br i1 %47, label %48, label %189
48: ; preds = %43
%49 = load i32*, i32** %7
%50 = load i32, i32* %16
%51 = sext i32 %50 to i64
%52 = getelementptr inbounds i32, i32* %49, i64 %51
%53 = load i32, i32* %52
%54 = load i32*, i32** %7
%55 = load i32, i32* %17
%56 = sext i32 %55 to i64
%57 = getelementptr inbounds i32, i32* %54, i64 %56
%58 = load i32, i32* %57
%59 = sub nsw i32 %53, %58
%60 = sitofp i32 %59 to float
store float %60, float* %12
%61 = load float, float* %11
%62 = load float, float* %12
%63 = fmul float %61, %62
store float %63, float* %11
%64 = load i32, i32* %16
%65 = load i32, i32* %5
%66 = icmp sle i32 %64, %65
br i1 %66, label %67, label %86
67: ; preds = %48
%68 = load float*, float** %8
%69 = load i32, i32* %9
%70 = add nsw i32 %69, 1
%71 = load i32, i32* %9
%72 = add nsw i32 %71, 1
%73 = mul nsw i32 %70, %72
%74 = load i32, i32* %16
%75 = mul nsw i32 %73, %74
%76 = load i32, i32* %9
%77 = add nsw i32 %76, 1
%78 = load i32, i32* %16
%79 = sub nsw i32 %78, 1
%80 = mul nsw i32 %77, %79
%81 = add nsw i32 %75, %80
%82 = load i32, i32* %17
%83 = add nsw i32 %81, %82
%84 = sext i32 %83 to i64
%85 = getelementptr inbounds float, float* %68, i64 %84
store float 0.000000e+00, float* %85
br label %86
86: ; preds = %67, %48
%87 = load i32, i32* %16
%88 = load i32, i32* %5
%89 = icmp slt i32 %87, %88
br i1 %89, label %90, label %92
90: ; preds = %86
%91 = load i32, i32* %16
store i32 %91, i32* %13
br label %94
92: ; preds = %86
%93 = load i32, i32* %5
store i32 %93, i32* %13
br label %94
94: ; preds = %92, %90
store i32 0, i32* %18
br label %95
95: ; preds = %182, %94
%96 = load i32, i32* %18
%97 = load i32, i32* %13
%98 = icmp sle i32 %96, %97
br i1 %98, label %99, label %185
99: ; preds = %95
%100 = load i32, i32* %18
%101 = icmp eq i32 %100, 0
br i1 %101, label %102, label %136
102: ; preds = %99
%103 = load i32*, i32** %7
%104 = load i32, i32* %16
%105 = sext i32 %104 to i64
%106 = getelementptr inbounds i32, i32* %103, i64 %105
%107 = load i32, i32* %106
%108 = load i32, i32* %9
%109 = load i32, i32* %18
%110 = load i32, i32* %16
%111 = sub nsw i32 %110, 1
%112 = load i32, i32* %17
%113 = load float*, float** %8
%114 = call i32 @get_weight(i32 %108, i32 %109, i32 %111, i32 %112, float* %113)
%115 = mul nsw i32 %107, %114
%116 = sitofp i32 %115 to float
%117 = load float, float* %12
%118 = fdiv float %116, %117
%119 = load float*, float** %8
%120 = load i32, i32* %9
%121 = add nsw i32 %120, 1
%122 = load i32, i32* %9
%123 = add nsw i32 %122, 1
%124 = mul nsw i32 %121, %123
%125 = load i32, i32* %18
%126 = mul nsw i32 %124, %125
%127 = load i32, i32* %9
%128 = add nsw i32 %127, 1
%129 = load i32, i32* %16
%130 = mul nsw i32 %128, %129
%131 = add nsw i32 %126, %130
%132 = load i32, i32* %17
%133 = add nsw i32 %131, %132
%134 = sext i32 %133 to i64
%135 = getelementptr inbounds float, float* %119, i64 %134
store float %118, float* %135
br label %181
136: ; preds = %99
%137 = load i32*, i32** %7
%138 = load i32, i32* %16
%139 = sext i32 %138 to i64
%140 = getelementptr inbounds i32, i32* %137, i64 %139
%141 = load i32, i32* %140
%142 = load i32, i32* %9
%143 = load i32, i32* %18
%144 = load i32, i32* %16
%145 = sub nsw i32 %144, 1
%146 = load i32, i32* %17
%147 = load float*, float** %8
%148 = call i32 @get_weight(i32 %142, i32 %143, i32 %145, i32 %146, float* %147)
%149 = mul nsw i32 %141, %148
%150 = load i32, i32* %18
%151 = load i32, i32* %9
%152 = load i32, i32* %18
%153 = sub nsw i32 %152, 1
%154 = load i32, i32* %16
%155 = sub nsw i32 %154, 1
%156 = load i32, i32* %17
%157 = load float*, float** %8
%158 = call i32 @get_weight(i32 %151, i32 %153, i32 %155, i32 %156, float* %157)
%159 = mul nsw i32 %150, %158
%160 = sub nsw i32 %149, %159
%161 = sitofp i32 %160 to float
%162 = load float, float* %12
%163 = fdiv float %161, %162
%164 = load float*, float** %8
%165 = load i32, i32* %9
%166 = add nsw i32 %165, 1
%167 = load i32, i32* %9
%168 = add nsw i32 %167, 1
%169 = mul nsw i32 %166, %168
%170 = load i32, i32* %18
%171 = mul nsw i32 %169, %170
%172 = load i32, i32* %9
%173 = add nsw i32 %172, 1
%174 = load i32, i32* %16
%175 = mul nsw i32 %173, %174
%176 = add nsw i32 %171, %175
%177 = load i32, i32* %17
%178 = add nsw i32 %176, %177
%179 = sext i32 %178 to i64
%180 = getelementptr inbounds float, float* %164, i64 %179
store float %163, float* %180
br label %181
181: ; preds = %136, %102
br label %182
182: ; preds = %181
%183 = load i32, i32* %18
%184 = add nsw i32 %183, 1
store i32 %184, i32* %18
br label %95
185: ; preds = %95
br label %186
186: ; preds = %185
%187 = load i32, i32* %17
%188 = add nsw i32 %187, 1
store i32 %188, i32* %17
br label %43
189: ; preds = %43
store i32 0, i32* %18
br label %190
190: ; preds = %287, %189
%191 = load i32, i32* %18
%192 = load i32, i32* %13
%193 = icmp sle i32 %191, %192
br i1 %193, label %194, label %290
194: ; preds = %190
%195 = load i32, i32* %18
%196 = icmp eq i32 %195, 0
br i1 %196, label %197, label %236
197: ; preds = %194
%198 = load float, float* %10
%199 = load float, float* %11
%200 = fdiv float %198, %199
%201 = load i32*, i32** %7
%202 = load i32, i32* %16
%203 = sub nsw i32 %202, 1
%204 = sext i32 %203 to i64
%205 = getelementptr inbounds i32, i32* %201, i64 %204
%206 = load i32, i32* %205
%207 = mul nsw i32 -1, %206
%208 = load i32, i32* %9
%209 = load i32, i32* %18
%210 = load i32, i32* %16
%211 = sub nsw i32 %210, 1
%212 = load i32, i32* %16
%213 = sub nsw i32 %212, 1
%214 = load float*, float** %8
%215 = call i32 @get_weight(i32 %208, i32 %209, i32 %211, i32 %213, float* %214)
%216 = mul nsw i32 %207, %215
%217 = sitofp i32 %216 to float
%218 = fmul float %200, %217
%219 = load float*, float** %8
%220 = load i32, i32* %9
%221 = add nsw i32 %220, 1
%222 = load i32, i32* %9
%223 = add nsw i32 %222, 1
%224 = mul nsw i32 %221, %223
%225 = load i32, i32* %18
%226 = mul nsw i32 %224, %225
%227 = load i32, i32* %9
%228 = add nsw i32 %227, 1
%229 = load i32, i32* %16
%230 = mul nsw i32 %228, %229
%231 = add nsw i32 %226, %230
%232 = load i32, i32* %16
%233 = add nsw i32 %231, %232
%234 = sext i32 %233 to i64
%235 = getelementptr inbounds float, float* %219, i64 %234
store float %218, float* %235
br label %286
236: ; preds = %194
%237 = load float, float* %10
%238 = load float, float* %11
%239 = fdiv float %237, %238
%240 = load i32, i32* %18
%241 = load i32, i32* %9
%242 = load i32, i32* %18
%243 = sub nsw i32 %242, 1
%244 = load i32, i32* %16
%245 = sub nsw i32 %244, 1
%246 = load i32, i32* %16
%247 = sub nsw i32 %246, 1
%248 = load float*, float** %8
%249 = call i32 @get_weight(i32 %241, i32 %243, i32 %245, i32 %247, float* %248)
%250 = mul nsw i32 %240, %249
%251 = load i32*, i32** %7
%252 = load i32, i32* %16
%253 = sub nsw i32 %252, 1
%254 = sext i32 %253 to i64
%255 = getelementptr inbounds i32, i32* %251, i64 %254
%256 = load i32, i32* %255
%257 = load i32, i32* %9
%258 = load i32, i32* %18
%259 = load i32, i32* %16
%260 = sub nsw i32 %259, 1
%261 = load i32, i32* %16
%262 = sub nsw i32 %261, 1
%263 = load float*, float** %8
%264 = call i32 @get_weight(i32 %257, i32 %258, i32 %260, i32 %262, float* %263)
%265 = mul nsw i32 %256, %264
%266 = sub nsw i32 %250, %265
%267 = sitofp i32 %266 to float
%268 = fmul float %239, %267
%269 = load float*, float** %8
%270 = load i32, i32* %9
%271 = add nsw i32 %270, 1
%272 = load i32, i32* %9
%273 = add nsw i32 %272, 1
%274 = mul nsw i32 %271, %273
%275 = load i32, i32* %18
%276 = mul nsw i32 %274, %275
%277 = load i32, i32* %9
%278 = add nsw i32 %277, 1
%279 = load i32, i32* %16
%280 = mul nsw i32 %278, %279
%281 = add nsw i32 %276, %280
%282 = load i32, i32* %16
%283 = add nsw i32 %281, %282
%284 = sext i32 %283 to i64
%285 = getelementptr inbounds float, float* %269, i64 %284
store float %268, float* %285
br label %286
286: ; preds = %236, %197
br label %287
287: ; preds = %286
%288 = load i32, i32* %18
%289 = add nsw i32 %288, 1
store i32 %289, i32* %18
br label %190
290: ; preds = %190
%291 = load float, float* %11
store float %291, float* %10
br label %292
292: ; preds = %290
%293 = load i32, i32* %16
%294 = add nsw i32 %293, 1
store i32 %294, i32* %16
br label %38
295: ; preds = %38
ret void
}

LLVM IR (-O2):
@.str = [19 x i8] c"delta %d %d %d=%f\0A\00"
define void @fornberg_weights(i32 %0, i32 %1, i32* nocapture %2, float* %3) {
%5 = add i32 %1, 1
store float 1.000000e+00, float* %3
%6 = tail call i32 @get_weight(i32 %1, i32 0, i32 0, i32 0, float* %3)
%7 = tail call i32 @printf(i8* getelementptr inbounds ([19 x i8], [19 x i8]* @.str, i64 0, i64 0), i32 0, i32 0, i32 0, i32 %6)
%8 = icmp slt i32 %1, 1
br i1 %8, label %132, label %9
9: ; preds = %4
%10 = mul nsw i32 %5, %5
%11 = zext i32 %10 to i64
%12 = sext i32 %5 to i64
%13 = sext i32 %0 to i64
%14 = zext i32 %5 to i64
br label %15
15: ; preds = %9, %129
%16 = phi i64 [ 1, %9 ], [ %130, %129 ]
%17 = phi float [ 1.000000e+00, %9 ], [ %64, %129 ]
%18 = trunc i64 %16 to i32
%19 = call i32 @llvm.smin.i32(i32 %0, i32 %18)
%20 = add i32 %19, 1
%21 = add nsw i64 %16, -1
%22 = getelementptr inbounds i32, i32* %2, i64 %16
%23 = icmp sgt i64 %16, %13
%24 = icmp slt i64 %16, %13
%25 = select i1 %24, i32 %18, i32 %0
%26 = icmp slt i32 %25, 0
%27 = mul nsw i64 %16, %11
%28 = mul nsw i64 %21, %12
%29 = add nsw i64 %28, %27
%30 = mul nsw i64 %16, %12
%31 = trunc i64 %30 to i32
%32 = zext i32 %20 to i64
%33 = trunc i64 %21 to i32
%34 = icmp eq i32 %19, 0
%35 = trunc i64 %21 to i32
br label %56
36: ; preds = %98
br i1 %26, label %129, label %37
37: ; preds = %103, %36
%38 = fdiv float %17, %64
%39 = getelementptr inbounds i32, i32* %2, i64 %21
%40 = mul nsw i32 %5, %18
%41 = add i32 %40, %18
%42 = zext i32 %20 to i64
%43 = sext i32 %40 to i64
%44 = add nsw i64 %16, %43
%45 = getelementptr inbounds float, float* %3, i64 %44
%46 = load i32, i32* %39
%47 = trunc i64 %21 to i32
%48 = tail call i32 @get_weight(i32 %1, i32 0, i32 %47, i32 %47, float* %3)
%49 = mul i32 %46, %48
%50 = sub i32 0, %49
%51 = sitofp i32 %50 to float
%52 = fmul float %38, %51
store float %52, float* %45
%53 = icmp eq i32 %19, 0
br i1 %53, label %129, label %54
54: ; preds = %37
%55 = trunc i64 %21 to i32
br label %109
56: ; preds = %101, %15
%57 = phi i64 [ 0, %15 ], [ %102, %101 ]
%58 = phi float [ 1.000000e+00, %15 ], [ %64, %101 ]
%59 = load i32, i32* %22
%60 = getelementptr inbounds i32, i32* %2, i64 %57
%61 = load i32, i32* %60
%62 = sub nsw i32 %59, %61
%63 = sitofp i32 %62 to float
%64 = fmul float %58, %63
br i1 %23, label %68, label %65
65: ; preds = %56
%66 = add nsw i64 %29, %57
%67 = getelementptr inbounds float, float* %3, i64 %66
store float 0.000000e+00, float* %67
br label %68
68: ; preds = %65, %56
br i1 %26, label %106, label %69
69: ; preds = %68
%70 = trunc i64 %57 to i32
%71 = add i32 %70, %31
%72 = add nsw i64 %57, %30
%73 = getelementptr inbounds float, float* %3, i64 %72
%74 = tail call i32 @get_weight(i32 %1, i32 0, i32 %33, i32 %70, float* %3)
%75 = mul nsw i32 %74, %59
%76 = sitofp i32 %75 to float
%77 = fdiv float %76, %63
store float %77, float* %73
br i1 %34, label %103, label %78
78: ; preds = %69, %78
%79 = phi i64 [ %96, %78 ], [ 1, %69 ]
%80 = load i32, i32* %22
%81 = trunc i64 %79 to i32
%82 = tail call i32 @get_weight(i32 %1, i32 %81, i32 %35, i32 %70, float* %3)
%83 = mul nsw i32 %82, %80
%84 = trunc i64 %79 to i32
%85 = add i32 %84, -1
%86 = tail call i32 @get_weight(i32 %1, i32 %85, i32 %35, i32 %70, float* %3)
%87 = mul nsw i32 %86, %81
%88 = sub nsw i32 %83, %87
%89 = sitofp i32 %88 to float
%90 = fdiv float %89, %63
%91 = trunc i64 %79 to i32
%92 = mul i32 %10, %91
%93 = add i32 %71, %92
%94 = sext i32 %93 to i64
%95 = getelementptr inbounds float, float* %3, i64 %94
store float %90, float* %95
%96 = add nuw nsw i64 %79, 1
%97 = icmp eq i64 %96, %32
br i1 %97, label %98, label %78
98: ; preds = %78
%99 = add nuw nsw i64 %57, 1
%100 = icmp eq i64 %99, %16
br i1 %100, label %36, label %101
101: ; preds = %98, %106, %103
%102 = phi i64 [ %99, %98 ], [ %107, %106 ], [ %104, %103 ]
br label %56
103: ; preds = %69
%104 = add nuw nsw i64 %57, 1
%105 = icmp eq i64 %104, %16
br i1 %105, label %37, label %101
106: ; preds = %68
%107 = add nuw nsw i64 %57, 1
%108 = icmp eq i64 %107, %16
br i1 %108, label %129, label %101
109: ; preds = %54, %109
%110 = phi i64 [ 1, %54 ], [ %127, %109 ]
%111 = trunc i64 %110 to i32
%112 = add i32 %111, -1
%113 = tail call i32 @get_weight(i32 %1, i32 %112, i32 %55, i32 %55, float* %3)
%114 = trunc i64 %110 to i32
%115 = mul nsw i32 %113, %114
%116 = load i32, i32* %39
%117 = tail call i32 @get_weight(i32 %1, i32 %114, i32 %55, i32 %55, float* %3)
%118 = mul nsw i32 %117, %116
%119 = sub nsw i32 %115, %118
%120 = sitofp i32 %119 to float
%121 = fmul float %38, %120
%122 = trunc i64 %110 to i32
%123 = mul i32 %10, %122
%124 = add i32 %41, %123
%125 = sext i32 %124 to i64
%126 = getelementptr inbounds float, float* %3, i64 %125
store float %121, float* %126
%127 = add nuw nsw i64 %110, 1
%128 = icmp eq i64 %127, %42
br i1 %128, label %129, label %109
129: ; preds = %106, %109, %37, %36
%130 = add nuw nsw i64 %16, 1
%131 = icmp eq i64 %130, %14
br i1 %131, label %132, label %15
132: ; preds = %129, %4
ret void
}