#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*/ free_vector(float *, int, int);
float sqrt(float);
float *vector(int, int);
void cntab1(int **nn, int ni, int nj, float *chisq, float *df, float *prob,
            float *cramrv, float *ccc) {
  float gammq(float a, float x);
  int nnj, nni, j, i, minij;
  float sum = 0.0, expctd, *sumi, *sumj, temp;
  sumi = vector(1, ni);
  sumj = vector(1, nj);
  nni = ni;
  nnj = nj;
  for (i = 1; i <= ni; i++) {
    sumi[i] = 0.0;
    for (j = 1; j <= nj; j++) {
      sumi[i] += nn[i][j];
      sum += nn[i][j];
    }
    if (sumi[i] == 0.0)
      --nni;
  }
  for (j = 1; j <= nj; j++) {
    sumj[j] = 0.0;
    for (i = 1; i <= ni; i++)
      sumj[j] += nn[i][j];
    if (sumj[j] == 0.0)
      --nnj;
  }
  *df = nni * nnj - nni - nnj + 1;
  *chisq = 0.0;
  for (i = 1; i <= ni; i++) {
    for (j = 1; j <= nj; j++) {
      expctd = sumj[j] * sumi[i] / sum;
      temp = nn[i][j] - expctd;
      *chisq += temp * temp / (expctd + 1.0e-30);
    }
  }
  *prob = gammq(0.5 * (*df), 0.5 * (*chisq));
  minij = nni < nnj ? nni - 1 : nnj - 1;
  *cramrv = sqrt(*chisq / (sum * minij));
  *ccc = sqrt(*chisq / (*chisq + sum));
  free_vector(sumj, 1, nj);
  free_vector(sumi, 1, ni);
}


x86 assembly (clang, -O0):
<cntab1>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%rax
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
mov    %rcx,-0x18(%rbp)
mov    %r8,-0x20(%rbp)
mov    %r9,-0x28(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x40(%rbp)
mov    -0xc(%rbp),%esi
mov    $0x1,%edi
call   1050 <vector@plt>
mov    %rax,-0x50(%rbp)
mov    -0x10(%rbp),%esi
mov    $0x1,%edi
call   1050 <vector@plt>
mov    %rax,-0x58(%rbp)
mov    -0xc(%rbp),%eax
mov    %eax,-0x30(%rbp)
mov    -0x10(%rbp),%eax
mov    %eax,-0x2c(%rbp)
movl   $0x1,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0xc(%rbp),%eax
jg     1266 <cntab1+0x126>
mov    -0x50(%rbp),%rax
movslq -0x38(%rbp),%rcx
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax,%rcx,4)
movl   $0x1,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x10(%rbp),%eax
jg     1226 <cntab1+0xe6>
mov    -0x8(%rbp),%rax
movslq -0x38(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x34(%rbp),%rcx
cvtsi2ssl (%rax,%rcx,4),%xmm0
mov    -0x50(%rbp),%rax
movslq -0x38(%rbp),%rcx
addss  (%rax,%rcx,4),%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x8(%rbp),%rax
movslq -0x38(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x34(%rbp),%rcx
cvtsi2ssl (%rax,%rcx,4),%xmm0
addss  -0x40(%rbp),%xmm0
movss  %xmm0,-0x40(%rbp)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    11c6 <cntab1+0x86>
mov    -0x50(%rbp),%rax
movslq -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1253 <cntab1+0x113>
jp     1253 <cntab1+0x113>
mov    -0x30(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x30(%rbp)
jmp    1258 <cntab1+0x118>
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    11a3 <cntab1+0x63>
movl   $0x1,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x10(%rbp),%eax
jg     1311 <cntab1+0x1d1>
mov    -0x58(%rbp),%rax
movslq -0x34(%rbp),%rcx
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax,%rcx,4)
movl   $0x1,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0xc(%rbp),%eax
jg     12d1 <cntab1+0x191>
mov    -0x8(%rbp),%rax
movslq -0x38(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x34(%rbp),%rcx
cvtsi2ssl (%rax,%rcx,4),%xmm0
mov    -0x58(%rbp),%rax
movslq -0x34(%rbp),%rcx
addss  (%rax,%rcx,4),%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    1290 <cntab1+0x150>
mov    -0x58(%rbp),%rax
movslq -0x34(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    12fe <cntab1+0x1be>
jp     12fe <cntab1+0x1be>
mov    -0x2c(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x2c(%rbp)
jmp    1303 <cntab1+0x1c3>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    126d <cntab1+0x12d>
mov    -0x30(%rbp),%eax
imul   -0x2c(%rbp),%eax
sub    -0x30(%rbp),%eax
sub    -0x2c(%rbp),%eax
add    $0x1,%eax
cvtsi2ss %eax,%xmm0
mov    -0x20(%rbp),%rax
movss  %xmm0,(%rax)
mov    -0x18(%rbp),%rax
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax)
movl   $0x1,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0xc(%rbp),%eax
jg     1401 <cntab1+0x2c1>
movl   $0x1,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x10(%rbp),%eax
jg     13ee <cntab1+0x2ae>
mov    -0x58(%rbp),%rax
movslq -0x34(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x50(%rbp),%rax
movslq -0x38(%rbp),%rcx
mulss  (%rax,%rcx,4),%xmm0
divss  -0x40(%rbp),%xmm0
movss  %xmm0,-0x44(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x38(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x34(%rbp),%rcx
cvtsi2ssl (%rax,%rcx,4),%xmm0
subss  -0x44(%rbp),%xmm0
movss  %xmm0,-0x5c(%rbp)
movss  -0x5c(%rbp),%xmm0
mulss  -0x5c(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
movss  -0x44(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xc48(%rip),%xmm2
addsd  %xmm2,%xmm0
divsd  %xmm0,%xmm1
mov    -0x18(%rbp),%rax
movss  (%rax),%xmm0
cvtss2sd %xmm0,%xmm0
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,(%rax)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1352 <cntab1+0x212>
jmp    13f3 <cntab1+0x2b3>
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    133f <cntab1+0x1ff>
mov    -0x20(%rbp),%rax
movss  (%rax),%xmm0
cvtss2sd %xmm0,%xmm1
movsd  0xbeb(%rip),%xmm0
mulsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
mov    -0x18(%rbp),%rax
movss  (%rax),%xmm1
cvtss2sd %xmm1,%xmm2
movsd  0xbcf(%rip),%xmm1
mulsd  %xmm2,%xmm1
cvtsd2ss %xmm1,%xmm1
call   1030 <gammq@plt>
mov    -0x28(%rbp),%rax
movss  %xmm0,(%rax)
mov    -0x30(%rbp),%eax
cmp    -0x2c(%rbp),%eax
jge    1460 <cntab1+0x320>
mov    -0x30(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x60(%rbp)
jmp    1469 <cntab1+0x329>
mov    -0x2c(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x60(%rbp)
mov    -0x60(%rbp),%eax
mov    %eax,-0x3c(%rbp)
mov    -0x18(%rbp),%rax
movss  (%rax),%xmm0
movss  -0x40(%rbp),%xmm1
cvtsi2ssl -0x3c(%rbp),%xmm2
mulss  %xmm2,%xmm1
divss  %xmm1,%xmm0
call   1060 <sqrt@plt>
mov    0x10(%rbp),%rax
movss  %xmm0,(%rax)
mov    -0x18(%rbp),%rax
movss  (%rax),%xmm0
mov    -0x18(%rbp),%rax
movss  (%rax),%xmm1
addss  -0x40(%rbp),%xmm1
divss  %xmm1,%xmm0
call   1060 <sqrt@plt>
mov    0x18(%rbp),%rax
movss  %xmm0,(%rax)
mov    -0x58(%rbp),%rdi
mov    -0x10(%rbp),%edx
mov    $0x1,%esi
call   1040 <free_vector@plt>
mov    -0x50(%rbp),%rdi
mov    -0xc(%rbp),%edx
mov    $0x1,%esi
call   1040 <free_vector@plt>
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<cntab1>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x38,%rsp
mov    %r9,0x30(%rsp)
mov    %r8,0x8(%rsp)
mov    %rcx,%r13
mov    %edx,%r15d
mov    %esi,%ebx
mov    %rdi,%rbp
mov    0x78(%rsp),%rax
mov    %rax,0x18(%rsp)
mov    0x70(%rsp),%rax
mov    %rax,0x10(%rsp)
mov    $0x1,%edi
call   1050 <vector@plt>
mov    %rax,%r12
mov    $0x1,%edi
mov    %r15d,%esi
call   1050 <vector@plt>
mov    %rbx,%r10
mov    %rax,%r14
test   %r10d,%r10d
jle    12b4 <cntab1+0x174>
lea    0x1(%r15),%r8d
lea    0x1(%r10),%r9d
lea    -0x1(%r8),%rdx
add    $0xfffffffffffffffe,%r8
mov    %edx,%esi
and    $0x3,%esi
and    $0xfffffffffffffffc,%rdx
xorps  %xmm0,%xmm0
mov    $0x1,%edi
mov    %r10d,%eax
xorps  %xmm4,%xmm4
mov    %r15,%r11
mov    %rax,%rbx
jmp    11eb <cntab1+0xab>
nop
xorps  %xmm1,%xmm1
cmpeqss %xmm0,%xmm1
movd   %xmm1,%eax
add    %eax,%ebx
add    $0x1,%rdi
cmp    %r9,%rdi
je     12c0 <cntab1+0x180>
movl   $0x0,(%r12,%rdi,4)
test   %r11d,%r11d
jle    11d0 <cntab1+0x90>
mov    0x0(%rbp,%rdi,8),%rax
xorps  %xmm1,%xmm1
cmp    $0x3,%r8
jae    1220 <cntab1+0xe0>
mov    $0x1,%ecx
test   %rsi,%rsi
jne    1286 <cntab1+0x146>
jmp    12a9 <cntab1+0x169>
cs nopw 0x0(%rax,%rax,1)
nop
xor    %ecx,%ecx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
xorps  %xmm2,%xmm2
cvtsi2ssl 0x4(%rax,%rcx,4),%xmm2
xorps  %xmm3,%xmm3
cvtsi2ssl 0x8(%rax,%rcx,4),%xmm3
addss  %xmm2,%xmm1
addss  %xmm2,%xmm4
addss  %xmm3,%xmm1
addss  %xmm3,%xmm4
xorps  %xmm2,%xmm2
cvtsi2ssl 0xc(%rax,%rcx,4),%xmm2
addss  %xmm2,%xmm1
addss  %xmm2,%xmm4
xorps  %xmm2,%xmm2
cvtsi2ssl 0x10(%rax,%rcx,4),%xmm2
addss  %xmm2,%xmm1
addss  %xmm2,%xmm4
add    $0x4,%rcx
cmp    %rcx,%rdx
jne    1230 <cntab1+0xf0>
add    $0x1,%rcx
test   %rsi,%rsi
je     12a9 <cntab1+0x169>
lea    (%rax,%rcx,4),%rax
xor    %ecx,%ecx
nopl   0x0(%rax)
xorps  %xmm2,%xmm2
cvtsi2ssl (%rax,%rcx,4),%xmm2
addss  %xmm2,%xmm1
addss  %xmm2,%xmm4
add    $0x1,%rcx
cmp    %rcx,%rsi
jne    1290 <cntab1+0x150>
movss  %xmm1,(%r12,%rdi,4)
jmp    11d3 <cntab1+0x93>
xorps  %xmm4,%xmm4
mov    %r10d,%eax
mov    %r15,%r11
mov    %rax,%rbx
mov    %rbx,0x20(%rsp)
mov    %r11d,%r15d
test   %r11d,%r11d
jle    13ee <cntab1+0x2ae>
lea    0x1(%r10),%r8d
lea    0x1(%r11),%r9d
lea    -0x1(%r8),%rdx
add    $0xfffffffffffffffe,%r8
mov    %edx,%esi
and    $0x3,%esi
and    $0xfffffffffffffffc,%rdx
mov    $0x1,%edi
xorps  %xmm0,%xmm0
mov    %r11d,%r15d
jmp    131f <cntab1+0x1df>
nopw   0x0(%rax,%rax,1)
movss  %xmm1,(%r14,%rdi,4)
cmpeqss %xmm0,%xmm1
movd   %xmm1,%eax
add    %eax,%r15d
add    $0x1,%rdi
cmp    %r9,%rdi
je     13ee <cntab1+0x2ae>
movl   $0x0,(%r14,%rdi,4)
xorps  %xmm1,%xmm1
test   %r10d,%r10d
jle    1306 <cntab1+0x1c6>
xorps  %xmm1,%xmm1
cmp    $0x3,%r8
jae    1350 <cntab1+0x210>
mov    $0x1,%eax
test   %rsi,%rsi
jne    13ba <cntab1+0x27a>
jmp    1300 <cntab1+0x1c0>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
xor    %eax,%eax
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
mov    0x8(%rbp,%rax,8),%rcx
mov    0x10(%rbp,%rax,8),%rbx
xorps  %xmm2,%xmm2
cvtsi2ssl (%rcx,%rdi,4),%xmm2
addss  %xmm1,%xmm2
xorps  %xmm3,%xmm3
cvtsi2ssl (%rbx,%rdi,4),%xmm3
addss  %xmm2,%xmm3
mov    0x18(%rbp,%rax,8),%rcx
xorps  %xmm2,%xmm2
cvtsi2ssl (%rcx,%rdi,4),%xmm2
mov    0x20(%rbp,%rax,8),%rcx
xorps  %xmm1,%xmm1
cvtsi2ssl (%rcx,%rdi,4),%xmm1
addss  %xmm3,%xmm2
addss  %xmm2,%xmm1
add    $0x4,%rax
cmp    %rax,%rdx
jne    1360 <cntab1+0x220>
add    $0x1,%rax
test   %rsi,%rsi
je     1300 <cntab1+0x1c0>
lea    0x0(,%rax,8),%rax
add    %rbp,%rax
xor    %ecx,%ecx
nopw   0x0(%rax,%rax,1)
mov    (%rax,%rcx,8),%rbx
xorps  %xmm2,%xmm2
cvtsi2ssl (%rbx,%rdi,4),%xmm2
addss  %xmm2,%xmm1
add    $0x1,%rcx
cmp    %rcx,%rsi
jne    13d0 <cntab1+0x290>
jmp    1300 <cntab1+0x1c0>
mov    %r15d,%eax
mov    0x20(%rsp),%rcx
imul   %ecx,%eax
add    %r15d,%ecx
sub    %ecx,%eax
add    $0x1,%eax
xorps  %xmm0,%xmm0
cvtsi2ss %eax,%xmm0
mov    0x8(%rsp),%rbx
movss  %xmm0,(%rbx)
movl   $0x0,0x0(%r13)
xorps  %xmm1,%xmm1
test   %r10d,%r10d
jle    14ad <cntab1+0x36d>
lea    0x1(%r11),%eax
lea    0x1(%r10),%ecx
mov    $0x1,%edx
xorps  %xmm1,%xmm1
movsd  0xbc3(%rip),%xmm0
jmp    1449 <cntab1+0x309>
nop
add    $0x1,%rdx
cmp    %rcx,%rdx
je     14ad <cntab1+0x36d>
test   %r11d,%r11d
jle    1440 <cntab1+0x300>
mov    0x0(%rbp,%rdx,8),%rsi
mov    $0x1,%edi
nopl   0x0(%rax,%rax,1)
movss  (%r14,%rdi,4),%xmm2
mulss  (%r12,%rdx,4),%xmm2
xorps  %xmm3,%xmm3
cvtsi2ssl (%rsi,%rdi,4),%xmm3
divss  %xmm4,%xmm2
subss  %xmm2,%xmm3
mulss  %xmm3,%xmm3
cvtss2sd %xmm3,%xmm3
cvtss2sd %xmm2,%xmm2
addsd  %xmm0,%xmm2
divsd  %xmm2,%xmm3
cvtss2sd %xmm1,%xmm1
addsd  %xmm3,%xmm1
cvtsd2ss %xmm1,%xmm1
movss  %xmm1,0x0(%r13)
add    $0x1,%rdi
cmp    %rdi,%rax
jne    1460 <cntab1+0x320>
jmp    1440 <cntab1+0x300>
mov    %r11,%rbp
mov    %r10,0x28(%rsp)
movss  0xb4b(%rip),%xmm2
movss  (%rbx),%xmm0
mulss  %xmm2,%xmm0
mulss  %xmm2,%xmm1
movss  %xmm4,0x8(%rsp)
call   1030 <gammq@plt>
movss  0x8(%rsp),%xmm2
mov    0x30(%rsp),%rax
movss  %xmm0,(%rax)
mov    0x20(%rsp),%rax
cmp    %r15d,%eax
cmovl  %eax,%r15d
add    $0xffffffff,%r15d
movss  0x0(%r13),%xmm0
xorps  %xmm1,%xmm1
cvtsi2ss %r15d,%xmm1
mulss  %xmm2,%xmm1
divss  %xmm1,%xmm0
ucomiss 0xafc(%rip),%xmm0
jb     1518 <cntab1+0x3d8>
sqrtss %xmm0,%xmm0
jmp    1523 <cntab1+0x3e3>
call   1060 <sqrt@plt>
movss  0x8(%rsp),%xmm2
mov    0x10(%rsp),%rax
movss  %xmm0,(%rax)
movss  0x0(%r13),%xmm0
addss  %xmm0,%xmm2
divss  %xmm2,%xmm0
xorps  %xmm1,%xmm1
ucomiss %xmm1,%xmm0
jb     1548 <cntab1+0x408>
sqrtss %xmm0,%xmm0
jmp    154d <cntab1+0x40d>
call   1060 <sqrt@plt>
mov    %rbp,%rdx
mov    0x18(%rsp),%rax
movss  %xmm0,(%rax)
mov    %r14,%rdi
mov    $0x1,%esi
call   1040 <free_vector@plt>
mov    %r12,%rdi
mov    $0x1,%esi
mov    0x28(%rsp),%rdx
add    $0x38,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
jmp    1040 <free_vector@plt>

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

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