void get_3x3_submatrix(float *dst, const struct matrix4 *m, int i, int j) {
  const float *mf = (const float *)m;
  int ti, tj, idst, jdst;
  for (ti = 0; ti < 4; ti++) {
    if (ti < i)
      idst = ti;
    else if (ti > i)
      idst = ti - 1;
    else
      continue;
    for (tj = 0; tj < 4; tj++) {
      if (tj < j)
        jdst = tj;
      else if (tj > j)
        jdst = tj - 1;
      else
        continue;
      dst[(idst * 3) + jdst] = mf[(ti * 4) + tj];
    }
  }
}


x86 assembly (clang, -O0):
<get_3x3_submatrix>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
mov    -0x10(%rbp),%rax
mov    %rax,-0x20(%rbp)
movl   $0x0,-0x24(%rbp)
cmpl   $0x4,-0x24(%rbp)
jge    11fb <get_3x3_submatrix+0xfb>
mov    -0x24(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1142 <get_3x3_submatrix+0x42>
mov    -0x24(%rbp),%eax
mov    %eax,-0x2c(%rbp)
jmp    1166 <get_3x3_submatrix+0x66>
mov    -0x24(%rbp),%eax
cmp    -0x14(%rbp),%eax
jle    115c <get_3x3_submatrix+0x5c>
mov    -0x24(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1161 <get_3x3_submatrix+0x61>
jmp    11ed <get_3x3_submatrix+0xed>
jmp    1166 <get_3x3_submatrix+0x66>
movl   $0x0,-0x28(%rbp)
cmpl   $0x4,-0x28(%rbp)
jge    11e8 <get_3x3_submatrix+0xe8>
mov    -0x28(%rbp),%eax
cmp    -0x18(%rbp),%eax
jge    118e <get_3x3_submatrix+0x8e>
mov    -0x28(%rbp),%eax
mov    %eax,-0x30(%rbp)
jmp    11b2 <get_3x3_submatrix+0xb2>
mov    -0x28(%rbp),%eax
cmp    -0x18(%rbp),%eax
jle    11a8 <get_3x3_submatrix+0xa8>
mov    -0x28(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    11ad <get_3x3_submatrix+0xad>
jmp    11da <get_3x3_submatrix+0xda>
jmp    11b2 <get_3x3_submatrix+0xb2>
mov    -0x20(%rbp),%rax
mov    -0x24(%rbp),%ecx
shl    $0x2,%ecx
add    -0x28(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
imul   $0x3,-0x2c(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    116d <get_3x3_submatrix+0x6d>
jmp    11ed <get_3x3_submatrix+0xed>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1121 <get_3x3_submatrix+0x21>
pop    %rbp
ret

x86 assembly (clang, -O2):
<get_3x3_submatrix>:
movslq %edx,%r8
xor    %edx,%edx
jmp    114f <get_3x3_submatrix+0x4f>
nopw   0x0(%rax,%rax,1)
js     11ad <get_3x3_submatrix+0xad>
movss  0x4(%rsi),%xmm0
movss  %xmm0,(%rdi,%r10,4)
movss  0x8(%rsi),%xmm0
movss  %xmm0,0x4(%rdi,%r10,4)
mov    $0x2,%eax
movss  0xc(%rsi),%xmm0
add    %r9d,%eax
cltq
movss  %xmm0,(%rdi,%rax,4)
add    $0x1,%rdx
add    $0x10,%rsi
cmp    $0x4,%rdx
je     11bd <get_3x3_submatrix+0xbd>
cmp    %r8,%rdx
jge    1160 <get_3x3_submatrix+0x60>
mov    %edx,%eax
jmp    1165 <get_3x3_submatrix+0x65>
nopl   0x0(%rax,%rax,1)
jle    1141 <get_3x3_submatrix+0x41>
lea    -0x1(%rdx),%eax
lea    (%rax,%rax,2),%r9d
movslq %r9d,%r10
test   %ecx,%ecx
jle    1110 <get_3x3_submatrix+0x10>
movss  (%rsi),%xmm0
movss  %xmm0,(%rdi,%r10,4)
cmp    $0x1,%ecx
jbe    1121 <get_3x3_submatrix+0x21>
movss  0x4(%rsi),%xmm0
movss  %xmm0,0x4(%rdi,%r10,4)
mov    $0x2,%eax
cmp    $0x2,%ecx
jbe    1132 <get_3x3_submatrix+0x32>
movss  0x8(%rsi),%xmm0
movss  %xmm0,0x8(%rdi,%r10,4)
mov    $0x3,%eax
cmp    $0x4,%ecx
jae    1132 <get_3x3_submatrix+0x32>
jmp    1141 <get_3x3_submatrix+0x41>
movss  (%rsi),%xmm0
movss  %xmm0,-0x4(%rdi,%r10,4)
jmp    1116 <get_3x3_submatrix+0x16>
ret

LLVM IR (-O0):
%struct.matrix4 = type opaque
define void @get_3x3_submatrix(float* %0, %struct.matrix4* %1, i32 %2, i32 %3) {
%5 = alloca float*
%6 = alloca %struct.matrix4*
%7 = alloca i32
%8 = alloca i32
%9 = alloca float*
%10 = alloca i32
%11 = alloca i32
%12 = alloca i32
%13 = alloca i32
store float* %0, float** %5
store %struct.matrix4* %1, %struct.matrix4** %6
store i32 %2, i32* %7
store i32 %3, i32* %8
%14 = load %struct.matrix4*, %struct.matrix4** %6
%15 = bitcast %struct.matrix4* %14 to float*
store float* %15, float** %9
store i32 0, i32* %10
br label %16
16: ; preds = %73, %4
%17 = load i32, i32* %10
%18 = icmp slt i32 %17, 4
br i1 %18, label %19, label %76
19: ; preds = %16
%20 = load i32, i32* %10
%21 = load i32, i32* %7
%22 = icmp slt i32 %20, %21
br i1 %22, label %23, label %25
23: ; preds = %19
%24 = load i32, i32* %10
store i32 %24, i32* %12
br label %34
25: ; preds = %19
%26 = load i32, i32* %10
%27 = load i32, i32* %7
%28 = icmp sgt i32 %26, %27
br i1 %28, label %29, label %32
29: ; preds = %25
%30 = load i32, i32* %10
%31 = sub nsw i32 %30, 1
store i32 %31, i32* %12
br label %33
32: ; preds = %25
br label %73
33: ; preds = %29
br label %34
34: ; preds = %33, %23
store i32 0, i32* %11
br label %35
35: ; preds = %69, %34
%36 = load i32, i32* %11
%37 = icmp slt i32 %36, 4
br i1 %37, label %38, label %72
38: ; preds = %35
%39 = load i32, i32* %11
%40 = load i32, i32* %8
%41 = icmp slt i32 %39, %40
br i1 %41, label %42, label %44
42: ; preds = %38
%43 = load i32, i32* %11
store i32 %43, i32* %13
br label %53
44: ; preds = %38
%45 = load i32, i32* %11
%46 = load i32, i32* %8
%47 = icmp sgt i32 %45, %46
br i1 %47, label %48, label %51
48: ; preds = %44
%49 = load i32, i32* %11
%50 = sub nsw i32 %49, 1
store i32 %50, i32* %13
br label %52
51: ; preds = %44
br label %69
52: ; preds = %48
br label %53
53: ; preds = %52, %42
%54 = load float*, float** %9
%55 = load i32, i32* %10
%56 = mul nsw i32 %55, 4
%57 = load i32, i32* %11
%58 = add nsw i32 %56, %57
%59 = sext i32 %58 to i64
%60 = getelementptr inbounds float, float* %54, i64 %59
%61 = load float, float* %60
%62 = load float*, float** %5
%63 = load i32, i32* %12
%64 = mul nsw i32 %63, 3
%65 = load i32, i32* %13
%66 = add nsw i32 %64, %65
%67 = sext i32 %66 to i64
%68 = getelementptr inbounds float, float* %62, i64 %67
store float %61, float* %68
br label %69
69: ; preds = %53, %51
%70 = load i32, i32* %11
%71 = add nsw i32 %70, 1
store i32 %71, i32* %11
br label %35
72: ; preds = %35
br label %73
73: ; preds = %72, %32
%74 = load i32, i32* %10
%75 = add nsw i32 %74, 1
store i32 %75, i32* %10
br label %16
76: ; preds = %16
ret void
}

LLVM IR (-O2):
%struct.matrix4 = type opaque
define void @get_3x3_submatrix(float* nocapture writeonly %0, %struct.matrix4* nocapture %1, i32 %2, i32 %3) {
%5 = bitcast %struct.matrix4* %1 to float*
%6 = sext i32 %2 to i64
%7 = icmp sgt i32 %3, 0
%8 = icmp slt i32 %3, 0
%9 = icmp ugt i32 %3, 1
%10 = icmp ugt i32 %3, 2
%11 = icmp ugt i32 %3, 3
br label %12
12: ; preds = %4, %71
%13 = phi i64 [ 0, %4 ], [ %72, %71 ]
%14 = trunc i64 %13 to i32
%15 = icmp slt i64 %13, %6
br i1 %15, label %20, label %16
16: ; preds = %12
%17 = icmp sgt i64 %13, %6
br i1 %17, label %18, label %71
18: ; preds = %16
%19 = add nsw i32 %14, -1
br label %20
20: ; preds = %12, %18
%21 = phi i32 [ %19, %18 ], [ %14, %12 ]
%22 = shl nsw i64 %13, 2
%23 = mul nsw i32 %21, 3
br i1 %7, label %31, label %24
24: ; preds = %20
br i1 %8, label %25, label %36
25: ; preds = %24
%26 = getelementptr inbounds float, float* %5, i64 %22
%27 = load float, float* %26
%28 = add nsw i32 %23, -1
%29 = sext i32 %28 to i64
%30 = getelementptr inbounds float, float* %0, i64 %29
store float %27, float* %30
br label %36
31: ; preds = %20
%32 = getelementptr inbounds float, float* %5, i64 %22
%33 = load float, float* %32
%34 = sext i32 %23 to i64
%35 = getelementptr inbounds float, float* %0, i64 %34
store float %33, float* %35
br i1 %9, label %42, label %49
36: ; preds = %25, %24
%37 = or i64 %22, 1
%38 = getelementptr inbounds float, float* %5, i64 %37
%39 = load float, float* %38
%40 = sext i32 %23 to i64
%41 = getelementptr inbounds float, float* %0, i64 %40
store float %39, float* %41
br label %49
42: ; preds = %31
%43 = or i64 %22, 1
%44 = getelementptr inbounds float, float* %5, i64 %43
%45 = load float, float* %44
%46 = add nsw i32 %23, 1
%47 = sext i32 %46 to i64
%48 = getelementptr inbounds float, float* %0, i64 %47
store float %45, float* %48
br i1 %10, label %56, label %63
49: ; preds = %36, %31
%50 = or i64 %22, 2
%51 = getelementptr inbounds float, float* %5, i64 %50
%52 = load float, float* %51
%53 = add nsw i32 %23, 1
%54 = sext i32 %53 to i64
%55 = getelementptr inbounds float, float* %0, i64 %54
store float %52, float* %55
br label %63
56: ; preds = %42
%57 = or i64 %22, 2
%58 = getelementptr inbounds float, float* %5, i64 %57
%59 = load float, float* %58
%60 = add nsw i32 %23, 2
%61 = sext i32 %60 to i64
%62 = getelementptr inbounds float, float* %0, i64 %61
store float %59, float* %62
br i1 %11, label %63, label %71
63: ; preds = %42, %49, %56
%64 = phi i32 [ 3, %56 ], [ 2, %49 ], [ 2, %42 ]
%65 = or i64 %22, 3
%66 = getelementptr inbounds float, float* %5, i64 %65
%67 = load float, float* %66
%68 = add nsw i32 %64, %23
%69 = sext i32 %68 to i64
%70 = getelementptr inbounds float, float* %0, i64 %69
store float %67, float* %70
br label %71
71: ; preds = %63, %56, %16
%72 = add nuw nsw i64 %13, 1
%73 = icmp eq i64 %72, 4
br i1 %73, label %74, label %12
74: ; preds = %71
ret void
}