#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*/ clean_sort(double **, int *);
int /*<<< orphan*/ free(double *);
scalar_t__ malloc(int);
int nb_cols_lfmm(char *);
int nb_lines(char *, int);
int /*<<< orphan*/ print_summary_tracyWidom(int, char *, char *);
int /*<<< orphan*/ printf(char *, char *);
int /*<<< orphan*/ read_data_double(char *, int, int, double *);
int /*<<< orphan*/ tw(double *, double *, double *, double *, int);
int /*<<< orphan*/ write_data_tracyWidom(char *, int, double *, double *,
                                         double *, double *, double *);
void tracyWidom(char *input_file, char *output_file) {
  double *values, *pvalues, *twstat, *effectn, *percentage;
  int i, N, M;
  double sum;
  M = nb_cols_lfmm(input_file);
  N = nb_lines(input_file, M);
  if (M != 1) {
    printf("Tracy-Widom: Error %s has more than one column\n", input_file);
  }
  print_summary_tracyWidom(N, input_file, output_file);
  values = (double *)malloc(N * sizeof(double));
  read_data_double(input_file, N, 1, values);
  clean_sort(&values, &N);
  pvalues = (double *)malloc(N * sizeof(double));
  twstat = (double *)malloc(N * sizeof(double));
  effectn = (double *)malloc(N * sizeof(double));
  percentage = (double *)malloc(N * sizeof(double));
  tw(values, pvalues, twstat, effectn, N);
  sum = 0.0;
  for (i = 0; i < N; i++)
    sum += values[i];
  for (i = 0; i < N; i++)
    percentage[i] = values[i] / sum;
  write_data_tracyWidom(output_file, N, values, pvalues, twstat, effectn,
                        percentage);
  free(values);
  free(pvalues);
  free(twstat);
  free(effectn);
  free(percentage);
}


x86 assembly (clang, -O0):
<tracyWidom>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    -0x8(%rbp),%rdi
call   10b0 <nb_cols_lfmm@plt>
mov    %eax,-0x44(%rbp)
mov    -0x8(%rbp),%rdi
mov    -0x44(%rbp),%esi
call   1050 <nb_lines@plt>
mov    %eax,-0x40(%rbp)
cmpl   $0x1,-0x44(%rbp)
je     11e5 <tracyWidom+0x45>
mov    -0x8(%rbp),%rsi
lea    0xe20(%rip),%rdi
call   1060 <printf@plt>
mov    -0x40(%rbp),%edi
mov    -0x8(%rbp),%rsi
mov    -0x10(%rbp),%rdx
call   10c0 <print_summary_tracyWidom@plt>
movslq -0x40(%rbp),%rax
shl    $0x3,%rax
mov    %eax,%edi
call   10a0 <malloc@plt>
mov    %rax,-0x18(%rbp)
mov    -0x8(%rbp),%rdi
mov    -0x40(%rbp),%esi
mov    -0x18(%rbp),%rcx
mov    $0x1,%edx
call   1070 <read_data_double@plt>
lea    -0x18(%rbp),%rdi
lea    -0x40(%rbp),%rsi
call   1090 <clean_sort@plt>
movslq -0x40(%rbp),%rax
shl    $0x3,%rax
mov    %eax,%edi
call   10a0 <malloc@plt>
mov    %rax,-0x20(%rbp)
movslq -0x40(%rbp),%rax
shl    $0x3,%rax
mov    %eax,%edi
call   10a0 <malloc@plt>
mov    %rax,-0x28(%rbp)
movslq -0x40(%rbp),%rax
shl    $0x3,%rax
mov    %eax,%edi
call   10a0 <malloc@plt>
mov    %rax,-0x30(%rbp)
movslq -0x40(%rbp),%rax
shl    $0x3,%rax
mov    %eax,%edi
call   10a0 <malloc@plt>
mov    %rax,-0x38(%rbp)
mov    -0x18(%rbp),%rdi
mov    -0x20(%rbp),%rsi
mov    -0x28(%rbp),%rdx
mov    -0x30(%rbp),%rcx
mov    -0x40(%rbp),%r8d
call   1080 <tw@plt>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x50(%rbp)
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x40(%rbp),%eax
jge    12cf <tracyWidom+0x12f>
mov    -0x18(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    129e <tracyWidom+0xfe>
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x40(%rbp),%eax
jge    130f <tracyWidom+0x16f>
mov    -0x18(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  -0x50(%rbp),%xmm0
mov    -0x38(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    12d6 <tracyWidom+0x136>
mov    -0x10(%rbp),%rdi
mov    -0x40(%rbp),%esi
mov    -0x18(%rbp),%rdx
mov    -0x20(%rbp),%rcx
mov    -0x28(%rbp),%r8
mov    -0x30(%rbp),%r9
mov    -0x38(%rbp),%rax
mov    %rax,(%rsp)
call   1040 <write_data_tracyWidom@plt>
mov    -0x18(%rbp),%rdi
call   1030 <free@plt>
mov    -0x20(%rbp),%rdi
call   1030 <free@plt>
mov    -0x28(%rbp),%rdi
call   1030 <free@plt>
mov    -0x30(%rbp),%rdi
call   1030 <free@plt>
mov    -0x38(%rbp),%rdi
call   1030 <free@plt>
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<tracyWidom>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x18,%rsp
mov    %rsi,%r14
mov    %rdi,%r15
call   10b0 <nb_cols_lfmm@plt>
mov    %eax,%ebx
mov    %r15,%rdi
mov    %eax,%esi
call   1050 <nb_lines@plt>
mov    %eax,%ebp
mov    %eax,0xc(%rsp)
cmp    $0x1,%ebx
je     11df <tracyWidom+0x3f>
lea    0xe29(%rip),%rdi
mov    %r15,%rsi
call   1060 <printf@plt>
mov    %ebp,%edi
mov    %r15,%rsi
mov    %r14,%rdx
call   10c0 <print_summary_tracyWidom@plt>
mov    %ebp,%edi
shl    $0x3,%edi
call   10a0 <malloc@plt>
mov    %rax,0x10(%rsp)
mov    %r15,%rdi
mov    %ebp,%esi
mov    $0x1,%edx
mov    %rax,%rcx
call   1070 <read_data_double@plt>
lea    0x10(%rsp),%rdi
lea    0xc(%rsp),%rsi
call   1090 <clean_sort@plt>
mov    0xc(%rsp),%edi
shl    $0x3,%edi
call   10a0 <malloc@plt>
mov    %rax,%r15
mov    0xc(%rsp),%edi
shl    $0x3,%edi
call   10a0 <malloc@plt>
mov    %rax,%r12
mov    0xc(%rsp),%edi
shl    $0x3,%edi
call   10a0 <malloc@plt>
mov    %rax,%r13
mov    0xc(%rsp),%edi
shl    $0x3,%edi
call   10a0 <malloc@plt>
mov    %rax,%rbx
mov    0x10(%rsp),%rdi
mov    0xc(%rsp),%r8d
mov    %r15,%rsi
mov    %r12,%rdx
mov    %r13,%rcx
call   1080 <tw@plt>
mov    0xc(%rsp),%esi
test   %esi,%esi
jle    1389 <tracyWidom+0x1e9>
mov    0x10(%rsp),%rcx
lea    -0x1(%rsi),%rdx
mov    %esi,%eax
and    $0x7,%eax
cmp    $0x7,%rdx
jae    1298 <tracyWidom+0xf8>
xorpd  %xmm0,%xmm0
xor    %edx,%edx
jmp    12e8 <tracyWidom+0x148>
mov    %esi,%edi
and    $0xfffffff8,%edi
xorpd  %xmm0,%xmm0
xor    %edx,%edx
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
addsd  (%rcx,%rdx,8),%xmm0
addsd  0x8(%rcx,%rdx,8),%xmm0
addsd  0x10(%rcx,%rdx,8),%xmm0
addsd  0x18(%rcx,%rdx,8),%xmm0
addsd  0x20(%rcx,%rdx,8),%xmm0
addsd  0x28(%rcx,%rdx,8),%xmm0
addsd  0x30(%rcx,%rdx,8),%xmm0
addsd  0x38(%rcx,%rdx,8),%xmm0
add    $0x8,%rdx
cmp    %rdx,%rdi
jne    12b0 <tracyWidom+0x110>
test   %rax,%rax
je     130e <tracyWidom+0x16e>
lea    (%rcx,%rdx,8),%rcx
xor    %edx,%edx
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
addsd  (%rcx,%rdx,8),%xmm0
add    $0x1,%rdx
cmp    %rdx,%rax
jne    1300 <tracyWidom+0x160>
test   %esi,%esi
jle    1389 <tracyWidom+0x1e9>
mov    0x10(%rsp),%rdx
cmp    $0x1,%esi
je     1336 <tracyWidom+0x196>
lea    (%rdx,%rsi,8),%rax
cmp    %rax,%rbx
jae    13dc <tracyWidom+0x23c>
lea    (%rbx,%rsi,8),%rax
cmp    %rax,%rdx
jae    13dc <tracyWidom+0x23c>
xor    %eax,%eax
mov    %rax,%rcx
not    %rcx
test   $0x1,%sil
je     1356 <tracyWidom+0x1b6>
movsd  (%rdx,%rax,8),%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,(%rbx,%rax,8)
or     $0x1,%rax
add    %rsi,%rcx
je     138e <tracyWidom+0x1ee>
nopl   0x0(%rax,%rax,1)
movsd  (%rdx,%rax,8),%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,(%rbx,%rax,8)
movsd  0x8(%rdx,%rax,8),%xmm1
divsd  %xmm0,%xmm1
movsd  %xmm1,0x8(%rbx,%rax,8)
add    $0x2,%rax
cmp    %rax,%rsi
jne    1360 <tracyWidom+0x1c0>
jmp    138e <tracyWidom+0x1ee>
mov    0x10(%rsp),%rdx
mov    %rbx,(%rsp)
mov    %r14,%rdi
mov    %r15,%rcx
mov    %r12,%r8
mov    %r13,%r9
call   1040 <write_data_tracyWidom@plt>
mov    0x10(%rsp),%rdi
call   1030 <free@plt>
mov    %r15,%rdi
call   1030 <free@plt>
mov    %r12,%rdi
call   1030 <free@plt>
mov    %r13,%rdi
call   1030 <free@plt>
mov    %rbx,%rdi
call   1030 <free@plt>
add    $0x18,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
mov    %esi,%eax
and    $0xfffffffe,%eax
movapd %xmm0,%xmm1
unpcklpd %xmm0,%xmm1
lea    -0x2(%rax),%rcx
mov    %rcx,%rbp
shr    %rbp
add    $0x1,%rbp
test   %rcx,%rcx
je     145a <tracyWidom+0x2ba>
mov    %rbp,%rcx
and    $0xfffffffffffffffe,%rcx
xor    %edi,%edi
cs nopw 0x0(%rax,%rax,1)
nop
movupd (%rdx,%rdi,8),%xmm2
divpd  %xmm1,%xmm2
movupd %xmm2,(%rbx,%rdi,8)
movupd 0x10(%rdx,%rdi,8),%xmm2
divpd  %xmm1,%xmm2
movupd %xmm2,0x10(%rbx,%rdi,8)
add    $0x4,%rdi
add    $0xfffffffffffffffe,%rcx
jne    1410 <tracyWidom+0x270>
test   $0x1,%bpl
je     144c <tracyWidom+0x2ac>
movupd (%rdx,%rdi,8),%xmm2
divpd  %xmm1,%xmm2
movupd %xmm2,(%rbx,%rdi,8)
cmp    %rsi,%rax
je     138e <tracyWidom+0x1ee>
jmp    1338 <tracyWidom+0x198>
xor    %edi,%edi
test   $0x1,%bpl
jne    143e <tracyWidom+0x29e>
jmp    144c <tracyWidom+0x2ac>

LLVM IR (-O0):
@.str = [48 x i8] c"Tracy-Widom: Error %s has more than one column\0A\00"
define void @tracyWidom(i8* %0, i8* %1) {
%3 = alloca i8*
%4 = alloca i8*
%5 = alloca double*
%6 = alloca double*
%7 = alloca double*
%8 = alloca double*
%9 = alloca double*
%10 = alloca i32
%11 = alloca i32
%12 = alloca i32
%13 = alloca double
store i8* %0, i8** %3
store i8* %1, i8** %4
%14 = load i8*, i8** %3
%15 = call i32 @nb_cols_lfmm(i8* %14)
store i32 %15, i32* %12
%16 = load i8*, i8** %3
%17 = load i32, i32* %12
%18 = call i32 @nb_lines(i8* %16, i32 %17)
store i32 %18, i32* %11
%19 = load i32, i32* %12
%20 = icmp ne i32 %19, 1
br i1 %20, label %21, label %24
21: ; preds = %2
%22 = load i8*, i8** %3
%23 = call i32 @printf(i8* getelementptr inbounds ([48 x i8], [48 x i8]* @.str, i64 0, i64 0), i8* %22)
br label %24
24: ; preds = %21, %2
%25 = load i32, i32* %11
%26 = load i8*, i8** %3
%27 = load i8*, i8** %4
%28 = call i32 @print_summary_tracyWidom(i32 %25, i8* %26, i8* %27)
%29 = load i32, i32* %11
%30 = sext i32 %29 to i64
%31 = mul i64 %30, 8
%32 = trunc i64 %31 to i32
%33 = call i64 @malloc(i32 %32)
%34 = inttoptr i64 %33 to double*
store double* %34, double** %5
%35 = load i8*, i8** %3
%36 = load i32, i32* %11
%37 = load double*, double** %5
%38 = call i32 @read_data_double(i8* %35, i32 %36, i32 1, double* %37)
%39 = call i32 @clean_sort(double** %5, i32* %11)
%40 = load i32, i32* %11
%41 = sext i32 %40 to i64
%42 = mul i64 %41, 8
%43 = trunc i64 %42 to i32
%44 = call i64 @malloc(i32 %43)
%45 = inttoptr i64 %44 to double*
store double* %45, double** %6
%46 = load i32, i32* %11
%47 = sext i32 %46 to i64
%48 = mul i64 %47, 8
%49 = trunc i64 %48 to i32
%50 = call i64 @malloc(i32 %49)
%51 = inttoptr i64 %50 to double*
store double* %51, double** %7
%52 = load i32, i32* %11
%53 = sext i32 %52 to i64
%54 = mul i64 %53, 8
%55 = trunc i64 %54 to i32
%56 = call i64 @malloc(i32 %55)
%57 = inttoptr i64 %56 to double*
store double* %57, double** %8
%58 = load i32, i32* %11
%59 = sext i32 %58 to i64
%60 = mul i64 %59, 8
%61 = trunc i64 %60 to i32
%62 = call i64 @malloc(i32 %61)
%63 = inttoptr i64 %62 to double*
store double* %63, double** %9
%64 = load double*, double** %5
%65 = load double*, double** %6
%66 = load double*, double** %7
%67 = load double*, double** %8
%68 = load i32, i32* %11
%69 = call i32 @tw(double* %64, double* %65, double* %66, double* %67, i32 %68)
store double 0.000000e+00, double* %13
store i32 0, i32* %10
br label %70
70: ; preds = %82, %24
%71 = load i32, i32* %10
%72 = load i32, i32* %11
%73 = icmp slt i32 %71, %72
br i1 %73, label %74, label %85
74: ; preds = %70
%75 = load double*, double** %5
%76 = load i32, i32* %10
%77 = sext i32 %76 to i64
%78 = getelementptr inbounds double, double* %75, i64 %77
%79 = load double, double* %78
%80 = load double, double* %13
%81 = fadd double %80, %79
store double %81, double* %13
br label %82
82: ; preds = %74
%83 = load i32, i32* %10
%84 = add nsw i32 %83, 1
store i32 %84, i32* %10
br label %70
85: ; preds = %70
store i32 0, i32* %10
br label %86
86: ; preds = %102, %85
%87 = load i32, i32* %10
%88 = load i32, i32* %11
%89 = icmp slt i32 %87, %88
br i1 %89, label %90, label %105
90: ; preds = %86
%91 = load double*, double** %5
%92 = load i32, i32* %10
%93 = sext i32 %92 to i64
%94 = getelementptr inbounds double, double* %91, i64 %93
%95 = load double, double* %94
%96 = load double, double* %13
%97 = fdiv double %95, %96
%98 = load double*, double** %9
%99 = load i32, i32* %10
%100 = sext i32 %99 to i64
%101 = getelementptr inbounds double, double* %98, i64 %100
store double %97, double* %101
br label %102
102: ; preds = %90
%103 = load i32, i32* %10
%104 = add nsw i32 %103, 1
store i32 %104, i32* %10
br label %86
105: ; preds = %86
%106 = load i8*, i8** %4
%107 = load i32, i32* %11
%108 = load double*, double** %5
%109 = load double*, double** %6
%110 = load double*, double** %7
%111 = load double*, double** %8
%112 = load double*, double** %9
%113 = call i32 @write_data_tracyWidom(i8* %106, i32 %107, double* %108, double* %109, double* %110, double* %111, double* %112)
%114 = load double*, double** %5
%115 = call i32 @free(double* %114)
%116 = load double*, double** %6
%117 = call i32 @free(double* %116)
%118 = load double*, double** %7
%119 = call i32 @free(double* %118)
%120 = load double*, double** %8
%121 = call i32 @free(double* %120)
%122 = load double*, double** %9
%123 = call i32 @free(double* %122)
ret void
}

LLVM IR (-O2):
@.str = [48 x i8] c"Tracy-Widom: Error %s has more than one column\0A\00"
define void @tracyWidom(i8* %0, i8* %1) {
%3 = alloca double*
%4 = alloca i32
%5 = bitcast double** %3 to i8*
%6 = bitcast i32* %4 to i8*
%7 = tail call i32 @nb_cols_lfmm(i8* %0)
%8 = tail call i32 @nb_lines(i8* %0, i32 %7)
store i32 %8, i32* %4
%9 = icmp eq i32 %7, 1
br i1 %9, label %12, label %10
10: ; preds = %2
%11 = tail call i32 @printf(i8* getelementptr inbounds ([48 x i8], [48 x i8]* @.str, i64 0, i64 0), i8* %0)
br label %12
12: ; preds = %10, %2
%13 = tail call i32 @print_summary_tracyWidom(i32 %8, i8* %0, i8* %1)
%14 = shl i32 %8, 3
%15 = tail call i64 @malloc(i32 %14)
%16 = inttoptr i64 %15 to double*
store double* %16, double** %3
%17 = tail call i32 @read_data_double(i8* %0, i32 %8, i32 1, double* %16)
%18 = call i32 @clean_sort(double** %3, i32* %4)
%19 = load i32, i32* %4
%20 = shl i32 %19, 3
%21 = call i64 @malloc(i32 %20)
%22 = inttoptr i64 %21 to double*
%23 = load i32, i32* %4
%24 = shl i32 %23, 3
%25 = call i64 @malloc(i32 %24)
%26 = inttoptr i64 %25 to double*
%27 = load i32, i32* %4
%28 = shl i32 %27, 3
%29 = call i64 @malloc(i32 %28)
%30 = inttoptr i64 %29 to double*
%31 = load i32, i32* %4
%32 = shl i32 %31, 3
%33 = call i64 @malloc(i32 %32)
%34 = inttoptr i64 %33 to double*
%35 = inttoptr i64 %33 to i8*
%36 = load double*, double** %3
%37 = load i32, i32* %4
%38 = call i32 @tw(double* %36, double* %22, double* %26, double* %30, i32 %37)
%39 = load i32, i32* %4
%40 = icmp sgt i32 %39, 0
br i1 %40, label %41, label %66
41: ; preds = %12
%42 = load double*, double** %3
%43 = zext i32 %39 to i64
%44 = add nsw i64 %43, -1
%45 = and i64 %43, 7
%46 = icmp ult i64 %44, 7
br i1 %46, label %49, label %47
47: ; preds = %41
%48 = and i64 %43, 4294967288
br label %136
49: ; preds = %136, %41
%50 = phi double [ undef, %41 ], [ %170, %136 ]
%51 = phi i64 [ 0, %41 ], [ %171, %136 ]
%52 = phi double [ 0.000000e+00, %41 ], [ %170, %136 ]
%53 = icmp eq i64 %45, 0
br i1 %53, label %64, label %54
54: ; preds = %49, %54
%55 = phi i64 [ %61, %54 ], [ %51, %49 ]
%56 = phi double [ %60, %54 ], [ %52, %49 ]
%57 = phi i64 [ %62, %54 ], [ 0, %49 ]
%58 = getelementptr inbounds double, double* %42, i64 %55
%59 = load double, double* %58
%60 = fadd double %56, %59
%61 = add nuw nsw i64 %55, 1
%62 = add i64 %57, 1
%63 = icmp eq i64 %62, %45
br i1 %63, label %64, label %54
64: ; preds = %54, %49
%65 = phi double [ %50, %49 ], [ %60, %54 ]
br i1 %40, label %68, label %66
66: ; preds = %12, %64
%67 = load double*, double** %3
br label %187
68: ; preds = %64
%69 = load double*, double** %3
%70 = zext i32 %39 to i64
%71 = icmp eq i32 %39, 1
br i1 %71, label %121, label %72
72: ; preds = %68
%73 = getelementptr double, double* %34, i64 %43
%74 = getelementptr double, double* %69, i64 %43
%75 = bitcast double* %74 to i8*
%76 = icmp ult i8* %35, %75
%77 = icmp ult double* %69, %73
%78 = and i1 %76, %77
br i1 %78, label %121, label %79
79: ; preds = %72
%80 = and i64 %43, 4294967294
%81 = insertelement <2 x double> poison, double %65, i64 0
%82 = shufflevector <2 x double> %81, <2 x double> poison, <2 x i32> zeroinitializer
%83 = add nsw i64 %80, -2
%84 = lshr exact i64 %83, 1
%85 = add nuw i64 %84, 1
%86 = and i64 %85, 1
%87 = icmp eq i64 %83, 0
br i1 %87, label %109, label %88
88: ; preds = %79
%89 = and i64 %85, -2
br label %90
90: ; preds = %90, %88
%91 = phi i64 [ 0, %88 ], [ %106, %90 ]
%92 = phi i64 [ 0, %88 ], [ %107, %90 ]
%93 = getelementptr inbounds double, double* %69, i64 %91
%94 = bitcast double* %93 to <2 x double>*
%95 = load <2 x double>, <2 x double>* %94
%96 = fdiv <2 x double> %95, %82
%97 = getelementptr inbounds double, double* %34, i64 %91
%98 = bitcast double* %97 to <2 x double>*
store <2 x double> %96, <2 x double>* %98
%99 = or i64 %91, 2
%100 = getelementptr inbounds double, double* %69, i64 %99
%101 = bitcast double* %100 to <2 x double>*
%102 = load <2 x double>, <2 x double>* %101
%103 = fdiv <2 x double> %102, %82
%104 = getelementptr inbounds double, double* %34, i64 %99
%105 = bitcast double* %104 to <2 x double>*
store <2 x double> %103, <2 x double>* %105
%106 = add nuw i64 %91, 4
%107 = add nuw i64 %92, 2
%108 = icmp eq i64 %107, %89
br i1 %108, label %109, label %90
109: ; preds = %90, %79
%110 = phi i64 [ 0, %79 ], [ %106, %90 ]
%111 = icmp eq i64 %86, 0
br i1 %111, label %119, label %112
112: ; preds = %109
%113 = getelementptr inbounds double, double* %69, i64 %110
%114 = bitcast double* %113 to <2 x double>*
%115 = load <2 x double>, <2 x double>* %114
%116 = fdiv <2 x double> %115, %82
%117 = getelementptr inbounds double, double* %34, i64 %110
%118 = bitcast double* %117 to <2 x double>*
store <2 x double> %116, <2 x double>* %118
br label %119
119: ; preds = %109, %112
%120 = icmp eq i64 %80, %43
br i1 %120, label %187, label %121
121: ; preds = %72, %68, %119
%122 = phi i64 [ 0, %72 ], [ 0, %68 ], [ %80, %119 ]
%123 = xor i64 %122, -1
%124 = and i64 %43, 1
%125 = icmp eq i64 %124, 0
br i1 %125, label %132, label %126
126: ; preds = %121
%127 = getelementptr inbounds double, double* %69, i64 %122
%128 = load double, double* %127
%129 = fdiv double %128, %65
%130 = getelementptr inbounds double, double* %34, i64 %122
store double %129, double* %130
%131 = or i64 %122, 1
br label %132
132: ; preds = %126, %121
%133 = phi i64 [ %122, %121 ], [ %131, %126 ]
%134 = sub nsw i64 0, %43
%135 = icmp eq i64 %123, %134
br i1 %135, label %187, label %174
136: ; preds = %136, %47
%137 = phi i64 [ 0, %47 ], [ %171, %136 ]
%138 = phi double [ 0.000000e+00, %47 ], [ %170, %136 ]
%139 = phi i64 [ 0, %47 ], [ %172, %136 ]
%140 = getelementptr inbounds double, double* %42, i64 %137
%141 = load double, double* %140
%142 = fadd double %138, %141
%143 = or i64 %137, 1
%144 = getelementptr inbounds double, double* %42, i64 %143
%145 = load double, double* %144
%146 = fadd double %142, %145
%147 = or i64 %137, 2
%148 = getelementptr inbounds double, double* %42, i64 %147
%149 = load double, double* %148
%150 = fadd double %146, %149
%151 = or i64 %137, 3
%152 = getelementptr inbounds double, double* %42, i64 %151
%153 = load double, double* %152
%154 = fadd double %150, %153
%155 = or i64 %137, 4
%156 = getelementptr inbounds double, double* %42, i64 %155
%157 = load double, double* %156
%158 = fadd double %154, %157
%159 = or i64 %137, 5
%160 = getelementptr inbounds double, double* %42, i64 %159
%161 = load double, double* %160
%162 = fadd double %158, %161
%163 = or i64 %137, 6
%164 = getelementptr inbounds double, double* %42, i64 %163
%165 = load double, double* %164
%166 = fadd double %162, %165
%167 = or i64 %137, 7
%168 = getelementptr inbounds double, double* %42, i64 %167
%169 = load double, double* %168
%170 = fadd double %166, %169
%171 = add nuw nsw i64 %137, 8
%172 = add i64 %139, 8
%173 = icmp eq i64 %172, %48
br i1 %173, label %49, label %136
174: ; preds = %132, %174
%175 = phi i64 [ %185, %174 ], [ %133, %132 ]
%176 = getelementptr inbounds double, double* %69, i64 %175
%177 = load double, double* %176
%178 = fdiv double %177, %65
%179 = getelementptr inbounds double, double* %34, i64 %175
store double %178, double* %179
%180 = add nuw nsw i64 %175, 1
%181 = getelementptr inbounds double, double* %69, i64 %180
%182 = load double, double* %181
%183 = fdiv double %182, %65
%184 = getelementptr inbounds double, double* %34, i64 %180
store double %183, double* %184
%185 = add nuw nsw i64 %175, 2
%186 = icmp eq i64 %185, %70
br i1 %186, label %187, label %174
187: ; preds = %132, %174, %119, %66
%188 = phi double* [ %67, %66 ], [ %69, %119 ], [ %69, %174 ], [ %69, %132 ]
%189 = call i32 @write_data_tracyWidom(i8* %1, i32 %39, double* %188, double* %22, double* %26, double* %30, double* %34)
%190 = load double*, double** %3
%191 = call i32 @free(double* %190)
%192 = call i32 @free(double* %22)
%193 = call i32 @free(double* %26)
%194 = call i32 @free(double* %30)
%195 = call i32 @free(double* %34)
ret void
}