#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 */
double *calloc(int, int);
int /*<<< orphan*/ free(double *);
double pow(double, int);
int sqrt(int);
void fast_spline(double *s, int N, double *b, int k_0) {
  int k;
  double z_1;
  double *c_plus, *c_minus;
  z_1 = -2. + sqrt(3.);
  c_plus = calloc(N, sizeof(double));
  c_plus[0] = 0.;
  for (k = 0; k < k_0; k++) {
    c_plus[0] += s[k] * pow(z_1, k);
  }
  for (k = 1; k < N; k++) {
    c_plus[k] = s[k] + z_1 * c_plus[k - 1];
  }
  b[N - 1] = (-z_1 / (1 - z_1 * z_1)) * (c_plus[N - 1] + z_1 * c_plus[N - 2]);
  for (k = N - 2; k >= 0; k--) {
    b[k] = z_1 * (b[k + 1] - c_plus[k]);
  };
  free(c_plus);
}


x86 assembly (clang, -O0):
<fast_spline>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    $0x3,%edi
call   1060 <sqrt@plt>
cvtsi2sd %eax,%xmm1
movsd  0xe94(%rip),%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0xc(%rbp),%edi
mov    $0x8,%esi
call   1040 <calloc@plt>
mov    %rax,-0x30(%rbp)
mov    -0x30(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax)
movl   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    11ed <fast_spline+0xad>
mov    -0x8(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x28(%rbp),%xmm0
mov    -0x20(%rbp),%edi
call   1050 <pow@plt>
movaps %xmm0,%xmm2
movsd  -0x40(%rbp),%xmm0
mov    -0x30(%rbp),%rax
movsd  (%rax),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    1198 <fast_spline+0x58>
movl   $0x1,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1247 <fast_spline+0x107>
mov    -0x8(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x28(%rbp),%xmm0
mov    -0x30(%rbp),%rax
mov    -0x20(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x30(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    11f4 <fast_spline+0xb4>
movsd  -0x28(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  -0x28(%rbp),%xmm1
movsd  -0x28(%rbp),%xmm2
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm1
movsd  0xd78(%rip),%xmm2
addsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
mov    -0x30(%rbp),%rax
mov    -0xc(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x28(%rbp),%xmm1
mov    -0x30(%rbp),%rax
mov    -0xc(%rbp),%ecx
sub    $0x2,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm3
mulsd  %xmm3,%xmm1
addsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
mov    -0xc(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0xc(%rbp),%eax
sub    $0x2,%eax
mov    %eax,-0x20(%rbp)
cmpl   $0x0,-0x20(%rbp)
jl     1335 <fast_spline+0x1f5>
movsd  -0x28(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x20(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x30(%rbp),%rax
movslq -0x20(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
movslq -0x20(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x20(%rbp)
jmp    12e8 <fast_spline+0x1a8>
mov    -0x30(%rbp),%rdi
call   1030 <free@plt>
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<fast_spline>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x28,%rsp
mov    %ecx,%ebx
mov    %rdx,%r14
mov    %esi,%r15d
mov    %rdi,%r13
mov    $0x3,%edi
call   1060 <sqrt@plt>
cvtsi2sd %eax,%xmm0
addsd  0xe91(%rip),%xmm0
movapd %xmm0,0x10(%rsp)
mov    %r15d,%edi
mov    $0x8,%esi
call   1040 <calloc@plt>
mov    %rax,%r12
movq   $0x0,(%rax)
test   %ebx,%ebx
jle    11b9 <fast_spline+0x79>
mov    %r15d,0x24(%rsp)
mov    %ebx,%r15d
lea    -0x1(%r15),%rax
mov    %r15d,%ebp
and    $0x3,%ebp
cmp    $0x3,%rax
movapd 0x10(%rsp),%xmm3
jae    11c8 <fast_spline+0x88>
xorpd  %xmm0,%xmm0
xor    %ebx,%ebx
jmp    12b4 <fast_spline+0x174>
xorpd  %xmm0,%xmm0
movapd 0x10(%rsp),%xmm3
jmp    130c <fast_spline+0x1cc>
and    $0xfffffffc,%r15d
xorpd  %xmm0,%xmm0
xor    %ebx,%ebx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movsd  %xmm0,(%rsp)
movsd  0x0(%r13,%rbx,8),%xmm1
movsd  %xmm1,0x8(%rsp)
movapd %xmm3,%xmm0
mov    %ebx,%edi
call   1050 <pow@plt>
mulsd  0x8(%rsp),%xmm0
addsd  (%rsp),%xmm0
movsd  %xmm0,(%rsp)
movsd  %xmm0,(%r12)
movsd  0x8(%r13,%rbx,8),%xmm0
movsd  %xmm0,0x8(%rsp)
lea    0x1(%rbx),%edi
movapd 0x10(%rsp),%xmm0
call   1050 <pow@plt>
mulsd  0x8(%rsp),%xmm0
addsd  (%rsp),%xmm0
movsd  %xmm0,(%rsp)
movsd  %xmm0,(%r12)
movsd  0x10(%r13,%rbx,8),%xmm0
movsd  %xmm0,0x8(%rsp)
lea    0x2(%rbx),%edi
movapd 0x10(%rsp),%xmm0
call   1050 <pow@plt>
mulsd  0x8(%rsp),%xmm0
addsd  (%rsp),%xmm0
movsd  %xmm0,(%rsp)
movsd  %xmm0,(%r12)
movsd  0x18(%r13,%rbx,8),%xmm0
movsd  %xmm0,0x8(%rsp)
lea    0x3(%rbx),%edi
movapd 0x10(%rsp),%xmm0
call   1050 <pow@plt>
movapd 0x10(%rsp),%xmm3
mulsd  0x8(%rsp),%xmm0
addsd  (%rsp),%xmm0
movsd  %xmm0,(%r12)
add    $0x4,%rbx
cmp    %r15,%rbx
jne    11e0 <fast_spline+0xa0>
test   %rbp,%rbp
mov    0x24(%rsp),%r15d
je     130c <fast_spline+0x1cc>
xchg   %ax,%ax
movsd  %xmm0,(%rsp)
movsd  0x0(%r13,%rbx,8),%xmm0
movsd  %xmm0,0x8(%rsp)
movapd %xmm3,%xmm0
mov    %ebx,%edi
call   1050 <pow@plt>
movapd 0x10(%rsp),%xmm3
mulsd  0x8(%rsp),%xmm0
movsd  (%rsp),%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,(%rsp)
movsd  (%rsp),%xmm0
movsd  %xmm0,(%r12)
add    $0x1,%rbx
add    $0xffffffffffffffff,%rbp
jne    12c0 <fast_spline+0x180>
cmp    $0x2,%r15d
jl     13c7 <fast_spline+0x287>
mov    %r15d,%esi
lea    -0x1(%rsi),%rdx
add    $0xfffffffffffffffe,%rsi
mov    %edx,%eax
and    $0x3,%eax
mov    $0x1,%ecx
cmp    $0x3,%rsi
jb     1395 <fast_spline+0x255>
and    $0xfffffffffffffffc,%rdx
xor    %ecx,%ecx
nopw   0x0(%rax,%rax,1)
mulsd  %xmm3,%xmm0
addsd  0x8(%r13,%rcx,8),%xmm0
movsd  %xmm0,0x8(%r12,%rcx,8)
mulsd  %xmm3,%xmm0
addsd  0x10(%r13,%rcx,8),%xmm0
movsd  %xmm0,0x10(%r12,%rcx,8)
mulsd  %xmm3,%xmm0
addsd  0x18(%r13,%rcx,8),%xmm0
movsd  %xmm0,0x18(%r12,%rcx,8)
mulsd  %xmm3,%xmm0
addsd  0x20(%r13,%rcx,8),%xmm0
movsd  %xmm0,0x20(%r12,%rcx,8)
add    $0x4,%rcx
cmp    %rcx,%rdx
jne    1340 <fast_spline+0x200>
add    $0x1,%rcx
test   %rax,%rax
je     13c7 <fast_spline+0x287>
lea    (%r12,%rcx,8),%rdx
lea    0x0(,%rcx,8),%rcx
add    %r13,%rcx
xor    %esi,%esi
nopl   0x0(%rax,%rax,1)
mulsd  %xmm3,%xmm0
addsd  (%rcx,%rsi,8),%xmm0
movsd  %xmm0,(%rdx,%rsi,8)
add    $0x1,%rsi
cmp    %rsi,%rax
jne    13b0 <fast_spline+0x270>
movapd 0xc41(%rip),%xmm0
xorpd  %xmm3,%xmm0
movapd %xmm3,%xmm1
mulsd  %xmm3,%xmm1
movsd  0xc25(%rip),%xmm2
subsd  %xmm1,%xmm2
movslq %r15d,%rax
movsd  -0x10(%r12,%rax,8),%xmm1
mulsd  %xmm3,%xmm1
addsd  -0x8(%r12,%rax,8),%xmm1
divsd  %xmm2,%xmm0
mulsd  %xmm0,%xmm1
movsd  %xmm1,-0x8(%r14,%rax,8)
cmp    $0x2,%r15d
jl     1473 <fast_spline+0x333>
add    $0xfffffffe,%r15d
mov    %r15d,%eax
test   $0x1,%r15b
jne    1439 <fast_spline+0x2f9>
movsd  0x8(%r14,%rax,8),%xmm0
subsd  (%r12,%rax,8),%xmm0
mulsd  %xmm3,%xmm0
movsd  %xmm0,(%r14,%rax,8)
add    $0xffffffffffffffff,%rax
test   %r15d,%r15d
je     1473 <fast_spline+0x333>
xchg   %ax,%ax
movsd  0x8(%r14,%rax,8),%xmm0
subsd  (%r12,%rax,8),%xmm0
mulsd  %xmm3,%xmm0
movsd  %xmm0,(%r14,%rax,8)
subsd  -0x8(%r12,%rax,8),%xmm0
mulsd  %xmm3,%xmm0
movsd  %xmm0,-0x8(%r14,%rax,8)
add    $0xfffffffffffffffe,%rax
cmp    $0xffffffffffffffff,%rax
jne    1440 <fast_spline+0x300>
mov    %r12,%rdi
add    $0x28,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
jmp    1030 <free@plt>

LLVM IR (-O0):
define void @fast_spline(double* %0, i32 %1, double* %2, i32 %3) {
%5 = alloca double*
%6 = alloca i32
%7 = alloca double*
%8 = alloca i32
%9 = alloca i32
%10 = alloca double
%11 = alloca double*
%12 = alloca double*
store double* %0, double** %5
store i32 %1, i32* %6
store double* %2, double** %7
store i32 %3, i32* %8
%13 = call i32 @sqrt(i32 3)
%14 = sitofp i32 %13 to double
%15 = fadd double -2.000000e+00, %14
store double %15, double* %10
%16 = load i32, i32* %6
%17 = call double* @calloc(i32 %16, i32 8)
store double* %17, double** %11
%18 = load double*, double** %11
%19 = getelementptr inbounds double, double* %18, i64 0
store double 0.000000e+00, double* %19
store i32 0, i32* %9
br label %20
20: ; preds = %37, %4
%21 = load i32, i32* %9
%22 = load i32, i32* %8
%23 = icmp slt i32 %21, %22
br i1 %23, label %24, label %40
24: ; preds = %20
%25 = load double*, double** %5
%26 = load i32, i32* %9
%27 = sext i32 %26 to i64
%28 = getelementptr inbounds double, double* %25, i64 %27
%29 = load double, double* %28
%30 = load double, double* %10
%31 = load i32, i32* %9
%32 = call double @pow(double %30, i32 %31)
%33 = load double*, double** %11
%34 = getelementptr inbounds double, double* %33, i64 0
%35 = load double, double* %34
%36 = call double @llvm.fmuladd.f64(double %29, double %32, double %35)
store double %36, double* %34
br label %37
37: ; preds = %24
%38 = load i32, i32* %9
%39 = add nsw i32 %38, 1
store i32 %39, i32* %9
br label %20
40: ; preds = %20
store i32 1, i32* %9
br label %41
41: ; preds = %63, %40
%42 = load i32, i32* %9
%43 = load i32, i32* %6
%44 = icmp slt i32 %42, %43
br i1 %44, label %45, label %66
45: ; preds = %41
%46 = load double*, double** %5
%47 = load i32, i32* %9
%48 = sext i32 %47 to i64
%49 = getelementptr inbounds double, double* %46, i64 %48
%50 = load double, double* %49
%51 = load double, double* %10
%52 = load double*, double** %11
%53 = load i32, i32* %9
%54 = sub nsw i32 %53, 1
%55 = sext i32 %54 to i64
%56 = getelementptr inbounds double, double* %52, i64 %55
%57 = load double, double* %56
%58 = call double @llvm.fmuladd.f64(double %51, double %57, double %50)
%59 = load double*, double** %11
%60 = load i32, i32* %9
%61 = sext i32 %60 to i64
%62 = getelementptr inbounds double, double* %59, i64 %61
store double %58, double* %62
br label %63
63: ; preds = %45
%64 = load i32, i32* %9
%65 = add nsw i32 %64, 1
store i32 %65, i32* %9
br label %41
66: ; preds = %41
%67 = load double, double* %10
%68 = fneg double %67
%69 = load double, double* %10
%70 = load double, double* %10
%71 = fneg double %69
%72 = call double @llvm.fmuladd.f64(double %71, double %70, double 1.000000e+00)
%73 = fdiv double %68, %72
%74 = load double*, double** %11
%75 = load i32, i32* %6
%76 = sub nsw i32 %75, 1
%77 = sext i32 %76 to i64
%78 = getelementptr inbounds double, double* %74, i64 %77
%79 = load double, double* %78
%80 = load double, double* %10
%81 = load double*, double** %11
%82 = load i32, i32* %6
%83 = sub nsw i32 %82, 2
%84 = sext i32 %83 to i64
%85 = getelementptr inbounds double, double* %81, i64 %84
%86 = load double, double* %85
%87 = call double @llvm.fmuladd.f64(double %80, double %86, double %79)
%88 = fmul double %73, %87
%89 = load double*, double** %7
%90 = load i32, i32* %6
%91 = sub nsw i32 %90, 1
%92 = sext i32 %91 to i64
%93 = getelementptr inbounds double, double* %89, i64 %92
store double %88, double* %93
%94 = load i32, i32* %6
%95 = sub nsw i32 %94, 2
store i32 %95, i32* %9
br label %96
96: ; preds = %118, %66
%97 = load i32, i32* %9
%98 = icmp sge i32 %97, 0
br i1 %98, label %99, label %121
99: ; preds = %96
%100 = load double, double* %10
%101 = load double*, double** %7
%102 = load i32, i32* %9
%103 = add nsw i32 %102, 1
%104 = sext i32 %103 to i64
%105 = getelementptr inbounds double, double* %101, i64 %104
%106 = load double, double* %105
%107 = load double*, double** %11
%108 = load i32, i32* %9
%109 = sext i32 %108 to i64
%110 = getelementptr inbounds double, double* %107, i64 %109
%111 = load double, double* %110
%112 = fsub double %106, %111
%113 = fmul double %100, %112
%114 = load double*, double** %7
%115 = load i32, i32* %9
%116 = sext i32 %115 to i64
%117 = getelementptr inbounds double, double* %114, i64 %116
store double %113, double* %117
br label %118
118: ; preds = %99
%119 = load i32, i32* %9
%120 = add nsw i32 %119, -1
store i32 %120, i32* %9
br label %96
121: ; preds = %96
%122 = load double*, double** %11
%123 = call i32 @free(double* %122)
ret void
}

LLVM IR (-O2):
define void @fast_spline(double* nocapture %0, i32 %1, double* nocapture %2, i32 %3) {
%5 = tail call i32 @sqrt(i32 3)
%6 = sitofp i32 %5 to double
%7 = fadd double %6, -2.000000e+00
%8 = tail call double* @calloc(i32 %1, i32 8)
store double 0.000000e+00, double* %8
%9 = icmp sgt i32 %3, 0
br i1 %9, label %10, label %34
10: ; preds = %4
%11 = zext i32 %3 to i64
%12 = add nsw i64 %11, -1
%13 = and i64 %11, 3
%14 = icmp ult i64 %12, 3
br i1 %14, label %17, label %15
15: ; preds = %10
%16 = and i64 %11, 4294967292
br label %45
17: ; preds = %45, %10
%18 = phi double [ undef, %10 ], [ %71, %45 ]
%19 = phi double [ 0.000000e+00, %10 ], [ %71, %45 ]
%20 = phi i64 [ 0, %10 ], [ %72, %45 ]
%21 = icmp eq i64 %13, 0
br i1 %21, label %34, label %22
22: ; preds = %17, %22
%23 = phi double [ %30, %22 ], [ %19, %17 ]
%24 = phi i64 [ %31, %22 ], [ %20, %17 ]
%25 = phi i64 [ %32, %22 ], [ 0, %17 ]
%26 = getelementptr inbounds double, double* %0, i64 %24
%27 = load double, double* %26
%28 = trunc i64 %24 to i32
%29 = tail call double @pow(double %7, i32 %28)
%30 = tail call double @llvm.fmuladd.f64(double %27, double %29, double %23)
store double %30, double* %8
%31 = add nuw nsw i64 %24, 1
%32 = add i64 %25, 1
%33 = icmp eq i64 %32, %13
br i1 %33, label %34, label %22
34: ; preds = %17, %22, %4
%35 = phi double [ 0.000000e+00, %4 ], [ %18, %17 ], [ %30, %22 ]
%36 = icmp sgt i32 %1, 1
br i1 %36, label %37, label %116
37: ; preds = %34
%38 = zext i32 %1 to i64
%39 = add nsw i64 %38, -1
%40 = add nsw i64 %38, -2
%41 = and i64 %39, 3
%42 = icmp ult i64 %40, 3
br i1 %42, label %101, label %43
43: ; preds = %37
%44 = and i64 %39, -4
br label %75
45: ; preds = %45, %15
%46 = phi double [ 0.000000e+00, %15 ], [ %71, %45 ]
%47 = phi i64 [ 0, %15 ], [ %72, %45 ]
%48 = phi i64 [ 0, %15 ], [ %73, %45 ]
%49 = getelementptr inbounds double, double* %0, i64 %47
%50 = load double, double* %49
%51 = trunc i64 %47 to i32
%52 = tail call double @pow(double %7, i32 %51)
%53 = tail call double @llvm.fmuladd.f64(double %50, double %52, double %46)
store double %53, double* %8
%54 = or i64 %47, 1
%55 = getelementptr inbounds double, double* %0, i64 %54
%56 = load double, double* %55
%57 = trunc i64 %54 to i32
%58 = tail call double @pow(double %7, i32 %57)
%59 = tail call double @llvm.fmuladd.f64(double %56, double %58, double %53)
store double %59, double* %8
%60 = or i64 %47, 2
%61 = getelementptr inbounds double, double* %0, i64 %60
%62 = load double, double* %61
%63 = trunc i64 %60 to i32
%64 = tail call double @pow(double %7, i32 %63)
%65 = tail call double @llvm.fmuladd.f64(double %62, double %64, double %59)
store double %65, double* %8
%66 = or i64 %47, 3
%67 = getelementptr inbounds double, double* %0, i64 %66
%68 = load double, double* %67
%69 = trunc i64 %66 to i32
%70 = tail call double @pow(double %7, i32 %69)
%71 = tail call double @llvm.fmuladd.f64(double %68, double %70, double %65)
store double %71, double* %8
%72 = add nuw nsw i64 %47, 4
%73 = add i64 %48, 4
%74 = icmp eq i64 %73, %16
br i1 %74, label %17, label %45
75: ; preds = %75, %43
%76 = phi double [ %35, %43 ], [ %96, %75 ]
%77 = phi i64 [ 1, %43 ], [ %98, %75 ]
%78 = phi i64 [ 0, %43 ], [ %99, %75 ]
%79 = getelementptr inbounds double, double* %0, i64 %77
%80 = load double, double* %79
%81 = tail call double @llvm.fmuladd.f64(double %7, double %76, double %80)
%82 = getelementptr inbounds double, double* %8, i64 %77
store double %81, double* %82
%83 = add nuw nsw i64 %77, 1
%84 = getelementptr inbounds double, double* %0, i64 %83
%85 = load double, double* %84
%86 = tail call double @llvm.fmuladd.f64(double %7, double %81, double %85)
%87 = getelementptr inbounds double, double* %8, i64 %83
store double %86, double* %87
%88 = add nuw nsw i64 %77, 2
%89 = getelementptr inbounds double, double* %0, i64 %88
%90 = load double, double* %89
%91 = tail call double @llvm.fmuladd.f64(double %7, double %86, double %90)
%92 = getelementptr inbounds double, double* %8, i64 %88
store double %91, double* %92
%93 = add nuw nsw i64 %77, 3
%94 = getelementptr inbounds double, double* %0, i64 %93
%95 = load double, double* %94
%96 = tail call double @llvm.fmuladd.f64(double %7, double %91, double %95)
%97 = getelementptr inbounds double, double* %8, i64 %93
store double %96, double* %97
%98 = add nuw nsw i64 %77, 4
%99 = add i64 %78, 4
%100 = icmp eq i64 %99, %44
br i1 %100, label %101, label %75
101: ; preds = %75, %37
%102 = phi double [ %35, %37 ], [ %96, %75 ]
%103 = phi i64 [ 1, %37 ], [ %98, %75 ]
%104 = icmp eq i64 %41, 0
br i1 %104, label %116, label %105
105: ; preds = %101, %105
%106 = phi double [ %111, %105 ], [ %102, %101 ]
%107 = phi i64 [ %113, %105 ], [ %103, %101 ]
%108 = phi i64 [ %114, %105 ], [ 0, %101 ]
%109 = getelementptr inbounds double, double* %0, i64 %107
%110 = load double, double* %109
%111 = tail call double @llvm.fmuladd.f64(double %7, double %106, double %110)
%112 = getelementptr inbounds double, double* %8, i64 %107
store double %111, double* %112
%113 = add nuw nsw i64 %107, 1
%114 = add i64 %108, 1
%115 = icmp eq i64 %114, %41
br i1 %115, label %116, label %105
116: ; preds = %101, %105, %34
%117 = fneg double %7
%118 = tail call double @llvm.fmuladd.f64(double %117, double %7, double 1.000000e+00)
%119 = fdiv double %117, %118
%120 = add nsw i32 %1, -1
%121 = sext i32 %120 to i64
%122 = getelementptr inbounds double, double* %8, i64 %121
%123 = load double, double* %122
%124 = add nsw i32 %1, -2
%125 = sext i32 %124 to i64
%126 = getelementptr inbounds double, double* %8, i64 %125
%127 = load double, double* %126
%128 = tail call double @llvm.fmuladd.f64(double %7, double %127, double %123)
%129 = fmul double %119, %128
%130 = getelementptr inbounds double, double* %2, i64 %121
store double %129, double* %130
br i1 %36, label %131, label %166
131: ; preds = %116
%132 = zext i32 %124 to i64
%133 = and i64 %132, 1
%134 = icmp eq i64 %133, 0
br i1 %134, label %135, label %145
135: ; preds = %131
%136 = add nuw nsw i64 %132, 1
%137 = getelementptr inbounds double, double* %2, i64 %136
%138 = load double, double* %137
%139 = getelementptr inbounds double, double* %8, i64 %132
%140 = load double, double* %139
%141 = fsub double %138, %140
%142 = fmul double %7, %141
%143 = getelementptr inbounds double, double* %2, i64 %132
store double %142, double* %143
%144 = add nsw i64 %132, -1
br label %145
145: ; preds = %135, %131
%146 = phi i64 [ %132, %131 ], [ %144, %135 ]
%147 = icmp eq i32 %124, 0
br i1 %147, label %166, label %148
148: ; preds = %145, %148
%149 = phi i64 [ %164, %148 ], [ %146, %145 ]
%150 = add nuw nsw i64 %149, 1
%151 = getelementptr inbounds double, double* %2, i64 %150
%152 = load double, double* %151
%153 = getelementptr inbounds double, double* %8, i64 %149
%154 = load double, double* %153
%155 = fsub double %152, %154
%156 = fmul double %7, %155
%157 = getelementptr inbounds double, double* %2, i64 %149
store double %156, double* %157
%158 = add nsw i64 %149, -1
%159 = getelementptr inbounds double, double* %8, i64 %158
%160 = load double, double* %159
%161 = fsub double %156, %160
%162 = fmul double %7, %161
%163 = getelementptr inbounds double, double* %2, i64 %158
store double %162, double* %163
%164 = add nsw i64 %149, -2
%165 = icmp eq i64 %158, 0
br i1 %165, label %166, label %148
166: ; preds = %145, %148, %116
%167 = tail call i32 @free(double* %8)
ret void
}