void extendMatrix(double **a, int m, int n, double *row, double *column) {
  int i;
  for (i = 0; i < n + 1; i++) {
    a[m][i] = row[i];
  }
  for (i = 0; i < m; i++) {
    a[i][n] = column[i];
  }
}


x86 assembly (clang, -O0):
<extendMatrix>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
mov    %rcx,-0x18(%rbp)
mov    %r8,-0x20(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
mov    -0x10(%rbp),%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jge    115e <extendMatrix+0x5e>
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
movslq -0xc(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    111d <extendMatrix+0x1d>
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    11a1 <extendMatrix+0xa1>
mov    -0x20(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x24(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x10(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1165 <extendMatrix+0x65>
pop    %rbp
ret

x86 assembly (clang, -O2):
<extendMatrix>:
push   %r15
push   %r14
push   %rbx
test   %edx,%edx
js     1288 <extendMatrix+0x188>
movslq %esi,%rax
mov    (%rdi,%rax,8),%r10
lea    0x1(%rdx),%r9d
cmp    $0x3,%edx
jb     1137 <extendMatrix+0x37>
lea    (%rcx,%r9,8),%rax
cmp    %rax,%r10
jae    11b0 <extendMatrix+0xb0>
lea    (%r10,%r9,8),%rax
cmp    %rcx,%rax
jbe    11b0 <extendMatrix+0xb0>
xor    %eax,%eax
mov    %rax,%r11
not    %r11
add    %r9,%r11
mov    %r9,%rbx
and    $0x3,%rbx
je     1165 <extendMatrix+0x65>
nopl   0x0(%rax,%rax,1)
movsd  (%rcx,%rax,8),%xmm0
movsd  %xmm0,(%r10,%rax,8)
add    $0x1,%rax
add    $0xffffffffffffffff,%rbx
jne    1150 <extendMatrix+0x50>
cmp    $0x3,%r11
jb     1288 <extendMatrix+0x188>
nop
movsd  (%rcx,%rax,8),%xmm0
movsd  %xmm0,(%r10,%rax,8)
movsd  0x8(%rcx,%rax,8),%xmm0
movsd  %xmm0,0x8(%r10,%rax,8)
movsd  0x10(%rcx,%rax,8),%xmm0
movsd  %xmm0,0x10(%r10,%rax,8)
movsd  0x18(%rcx,%rax,8),%xmm0
movsd  %xmm0,0x18(%r10,%rax,8)
add    $0x4,%rax
cmp    %rax,%r9
jne    1170 <extendMatrix+0x70>
jmp    1288 <extendMatrix+0x188>
mov    %r9d,%eax
and    $0xfffffffc,%eax
lea    -0x4(%rax),%rbx
mov    %rbx,%r14
shr    $0x2,%r14
add    $0x1,%r14
mov    %r14d,%r11d
and    $0x3,%r11d
cmp    $0xc,%rbx
jae    11d6 <extendMatrix+0xd6>
xor    %ebx,%ebx
jmp    1240 <extendMatrix+0x140>
and    $0xfffffffffffffffc,%r14
xor    %ebx,%ebx
nopl   0x0(%rax)
movups (%rcx,%rbx,8),%xmm0
movups 0x10(%rcx,%rbx,8),%xmm1
movups %xmm0,(%r10,%rbx,8)
movups %xmm1,0x10(%r10,%rbx,8)
movups 0x20(%rcx,%rbx,8),%xmm0
movups 0x30(%rcx,%rbx,8),%xmm1
movups %xmm0,0x20(%r10,%rbx,8)
movups %xmm1,0x30(%r10,%rbx,8)
movups 0x40(%rcx,%rbx,8),%xmm0
movups 0x50(%rcx,%rbx,8),%xmm1
movups %xmm0,0x40(%r10,%rbx,8)
movups %xmm1,0x50(%r10,%rbx,8)
movups 0x60(%rcx,%rbx,8),%xmm0
movups 0x70(%rcx,%rbx,8),%xmm1
movups %xmm0,0x60(%r10,%rbx,8)
movups %xmm1,0x70(%r10,%rbx,8)
add    $0x10,%rbx
add    $0xfffffffffffffffc,%r14
jne    11e0 <extendMatrix+0xe0>
test   %r11,%r11
je     127f <extendMatrix+0x17f>
lea    (%r10,%rbx,8),%r14
add    $0x10,%r14
lea    (%rcx,%rbx,8),%r15
add    $0x10,%r15
shl    $0x5,%r11
xor    %ebx,%ebx
nopl   0x0(%rax,%rax,1)
movups -0x10(%r15,%rbx,1),%xmm0
movups (%r15,%rbx,1),%xmm1
movups %xmm0,-0x10(%r14,%rbx,1)
movups %xmm1,(%r14,%rbx,1)
add    $0x20,%rbx
cmp    %rbx,%r11
jne    1260 <extendMatrix+0x160>
cmp    %r9,%rax
jne    1139 <extendMatrix+0x39>
test   %esi,%esi
jle    1327 <extendMatrix+0x227>
movslq %edx,%rax
mov    %esi,%edx
lea    -0x1(%rdx),%rsi
mov    %edx,%ecx
and    $0x3,%ecx
cmp    $0x3,%rsi
jae    12a8 <extendMatrix+0x1a8>
xor    %esi,%esi
jmp    12fb <extendMatrix+0x1fb>
and    $0xfffffffc,%edx
xor    %esi,%esi
nopl   (%rax)
movsd  (%r8,%rsi,8),%xmm0
mov    (%rdi,%rsi,8),%rbx
movsd  %xmm0,(%rbx,%rax,8)
movsd  0x8(%r8,%rsi,8),%xmm0
mov    0x8(%rdi,%rsi,8),%rbx
movsd  %xmm0,(%rbx,%rax,8)
movsd  0x10(%r8,%rsi,8),%xmm0
mov    0x10(%rdi,%rsi,8),%rbx
movsd  %xmm0,(%rbx,%rax,8)
movsd  0x18(%r8,%rsi,8),%xmm0
mov    0x18(%rdi,%rsi,8),%rbx
movsd  %xmm0,(%rbx,%rax,8)
add    $0x4,%rsi
cmp    %rsi,%rdx
jne    12b0 <extendMatrix+0x1b0>
test   %rcx,%rcx
je     1327 <extendMatrix+0x227>
lea    (%rdi,%rsi,8),%rdx
lea    (%r8,%rsi,8),%rsi
xor    %edi,%edi
nopw   0x0(%rax,%rax,1)
movsd  (%rsi,%rdi,8),%xmm0
mov    (%rdx,%rdi,8),%rbx
movsd  %xmm0,(%rbx,%rax,8)
add    $0x1,%rdi
cmp    %rdi,%rcx
jne    1310 <extendMatrix+0x210>
pop    %rbx
pop    %r14
pop    %r15
ret

LLVM IR (-O0):
define void @extendMatrix(double** %0, i32 %1, i32 %2, double* %3, double* %4) {
%6 = alloca double**
%7 = alloca i32
%8 = alloca i32
%9 = alloca double*
%10 = alloca double*
%11 = alloca i32
store double** %0, double*** %6
store i32 %1, i32* %7
store i32 %2, i32* %8
store double* %3, double** %9
store double* %4, double** %10
store i32 0, i32* %11
br label %12
12: ; preds = %31, %5
%13 = load i32, i32* %11
%14 = load i32, i32* %8
%15 = add nsw i32 %14, 1
%16 = icmp slt i32 %13, %15
br i1 %16, label %17, label %34
17: ; preds = %12
%18 = load double*, double** %9
%19 = load i32, i32* %11
%20 = sext i32 %19 to i64
%21 = getelementptr inbounds double, double* %18, i64 %20
%22 = load double, double* %21
%23 = load double**, double*** %6
%24 = load i32, i32* %7
%25 = sext i32 %24 to i64
%26 = getelementptr inbounds double*, double** %23, i64 %25
%27 = load double*, double** %26
%28 = load i32, i32* %11
%29 = sext i32 %28 to i64
%30 = getelementptr inbounds double, double* %27, i64 %29
store double %22, double* %30
br label %31
31: ; preds = %17
%32 = load i32, i32* %11
%33 = add nsw i32 %32, 1
store i32 %33, i32* %11
br label %12
34: ; preds = %12
store i32 0, i32* %11
br label %35
35: ; preds = %53, %34
%36 = load i32, i32* %11
%37 = load i32, i32* %7
%38 = icmp slt i32 %36, %37
br i1 %38, label %39, label %56
39: ; preds = %35
%40 = load double*, double** %10
%41 = load i32, i32* %11
%42 = sext i32 %41 to i64
%43 = getelementptr inbounds double, double* %40, i64 %42
%44 = load double, double* %43
%45 = load double**, double*** %6
%46 = load i32, i32* %11
%47 = sext i32 %46 to i64
%48 = getelementptr inbounds double*, double** %45, i64 %47
%49 = load double*, double** %48
%50 = load i32, i32* %8
%51 = sext i32 %50 to i64
%52 = getelementptr inbounds double, double* %49, i64 %51
store double %44, double* %52
br label %53
53: ; preds = %39
%54 = load i32, i32* %11
%55 = add nsw i32 %54, 1
store i32 %55, i32* %11
br label %35
56: ; preds = %35
ret void
}

LLVM IR (-O2):
define void @extendMatrix(double** nocapture %0, i32 %1, i32 %2, double* nocapture %3, double* nocapture %4) {
%6 = icmp slt i32 %2, 0
br i1 %6, label %117, label %7
7: ; preds = %5
%8 = sext i32 %1 to i64
%9 = getelementptr inbounds double*, double** %0, i64 %8
%10 = load double*, double** %9
%11 = add nuw i32 %2, 1
%12 = zext i32 %11 to i64
%13 = icmp ult i32 %2, 3
br i1 %13, label %99, label %14
14: ; preds = %7
%15 = getelementptr double, double* %10, i64 %12
%16 = getelementptr double, double* %3, i64 %12
%17 = icmp ult double* %10, %16
%18 = icmp ugt double* %15, %3
%19 = and i1 %17, %18
br i1 %19, label %99, label %20
20: ; preds = %14
%21 = and i64 %12, 4294967292
%22 = add nsw i64 %21, -4
%23 = lshr exact i64 %22, 2
%24 = add nuw nsw i64 %23, 1
%25 = and i64 %24, 3
%26 = icmp ult i64 %22, 12
br i1 %26, label %78, label %27
27: ; preds = %20
%28 = and i64 %24, 9223372036854775804
br label %29
29: ; preds = %29, %27
%30 = phi i64 [ 0, %27 ], [ %75, %29 ]
%31 = phi i64 [ 0, %27 ], [ %76, %29 ]
%32 = getelementptr inbounds double, double* %3, i64 %30
%33 = bitcast double* %32 to <2 x double>*
%34 = load <2 x double>, <2 x double>* %33
%35 = getelementptr inbounds double, double* %32, i64 2
%36 = bitcast double* %35 to <2 x double>*
%37 = load <2 x double>, <2 x double>* %36
%38 = getelementptr inbounds double, double* %10, i64 %30
%39 = bitcast double* %38 to <2 x double>*
store <2 x double> %34, <2 x double>* %39
%40 = getelementptr inbounds double, double* %38, i64 2
%41 = bitcast double* %40 to <2 x double>*
store <2 x double> %37, <2 x double>* %41
%42 = or i64 %30, 4
%43 = getelementptr inbounds double, double* %3, i64 %42
%44 = bitcast double* %43 to <2 x double>*
%45 = load <2 x double>, <2 x double>* %44
%46 = getelementptr inbounds double, double* %43, i64 2
%47 = bitcast double* %46 to <2 x double>*
%48 = load <2 x double>, <2 x double>* %47
%49 = getelementptr inbounds double, double* %10, i64 %42
%50 = bitcast double* %49 to <2 x double>*
store <2 x double> %45, <2 x double>* %50
%51 = getelementptr inbounds double, double* %49, i64 2
%52 = bitcast double* %51 to <2 x double>*
store <2 x double> %48, <2 x double>* %52
%53 = or i64 %30, 8
%54 = getelementptr inbounds double, double* %3, i64 %53
%55 = bitcast double* %54 to <2 x double>*
%56 = load <2 x double>, <2 x double>* %55
%57 = getelementptr inbounds double, double* %54, i64 2
%58 = bitcast double* %57 to <2 x double>*
%59 = load <2 x double>, <2 x double>* %58
%60 = getelementptr inbounds double, double* %10, i64 %53
%61 = bitcast double* %60 to <2 x double>*
store <2 x double> %56, <2 x double>* %61
%62 = getelementptr inbounds double, double* %60, i64 2
%63 = bitcast double* %62 to <2 x double>*
store <2 x double> %59, <2 x double>* %63
%64 = or i64 %30, 12
%65 = getelementptr inbounds double, double* %3, i64 %64
%66 = bitcast double* %65 to <2 x double>*
%67 = load <2 x double>, <2 x double>* %66
%68 = getelementptr inbounds double, double* %65, i64 2
%69 = bitcast double* %68 to <2 x double>*
%70 = load <2 x double>, <2 x double>* %69
%71 = getelementptr inbounds double, double* %10, i64 %64
%72 = bitcast double* %71 to <2 x double>*
store <2 x double> %67, <2 x double>* %72
%73 = getelementptr inbounds double, double* %71, i64 2
%74 = bitcast double* %73 to <2 x double>*
store <2 x double> %70, <2 x double>* %74
%75 = add nuw i64 %30, 16
%76 = add nuw i64 %31, 4
%77 = icmp eq i64 %76, %28
br i1 %77, label %78, label %29
78: ; preds = %29, %20
%79 = phi i64 [ 0, %20 ], [ %75, %29 ]
%80 = icmp eq i64 %25, 0
br i1 %80, label %97, label %81
81: ; preds = %78, %81
%82 = phi i64 [ %94, %81 ], [ %79, %78 ]
%83 = phi i64 [ %95, %81 ], [ 0, %78 ]
%84 = getelementptr inbounds double, double* %3, i64 %82
%85 = bitcast double* %84 to <2 x double>*
%86 = load <2 x double>, <2 x double>* %85
%87 = getelementptr inbounds double, double* %84, i64 2
%88 = bitcast double* %87 to <2 x double>*
%89 = load <2 x double>, <2 x double>* %88
%90 = getelementptr inbounds double, double* %10, i64 %82
%91 = bitcast double* %90 to <2 x double>*
store <2 x double> %86, <2 x double>* %91
%92 = getelementptr inbounds double, double* %90, i64 2
%93 = bitcast double* %92 to <2 x double>*
store <2 x double> %89, <2 x double>* %93
%94 = add nuw i64 %82, 4
%95 = add i64 %83, 1
%96 = icmp eq i64 %95, %25
br i1 %96, label %97, label %81
97: ; preds = %81, %78
%98 = icmp eq i64 %21, %12
br i1 %98, label %117, label %99
99: ; preds = %14, %7, %97
%100 = phi i64 [ 0, %14 ], [ 0, %7 ], [ %21, %97 ]
%101 = xor i64 %100, -1
%102 = add nsw i64 %101, %12
%103 = and i64 %12, 3
%104 = icmp eq i64 %103, 0
br i1 %104, label %114, label %105
105: ; preds = %99, %105
%106 = phi i64 [ %111, %105 ], [ %100, %99 ]
%107 = phi i64 [ %112, %105 ], [ 0, %99 ]
%108 = getelementptr inbounds double, double* %3, i64 %106
%109 = load double, double* %108
%110 = getelementptr inbounds double, double* %10, i64 %106
store double %109, double* %110
%111 = add nuw nsw i64 %106, 1
%112 = add i64 %107, 1
%113 = icmp eq i64 %112, %103
br i1 %113, label %114, label %105
114: ; preds = %105, %99
%115 = phi i64 [ %100, %99 ], [ %111, %105 ]
%116 = icmp ult i64 %102, 3
br i1 %116, label %117, label %127
117: ; preds = %114, %127, %97, %5
%118 = icmp sgt i32 %1, 0
br i1 %118, label %119, label %189
119: ; preds = %117
%120 = sext i32 %2 to i64
%121 = zext i32 %1 to i64
%122 = add nsw i64 %121, -1
%123 = and i64 %121, 3
%124 = icmp ult i64 %122, 3
br i1 %124, label %175, label %125
125: ; preds = %119
%126 = and i64 %121, 4294967292
br label %146
127: ; preds = %114, %127
%128 = phi i64 [ %144, %127 ], [ %115, %114 ]
%129 = getelementptr inbounds double, double* %3, i64 %128
%130 = load double, double* %129
%131 = getelementptr inbounds double, double* %10, i64 %128
store double %130, double* %131
%132 = add nuw nsw i64 %128, 1
%133 = getelementptr inbounds double, double* %3, i64 %132
%134 = load double, double* %133
%135 = getelementptr inbounds double, double* %10, i64 %132
store double %134, double* %135
%136 = add nuw nsw i64 %128, 2
%137 = getelementptr inbounds double, double* %3, i64 %136
%138 = load double, double* %137
%139 = getelementptr inbounds double, double* %10, i64 %136
store double %138, double* %139
%140 = add nuw nsw i64 %128, 3
%141 = getelementptr inbounds double, double* %3, i64 %140
%142 = load double, double* %141
%143 = getelementptr inbounds double, double* %10, i64 %140
store double %142, double* %143
%144 = add nuw nsw i64 %128, 4
%145 = icmp eq i64 %144, %12
br i1 %145, label %117, label %127
146: ; preds = %146, %125
%147 = phi i64 [ 0, %125 ], [ %172, %146 ]
%148 = phi i64 [ 0, %125 ], [ %173, %146 ]
%149 = getelementptr inbounds double, double* %4, i64 %147
%150 = load double, double* %149
%151 = getelementptr inbounds double*, double** %0, i64 %147
%152 = load double*, double** %151
%153 = getelementptr inbounds double, double* %152, i64 %120
store double %150, double* %153
%154 = or i64 %147, 1
%155 = getelementptr inbounds double, double* %4, i64 %154
%156 = load double, double* %155
%157 = getelementptr inbounds double*, double** %0, i64 %154
%158 = load double*, double** %157
%159 = getelementptr inbounds double, double* %158, i64 %120
store double %156, double* %159
%160 = or i64 %147, 2
%161 = getelementptr inbounds double, double* %4, i64 %160
%162 = load double, double* %161
%163 = getelementptr inbounds double*, double** %0, i64 %160
%164 = load double*, double** %163
%165 = getelementptr inbounds double, double* %164, i64 %120
store double %162, double* %165
%166 = or i64 %147, 3
%167 = getelementptr inbounds double, double* %4, i64 %166
%168 = load double, double* %167
%169 = getelementptr inbounds double*, double** %0, i64 %166
%170 = load double*, double** %169
%171 = getelementptr inbounds double, double* %170, i64 %120
store double %168, double* %171
%172 = add nuw nsw i64 %147, 4
%173 = add i64 %148, 4
%174 = icmp eq i64 %173, %126
br i1 %174, label %175, label %146
175: ; preds = %146, %119
%176 = phi i64 [ 0, %119 ], [ %172, %146 ]
%177 = icmp eq i64 %123, 0
br i1 %177, label %189, label %178
178: ; preds = %175, %178
%179 = phi i64 [ %186, %178 ], [ %176, %175 ]
%180 = phi i64 [ %187, %178 ], [ 0, %175 ]
%181 = getelementptr inbounds double, double* %4, i64 %179
%182 = load double, double* %181
%183 = getelementptr inbounds double*, double** %0, i64 %179
%184 = load double*, double** %183
%185 = getelementptr inbounds double, double* %184, i64 %120
store double %182, double* %185
%186 = add nuw nsw i64 %179, 1
%187 = add i64 %180, 1
%188 = icmp eq i64 %187, %123
br i1 %188, label %189, label %178
189: ; preds = %175, %178, %117
ret void
}