#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*/ tf_dtx(double *, int, int, double *, double *);
void tf_dtxtil(double *x, int n, int k, double *a, double *b) {
  int i;
  if (k > 0)
    for (i = 0; i < n - k; i++) {
      a[i] = a[i] * k / (x[k + i] - x[i]);
    }
  tf_dtx(x, n, k, a, b);
}


x86 assembly (clang, -O0):
<tf_dtxtil>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
mov    %rcx,-0x18(%rbp)
mov    %r8,-0x20(%rbp)
cmpl   $0x0,-0x10(%rbp)
jle    11a5 <tf_dtxtil+0x95>
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
mov    -0xc(%rbp),%ecx
sub    -0x10(%rbp),%ecx
cmp    %ecx,%eax
jge    11a0 <tf_dtxtil+0x90>
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
cvtsi2sdl -0x10(%rbp),%xmm1
mulsd  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
mov    -0x10(%rbp),%ecx
add    -0x24(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm1
divsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    113b <tf_dtxtil+0x2b>
jmp    11a5 <tf_dtxtil+0x95>
mov    -0x8(%rbp),%rdi
mov    -0xc(%rbp),%esi
mov    -0x10(%rbp),%edx
mov    -0x18(%rbp),%rcx
mov    -0x20(%rbp),%r8
call   1030 <tf_dtx@plt>
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<tf_dtxtil>:
push   %rbp
push   %r15
push   %r14
push   %rbx
test   %edx,%edx
jle    12dc <tf_dtxtil+0x1cc>
mov    %esi,%eax
sub    %edx,%eax
test   %eax,%eax
jle    12dc <tf_dtxtil+0x1cc>
cvtsi2sd %edx,%xmm0
mov    %edx,%r10d
mov    %eax,%r9d
cmp    $0x1,%eax
jne    11ca <tf_dtxtil+0xba>
xor    %eax,%eax
mov    %rax,%rbx
not    %rbx
test   $0x1,%r9b
je     116f <tf_dtxtil+0x5f>
movsd  (%rcx,%rax,8),%xmm1
lea    (%rax,%r10,1),%rbp
movsd  (%rdi,%rbp,8),%xmm2
subsd  (%rdi,%rax,8),%xmm2
mulsd  %xmm0,%xmm1
divsd  %xmm2,%xmm1
movsd  %xmm1,(%rcx,%rax,8)
or     $0x1,%rax
add    %r9,%rbx
je     12dc <tf_dtxtil+0x1cc>
lea    (%rdi,%r10,8),%rbx
add    $0x8,%rbx
movsd  (%rcx,%rax,8),%xmm1
movsd  -0x8(%rbx,%rax,8),%xmm2
subsd  (%rdi,%rax,8),%xmm2
mulsd  %xmm0,%xmm1
divsd  %xmm2,%xmm1
movsd  %xmm1,(%rcx,%rax,8)
movsd  0x8(%rcx,%rax,8),%xmm1
mulsd  %xmm0,%xmm1
movsd  (%rbx,%rax,8),%xmm2
subsd  0x8(%rdi,%rax,8),%xmm2
divsd  %xmm2,%xmm1
movsd  %xmm1,0x8(%rcx,%rax,8)
add    $0x2,%rax
cmp    %rax,%r9
jne    1180 <tf_dtxtil+0x70>
jmp    12dc <tf_dtxtil+0x1cc>
lea    (%rcx,%r9,8),%rax
lea    (%rdi,%r9,8),%rbx
lea    (%rdi,%r10,8),%r11
cmp    %rcx,%rbx
seta   %bpl
lea    (%r10,%r9,1),%rbx
cmp    %rdi,%rax
seta   %r14b
lea    (%rdi,%rbx,8),%rbx
cmp    %rcx,%rbx
seta   %bl
cmp    %rax,%r11
setb   %r11b
xor    %eax,%eax
test   %r14b,%bpl
jne    113f <tf_dtxtil+0x2f>
and    %r11b,%bl
jne    113f <tf_dtxtil+0x2f>
mov    %r9d,%eax
and    $0xfffffffe,%eax
movapd %xmm0,%xmm1
unpcklpd %xmm0,%xmm1
lea    -0x2(%rax),%rbp
mov    %rbp,%r14
shr    %r14
add    $0x1,%r14
test   %rbp,%rbp
je     12e7 <tf_dtxtil+0x1d7>
mov    %r14,%rbx
and    $0xfffffffffffffffe,%rbx
lea    (%rdi,%r10,8),%r15
add    $0x10,%r15
xor    %r11d,%r11d
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movupd (%rcx,%r11,8),%xmm2
movupd 0x10(%rcx,%r11,8),%xmm3
mulpd  %xmm1,%xmm2
movupd -0x10(%r15,%r11,8),%xmm4
movupd (%rdi,%r11,8),%xmm5
subpd  %xmm5,%xmm4
divpd  %xmm4,%xmm2
movupd %xmm2,(%rcx,%r11,8)
mulpd  %xmm1,%xmm3
movupd (%r15,%r11,8),%xmm2
movupd 0x10(%rdi,%r11,8),%xmm4
subpd  %xmm4,%xmm2
divpd  %xmm2,%xmm3
movupd %xmm3,0x10(%rcx,%r11,8)
add    $0x4,%r11
add    $0xfffffffffffffffe,%rbx
jne    1250 <tf_dtxtil+0x140>
test   $0x1,%r14b
je     12d3 <tf_dtxtil+0x1c3>
movupd (%rcx,%r11,8),%xmm2
mulpd  %xmm1,%xmm2
lea    (%r11,%r10,1),%rbx
movupd (%rdi,%rbx,8),%xmm1
movupd (%rdi,%r11,8),%xmm3
subpd  %xmm3,%xmm1
divpd  %xmm1,%xmm2
movupd %xmm2,(%rcx,%r11,8)
cmp    %r9,%rax
jne    113f <tf_dtxtil+0x2f>
pop    %rbx
pop    %r14
pop    %r15
pop    %rbp
jmp    1030 <tf_dtx@plt>
xor    %r11d,%r11d
test   $0x1,%r14b
jne    12ac <tf_dtxtil+0x19c>
jmp    12d3 <tf_dtxtil+0x1c3>

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

LLVM IR (-O2):
define void @tf_dtxtil(double* %0, i32 %1, i32 %2, double* %3, double* %4) {
%6 = icmp sgt i32 %2, 0
br i1 %6, label %7, label %140
7: ; preds = %5
%8 = sub i32 %1, %2
%9 = icmp sgt i32 %8, 0
br i1 %9, label %10, label %140
10: ; preds = %7
%11 = sitofp i32 %2 to double
%12 = zext i32 %2 to i64
%13 = zext i32 %8 to i64
%14 = icmp eq i32 %8, 1
br i1 %14, label %94, label %15
15: ; preds = %10
%16 = getelementptr double, double* %3, i64 %13
%17 = getelementptr double, double* %0, i64 %13
%18 = getelementptr double, double* %0, i64 %12
%19 = add nuw nsw i64 %12, %13
%20 = getelementptr double, double* %0, i64 %19
%21 = icmp ugt double* %17, %3
%22 = icmp ugt double* %16, %0
%23 = and i1 %21, %22
%24 = icmp ugt double* %20, %3
%25 = icmp ult double* %18, %16
%26 = and i1 %24, %25
%27 = or i1 %23, %26
br i1 %27, label %94, label %28
28: ; preds = %15
%29 = and i64 %13, 4294967294
%30 = insertelement <2 x double> poison, double %11, i64 0
%31 = shufflevector <2 x double> %30, <2 x double> poison, <2 x i32> zeroinitializer
%32 = add nsw i64 %29, -2
%33 = lshr exact i64 %32, 1
%34 = add nuw i64 %33, 1
%35 = and i64 %34, 1
%36 = icmp eq i64 %32, 0
br i1 %36, label %74, label %37
37: ; preds = %28
%38 = and i64 %34, -2
br label %39
39: ; preds = %39, %37
%40 = phi i64 [ 0, %37 ], [ %71, %39 ]
%41 = phi i64 [ 0, %37 ], [ %72, %39 ]
%42 = getelementptr inbounds double, double* %3, i64 %40
%43 = bitcast double* %42 to <2 x double>*
%44 = load <2 x double>, <2 x double>* %43
%45 = fmul <2 x double> %44, %31
%46 = add nuw nsw i64 %40, %12
%47 = getelementptr inbounds double, double* %0, i64 %46
%48 = bitcast double* %47 to <2 x double>*
%49 = load <2 x double>, <2 x double>* %48
%50 = getelementptr inbounds double, double* %0, i64 %40
%51 = bitcast double* %50 to <2 x double>*
%52 = load <2 x double>, <2 x double>* %51
%53 = fsub <2 x double> %49, %52
%54 = fdiv <2 x double> %45, %53
%55 = bitcast double* %42 to <2 x double>*
store <2 x double> %54, <2 x double>* %55
%56 = or i64 %40, 2
%57 = getelementptr inbounds double, double* %3, i64 %56
%58 = bitcast double* %57 to <2 x double>*
%59 = load <2 x double>, <2 x double>* %58
%60 = fmul <2 x double> %59, %31
%61 = add nuw nsw i64 %56, %12
%62 = getelementptr inbounds double, double* %0, i64 %61
%63 = bitcast double* %62 to <2 x double>*
%64 = load <2 x double>, <2 x double>* %63
%65 = getelementptr inbounds double, double* %0, i64 %56
%66 = bitcast double* %65 to <2 x double>*
%67 = load <2 x double>, <2 x double>* %66
%68 = fsub <2 x double> %64, %67
%69 = fdiv <2 x double> %60, %68
%70 = bitcast double* %57 to <2 x double>*
store <2 x double> %69, <2 x double>* %70
%71 = add nuw i64 %40, 4
%72 = add nuw i64 %41, 2
%73 = icmp eq i64 %72, %38
br i1 %73, label %74, label %39
74: ; preds = %39, %28
%75 = phi i64 [ 0, %28 ], [ %71, %39 ]
%76 = icmp eq i64 %35, 0
br i1 %76, label %92, label %77
77: ; preds = %74
%78 = getelementptr inbounds double, double* %3, i64 %75
%79 = bitcast double* %78 to <2 x double>*
%80 = load <2 x double>, <2 x double>* %79
%81 = fmul <2 x double> %80, %31
%82 = add nuw nsw i64 %75, %12
%83 = getelementptr inbounds double, double* %0, i64 %82
%84 = bitcast double* %83 to <2 x double>*
%85 = load <2 x double>, <2 x double>* %84
%86 = getelementptr inbounds double, double* %0, i64 %75
%87 = bitcast double* %86 to <2 x double>*
%88 = load <2 x double>, <2 x double>* %87
%89 = fsub <2 x double> %85, %88
%90 = fdiv <2 x double> %81, %89
%91 = bitcast double* %78 to <2 x double>*
store <2 x double> %90, <2 x double>* %91
br label %92
92: ; preds = %74, %77
%93 = icmp eq i64 %29, %13
br i1 %93, label %140, label %94
94: ; preds = %15, %10, %92
%95 = phi i64 [ 0, %15 ], [ 0, %10 ], [ %29, %92 ]
%96 = xor i64 %95, -1
%97 = and i64 %13, 1
%98 = icmp eq i64 %97, 0
br i1 %98, label %111, label %99
99: ; preds = %94
%100 = getelementptr inbounds double, double* %3, i64 %95
%101 = load double, double* %100
%102 = fmul double %101, %11
%103 = add nuw nsw i64 %95, %12
%104 = getelementptr inbounds double, double* %0, i64 %103
%105 = load double, double* %104
%106 = getelementptr inbounds double, double* %0, i64 %95
%107 = load double, double* %106
%108 = fsub double %105, %107
%109 = fdiv double %102, %108
store double %109, double* %100
%110 = or i64 %95, 1
br label %111
111: ; preds = %99, %94
%112 = phi i64 [ %95, %94 ], [ %110, %99 ]
%113 = sub nsw i64 0, %13
%114 = icmp eq i64 %96, %113
br i1 %114, label %140, label %115
115: ; preds = %111, %115
%116 = phi i64 [ %138, %115 ], [ %112, %111 ]
%117 = getelementptr inbounds double, double* %3, i64 %116
%118 = load double, double* %117
%119 = fmul double %118, %11
%120 = add nuw nsw i64 %116, %12
%121 = getelementptr inbounds double, double* %0, i64 %120
%122 = load double, double* %121
%123 = getelementptr inbounds double, double* %0, i64 %116
%124 = load double, double* %123
%125 = fsub double %122, %124
%126 = fdiv double %119, %125
store double %126, double* %117
%127 = add nuw nsw i64 %116, 1
%128 = getelementptr inbounds double, double* %3, i64 %127
%129 = load double, double* %128
%130 = fmul double %129, %11
%131 = add nuw nsw i64 %127, %12
%132 = getelementptr inbounds double, double* %0, i64 %131
%133 = load double, double* %132
%134 = getelementptr inbounds double, double* %0, i64 %127
%135 = load double, double* %134
%136 = fsub double %133, %135
%137 = fdiv double %130, %136
store double %137, double* %128
%138 = add nuw nsw i64 %116, 2
%139 = icmp eq i64 %138, %13
br i1 %139, label %140, label %115
140: ; preds = %111, %115, %92, %7, %5
%141 = tail call i32 @tf_dtx(double* %0, i32 %1, i32 %2, double* %3, double* %4)
ret void
}