#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*/ ENOMEM;
int /*<<< orphan*/ assert(int /*<<< orphan*/);
double *calloc(size_t, int);
int /*<<< orphan*/ error(int /*<<< orphan*/, int /*<<< orphan*/, char *);
int /*<<< orphan*/ isgreater(double, double);
double **malloc(size_t);
double **rmtable(double *freq, size_t geno) {
  double **table = malloc(geno * sizeof(double *));
  if (table == NULL)
    error(0, ENOMEM, "Null pointer\n");
  double denom = 0.0F;
  for (int i = 0; i < geno; i++) {
    table[i] = calloc(geno, sizeof(double));
    if (table[i] == NULL)
      error(0, ENOMEM, "Null pointer\n");
    for (int j = 0; j < geno; j++) {
      table[i][j] = freq[i] * freq[j];
      denom += table[i][j];
    }
  }
  assert(isgreater(denom, 0.0));
  for (int i = 0; i < geno; i++)
    for (int j = 0; j < geno; j++)
      table[i][j] /= denom;
  return table;
}


x86 assembly (clang, -O0):
<rmtable>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    -0x10(%rbp),%rdi
shl    $0x3,%rdi
call   1050 <malloc@plt>
mov    %rax,-0x18(%rbp)
cmpq   $0x0,-0x18(%rbp)
jne    1193 <rmtable+0x43>
mov    0x2e5d(%rip),%rax
mov    (%rax),%esi
xor    %edi,%edi
lea    0xe72(%rip),%rdx
call   1060 <error@plt>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x20(%rbp)
movl   $0x0,-0x24(%rbp)
movslq -0x24(%rbp),%rax
cmp    -0x10(%rbp),%rax
jae    127b <rmtable+0x12b>
mov    -0x10(%rbp),%rdi
mov    $0x8,%esi
call   1030 <calloc@plt>
mov    %rax,%rdx
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
mov    %rdx,(%rax,%rcx,8)
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
cmpq   $0x0,(%rax,%rcx,8)
jne    11f7 <rmtable+0xa7>
mov    0x2df9(%rip),%rax
mov    (%rax),%esi
xor    %edi,%edi
lea    0xe0e(%rip),%rdx
call   1060 <error@plt>
movl   $0x0,-0x28(%rbp)
movslq -0x28(%rbp),%rax
cmp    -0x10(%rbp),%rax
jae    1268 <rmtable+0x118>
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x28(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x28(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    11fe <rmtable+0xae>
jmp    126d <rmtable+0x11d>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    11a2 <rmtable+0x52>
movsd  -0x20(%rbp),%xmm0
xorps  %xmm1,%xmm1
call   1070 <isgreater@plt>
mov    %eax,%edi
call   1040 <assert@plt>
movl   $0x0,-0x2c(%rbp)
movslq -0x2c(%rbp),%rax
cmp    -0x10(%rbp),%rax
jae    12fd <rmtable+0x1ad>
movl   $0x0,-0x30(%rbp)
movslq -0x30(%rbp),%rax
cmp    -0x10(%rbp),%rax
jae    12ea <rmtable+0x19a>
movsd  -0x20(%rbp),%xmm1
mov    -0x18(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    12ab <rmtable+0x15b>
jmp    12ef <rmtable+0x19f>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1296 <rmtable+0x146>
mov    -0x18(%rbp),%rax
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<rmtable>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x18,%rsp
mov    %rsi,%r15
mov    %rdi,%rbx
lea    0x0(,%rsi,8),%rdi
call   1050 <malloc@plt>
mov    %rax,%r14
test   %rax,%rax
je     126b <rmtable+0x11b>
test   %r15,%r15
je     128b <rmtable+0x13b>
mov    %r15,%rbp
and    $0xfffffffffffffffe,%rbp
xorpd  %xmm0,%xmm0
xor    %r12d,%r12d
jmp    11ad <rmtable+0x5d>
cs nopw 0x0(%rax,%rax,1)
add    $0x1,%r12
cmp    %r15,%r12
je     128f <rmtable+0x13f>
movapd %xmm0,(%rsp)
mov    %r15,%rdi
mov    $0x8,%esi
call   1030 <calloc@plt>
mov    %rax,%r13
mov    %rax,(%r14,%r12,8)
test   %rax,%rax
jne    11e2 <rmtable+0x92>
mov    0x2e0e(%rip),%rax
mov    (%rax),%esi
xor    %edi,%edi
lea    0xe23(%rip),%rdx
call   1060 <error@plt>
xor    %eax,%eax
cmp    $0x1,%r15
jne    1200 <rmtable+0xb0>
movapd (%rsp),%xmm0
jmp    1246 <rmtable+0xf6>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movapd (%rsp),%xmm0
cs nopw 0x0(%rax,%rax,1)
nop
movsd  (%rbx,%r12,8),%xmm1
mulsd  (%rbx,%rax,8),%xmm1
movsd  %xmm1,0x0(%r13,%rax,8)
addsd  %xmm0,%xmm1
movsd  (%rbx,%r12,8),%xmm0
mulsd  0x8(%rbx,%rax,8),%xmm0
movsd  %xmm0,0x8(%r13,%rax,8)
addsd  %xmm1,%xmm0
add    $0x2,%rax
cmp    %rax,%rbp
jne    1210 <rmtable+0xc0>
test   $0x1,%r15b
je     11a0 <rmtable+0x50>
movsd  (%rbx,%r12,8),%xmm1
mulsd  (%rbx,%rax,8),%xmm1
movsd  %xmm1,0x0(%r13,%rax,8)
addsd  %xmm1,%xmm0
jmp    11a0 <rmtable+0x50>
mov    0x2d6e(%rip),%rax
mov    (%rax),%esi
lea    0xd85(%rip),%rdx
xor    %edi,%edi
call   1060 <error@plt>
test   %r15,%r15
jne    1186 <rmtable+0x36>
xorpd  %xmm0,%xmm0
movapd %xmm0,(%rsp)
xorpd  %xmm1,%xmm1
call   1070 <isgreater@plt>
mov    %eax,%edi
call   1040 <assert@plt>
movapd (%rsp),%xmm3
test   %r15,%r15
je     138e <rmtable+0x23e>
lea    -0x2(%r15),%rax
mov    %rax,%rcx
shr    %rcx
add    $0x1,%rcx
mov    %r15,%rdx
and    $0xfffffffffffffffe,%rdx
movapd %xmm3,%xmm0
unpcklpd %xmm3,%xmm0
mov    %rcx,%r8
and    $0xfffffffffffffffe,%r8
xor    %edi,%edi
jmp    12ed <rmtable+0x19d>
nopw   0x0(%rax,%rax,1)
add    $0x1,%rdi
cmp    %r15,%rdi
je     138e <rmtable+0x23e>
mov    (%r14,%rdi,8),%rbp
cmp    $0x2,%r15
jae    1300 <rmtable+0x1b0>
xor    %esi,%esi
jmp    1370 <rmtable+0x220>
nopl   0x0(%rax,%rax,1)
cmp    $0x2,%rax
jae    1310 <rmtable+0x1c0>
xor    %ebx,%ebx
jmp    134a <rmtable+0x1fa>
nopw   0x0(%rax,%rax,1)
mov    %r8,%rsi
xor    %ebx,%ebx
cs nopw 0x0(%rax,%rax,1)
nop
movupd 0x0(%rbp,%rbx,8),%xmm1
movupd 0x10(%rbp,%rbx,8),%xmm2
divpd  %xmm0,%xmm1
movupd %xmm1,0x0(%rbp,%rbx,8)
divpd  %xmm0,%xmm2
movupd %xmm2,0x10(%rbp,%rbx,8)
add    $0x4,%rbx
add    $0xfffffffffffffffe,%rsi
jne    1320 <rmtable+0x1d0>
test   $0x1,%cl
je     135f <rmtable+0x20f>
movupd 0x0(%rbp,%rbx,8),%xmm1
divpd  %xmm0,%xmm1
movupd %xmm1,0x0(%rbp,%rbx,8)
mov    %rdx,%rsi
cmp    %r15,%rdx
je     12e0 <rmtable+0x190>
nopl   0x0(%rax,%rax,1)
movsd  0x0(%rbp,%rsi,8),%xmm1
divsd  %xmm3,%xmm1
movsd  %xmm1,0x0(%rbp,%rsi,8)
add    $0x1,%rsi
cmp    %rsi,%r15
jne    1370 <rmtable+0x220>
jmp    12e0 <rmtable+0x190>
mov    %r14,%rax
add    $0x18,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
@ENOMEM = global i32 0
@.str = [14 x i8] c"Null pointer\0A\00"
define noalias double** @rmtable(double* nocapture %0, i64 %1) {
%3 = shl i64 %1, 3
%4 = tail call double** @malloc(i64 %3)
%5 = icmp eq double** %4, null
br i1 %5, label %6, label %9
6: ; preds = %2
%7 = load i32, i32* @ENOMEM
%8 = tail call i32 @error(i32 0, i32 %7, i8* getelementptr inbounds ([14 x i8], [14 x i8]* @.str, i64 0, i64 0))
br label %9
9: ; preds = %6, %2
%10 = icmp eq i64 %1, 0
br i1 %10, label %16, label %11
11: ; preds = %9
%12 = and i64 %1, 1
%13 = icmp eq i64 %1, 1
%14 = and i64 %1, -2
%15 = icmp eq i64 %12, 0
br label %33
16: ; preds = %55, %9
%17 = phi double [ 0.000000e+00, %9 ], [ %56, %55 ]
%18 = tail call i32 @isgreater(double %17, double 0.000000e+00)
%19 = tail call i32 @assert(i32 %18)
br i1 %10, label %112, label %20
20: ; preds = %16
%21 = add i64 %1, -2
%22 = lshr i64 %21, 1
%23 = add nuw i64 %22, 1
%24 = icmp ult i64 %1, 2
%25 = and i64 %1, -2
%26 = insertelement <2 x double> poison, double %17, i64 0
%27 = shufflevector <2 x double> %26, <2 x double> poison, <2 x i32> zeroinitializer
%28 = and i64 %23, 1
%29 = icmp ult i64 %21, 2
%30 = and i64 %23, -2
%31 = icmp eq i64 %28, 0
%32 = icmp eq i64 %25, %1
br label %79
33: ; preds = %11, %55
%34 = phi i64 [ %57, %55 ], [ 0, %11 ]
%35 = phi double [ %56, %55 ], [ 0.000000e+00, %11 ]
%36 = tail call double* @calloc(i64 %1, i32 8)
%37 = getelementptr inbounds double*, double** %4, i64 %34
store double* %36, double** %37
%38 = icmp eq double* %36, null
br i1 %38, label %39, label %42
39: ; preds = %33
%40 = load i32, i32* @ENOMEM
%41 = tail call i32 @error(i32 0, i32 %40, i8* getelementptr inbounds ([14 x i8], [14 x i8]* @.str, i64 0, i64 0))
br label %42
42: ; preds = %39, %33
%43 = getelementptr inbounds double, double* %0, i64 %34
br i1 %13, label %44, label %59
44: ; preds = %59, %42
%45 = phi double [ undef, %42 ], [ %75, %59 ]
%46 = phi i64 [ 0, %42 ], [ %76, %59 ]
%47 = phi double [ %35, %42 ], [ %75, %59 ]
br i1 %15, label %55, label %48
48: ; preds = %44
%49 = load double, double* %43
%50 = getelementptr inbounds double, double* %0, i64 %46
%51 = load double, double* %50
%52 = fmul double %49, %51
%53 = getelementptr inbounds double, double* %36, i64 %46
store double %52, double* %53
%54 = fadd double %47, %52
br label %55
55: ; preds = %44, %48
%56 = phi double [ %45, %44 ], [ %54, %48 ]
%57 = add nuw i64 %34, 1
%58 = icmp eq i64 %57, %1
br i1 %58, label %16, label %33
59: ; preds = %42, %59
%60 = phi i64 [ %76, %59 ], [ 0, %42 ]
%61 = phi double [ %75, %59 ], [ %35, %42 ]
%62 = phi i64 [ %77, %59 ], [ 0, %42 ]
%63 = load double, double* %43
%64 = getelementptr inbounds double, double* %0, i64 %60
%65 = load double, double* %64
%66 = fmul double %63, %65
%67 = getelementptr inbounds double, double* %36, i64 %60
store double %66, double* %67
%68 = fadd double %61, %66
%69 = or i64 %60, 1
%70 = load double, double* %43
%71 = getelementptr inbounds double, double* %0, i64 %69
%72 = load double, double* %71
%73 = fmul double %70, %72
%74 = getelementptr inbounds double, double* %36, i64 %69
store double %73, double* %74
%75 = fadd double %68, %73
%76 = add nuw nsw i64 %60, 2
%77 = add i64 %62, 2
%78 = icmp eq i64 %77, %14
br i1 %78, label %44, label %59
79: ; preds = %20, %113
%80 = phi i64 [ %114, %113 ], [ 0, %20 ]
%81 = getelementptr inbounds double*, double** %4, i64 %80
%82 = load double*, double** %81
br i1 %24, label %110, label %83
83: ; preds = %79
br i1 %29, label %101, label %84
84: ; preds = %83, %84
%85 = phi i64 [ %98, %84 ], [ 0, %83 ]
%86 = phi i64 [ %99, %84 ], [ 0, %83 ]
%87 = getelementptr inbounds double, double* %82, i64 %85
%88 = bitcast double* %87 to <2 x double>*
%89 = load <2 x double>, <2 x double>* %88
%90 = fdiv <2 x double> %89, %27
%91 = bitcast double* %87 to <2 x double>*
store <2 x double> %90, <2 x double>* %91
%92 = or i64 %85, 2
%93 = getelementptr inbounds double, double* %82, i64 %92
%94 = bitcast double* %93 to <2 x double>*
%95 = load <2 x double>, <2 x double>* %94
%96 = fdiv <2 x double> %95, %27
%97 = bitcast double* %93 to <2 x double>*
store <2 x double> %96, <2 x double>* %97
%98 = add nuw i64 %85, 4
%99 = add nuw i64 %86, 2
%100 = icmp eq i64 %99, %30
br i1 %100, label %101, label %84
101: ; preds = %84, %83
%102 = phi i64 [ 0, %83 ], [ %98, %84 ]
br i1 %31, label %109, label %103
103: ; preds = %101
%104 = getelementptr inbounds double, double* %82, i64 %102
%105 = bitcast double* %104 to <2 x double>*
%106 = load <2 x double>, <2 x double>* %105
%107 = fdiv <2 x double> %106, %27
%108 = bitcast double* %104 to <2 x double>*
store <2 x double> %107, <2 x double>* %108
br label %109
109: ; preds = %101, %103
br i1 %32, label %113, label %110
110: ; preds = %79, %109
%111 = phi i64 [ 0, %79 ], [ %25, %109 ]
br label %116
112: ; preds = %113, %16
ret double** %4
113: ; preds = %116, %109
%114 = add nuw nsw i64 %80, 1
%115 = icmp eq i64 %114, %1
br i1 %115, label %112, label %79
116: ; preds = %110, %116
%117 = phi i64 [ %121, %116 ], [ %111, %110 ]
%118 = getelementptr inbounds double, double* %82, i64 %117
%119 = load double, double* %118
%120 = fdiv double %119, %17
store double %120, double* %118
%121 = add nuw nsw i64 %117, 1
%122 = icmp eq i64 %121, %1
br i1 %122, label %113, label %116
}