double **allocateDblMatrix(unsigned n_rows, unsigned n_cols) {
  double **mat_p;
  double *arr;
  const int arr_len = n_rows * n_cols * (sizeof(double));
  unsigned i, j;
  ;
  assert(n_rows > 0);
  assert(n_cols > 0);
  mat_p = (double **)malloc(n_rows * sizeof(double *));
  if (mat_p == 0L) {
    PyErr_NoMemory();
    return 0L;
  }
  arr = (double *)malloc(arr_len);
  if (arr == 0L) {
    free(mat_p);
    PyErr_NoMemory();
    return 0L;
  }
  for (i = 0; i < n_rows; ++i) {
    mat_p[i] = arr;
    for (j = 0; j < n_cols; ++j)
      *arr++ = 0.0;
  }
  return mat_p;
}


x86 assembly (clang, -O0):
<allocateDblMatrix>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%rsp
mov    %edi,-0xc(%rbp)
mov    %esi,-0x10(%rbp)
mov    -0xc(%rbp),%eax
imul   -0x10(%rbp),%eax
mov    %eax,%eax
shl    $0x3,%rax
mov    %eax,-0x24(%rbp)
cmpl   $0x0,-0xc(%rbp)
seta   %al
and    $0x1,%al
movzbl %al,%edi
mov    $0x0,%al
call   1050 <assert@plt>
cmpl   $0x0,-0x10(%rbp)
seta   %al
and    $0x1,%al
movzbl %al,%edi
mov    $0x0,%al
call   1050 <assert@plt>
mov    -0xc(%rbp),%eax
mov    %eax,%edi
shl    $0x3,%rdi
call   1060 <malloc@plt>
mov    %rax,-0x18(%rbp)
cmpq   $0x0,-0x18(%rbp)
jne    11b5 <allocateDblMatrix+0x75>
mov    $0x0,%al
call   1040 <PyErr_NoMemory@plt>
movq   $0x0,-0x8(%rbp)
jmp    125e <allocateDblMatrix+0x11e>
movslq -0x24(%rbp),%rdi
call   1060 <malloc@plt>
mov    %rax,-0x20(%rbp)
cmpq   $0x0,-0x20(%rbp)
jne    11ea <allocateDblMatrix+0xaa>
mov    -0x18(%rbp),%rdi
call   1030 <free@plt>
mov    $0x0,%al
call   1040 <PyErr_NoMemory@plt>
movq   $0x0,-0x8(%rbp)
jmp    125e <allocateDblMatrix+0x11e>
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0xc(%rbp),%eax
jae    1256 <allocateDblMatrix+0x116>
mov    -0x20(%rbp),%rdx
mov    -0x18(%rbp),%rax
mov    -0x28(%rbp),%ecx
mov    %rdx,(%rax,%rcx,8)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x10(%rbp),%eax
jae    1243 <allocateDblMatrix+0x103>
mov    -0x20(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x20(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1213 <allocateDblMatrix+0xd3>
jmp    1248 <allocateDblMatrix+0x108>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    11f1 <allocateDblMatrix+0xb1>
mov    -0x18(%rbp),%rax
mov    %rax,-0x8(%rbp)
mov    -0x8(%rbp),%rax
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<allocateDblMatrix>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x18,%rsp
mov    %esi,%r15d
mov    %edi,%ebp
xor    %edi,%edi
test   %ebp,%ebp
setne  %dil
xor    %eax,%eax
call   1060 <assert@plt>
xor    %edi,%edi
test   %r15d,%r15d
setne  %dil
xor    %eax,%eax
call   1060 <assert@plt>
mov    %ebp,%r13d
lea    0x0(,%r13,8),%rdi
call   1070 <malloc@plt>
test   %rax,%rax
je     1232 <allocateDblMatrix+0xe2>
mov    %rax,%r14
mov    %ebp,%eax
imul   %r15d,%eax
shl    $0x3,%eax
movslq %eax,%rdi
call   1070 <malloc@plt>
test   %rax,%rax
je     122a <allocateDblMatrix+0xda>
test   %ebp,%ebp
je     123c <allocateDblMatrix+0xec>
mov    %rax,%rbx
mov    %r15d,%ebp
shl    $0x3,%rbp
lea    -0x1(%r15),%r12d
add    $0x1,%r12
lea    -0x1(%r13),%rax
mov    %r13d,%ecx
and    $0x3,%ecx
mov    %rcx,0x10(%rsp)
cmp    $0x3,%rax
jae    124e <allocateDblMatrix+0xfe>
xor    %r13d,%r13d
cmpq   $0x0,0x10(%rsp)
je     123c <allocateDblMatrix+0xec>
mov    %rbp,0x8(%rsp)
lea    (%r14,%r13,8),%r13
xor    %ebp,%ebp
jmp    120b <allocateDblMatrix+0xbb>
nopl   0x0(%rax,%rax,1)
add    $0x1,%rbp
cmp    %rbp,0x10(%rsp)
je     123c <allocateDblMatrix+0xec>
mov    %rbx,0x0(%r13,%rbp,8)
test   %r15d,%r15d
je     1200 <allocateDblMatrix+0xb0>
mov    %rbx,%rdi
xor    %esi,%esi
mov    0x8(%rsp),%rdx
call   1050 <memset@plt>
lea    (%rbx,%r12,8),%rbx
jmp    1200 <allocateDblMatrix+0xb0>
mov    %r14,%rdi
call   1030 <free@plt>
xor    %r14d,%r14d
xor    %eax,%eax
call   1040 <PyErr_NoMemory@plt>
mov    %r14,%rax
add    $0x18,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
and    $0xfffffffc,%r13d
mov    %r13,0x8(%rsp)
xor    %r13d,%r13d
jmp    126f <allocateDblMatrix+0x11f>
nopl   0x0(%rax)
add    $0x4,%r13
cmp    %r13,0x8(%rsp)
je     11e6 <allocateDblMatrix+0x96>
mov    %rbx,(%r14,%r13,8)
test   %r15d,%r15d
je     1289 <allocateDblMatrix+0x139>
mov    %rbx,%rdi
xor    %esi,%esi
mov    %rbp,%rdx
call   1050 <memset@plt>
lea    (%rbx,%r12,8),%rbx
mov    %rbx,0x8(%r14,%r13,8)
test   %r15d,%r15d
je     12a4 <allocateDblMatrix+0x154>
mov    %rbx,%rdi
xor    %esi,%esi
mov    %rbp,%rdx
call   1050 <memset@plt>
lea    (%rbx,%r12,8),%rbx
mov    %rbx,0x10(%r14,%r13,8)
test   %r15d,%r15d
je     12bf <allocateDblMatrix+0x16f>
mov    %rbx,%rdi
xor    %esi,%esi
mov    %rbp,%rdx
call   1050 <memset@plt>
lea    (%rbx,%r12,8),%rbx
mov    %rbx,0x18(%r14,%r13,8)
test   %r15d,%r15d
je     1260 <allocateDblMatrix+0x110>
mov    %rbx,%rdi
xor    %esi,%esi
mov    %rbp,%rdx
call   1050 <memset@plt>
lea    (%rbx,%r12,8),%rbx
jmp    1260 <allocateDblMatrix+0x110>

LLVM IR (-O0):
define double** @allocateDblMatrix(i32 %0, i32 %1) {
%3 = alloca double**
%4 = alloca i32
%5 = alloca i32
%6 = alloca double**
%7 = alloca double*
%8 = alloca i32
%9 = alloca i32
%10 = alloca i32
store i32 %0, i32* %4
store i32 %1, i32* %5
%11 = load i32, i32* %4
%12 = load i32, i32* %5
%13 = mul i32 %11, %12
%14 = zext i32 %13 to i64
%15 = mul i64 %14, 8
%16 = trunc i64 %15 to i32
store i32 %16, i32* %8
%17 = load i32, i32* %4
%18 = icmp ugt i32 %17, 0
%19 = zext i1 %18 to i32
%20 = call i32 (i32, ...) bitcast (i32 (...)* @assert to i32 (i32, ...)*)(i32 %19)
%21 = load i32, i32* %5
%22 = icmp ugt i32 %21, 0
%23 = zext i1 %22 to i32
%24 = call i32 (i32, ...) bitcast (i32 (...)* @assert to i32 (i32, ...)*)(i32 %23)
%25 = load i32, i32* %4
%26 = zext i32 %25 to i64
%27 = mul i64 %26, 8
%28 = call i8* @malloc(i64 %27)
%29 = bitcast i8* %28 to double**
store double** %29, double*** %6
%30 = load double**, double*** %6
%31 = icmp eq double** %30, null
br i1 %31, label %32, label %34
32: ; preds = %2
%33 = call i32 (...) @PyErr_NoMemory()
store double** null, double*** %3
br label %72
34: ; preds = %2
%35 = load i32, i32* %8
%36 = sext i32 %35 to i64
%37 = call i8* @malloc(i64 %36)
%38 = bitcast i8* %37 to double*
store double* %38, double** %7
%39 = load double*, double** %7
%40 = icmp eq double* %39, null
br i1 %40, label %41, label %45
41: ; preds = %34
%42 = load double**, double*** %6
%43 = bitcast double** %42 to i8*
%44 = call i32 (...) @PyErr_NoMemory()
store double** null, double*** %3
br label %72
45: ; preds = %34
store i32 0, i32* %9
br label %46
46: ; preds = %67, %45
%47 = load i32, i32* %9
%48 = load i32, i32* %4
%49 = icmp ult i32 %47, %48
br i1 %49, label %50, label %70
50: ; preds = %46
%51 = load double*, double** %7
%52 = load double**, double*** %6
%53 = load i32, i32* %9
%54 = zext i32 %53 to i64
%55 = getelementptr inbounds double*, double** %52, i64 %54
store double* %51, double** %55
store i32 0, i32* %10
br label %56
56: ; preds = %63, %50
%57 = load i32, i32* %10
%58 = load i32, i32* %5
%59 = icmp ult i32 %57, %58
br i1 %59, label %60, label %66
60: ; preds = %56
%61 = load double*, double** %7
%62 = getelementptr inbounds double, double* %61, i32 1
store double* %62, double** %7
store double 0.000000e+00, double* %61
br label %63
63: ; preds = %60
%64 = load i32, i32* %10
%65 = add i32 %64, 1
store i32 %65, i32* %10
br label %56
66: ; preds = %56
br label %67
67: ; preds = %66
%68 = load i32, i32* %9
%69 = add i32 %68, 1
store i32 %69, i32* %9
br label %46
70: ; preds = %46
%71 = load double**, double*** %6
store double** %71, double*** %3
br label %72
72: ; preds = %70, %41, %32
%73 = load double**, double*** %3
ret double** %73
}

LLVM IR (-O2):
define double** @allocateDblMatrix(i32 %0, i32 %1) {
%3 = icmp ne i32 %0, 0
%4 = zext i1 %3 to i32
%5 = tail call i32 (i32, ...) bitcast (i32 (...)* @assert to i32 (i32, ...)*)(i32 %4)
%6 = icmp ne i32 %1, 0
%7 = zext i1 %6 to i32
%8 = tail call i32 (i32, ...) bitcast (i32 (...)* @assert to i32 (i32, ...)*)(i32 %7)
%9 = zext i32 %0 to i64
%10 = shl nuw nsw i64 %9, 3
%11 = tail call i8* @malloc(i64 %10)
%12 = bitcast i8* %11 to double**
%13 = icmp eq i8* %11, null
br i1 %13, label %70, label %14
14: ; preds = %2
%15 = shl i32 %0, 3
%16 = mul i32 %15, %1
%17 = sext i32 %16 to i64
%18 = tail call i8* @malloc(i64 %17)
%19 = icmp eq i8* %18, null
br i1 %19, label %20, label %21
20: ; preds = %14
tail call void @free(i8* %11)
br label %70
21: ; preds = %14
%22 = icmp eq i32 %0, 0
br i1 %22, label %89, label %23
23: ; preds = %21
%24 = bitcast i8* %18 to double*
%25 = icmp eq i32 %1, 0
%26 = zext i32 %1 to i64
%27 = shl nuw nsw i64 %26, 3
%28 = add i32 %1, -1
%29 = zext i32 %28 to i64
%30 = add nuw nsw i64 %29, 1
%31 = add nsw i64 %9, -1
%32 = and i64 %9, 3
%33 = icmp ult i64 %31, 3
br i1 %33, label %72, label %34
34: ; preds = %23
%35 = and i64 %9, 4294967292
br label %36
36: ; preds = %65, %34
%37 = phi i64 [ 0, %34 ], [ %67, %65 ]
%38 = phi double* [ %24, %34 ], [ %66, %65 ]
%39 = phi i64 [ 0, %34 ], [ %68, %65 ]
%40 = getelementptr inbounds double*, double** %12, i64 %37
store double* %38, double** %40
br i1 %25, label %44, label %41
41: ; preds = %36
%42 = bitcast double* %38 to i8*
%43 = getelementptr double, double* %38, i64 %30
br label %44
44: ; preds = %41, %36
%45 = phi double* [ %38, %36 ], [ %43, %41 ]
%46 = or i64 %37, 1
%47 = getelementptr inbounds double*, double** %12, i64 %46
store double* %45, double** %47
br i1 %25, label %51, label %48
48: ; preds = %44
%49 = bitcast double* %45 to i8*
%50 = getelementptr double, double* %45, i64 %30
br label %51
51: ; preds = %48, %44
%52 = phi double* [ %45, %44 ], [ %50, %48 ]
%53 = or i64 %37, 2
%54 = getelementptr inbounds double*, double** %12, i64 %53
store double* %52, double** %54
br i1 %25, label %58, label %55
55: ; preds = %51
%56 = bitcast double* %52 to i8*
%57 = getelementptr double, double* %52, i64 %30
br label %58
58: ; preds = %55, %51
%59 = phi double* [ %52, %51 ], [ %57, %55 ]
%60 = or i64 %37, 3
%61 = getelementptr inbounds double*, double** %12, i64 %60
store double* %59, double** %61
br i1 %25, label %65, label %62
62: ; preds = %58
%63 = bitcast double* %59 to i8*
%64 = getelementptr double, double* %59, i64 %30
br label %65
65: ; preds = %62, %58
%66 = phi double* [ %59, %58 ], [ %64, %62 ]
%67 = add nuw nsw i64 %37, 4
%68 = add i64 %39, 4
%69 = icmp eq i64 %68, %35
br i1 %69, label %72, label %36
70: ; preds = %2, %20
%71 = tail call i32 (...) @PyErr_NoMemory()
br label %89
72: ; preds = %65, %23
%73 = phi i64 [ 0, %23 ], [ %67, %65 ]
%74 = phi double* [ %24, %23 ], [ %66, %65 ]
%75 = icmp eq i64 %32, 0
br i1 %75, label %89, label %76
76: ; preds = %72, %84
%77 = phi i64 [ %86, %84 ], [ %73, %72 ]
%78 = phi double* [ %85, %84 ], [ %74, %72 ]
%79 = phi i64 [ %87, %84 ], [ 0, %72 ]
%80 = getelementptr inbounds double*, double** %12, i64 %77
store double* %78, double** %80
br i1 %25, label %84, label %81
81: ; preds = %76
%82 = bitcast double* %78 to i8*
%83 = getelementptr double, double* %78, i64 %30
br label %84
84: ; preds = %81, %76
%85 = phi double* [ %78, %76 ], [ %83, %81 ]
%86 = add nuw nsw i64 %77, 1
%87 = add i64 %79, 1
%88 = icmp eq i64 %87, %32
br i1 %88, label %89, label %76
89: ; preds = %72, %84, %70, %21
%90 = phi double** [ %12, %21 ], [ null, %70 ], [ %12, %84 ], [ %12, %72 ]
ret double** %90
}