#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 galois_single_divide(int, unsigned int, unsigned int);
unsigned int jerasure_invert_bitmatrix(int *, int *, unsigned int);
unsigned int jerasure_invertible_bitmatrix(int *, unsigned int);
int *jerasure_matrix_multiply(int *, int *, unsigned int, unsigned int,
                              unsigned int, unsigned int, int);
int *jerasure_matrix_to_bitmatrix(unsigned int, unsigned int, unsigned int,
                                  int *);
int /*<<< orphan*/ jerasure_print_bitmatrix(int *, unsigned int, unsigned int,
                                            unsigned int);
scalar_t__ malloc(int);
int /*<<< orphan*/ memcpy(int *, int *, int);
int /*<<< orphan*/ printf(char *, ...);
scalar_t__ sscanf(char *, char *, unsigned int *);
int /*<<< orphan*/ usage(char *);
int main(int argc, char **argv) {
  unsigned int k, w, i, j, n;
  int *matrix;
  int *bitmatrix;
  int *bitmatrix_copy;
  int *inverse;
  int *identity;
  if (argc != 3)
    usage(NULL);
  if (sscanf(argv[1], "%d", &k) == 0 || k <= 0)
    usage("Bad k");
  if (sscanf(argv[2], "%d", &w) == 0 || w <= 0 || w > 31)
    usage("Bad w");
  if (k >= (1 << w))
    usage("K too big");
  matrix = (int *)malloc(sizeof(int) * (k * k));
  bitmatrix_copy = (int *)malloc(sizeof(int) * (k * w * k * w));
  inverse = (int *)malloc(sizeof(int) * (k * w * k * w));
  for (i = 0; i < k; i++) {
    for (j = 0; j < k; j++) {
      n = i ^ ((1 << w) - 1 - j);
      matrix[i * k + j] = (n == 0) ? 0 : galois_single_divide(1, n, w);
    }
  }
  bitmatrix = jerasure_matrix_to_bitmatrix(k, k, w, matrix);
  printf("The Cauchy Bit-Matrix:\n");
  jerasure_print_bitmatrix(bitmatrix, k * w, k * w, w);
  memcpy(bitmatrix_copy, bitmatrix, sizeof(int) * k * w * k * w);
  i = jerasure_invertible_bitmatrix(bitmatrix_copy, k * w);
  printf("\nInvertible: %s\n", (i == 1) ? "Yes" : "No");
  if (i == 1) {
    printf("\nInverse:\n");
    memcpy(bitmatrix_copy, bitmatrix, sizeof(int) * k * w * k * w);
    i = jerasure_invert_bitmatrix(bitmatrix_copy, inverse, k * w);
    jerasure_print_bitmatrix(inverse, k * w, k * w, w);
    identity = jerasure_matrix_multiply(inverse, bitmatrix, k * w, k * w, k * w,
                                        k * w, 2);
    printf("\nInverse times matrix (should be identity):\n");
    jerasure_print_bitmatrix(identity, k * w, k * w, w);
  }
  return 0;
}


x86 assembly (clang, -O0):
<main>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
movl   $0x0,-0x4(%rbp)
mov    %edi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
cmpl   $0x3,-0x8(%rbp)
je     11d9 <main+0x29>
xor    %eax,%eax
mov    %eax,%edi
call   10b0 <usage@plt>
mov    -0x10(%rbp),%rax
mov    0x8(%rax),%rdi
lea    0xe18(%rip),%rsi
lea    -0x14(%rbp),%rdx
call   10a0 <sscanf@plt>
cmp    $0x0,%rax
je     1205 <main+0x55>
cmpl   $0x0,-0x14(%rbp)
ja     1211 <main+0x61>
lea    0xdf7(%rip),%rdi
call   10b0 <usage@plt>
mov    -0x10(%rbp),%rax
mov    0x10(%rax),%rdi
lea    0xde0(%rip),%rsi
lea    -0x18(%rbp),%rdx
call   10a0 <sscanf@plt>
cmp    $0x0,%rax
je     1247 <main+0x97>
cmpl   $0x0,-0x18(%rbp)
jbe    1247 <main+0x97>
cmpl   $0x1f,-0x18(%rbp)
jbe    1253 <main+0xa3>
lea    0xdbb(%rip),%rdi
call   10b0 <usage@plt>
mov    -0x14(%rbp),%eax
mov    -0x18(%rbp),%ecx
mov    $0x1,%edx
shl    %cl,%edx
mov    %edx,%ecx
cmp    %ecx,%eax
jb     1276 <main+0xc6>
lea    0xd9e(%rip),%rdi
call   10b0 <usage@plt>
mov    -0x14(%rbp),%eax
imul   -0x14(%rbp),%eax
mov    %eax,%eax
shl    $0x2,%rax
mov    %eax,%edi
call   1090 <malloc@plt>
mov    %rax,-0x30(%rbp)
mov    -0x14(%rbp),%eax
imul   -0x18(%rbp),%eax
imul   -0x14(%rbp),%eax
imul   -0x18(%rbp),%eax
mov    %eax,%eax
shl    $0x2,%rax
mov    %eax,%edi
call   1090 <malloc@plt>
mov    %rax,-0x40(%rbp)
mov    -0x14(%rbp),%eax
imul   -0x18(%rbp),%eax
imul   -0x14(%rbp),%eax
imul   -0x18(%rbp),%eax
mov    %eax,%eax
shl    $0x2,%rax
mov    %eax,%edi
call   1090 <malloc@plt>
mov    %rax,-0x48(%rbp)
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jae    136c <main+0x1bc>
movl   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x14(%rbp),%eax
jae    1359 <main+0x1a9>
mov    -0x1c(%rbp),%eax
mov    -0x18(%rbp),%ecx
mov    $0x1,%edx
shl    %cl,%edx
mov    %edx,%ecx
sub    $0x1,%ecx
sub    -0x20(%rbp),%ecx
xor    %ecx,%eax
mov    %eax,-0x24(%rbp)
cmpl   $0x0,-0x24(%rbp)
jne    1322 <main+0x172>
xor    %eax,%eax
mov    %eax,-0x54(%rbp)
jmp    1335 <main+0x185>
mov    -0x24(%rbp),%esi
mov    -0x18(%rbp),%edx
mov    $0x1,%edi
call   1050 <galois_single_divide@plt>
mov    %eax,-0x54(%rbp)
mov    -0x54(%rbp),%edx
mov    -0x30(%rbp),%rax
mov    -0x1c(%rbp),%ecx
imul   -0x14(%rbp),%ecx
add    -0x20(%rbp),%ecx
mov    %ecx,%ecx
mov    %edx,(%rax,%rcx,4)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    12e8 <main+0x138>
jmp    135e <main+0x1ae>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    12d5 <main+0x125>
mov    -0x14(%rbp),%edi
mov    -0x14(%rbp),%esi
mov    -0x18(%rbp),%edx
mov    -0x30(%rbp),%rcx
call   1030 <jerasure_matrix_to_bitmatrix@plt>
mov    %rax,-0x38(%rbp)
lea    0xc90(%rip),%rdi
mov    $0x0,%al
call   1060 <printf@plt>
mov    -0x38(%rbp),%rdi
mov    -0x14(%rbp),%esi
imul   -0x18(%rbp),%esi
mov    -0x14(%rbp),%edx
imul   -0x18(%rbp),%edx
mov    -0x18(%rbp),%ecx
call   10d0 <jerasure_print_bitmatrix@plt>
mov    -0x40(%rbp),%rdi
mov    -0x38(%rbp),%rsi
mov    -0x14(%rbp),%eax
shl    $0x2,%rax
mov    -0x18(%rbp),%ecx
imul   %rcx,%rax
mov    -0x14(%rbp),%ecx
imul   %rcx,%rax
mov    -0x18(%rbp),%ecx
imul   %rcx,%rax
mov    %eax,%edx
call   1080 <memcpy@plt>
mov    -0x40(%rbp),%rdi
mov    -0x14(%rbp),%esi
imul   -0x18(%rbp),%esi
call   1040 <jerasure_invertible_bitmatrix@plt>
mov    %eax,-0x1c(%rbp)
mov    -0x1c(%rbp),%ecx
lea    0xc54(%rip),%rsi
lea    0xc49(%rip),%rax
cmp    $0x1,%ecx
cmove  %rax,%rsi
lea    0xc2a(%rip),%rdi
mov    $0x0,%al
call   1060 <printf@plt>
cmpl   $0x1,-0x1c(%rbp)
jne    14e2 <main+0x332>
lea    0xc2a(%rip),%rdi
mov    $0x0,%al
call   1060 <printf@plt>
mov    -0x40(%rbp),%rdi
mov    -0x38(%rbp),%rsi
mov    -0x14(%rbp),%eax
shl    $0x2,%rax
mov    -0x18(%rbp),%ecx
imul   %rcx,%rax
mov    -0x14(%rbp),%ecx
imul   %rcx,%rax
mov    -0x18(%rbp),%ecx
imul   %rcx,%rax
mov    %eax,%edx
call   1080 <memcpy@plt>
mov    -0x40(%rbp),%rdi
mov    -0x48(%rbp),%rsi
mov    -0x14(%rbp),%edx
imul   -0x18(%rbp),%edx
call   10c0 <jerasure_invert_bitmatrix@plt>
mov    %eax,-0x1c(%rbp)
mov    -0x48(%rbp),%rdi
mov    -0x14(%rbp),%esi
imul   -0x18(%rbp),%esi
mov    -0x14(%rbp),%edx
imul   -0x18(%rbp),%edx
mov    -0x18(%rbp),%ecx
call   10d0 <jerasure_print_bitmatrix@plt>
mov    -0x48(%rbp),%rdi
mov    -0x38(%rbp),%rsi
mov    -0x14(%rbp),%edx
imul   -0x18(%rbp),%edx
mov    -0x14(%rbp),%ecx
imul   -0x18(%rbp),%ecx
mov    -0x14(%rbp),%r8d
imul   -0x18(%rbp),%r8d
mov    -0x14(%rbp),%r9d
imul   -0x18(%rbp),%r9d
movl   $0x2,(%rsp)
call   1070 <jerasure_matrix_multiply@plt>
mov    %rax,-0x50(%rbp)
lea    0xb93(%rip),%rdi
mov    $0x0,%al
call   1060 <printf@plt>
mov    -0x50(%rbp),%rdi
mov    -0x14(%rbp),%esi
imul   -0x18(%rbp),%esi
mov    -0x14(%rbp),%edx
imul   -0x18(%rbp),%edx
mov    -0x18(%rbp),%ecx
call   10d0 <jerasure_print_bitmatrix@plt>
xor    %eax,%eax
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<main>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x18,%rsp
mov    %rsi,%rbx
cmp    $0x3,%edi
je     11dd <main+0x1d>
xor    %edi,%edi
call   10c0 <usage@plt>
mov    0x8(%rbx),%rdi
lea    0xe18(%rip),%rsi
lea    0x10(%rsp),%rdx
call   10b0 <sscanf@plt>
test   %rax,%rax
je     11fe <main+0x3e>
cmpl   $0x0,0x10(%rsp)
jne    120a <main+0x4a>
lea    0xdfe(%rip),%rdi
call   10c0 <usage@plt>
mov    0x10(%rbx),%rdi
lea    0xdeb(%rip),%rsi
lea    0x14(%rsp),%rdx
call   10b0 <sscanf@plt>
test   %rax,%rax
je     1231 <main+0x71>
mov    0x14(%rsp),%ecx
test   %ecx,%ecx
je     1231 <main+0x71>
cmp    $0x20,%ecx
jb     1241 <main+0x81>
lea    0xdd1(%rip),%rdi
call   10c0 <usage@plt>
mov    0x14(%rsp),%ecx
mov    0x10(%rsp),%edi
mov    %edi,%eax
shr    %cl,%eax
test   %eax,%eax
je     125d <main+0x9d>
lea    0xdbb(%rip),%rdi
call   10c0 <usage@plt>
mov    0x10(%rsp),%edi
imul   %edi,%edi
shl    $0x2,%edi
call   10a0 <malloc@plt>
mov    %rax,%r12
mov    0x14(%rsp),%edi
imul   0x10(%rsp),%edi
imul   %edi,%edi
shl    $0x2,%edi
call   10a0 <malloc@plt>
mov    %rax,%r15
mov    0x14(%rsp),%edi
imul   0x10(%rsp),%edi
imul   %edi,%edi
shl    $0x2,%edi
call   10a0 <malloc@plt>
mov    %rax,%r14
mov    0x10(%rsp),%edi
test   %edi,%edi
je     1310 <main+0x150>
xor    %r13d,%r13d
jmp    12b9 <main+0xf9>
cs nopw 0x0(%rax,%rax,1)
add    $0x1,%r13d
cmp    %edi,%r13d
jae    1312 <main+0x152>
mov    $0xffffffff,%ebx
xor    %ebp,%ebp
jmp    12f8 <main+0x138>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
mov    $0x1,%edi
mov    %ecx,%edx
call   1060 <galois_single_divide@plt>
mov    0x10(%rsp),%edi
imul   %r13d,%edi
add    %ebp,%edi
mov    %eax,(%r12,%rdi,4)
add    $0x1,%ebp
mov    0x10(%rsp),%edi
add    $0xffffffff,%ebx
cmp    %edi,%ebp
jae    12b0 <main+0xf0>
mov    0x14(%rsp),%ecx
mov    $0xffffffff,%eax
shl    %cl,%eax
mov    %ebx,%esi
sub    %eax,%esi
xor    %r13d,%esi
jne    12d0 <main+0x110>
xor    %eax,%eax
jmp    12e0 <main+0x120>
xor    %edi,%edi
mov    0x14(%rsp),%edx
mov    %edi,%esi
mov    %r12,%rcx
call   1040 <jerasure_matrix_to_bitmatrix@plt>
mov    %rax,%rbx
lea    0xd07(%rip),%rdi
call   1030 <puts@plt>
mov    0x14(%rsp),%ecx
mov    0x10(%rsp),%edx
imul   %ecx,%edx
mov    %rbx,%rdi
mov    %edx,%esi
call   10e0 <jerasure_print_bitmatrix@plt>
mov    0x14(%rsp),%edx
imul   0x10(%rsp),%edx
imul   %edx,%edx
shl    $0x2,%edx
mov    %r15,%rdi
mov    %rbx,%rsi
call   1090 <memcpy@plt>
mov    0x14(%rsp),%esi
imul   0x10(%rsp),%esi
mov    %r15,%rdi
call   1050 <jerasure_invertible_bitmatrix@plt>
mov    %eax,%ebp
cmp    $0x1,%eax
lea    0xcaf(%rip),%rax
lea    0xcac(%rip),%rsi
cmove  %rax,%rsi
lea    0xc8c(%rip),%rdi
xor    %eax,%eax
call   1070 <printf@plt>
cmp    $0x1,%ebp
jne    1433 <main+0x273>
lea    0xca4(%rip),%rdi
call   1030 <puts@plt>
mov    0x14(%rsp),%edx
imul   0x10(%rsp),%edx
imul   %edx,%edx
shl    $0x2,%edx
mov    %r15,%rdi
mov    %rbx,%rsi
call   1090 <memcpy@plt>
mov    0x14(%rsp),%edx
imul   0x10(%rsp),%edx
mov    %r15,%rdi
mov    %r14,%rsi
call   10d0 <jerasure_invert_bitmatrix@plt>
mov    0x14(%rsp),%ecx
mov    0x10(%rsp),%edx
imul   %ecx,%edx
mov    %r14,%rdi
mov    %edx,%esi
call   10e0 <jerasure_print_bitmatrix@plt>
mov    0x14(%rsp),%ecx
imul   0x10(%rsp),%ecx
movl   $0x2,(%rsp)
mov    %r14,%rdi
mov    %rbx,%rsi
mov    %ecx,%edx
mov    %ecx,%r8d
mov    %ecx,%r9d
call   1080 <jerasure_matrix_multiply@plt>
mov    %rax,%rbx
lea    0xc39(%rip),%rdi
call   1030 <puts@plt>
mov    0x14(%rsp),%ecx
mov    0x10(%rsp),%edx
imul   %ecx,%edx
mov    %rbx,%rdi
mov    %edx,%esi
call   10e0 <jerasure_print_bitmatrix@plt>
xor    %eax,%eax
add    $0x18,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
@.str = [3 x i8] c"%d\00"
@.str.1 = [6 x i8] c"Bad k\00"
@.str.2 = [6 x i8] c"Bad w\00"
@.str.3 = [10 x i8] c"K too big\00"
@.str.4 = [24 x i8] c"The Cauchy Bit-Matrix:\0A\00"
@.str.5 = [17 x i8] c"\0AInvertible: %s\0A\00"
@.str.6 = [4 x i8] c"Yes\00"
@.str.7 = [3 x i8] c"No\00"
@.str.8 = [11 x i8] c"\0AInverse:\0A\00"
@.str.9 = [45 x i8] c"\0AInverse times matrix (should be identity):\0A\00"
define i32 @main(i32 %0, i8** %1) {
%3 = alloca i32
%4 = alloca i32
%5 = alloca i8**
%6 = alloca i32
%7 = alloca i32
%8 = alloca i32
%9 = alloca i32
%10 = alloca i32
%11 = alloca i32*
%12 = alloca i32*
%13 = alloca i32*
%14 = alloca i32*
%15 = alloca i32*
store i32 0, i32* %3
store i32 %0, i32* %4
store i8** %1, i8*** %5
%16 = load i32, i32* %4
%17 = icmp ne i32 %16, 3
br i1 %17, label %18, label %20
18: ; preds = %2
%19 = call i32 @usage(i8* null)
br label %20
20: ; preds = %18, %2
%21 = load i8**, i8*** %5
%22 = getelementptr inbounds i8*, i8** %21, i64 1
%23 = load i8*, i8** %22
%24 = call i64 @sscanf(i8* %23, i8* getelementptr inbounds ([3 x i8], [3 x i8]* @.str, i64 0, i64 0), i32* %6)
%25 = icmp eq i64 %24, 0
br i1 %25, label %29, label %26
26: ; preds = %20
%27 = load i32, i32* %6
%28 = icmp ule i32 %27, 0
br i1 %28, label %29, label %31
29: ; preds = %26, %20
%30 = call i32 @usage(i8* getelementptr inbounds ([6 x i8], [6 x i8]* @.str.1, i64 0, i64 0))
br label %31
31: ; preds = %29, %26
%32 = load i8**, i8*** %5
%33 = getelementptr inbounds i8*, i8** %32, i64 2
%34 = load i8*, i8** %33
%35 = call i64 @sscanf(i8* %34, i8* getelementptr inbounds ([3 x i8], [3 x i8]* @.str, i64 0, i64 0), i32* %7)
%36 = icmp eq i64 %35, 0
br i1 %36, label %43, label %37
37: ; preds = %31
%38 = load i32, i32* %7
%39 = icmp ule i32 %38, 0
br i1 %39, label %43, label %40
40: ; preds = %37
%41 = load i32, i32* %7
%42 = icmp ugt i32 %41, 31
br i1 %42, label %43, label %45
43: ; preds = %40, %37, %31
%44 = call i32 @usage(i8* getelementptr inbounds ([6 x i8], [6 x i8]* @.str.2, i64 0, i64 0))
br label %45
45: ; preds = %43, %40
%46 = load i32, i32* %6
%47 = load i32, i32* %7
%48 = shl i32 1, %47
%49 = icmp uge i32 %46, %48
br i1 %49, label %50, label %52
50: ; preds = %45
%51 = call i32 @usage(i8* getelementptr inbounds ([10 x i8], [10 x i8]* @.str.3, i64 0, i64 0))
br label %52
52: ; preds = %50, %45
%53 = load i32, i32* %6
%54 = load i32, i32* %6
%55 = mul i32 %53, %54
%56 = zext i32 %55 to i64
%57 = mul i64 4, %56
%58 = trunc i64 %57 to i32
%59 = call i64 @malloc(i32 %58)
%60 = inttoptr i64 %59 to i32*
store i32* %60, i32** %11
%61 = load i32, i32* %6
%62 = load i32, i32* %7
%63 = mul i32 %61, %62
%64 = load i32, i32* %6
%65 = mul i32 %63, %64
%66 = load i32, i32* %7
%67 = mul i32 %65, %66
%68 = zext i32 %67 to i64
%69 = mul i64 4, %68
%70 = trunc i64 %69 to i32
%71 = call i64 @malloc(i32 %70)
%72 = inttoptr i64 %71 to i32*
store i32* %72, i32** %13
%73 = load i32, i32* %6
%74 = load i32, i32* %7
%75 = mul i32 %73, %74
%76 = load i32, i32* %6
%77 = mul i32 %75, %76
%78 = load i32, i32* %7
%79 = mul i32 %77, %78
%80 = zext i32 %79 to i64
%81 = mul i64 4, %80
%82 = trunc i64 %81 to i32
%83 = call i64 @malloc(i32 %82)
%84 = inttoptr i64 %83 to i32*
store i32* %84, i32** %14
store i32 0, i32* %8
br label %85
85: ; preds = %123, %52
%86 = load i32, i32* %8
%87 = load i32, i32* %6
%88 = icmp ult i32 %86, %87
br i1 %88, label %89, label %126
89: ; preds = %85
store i32 0, i32* %9
br label %90
90: ; preds = %119, %89
%91 = load i32, i32* %9
%92 = load i32, i32* %6
%93 = icmp ult i32 %91, %92
br i1 %93, label %94, label %122
94: ; preds = %90
%95 = load i32, i32* %8
%96 = load i32, i32* %7
%97 = shl i32 1, %96
%98 = sub nsw i32 %97, 1
%99 = load i32, i32* %9
%100 = sub i32 %98, %99
%101 = xor i32 %95, %100
store i32 %101, i32* %10
%102 = load i32, i32* %10
%103 = icmp eq i32 %102, 0
br i1 %103, label %104, label %105
104: ; preds = %94
br label %109
105: ; preds = %94
%106 = load i32, i32* %10
%107 = load i32, i32* %7
%108 = call i32 @galois_single_divide(i32 1, i32 %106, i32 %107)
br label %109
109: ; preds = %105, %104
%110 = phi i32 [ 0, %104 ], [ %108, %105 ]
%111 = load i32*, i32** %11
%112 = load i32, i32* %8
%113 = load i32, i32* %6
%114 = mul i32 %112, %113
%115 = load i32, i32* %9
%116 = add i32 %114, %115
%117 = zext i32 %116 to i64
%118 = getelementptr inbounds i32, i32* %111, i64 %117
store i32 %110, i32* %118
br label %119
119: ; preds = %109
%120 = load i32, i32* %9
%121 = add i32 %120, 1
store i32 %121, i32* %9
br label %90
122: ; preds = %90
br label %123
123: ; preds = %122
%124 = load i32, i32* %8
%125 = add i32 %124, 1
store i32 %125, i32* %8
br label %85
126: ; preds = %85
%127 = load i32, i32* %6
%128 = load i32, i32* %6
%129 = load i32, i32* %7
%130 = load i32*, i32** %11
%131 = call i32* @jerasure_matrix_to_bitmatrix(i32 %127, i32 %128, i32 %129, i32* %130)
store i32* %131, i32** %12
%132 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([24 x i8], [24 x i8]* @.str.4, i64 0, i64 0))
%133 = load i32*, i32** %12
%134 = load i32, i32* %6
%135 = load i32, i32* %7
%136 = mul i32 %134, %135
%137 = load i32, i32* %6
%138 = load i32, i32* %7
%139 = mul i32 %137, %138
%140 = load i32, i32* %7
%141 = call i32 @jerasure_print_bitmatrix(i32* %133, i32 %136, i32 %139, i32 %140)
%142 = load i32*, i32** %13
%143 = load i32*, i32** %12
%144 = load i32, i32* %6
%145 = zext i32 %144 to i64
%146 = mul i64 4, %145
%147 = load i32, i32* %7
%148 = zext i32 %147 to i64
%149 = mul i64 %146, %148
%150 = load i32, i32* %6
%151 = zext i32 %150 to i64
%152 = mul i64 %149, %151
%153 = load i32, i32* %7
%154 = zext i32 %153 to i64
%155 = mul i64 %152, %154
%156 = trunc i64 %155 to i32
%157 = call i32 @memcpy(i32* %142, i32* %143, i32 %156)
%158 = load i32*, i32** %13
%159 = load i32, i32* %6
%160 = load i32, i32* %7
%161 = mul i32 %159, %160
%162 = call i32 @jerasure_invertible_bitmatrix(i32* %158, i32 %161)
store i32 %162, i32* %8
%163 = load i32, i32* %8
%164 = icmp eq i32 %163, 1
%165 = zext i1 %164 to i64
%166 = select i1 %164, i8* getelementptr inbounds ([4 x i8], [4 x i8]* @.str.6, i64 0, i64 0), i8* getelementptr inbounds ([3 x i8], [3 x i8]* @.str.7, i64 0, i64 0)
%167 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([17 x i8], [17 x i8]* @.str.5, i64 0, i64 0), i8* %166)
%168 = load i32, i32* %8
%169 = icmp eq i32 %168, 1
br i1 %169, label %170, label %228
170: ; preds = %126
%171 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([11 x i8], [11 x i8]* @.str.8, i64 0, i64 0))
%172 = load i32*, i32** %13
%173 = load i32*, i32** %12
%174 = load i32, i32* %6
%175 = zext i32 %174 to i64
%176 = mul i64 4, %175
%177 = load i32, i32* %7
%178 = zext i32 %177 to i64
%179 = mul i64 %176, %178
%180 = load i32, i32* %6
%181 = zext i32 %180 to i64
%182 = mul i64 %179, %181
%183 = load i32, i32* %7
%184 = zext i32 %183 to i64
%185 = mul i64 %182, %184
%186 = trunc i64 %185 to i32
%187 = call i32 @memcpy(i32* %172, i32* %173, i32 %186)
%188 = load i32*, i32** %13
%189 = load i32*, i32** %14
%190 = load i32, i32* %6
%191 = load i32, i32* %7
%192 = mul i32 %190, %191
%193 = call i32 @jerasure_invert_bitmatrix(i32* %188, i32* %189, i32 %192)
store i32 %193, i32* %8
%194 = load i32*, i32** %14
%195 = load i32, i32* %6
%196 = load i32, i32* %7
%197 = mul i32 %195, %196
%198 = load i32, i32* %6
%199 = load i32, i32* %7
%200 = mul i32 %198, %199
%201 = load i32, i32* %7
%202 = call i32 @jerasure_print_bitmatrix(i32* %194, i32 %197, i32 %200, i32 %201)
%203 = load i32*, i32** %14
%204 = load i32*, i32** %12
%205 = load i32, i32* %6
%206 = load i32, i32* %7
%207 = mul i32 %205, %206
%208 = load i32, i32* %6
%209 = load i32, i32* %7
%210 = mul i32 %208, %209
%211 = load i32, i32* %6
%212 = load i32, i32* %7
%213 = mul i32 %211, %212
%214 = load i32, i32* %6
%215 = load i32, i32* %7
%216 = mul i32 %214, %215
%217 = call i32* @jerasure_matrix_multiply(i32* %203, i32* %204, i32 %207, i32 %210, i32 %213, i32 %216, i32 2)
store i32* %217, i32** %15
%218 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([45 x i8], [45 x i8]* @.str.9, i64 0, i64 0))
%219 = load i32*, i32** %15
%220 = load i32, i32* %6
%221 = load i32, i32* %7
%222 = mul i32 %220, %221
%223 = load i32, i32* %6
%224 = load i32, i32* %7
%225 = mul i32 %223, %224
%226 = load i32, i32* %7
%227 = call i32 @jerasure_print_bitmatrix(i32* %219, i32 %222, i32 %225, i32 %226)
br label %228
228: ; preds = %170, %126
ret i32 0
}

LLVM IR (-O2):
@.str = [3 x i8] c"%d\00"
@.str.1 = [6 x i8] c"Bad k\00"
@.str.2 = [6 x i8] c"Bad w\00"
@.str.3 = [10 x i8] c"K too big\00"
@.str.5 = [17 x i8] c"\0AInvertible: %s\0A\00"
@.str.6 = [4 x i8] c"Yes\00"
@.str.7 = [3 x i8] c"No\00"
@str = [23 x i8] c"The Cauchy Bit-Matrix:\00"
@str.10 = [10 x i8] c"\0AInverse:\00"
@str.11 = [44 x i8] c"\0AInverse times matrix (should be identity):\00"
define i32 @main(i32 %0, i8** nocapture %1) {
%3 = alloca i32
%4 = alloca i32
%5 = bitcast i32* %3 to i8*
%6 = bitcast i32* %4 to i8*
%7 = icmp eq i32 %0, 3
br i1 %7, label %10, label %8
8: ; preds = %2
%9 = tail call i32 @usage(i8* null)
br label %10
10: ; preds = %8, %2
%11 = getelementptr inbounds i8*, i8** %1, i64 1
%12 = load i8*, i8** %11
%13 = call i64 @sscanf(i8* %12, i8* getelementptr inbounds ([3 x i8], [3 x i8]* @.str, i64 0, i64 0), i32* %3)
%14 = icmp eq i64 %13, 0
%15 = load i32, i32* %3
%16 = icmp eq i32 %15, 0
%17 = select i1 %14, i1 true, i1 %16
br i1 %17, label %18, label %20
18: ; preds = %10
%19 = call i32 @usage(i8* getelementptr inbounds ([6 x i8], [6 x i8]* @.str.1, i64 0, i64 0))
br label %20
20: ; preds = %10, %18
%21 = getelementptr inbounds i8*, i8** %1, i64 2
%22 = load i8*, i8** %21
%23 = call i64 @sscanf(i8* %22, i8* getelementptr inbounds ([3 x i8], [3 x i8]* @.str, i64 0, i64 0), i32* %4)
%24 = icmp eq i64 %23, 0
%25 = load i32, i32* %4
%26 = icmp eq i32 %25, 0
%27 = select i1 %24, i1 true, i1 %26
%28 = icmp ugt i32 %25, 31
%29 = select i1 %27, i1 true, i1 %28
br i1 %29, label %30, label %33
30: ; preds = %20
%31 = call i32 @usage(i8* getelementptr inbounds ([6 x i8], [6 x i8]* @.str.2, i64 0, i64 0))
%32 = load i32, i32* %4
br label %33
33: ; preds = %20, %30
%34 = phi i32 [ %25, %20 ], [ %32, %30 ]
%35 = load i32, i32* %3
%36 = lshr i32 %35, %34
%37 = icmp eq i32 %36, 0
br i1 %37, label %41, label %38
38: ; preds = %33
%39 = call i32 @usage(i8* getelementptr inbounds ([10 x i8], [10 x i8]* @.str.3, i64 0, i64 0))
%40 = load i32, i32* %3
br label %41
41: ; preds = %38, %33
%42 = phi i32 [ %40, %38 ], [ %35, %33 ]
%43 = shl i32 %42, 2
%44 = mul i32 %43, %42
%45 = call i64 @malloc(i32 %44)
%46 = inttoptr i64 %45 to i32*
%47 = load i32, i32* %3
%48 = load i32, i32* %4
%49 = mul i32 %48, %47
%50 = shl i32 %49, 2
%51 = mul i32 %50, %49
%52 = call i64 @malloc(i32 %51)
%53 = inttoptr i64 %52 to i32*
%54 = load i32, i32* %3
%55 = load i32, i32* %4
%56 = mul i32 %55, %54
%57 = shl i32 %56, 2
%58 = mul i32 %57, %56
%59 = call i64 @malloc(i32 %58)
%60 = inttoptr i64 %59 to i32*
%61 = load i32, i32* %3
%62 = icmp eq i32 %61, 0
br i1 %62, label %91, label %63
63: ; preds = %41, %88
%64 = phi i32 [ %86, %88 ], [ %61, %41 ]
%65 = phi i32 [ %89, %88 ], [ 0, %41 ]
br label %66
66: ; preds = %63, %78
%67 = phi i32 [ %86, %78 ], [ %64, %63 ]
%68 = phi i32 [ %85, %78 ], [ 0, %63 ]
%69 = load i32, i32* %4
%70 = shl nsw i32 -1, %69
%71 = xor i32 %70, -1
%72 = sub i32 %71, %68
%73 = xor i32 %72, %65
%74 = icmp eq i32 %73, 0
br i1 %74, label %78, label %75
75: ; preds = %66
%76 = call i32 @galois_single_divide(i32 1, i32 %73, i32 %69)
%77 = load i32, i32* %3
br label %78
78: ; preds = %66, %75
%79 = phi i32 [ %77, %75 ], [ %67, %66 ]
%80 = phi i32 [ %76, %75 ], [ 0, %66 ]
%81 = mul i32 %79, %65
%82 = add i32 %81, %68
%83 = zext i32 %82 to i64
%84 = getelementptr inbounds i32, i32* %46, i64 %83
store i32 %80, i32* %84
%85 = add nuw i32 %68, 1
%86 = load i32, i32* %3
%87 = icmp ult i32 %85, %86
br i1 %87, label %66, label %88
88: ; preds = %78
%89 = add nuw i32 %65, 1
%90 = icmp ult i32 %89, %86
br i1 %90, label %63, label %91
91: ; preds = %88, %41
%92 = phi i32 [ 0, %41 ], [ %86, %88 ]
%93 = load i32, i32* %4
%94 = call i32* @jerasure_matrix_to_bitmatrix(i32 %92, i32 %92, i32 %93, i32* %46)
%95 = call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([23 x i8], [23 x i8]* @str, i64 0, i64 0))
%96 = load i32, i32* %3
%97 = load i32, i32* %4
%98 = mul i32 %97, %96
%99 = call i32 @jerasure_print_bitmatrix(i32* %94, i32 %98, i32 %98, i32 %97)
%100 = load i32, i32* %3
%101 = load i32, i32* %4
%102 = mul i32 %101, %100
%103 = shl i32 %102, 2
%104 = mul i32 %103, %102
%105 = call i32 @memcpy(i32* %53, i32* %94, i32 %104)
%106 = load i32, i32* %3
%107 = load i32, i32* %4
%108 = mul i32 %107, %106
%109 = call i32 @jerasure_invertible_bitmatrix(i32* %53, i32 %108)
%110 = icmp eq i32 %109, 1
%111 = select i1 %110, i8* getelementptr inbounds ([4 x i8], [4 x i8]* @.str.6, i64 0, i64 0), i8* getelementptr inbounds ([3 x i8], [3 x i8]* @.str.7, i64 0, i64 0)
%112 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([17 x i8], [17 x i8]* @.str.5, i64 0, i64 0), i8* %111)
br i1 %110, label %113, label %138
113: ; preds = %91
%114 = call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([10 x i8], [10 x i8]* @str.10, i64 0, i64 0))
%115 = load i32, i32* %3
%116 = load i32, i32* %4
%117 = mul i32 %116, %115
%118 = shl i32 %117, 2
%119 = mul i32 %118, %117
%120 = call i32 @memcpy(i32* %53, i32* %94, i32 %119)
%121 = load i32, i32* %3
%122 = load i32, i32* %4
%123 = mul i32 %122, %121
%124 = call i32 @jerasure_invert_bitmatrix(i32* %53, i32* %60, i32 %123)
%125 = load i32, i32* %3
%126 = load i32, i32* %4
%127 = mul i32 %126, %125
%128 = call i32 @jerasure_print_bitmatrix(i32* %60, i32 %127, i32 %127, i32 %126)
%129 = load i32, i32* %3
%130 = load i32, i32* %4
%131 = mul i32 %130, %129
%132 = call i32* @jerasure_matrix_multiply(i32* %60, i32* %94, i32 %131, i32 %131, i32 %131, i32 %131, i32 2)
%133 = call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([44 x i8], [44 x i8]* @str.11, i64 0, i64 0))
%134 = load i32, i32* %3
%135 = load i32, i32* %4
%136 = mul i32 %135, %134
%137 = call i32 @jerasure_print_bitmatrix(i32* %132, i32 %136, i32 %136, i32 %135)
br label %138
138: ; preds = %113, %91
ret i32 0
}