#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 betai(double, double, double);
int /*<<< orphan*/ fprintf(int /*<<< orphan*/, char *);
double sqrt(double);
int /*<<< orphan*/ stderr;
int ttest(double *data1, int n1, double *data2, int n2, double *t,
          double *prob) {
  double var1, var2, df, ave1, ave2;
  int i;
  double dev;
  if (n1 <= 1) {
    fprintf(stderr, "ttest: n1 must be > 1");
    return 0;
  } else {
    ave1 = .0;
    for (i = 0; i < n1; i++)
      ave1 += data1[i];
    ave1 /= n1;
    var1 = .0;
    for (i = 0; i < n1; i++) {
      dev = data1[i] - ave1;
      var1 += dev * dev;
    }
    var1 /= (n1 - 1);
  }
  if (n2 <= 1) {
    fprintf(stderr, "ttest: n2 must be > 1");
    return 0;
  } else {
    ave2 = .0;
    for (i = 0; i < n2; i++)
      ave2 += data2[i];
    ave2 /= n2;
    var2 = .0;
    for (i = 0; i < n2; i++) {
      dev = data2[i] - ave2;
      var2 += dev * dev;
    }
    var2 /= (n2 - 1);
  }
  *t = (ave1 - ave2) / sqrt(var1 / n1 + var2 / n2);
  df = (var1 / n1 + var2 / n2) * (var1 / n1 + var2 / n2) /
       ((var1 / n1) * (var1 / n1) / (n1 - 1) +
        (var2 / n2) * (var2 / n2) / (n2 - 1));
  *prob = betai(0.5 * df, 0.5, df / (df + (*t) * (*t)));
  return 1;
}


x86 assembly (clang, -O0):
<ttest>:
push   %rbp
mov    %rsp,%rbp
sub    $0x80,%rsp
mov    %rdi,-0x10(%rbp)
mov    %esi,-0x14(%rbp)
mov    %rdx,-0x20(%rbp)
mov    %ecx,-0x24(%rbp)
mov    %r8,-0x30(%rbp)
mov    %r9,-0x38(%rbp)
cmpl   $0x1,-0x14(%rbp)
jg     117c <ttest+0x4c>
mov    0x2e7e(%rip),%rax
mov    (%rax),%edi
lea    0xe9d(%rip),%rsi
call   1030 <fprintf@plt>
movl   $0x0,-0x4(%rbp)
jmp    1479 <ttest+0x349>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x58(%rbp)
movl   $0x0,-0x64(%rbp)
mov    -0x64(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    11bc <ttest+0x8c>
mov    -0x10(%rbp),%rax
movslq -0x64(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x58(%rbp),%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x64(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x64(%rbp)
jmp    118b <ttest+0x5b>
cvtsi2sdl -0x14(%rbp),%xmm1
movsd  -0x58(%rbp),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x40(%rbp)
movl   $0x0,-0x64(%rbp)
mov    -0x64(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    122b <ttest+0xfb>
mov    -0x10(%rbp),%rax
movslq -0x64(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x58(%rbp),%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x70(%rbp),%xmm0
movsd  -0x70(%rbp),%xmm2
movsd  -0x40(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x64(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x64(%rbp)
jmp    11de <ttest+0xae>
mov    -0x14(%rbp),%eax
sub    $0x1,%eax
cvtsi2sd %eax,%xmm1
movsd  -0x40(%rbp),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
cmpl   $0x1,-0x24(%rbp)
jg     126e <ttest+0x13e>
mov    0x2d8c(%rip),%rax
mov    (%rax),%edi
lea    0xdc1(%rip),%rsi
call   1030 <fprintf@plt>
movl   $0x0,-0x4(%rbp)
jmp    1479 <ttest+0x349>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x60(%rbp)
movl   $0x0,-0x64(%rbp)
mov    -0x64(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    12ae <ttest+0x17e>
mov    -0x20(%rbp),%rax
movslq -0x64(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x60(%rbp),%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    -0x64(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x64(%rbp)
jmp    127d <ttest+0x14d>
cvtsi2sdl -0x24(%rbp),%xmm1
movsd  -0x60(%rbp),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x48(%rbp)
movl   $0x0,-0x64(%rbp)
mov    -0x64(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    131d <ttest+0x1ed>
mov    -0x20(%rbp),%rax
movslq -0x64(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x60(%rbp),%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x70(%rbp),%xmm0
movsd  -0x70(%rbp),%xmm2
movsd  -0x48(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x64(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x64(%rbp)
jmp    12d0 <ttest+0x1a0>
mov    -0x24(%rbp),%eax
sub    $0x1,%eax
cvtsi2sd %eax,%xmm1
movsd  -0x48(%rbp),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x58(%rbp),%xmm0
subsd  -0x60(%rbp),%xmm0
movsd  %xmm0,-0x78(%rbp)
movsd  -0x40(%rbp),%xmm0
cvtsi2sdl -0x14(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  -0x48(%rbp),%xmm1
cvtsi2sdl -0x24(%rbp),%xmm2
divsd  %xmm2,%xmm1
addsd  %xmm1,%xmm0
call   1050 <sqrt@plt>
movaps %xmm0,%xmm1
movsd  -0x78(%rbp),%xmm0
divsd  %xmm1,%xmm0
mov    -0x30(%rbp),%rax
movsd  %xmm0,(%rax)
movsd  -0x40(%rbp),%xmm0
cvtsi2sdl -0x14(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  -0x48(%rbp),%xmm1
cvtsi2sdl -0x24(%rbp),%xmm2
divsd  %xmm2,%xmm1
addsd  %xmm1,%xmm0
movsd  -0x40(%rbp),%xmm1
cvtsi2sdl -0x14(%rbp),%xmm2
divsd  %xmm2,%xmm1
movsd  -0x48(%rbp),%xmm2
cvtsi2sdl -0x24(%rbp),%xmm3
divsd  %xmm3,%xmm2
addsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
movsd  -0x40(%rbp),%xmm1
cvtsi2sdl -0x14(%rbp),%xmm2
divsd  %xmm2,%xmm1
movsd  -0x40(%rbp),%xmm2
cvtsi2sdl -0x14(%rbp),%xmm3
divsd  %xmm3,%xmm2
mulsd  %xmm2,%xmm1
mov    -0x14(%rbp),%eax
sub    $0x1,%eax
cvtsi2sd %eax,%xmm2
divsd  %xmm2,%xmm1
movsd  -0x48(%rbp),%xmm2
cvtsi2sdl -0x24(%rbp),%xmm3
divsd  %xmm3,%xmm2
movsd  -0x48(%rbp),%xmm3
cvtsi2sdl -0x24(%rbp),%xmm4
divsd  %xmm4,%xmm3
mulsd  %xmm3,%xmm2
mov    -0x24(%rbp),%eax
sub    $0x1,%eax
cvtsi2sd %eax,%xmm3
divsd  %xmm3,%xmm2
addsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  0xbce(%rip),%xmm0
mulsd  -0x50(%rbp),%xmm0
movsd  -0x50(%rbp),%xmm2
movsd  -0x50(%rbp),%xmm3
mov    -0x30(%rbp),%rax
movsd  (%rax),%xmm1
mov    -0x30(%rbp),%rax
movsd  (%rax),%xmm4
mulsd  %xmm4,%xmm1
addsd  %xmm3,%xmm1
divsd  %xmm1,%xmm2
movsd  0xb9b(%rip),%xmm1
call   1040 <betai@plt>
mov    -0x38(%rbp),%rax
movsd  %xmm0,(%rax)
movl   $0x1,-0x4(%rbp)
mov    -0x4(%rbp),%eax
add    $0x80,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<ttest>:
push   %r15
push   %r14
push   %rbx
sub    $0x30,%rsp
cmp    $0x2,%esi
jge    1153 <ttest+0x23>
mov    0x2e9b(%rip),%rax
mov    (%rax),%edi
lea    0xeba(%rip),%rsi
jmp    12bf <ttest+0x18f>
mov    %r9,%r14
mov    %r8,%r15
mov    %esi,%r9d
lea    -0x1(%r9),%r8
mov    %r9d,%r10d
and    $0x7,%r10d
cmp    $0x7,%r8
jae    1175 <ttest+0x45>
xorpd  %xmm6,%xmm6
xor    %ebx,%ebx
jmp    11c8 <ttest+0x98>
mov    %r9d,%eax
and    $0xfffffff8,%eax
xorpd  %xmm6,%xmm6
xor    %ebx,%ebx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
addsd  (%rdi,%rbx,8),%xmm6
addsd  0x8(%rdi,%rbx,8),%xmm6
addsd  0x10(%rdi,%rbx,8),%xmm6
addsd  0x18(%rdi,%rbx,8),%xmm6
addsd  0x20(%rdi,%rbx,8),%xmm6
addsd  0x28(%rdi,%rbx,8),%xmm6
addsd  0x30(%rdi,%rbx,8),%xmm6
addsd  0x38(%rdi,%rbx,8),%xmm6
add    $0x8,%rbx
cmp    %rbx,%rax
jne    1190 <ttest+0x60>
test   %r10,%r10
je     11ee <ttest+0xbe>
lea    (%rdi,%rbx,8),%rax
xor    %ebx,%ebx
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
addsd  (%rax,%rbx,8),%xmm6
add    $0x1,%rbx
cmp    %rbx,%r10
jne    11e0 <ttest+0xb0>
cvtsi2sd %esi,%xmm0
divsd  %xmm0,%xmm6
test   %esi,%esi
jle    120e <ttest+0xde>
mov    %r9d,%eax
and    $0x3,%eax
cmp    $0x3,%r8
jae    1217 <ttest+0xe7>
xorpd  %xmm7,%xmm7
xor    %ebx,%ebx
jmp    1280 <ttest+0x150>
xorpd  %xmm7,%xmm7
jmp    12aa <ttest+0x17a>
and    $0xfffffffc,%r9d
xorpd  %xmm7,%xmm7
xor    %ebx,%ebx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movsd  (%rdi,%rbx,8),%xmm1
movsd  0x8(%rdi,%rbx,8),%xmm2
subsd  %xmm6,%xmm1
mulsd  %xmm1,%xmm1
addsd  %xmm7,%xmm1
subsd  %xmm6,%xmm2
mulsd  %xmm2,%xmm2
addsd  %xmm1,%xmm2
movsd  0x10(%rdi,%rbx,8),%xmm1
subsd  %xmm6,%xmm1
mulsd  %xmm1,%xmm1
addsd  %xmm2,%xmm1
movsd  0x18(%rdi,%rbx,8),%xmm7
subsd  %xmm6,%xmm7
mulsd  %xmm7,%xmm7
addsd  %xmm1,%xmm7
add    $0x4,%rbx
cmp    %rbx,%r9
jne    1230 <ttest+0x100>
test   %rax,%rax
je     12aa <ttest+0x17a>
lea    (%rdi,%rbx,8),%rdi
xor    %ebx,%ebx
nopl   0x0(%rax,%rax,1)
movsd  (%rdi,%rbx,8),%xmm1
subsd  %xmm6,%xmm1
mulsd  %xmm1,%xmm1
addsd  %xmm1,%xmm7
add    $0x1,%rbx
cmp    %rbx,%rax
jne    1290 <ttest+0x160>
cmp    $0x2,%ecx
jge    12cb <ttest+0x19b>
mov    0x2d2a(%rip),%rax
mov    (%rax),%edi
lea    0xd5f(%rip),%rsi
call   1030 <fprintf@plt>
xor    %eax,%eax
jmp    14ee <ttest+0x3be>
add    $0xffffffff,%esi
mov    %ecx,%r9d
lea    -0x1(%r9),%r8
mov    %r9d,%eax
and    $0x7,%eax
cmp    $0x7,%r8
jae    12e9 <ttest+0x1b9>
xorpd  %xmm1,%xmm1
xor    %edi,%edi
jmp    1338 <ttest+0x208>
mov    %r9d,%ebx
and    $0xfffffff8,%ebx
xorpd  %xmm1,%xmm1
xor    %edi,%edi
cs nopw 0x0(%rax,%rax,1)
nop
addsd  (%rdx,%rdi,8),%xmm1
addsd  0x8(%rdx,%rdi,8),%xmm1
addsd  0x10(%rdx,%rdi,8),%xmm1
addsd  0x18(%rdx,%rdi,8),%xmm1
addsd  0x20(%rdx,%rdi,8),%xmm1
addsd  0x28(%rdx,%rdi,8),%xmm1
addsd  0x30(%rdx,%rdi,8),%xmm1
addsd  0x38(%rdx,%rdi,8),%xmm1
add    $0x8,%rdi
cmp    %rdi,%rbx
jne    1300 <ttest+0x1d0>
cvtsi2sd %esi,%xmm9
test   %rax,%rax
je     135e <ttest+0x22e>
lea    (%rdx,%rdi,8),%rsi
xor    %edi,%edi
nopl   0x0(%rax,%rax,1)
addsd  (%rsi,%rdi,8),%xmm1
add    $0x1,%rdi
cmp    %rdi,%rax
jne    1350 <ttest+0x220>
divsd  %xmm9,%xmm7
xorps  %xmm2,%xmm2
cvtsi2sd %ecx,%xmm2
divsd  %xmm2,%xmm1
xorpd  %xmm8,%xmm8
xorpd  %xmm3,%xmm3
test   %ecx,%ecx
jle    141a <ttest+0x2ea>
mov    %r9d,%esi
and    $0x3,%esi
cmp    $0x3,%r8
jae    1393 <ttest+0x263>
xorpd  %xmm3,%xmm3
xor    %eax,%eax
jmp    13f0 <ttest+0x2c0>
and    $0xfffffffc,%r9d
xorpd  %xmm3,%xmm3
xor    %eax,%eax
nopl   (%rax)
movsd  (%rdx,%rax,8),%xmm4
movsd  0x8(%rdx,%rax,8),%xmm5
subsd  %xmm1,%xmm4
mulsd  %xmm4,%xmm4
addsd  %xmm3,%xmm4
subsd  %xmm1,%xmm5
mulsd  %xmm5,%xmm5
addsd  %xmm4,%xmm5
movsd  0x10(%rdx,%rax,8),%xmm4
subsd  %xmm1,%xmm4
mulsd  %xmm4,%xmm4
addsd  %xmm5,%xmm4
movsd  0x18(%rdx,%rax,8),%xmm3
subsd  %xmm1,%xmm3
mulsd  %xmm3,%xmm3
addsd  %xmm4,%xmm3
add    $0x4,%rax
cmp    %rax,%r9
jne    13a0 <ttest+0x270>
test   %rsi,%rsi
je     141a <ttest+0x2ea>
lea    (%rdx,%rax,8),%rax
xor    %edx,%edx
nopl   0x0(%rax,%rax,1)
movsd  (%rax,%rdx,8),%xmm4
subsd  %xmm1,%xmm4
mulsd  %xmm4,%xmm4
addsd  %xmm4,%xmm3
add    $0x1,%rdx
cmp    %rdx,%rsi
jne    1400 <ttest+0x2d0>
add    $0xffffffff,%ecx
xorps  %xmm4,%xmm4
cvtsi2sd %ecx,%xmm4
divsd  %xmm4,%xmm3
subsd  %xmm1,%xmm6
divsd  %xmm0,%xmm7
divsd  %xmm2,%xmm3
movapd %xmm7,%xmm2
addsd  %xmm3,%xmm2
ucomisd %xmm8,%xmm2
jb     144c <ttest+0x31c>
xorps  %xmm0,%xmm0
sqrtsd %xmm2,%xmm0
jmp    149d <ttest+0x36d>
movapd %xmm2,%xmm0
movsd  %xmm6,0x28(%rsp)
movsd  %xmm7,0x20(%rsp)
movsd  %xmm9,0x18(%rsp)
movsd  %xmm3,0x10(%rsp)
movsd  %xmm2,0x8(%rsp)
movsd  %xmm4,(%rsp)
call   1050 <sqrt@plt>
movsd  (%rsp),%xmm4
movsd  0x8(%rsp),%xmm2
movsd  0x10(%rsp),%xmm3
movsd  0x18(%rsp),%xmm9
movsd  0x20(%rsp),%xmm7
movsd  0x28(%rsp),%xmm6
divsd  %xmm0,%xmm6
movsd  %xmm6,(%r15)
mulsd  %xmm2,%xmm2
mulsd  %xmm7,%xmm7
divsd  %xmm9,%xmm7
mulsd  %xmm3,%xmm3
divsd  %xmm4,%xmm3
addsd  %xmm7,%xmm3
divsd  %xmm3,%xmm2
movsd  0xb35(%rip),%xmm1
movapd %xmm2,%xmm0
mulsd  %xmm1,%xmm0
mulsd  %xmm6,%xmm6
addsd  %xmm2,%xmm6
divsd  %xmm6,%xmm2
call   1040 <betai@plt>
movsd  %xmm0,(%r14)
mov    $0x1,%eax
add    $0x30,%rsp
pop    %rbx
pop    %r14
pop    %r15
ret

LLVM IR (-O0):
@stderr = global i32 0
@.str = [22 x i8] c"ttest: n1 must be > 1\00"
@.str.1 = [22 x i8] c"ttest: n2 must be > 1\00"
define i32 @ttest(double* %0, i32 %1, double* %2, i32 %3, double* %4, double* %5) {
%7 = alloca i32
%8 = alloca double*
%9 = alloca i32
%10 = alloca double*
%11 = alloca i32
%12 = alloca double*
%13 = alloca double*
%14 = alloca double
%15 = alloca double
%16 = alloca double
%17 = alloca double
%18 = alloca double
%19 = alloca i32
%20 = alloca double
store double* %0, double** %8
store i32 %1, i32* %9
store double* %2, double** %10
store i32 %3, i32* %11
store double* %4, double** %12
store double* %5, double** %13
%21 = load i32, i32* %9
%22 = icmp sle i32 %21, 1
br i1 %22, label %23, label %26
23: ; preds = %6
%24 = load i32, i32* @stderr
%25 = call i32 @fprintf(i32 %24, i8* getelementptr inbounds ([22 x i8], [22 x i8]* @.str, i64 0, i64 0))
store i32 0, i32* %7
br label %199
26: ; preds = %6
store double 0.000000e+00, double* %17
store i32 0, i32* %19
br label %27
27: ; preds = %39, %26
%28 = load i32, i32* %19
%29 = load i32, i32* %9
%30 = icmp slt i32 %28, %29
br i1 %30, label %31, label %42
31: ; preds = %27
%32 = load double*, double** %8
%33 = load i32, i32* %19
%34 = sext i32 %33 to i64
%35 = getelementptr inbounds double, double* %32, i64 %34
%36 = load double, double* %35
%37 = load double, double* %17
%38 = fadd double %37, %36
store double %38, double* %17
br label %39
39: ; preds = %31
%40 = load i32, i32* %19
%41 = add nsw i32 %40, 1
store i32 %41, i32* %19
br label %27
42: ; preds = %27
%43 = load i32, i32* %9
%44 = sitofp i32 %43 to double
%45 = load double, double* %17
%46 = fdiv double %45, %44
store double %46, double* %17
store double 0.000000e+00, double* %14
store i32 0, i32* %19
br label %47
47: ; preds = %63, %42
%48 = load i32, i32* %19
%49 = load i32, i32* %9
%50 = icmp slt i32 %48, %49
br i1 %50, label %51, label %66
51: ; preds = %47
%52 = load double*, double** %8
%53 = load i32, i32* %19
%54 = sext i32 %53 to i64
%55 = getelementptr inbounds double, double* %52, i64 %54
%56 = load double, double* %55
%57 = load double, double* %17
%58 = fsub double %56, %57
store double %58, double* %20
%59 = load double, double* %20
%60 = load double, double* %20
%61 = load double, double* %14
%62 = call double @llvm.fmuladd.f64(double %59, double %60, double %61)
store double %62, double* %14
br label %63
63: ; preds = %51
%64 = load i32, i32* %19
%65 = add nsw i32 %64, 1
store i32 %65, i32* %19
br label %47
66: ; preds = %47
%67 = load i32, i32* %9
%68 = sub nsw i32 %67, 1
%69 = sitofp i32 %68 to double
%70 = load double, double* %14
%71 = fdiv double %70, %69
store double %71, double* %14
br label %72
72: ; preds = %66
%73 = load i32, i32* %11
%74 = icmp sle i32 %73, 1
br i1 %74, label %75, label %78
75: ; preds = %72
%76 = load i32, i32* @stderr
%77 = call i32 @fprintf(i32 %76, i8* getelementptr inbounds ([22 x i8], [22 x i8]* @.str.1, i64 0, i64 0))
store i32 0, i32* %7
br label %199
78: ; preds = %72
store double 0.000000e+00, double* %18
store i32 0, i32* %19
br label %79
79: ; preds = %91, %78
%80 = load i32, i32* %19
%81 = load i32, i32* %11
%82 = icmp slt i32 %80, %81
br i1 %82, label %83, label %94
83: ; preds = %79
%84 = load double*, double** %10
%85 = load i32, i32* %19
%86 = sext i32 %85 to i64
%87 = getelementptr inbounds double, double* %84, i64 %86
%88 = load double, double* %87
%89 = load double, double* %18
%90 = fadd double %89, %88
store double %90, double* %18
br label %91
91: ; preds = %83
%92 = load i32, i32* %19
%93 = add nsw i32 %92, 1
store i32 %93, i32* %19
br label %79
94: ; preds = %79
%95 = load i32, i32* %11
%96 = sitofp i32 %95 to double
%97 = load double, double* %18
%98 = fdiv double %97, %96
store double %98, double* %18
store double 0.000000e+00, double* %15
store i32 0, i32* %19
br label %99
99: ; preds = %115, %94
%100 = load i32, i32* %19
%101 = load i32, i32* %11
%102 = icmp slt i32 %100, %101
br i1 %102, label %103, label %118
103: ; preds = %99
%104 = load double*, double** %10
%105 = load i32, i32* %19
%106 = sext i32 %105 to i64
%107 = getelementptr inbounds double, double* %104, i64 %106
%108 = load double, double* %107
%109 = load double, double* %18
%110 = fsub double %108, %109
store double %110, double* %20
%111 = load double, double* %20
%112 = load double, double* %20
%113 = load double, double* %15
%114 = call double @llvm.fmuladd.f64(double %111, double %112, double %113)
store double %114, double* %15
br label %115
115: ; preds = %103
%116 = load i32, i32* %19
%117 = add nsw i32 %116, 1
store i32 %117, i32* %19
br label %99
118: ; preds = %99
%119 = load i32, i32* %11
%120 = sub nsw i32 %119, 1
%121 = sitofp i32 %120 to double
%122 = load double, double* %15
%123 = fdiv double %122, %121
store double %123, double* %15
br label %124
124: ; preds = %118
%125 = load double, double* %17
%126 = load double, double* %18
%127 = fsub double %125, %126
%128 = load double, double* %14
%129 = load i32, i32* %9
%130 = sitofp i32 %129 to double
%131 = fdiv double %128, %130
%132 = load double, double* %15
%133 = load i32, i32* %11
%134 = sitofp i32 %133 to double
%135 = fdiv double %132, %134
%136 = fadd double %131, %135
%137 = call double @sqrt(double %136)
%138 = fdiv double %127, %137
%139 = load double*, double** %12
store double %138, double* %139
%140 = load double, double* %14
%141 = load i32, i32* %9
%142 = sitofp i32 %141 to double
%143 = fdiv double %140, %142
%144 = load double, double* %15
%145 = load i32, i32* %11
%146 = sitofp i32 %145 to double
%147 = fdiv double %144, %146
%148 = fadd double %143, %147
%149 = load double, double* %14
%150 = load i32, i32* %9
%151 = sitofp i32 %150 to double
%152 = fdiv double %149, %151
%153 = load double, double* %15
%154 = load i32, i32* %11
%155 = sitofp i32 %154 to double
%156 = fdiv double %153, %155
%157 = fadd double %152, %156
%158 = fmul double %148, %157
%159 = load double, double* %14
%160 = load i32, i32* %9
%161 = sitofp i32 %160 to double
%162 = fdiv double %159, %161
%163 = load double, double* %14
%164 = load i32, i32* %9
%165 = sitofp i32 %164 to double
%166 = fdiv double %163, %165
%167 = fmul double %162, %166
%168 = load i32, i32* %9
%169 = sub nsw i32 %168, 1
%170 = sitofp i32 %169 to double
%171 = fdiv double %167, %170
%172 = load double, double* %15
%173 = load i32, i32* %11
%174 = sitofp i32 %173 to double
%175 = fdiv double %172, %174
%176 = load double, double* %15
%177 = load i32, i32* %11
%178 = sitofp i32 %177 to double
%179 = fdiv double %176, %178
%180 = fmul double %175, %179
%181 = load i32, i32* %11
%182 = sub nsw i32 %181, 1
%183 = sitofp i32 %182 to double
%184 = fdiv double %180, %183
%185 = fadd double %171, %184
%186 = fdiv double %158, %185
store double %186, double* %16
%187 = load double, double* %16
%188 = fmul double 5.000000e-01, %187
%189 = load double, double* %16
%190 = load double, double* %16
%191 = load double*, double** %12
%192 = load double, double* %191
%193 = load double*, double** %12
%194 = load double, double* %193
%195 = call double @llvm.fmuladd.f64(double %192, double %194, double %190)
%196 = fdiv double %189, %195
%197 = call double @betai(double %188, double 5.000000e-01, double %196)
%198 = load double*, double** %13
store double %197, double* %198
store i32 1, i32* %7
br label %199
199: ; preds = %124, %75, %23
%200 = load i32, i32* %7
ret i32 %200
}

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