void add_zero_pad(double *X, int rows, int cols, int zrow, int zcol,
                  double *Y) {
  int r, c, i, j, u, v;
  r = rows + zrow;
  c = cols + zcol;
  for (i = 0; i < rows; ++i) {
    u = i * c;
    v = i * cols;
    for (j = 0; j < cols; ++j) {
      Y[u + j] = X[v + j];
    }
    for (j = cols; j < c; ++j) {
      Y[u + j] = 0.;
    }
  }
  for (i = rows; i < r; ++i) {
    u = i * c;
    for (j = 0; j < c; ++j) {
      Y[u + j] = 0.;
    }
  }
}


x86 assembly (clang, -O0):
<add_zero_pad>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
mov    %ecx,-0x14(%rbp)
mov    %r8d,-0x18(%rbp)
mov    %r9,-0x20(%rbp)
mov    -0xc(%rbp),%eax
add    -0x14(%rbp),%eax
mov    %eax,-0x24(%rbp)
mov    -0x10(%rbp),%eax
add    -0x18(%rbp),%eax
mov    %eax,-0x28(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    11df <add_zero_pad+0xdf>
mov    -0x2c(%rbp),%eax
imul   -0x28(%rbp),%eax
mov    %eax,-0x34(%rbp)
mov    -0x2c(%rbp),%eax
imul   -0x10(%rbp),%eax
mov    %eax,-0x38(%rbp)
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    1197 <add_zero_pad+0x97>
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x34(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1159 <add_zero_pad+0x59>
mov    -0x10(%rbp),%eax
mov    %eax,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x28(%rbp),%eax
jge    11cc <add_zero_pad+0xcc>
mov    -0x20(%rbp),%rax
mov    -0x34(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    119d <add_zero_pad+0x9d>
jmp    11d1 <add_zero_pad+0xd1>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1132 <add_zero_pad+0x32>
mov    -0xc(%rbp),%eax
mov    %eax,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    1244 <add_zero_pad+0x144>
mov    -0x2c(%rbp),%eax
imul   -0x28(%rbp),%eax
mov    %eax,-0x34(%rbp)
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x28(%rbp),%eax
jge    1231 <add_zero_pad+0x131>
mov    -0x20(%rbp),%rax
mov    -0x34(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1202 <add_zero_pad+0x102>
jmp    1236 <add_zero_pad+0x136>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    11e5 <add_zero_pad+0xe5>
pop    %rbp
ret

x86 assembly (clang, -O2):
<add_zero_pad>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x98,%rsp
mov    %ecx,0x1c(%rsp)
mov    %edx,%r15d
mov    %r8,%rax
mov    %r8,0x88(%rsp)
lea    (%r8,%r15,1),%eax
mov    %rax,0x10(%rsp)
mov    %esi,0xc(%rsp)
test   %esi,%esi
mov    %r9,0x20(%rsp)
jle    13ef <add_zero_pad+0x2df>
lea    0x1(%r15),%eax
mov    0x10(%rsp),%rdx
cmp    %eax,%edx
cmovg  %edx,%eax
mov    %r15d,%ecx
not    %ecx
add    %eax,%ecx
lea    0x8(,%rcx,8),%rax
mov    %rax,0x60(%rsp)
movslq %edx,%rax
movslq %r15d,%rcx
mov    0xc(%rsp),%edx
mov    %rdx,0x80(%rsp)
mov    %ecx,%ebx
mov    %ebx,%edx
and    $0xfffffffc,%edx
mov    %rdx,0x40(%rsp)
add    $0xfffffffffffffffc,%rdx
mov    %rdx,0x38(%rsp)
mov    %rdx,%rsi
shr    $0x2,%rsi
add    $0x1,%rsi
mov    %rsi,%rdx
mov    %rsi,0x30(%rsp)
and    $0xfffffffffffffffe,%rsi
mov    %rsi,0x28(%rsp)
mov    %ebx,%r8d
and    $0x3,%r8d
lea    0x30(%rdi),%r14
mov    %rcx,0x48(%rsp)
lea    0x0(,%rcx,8),%rcx
mov    %rcx,0x78(%rsp)
lea    0x30(%r9),%r12
mov    %rax,0x50(%rsp)
lea    0x0(,%rax,8),%rax
mov    %rax,0x70(%rsp)
xor    %r11d,%r11d
mov    %rdi,%r13
mov    %r9,%rbp
mov    %rdi,0x68(%rsp)
mov    %r8,0x58(%rsp)
jmp    1228 <add_zero_pad+0x118>
nopw   0x0(%rax,%rax,1)
add    $0x1,%r11
mov    0x78(%rsp),%rax
add    %rax,%r14
mov    0x70(%rsp),%rcx
add    %rcx,%r12
add    %rcx,%rbp
add    %rax,%r13
cmp    0x80(%rsp),%r11
je     13ef <add_zero_pad+0x2df>
test   %r15d,%r15d
jle    1390 <add_zero_pad+0x280>
cmp    $0x4,%r15d
jb     1273 <add_zero_pad+0x163>
mov    %r11,%rax
imul   0x50(%rsp),%rax
mov    %r11,%rcx
imul   0x48(%rsp),%rcx
lea    (%rcx,%rbx,1),%rdx
lea    (%r9,%rax,8),%rsi
lea    (%rdi,%rdx,8),%rdx
cmp    %rdx,%rsi
jae    1302 <add_zero_pad+0x1f2>
lea    (%rax,%rbx,1),%rdx
lea    (%r9,%rdx,8),%rdx
lea    (%rdi,%rcx,8),%rsi
cmp    %rdx,%rsi
jae    1302 <add_zero_pad+0x1f2>
xor    %eax,%eax
mov    %rax,%rcx
not    %rcx
add    %rbx,%rcx
test   %r8,%r8
je     12a7 <add_zero_pad+0x197>
mov    %r8,%rdx
cs nopw 0x0(%rax,%rax,1)
movsd  0x0(%r13,%rax,8),%xmm0
movsd  %xmm0,0x0(%rbp,%rax,8)
add    $0x1,%rax
add    $0xffffffffffffffff,%rdx
jne    1290 <add_zero_pad+0x180>
cmp    $0x3,%rcx
jb     1390 <add_zero_pad+0x280>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movsd  0x0(%r13,%rax,8),%xmm0
movsd  %xmm0,0x0(%rbp,%rax,8)
movsd  0x8(%r13,%rax,8),%xmm0
movsd  %xmm0,0x8(%rbp,%rax,8)
movsd  0x10(%r13,%rax,8),%xmm0
movsd  %xmm0,0x10(%rbp,%rax,8)
movsd  0x18(%r13,%rax,8),%xmm0
movsd  %xmm0,0x18(%rbp,%rax,8)
add    $0x4,%rax
cmp    %rax,%rbx
jne    12c0 <add_zero_pad+0x1b0>
jmp    1390 <add_zero_pad+0x280>
cmpq   $0x0,0x38(%rsp)
je     13e0 <add_zero_pad+0x2d0>
mov    0x28(%rsp),%rsi
xor    %edx,%edx
cs nopw 0x0(%rax,%rax,1)
nop
movups -0x30(%r14,%rdx,8),%xmm0
movups -0x20(%r14,%rdx,8),%xmm1
movups %xmm0,-0x30(%r12,%rdx,8)
movups %xmm1,-0x20(%r12,%rdx,8)
movups -0x10(%r14,%rdx,8),%xmm0
movups (%r14,%rdx,8),%xmm1
movups %xmm0,-0x10(%r12,%rdx,8)
movups %xmm1,(%r12,%rdx,8)
add    $0x8,%rdx
add    $0xfffffffffffffffe,%rsi
jne    1320 <add_zero_pad+0x210>
testb  $0x1,0x30(%rsp)
je     1379 <add_zero_pad+0x269>
add    %rdx,%rcx
movups (%rdi,%rcx,8),%xmm0
movups 0x10(%rdi,%rcx,8),%xmm1
add    %rdx,%rax
movups %xmm0,(%r9,%rax,8)
movups %xmm1,0x10(%r9,%rax,8)
mov    0x40(%rsp),%rcx
mov    %rcx,%rax
cmp    %rbx,%rcx
jne    1275 <add_zero_pad+0x165>
nopw   0x0(%rax,%rax,1)
cmpl   $0x0,0x88(%rsp)
jle    1200 <add_zero_pad+0xf0>
mov    0x10(%rsp),%rax
imul   %r11d,%eax
add    %r15d,%eax
cltq
lea    (%r9,%rax,8),%rdi
xor    %esi,%esi
mov    0x60(%rsp),%rdx
mov    %r11,0x90(%rsp)
call   1030 <memset@plt>
mov    0x90(%rsp),%r11
mov    0x58(%rsp),%r8
mov    0x68(%rsp),%rdi
mov    0x20(%rsp),%r9
jmp    1200 <add_zero_pad+0xf0>
xor    %edx,%edx
testb  $0x1,0x30(%rsp)
jne    135f <add_zero_pad+0x24f>
jmp    1379 <add_zero_pad+0x269>
mov    0x1c(%rsp),%r15d
test   %r15d,%r15d
mov    0xc(%rsp),%ebx
mov    0x10(%rsp),%r12
jle    1448 <add_zero_pad+0x338>
add    %ebx,%r15d
lea    -0x1(%r12),%eax
lea    0x8(,%rax,8),%r14
mov    %ebx,%ebp
imul   %r12d,%ebp
jmp    142b <add_zero_pad+0x31b>
nopw   0x0(%rax,%rax,1)
add    $0x1,%ebx
add    %r12d,%ebp
cmp    %r15d,%ebx
jge    1448 <add_zero_pad+0x338>
test   %r12d,%r12d
jle    1420 <add_zero_pad+0x310>
movslq %ebp,%rax
lea    (%r9,%rax,8),%rdi
xor    %esi,%esi
mov    %r14,%rdx
call   1030 <memset@plt>
mov    0x20(%rsp),%r9
jmp    1420 <add_zero_pad+0x310>
add    $0x98,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @add_zero_pad(double* %0, i32 %1, i32 %2, i32 %3, i32 %4, double* %5) {
%7 = alloca double*
%8 = alloca i32
%9 = alloca i32
%10 = alloca i32
%11 = alloca i32
%12 = alloca double*
%13 = alloca i32
%14 = alloca i32
%15 = alloca i32
%16 = alloca i32
%17 = alloca i32
%18 = alloca i32
store double* %0, double** %7
store i32 %1, i32* %8
store i32 %2, i32* %9
store i32 %3, i32* %10
store i32 %4, i32* %11
store double* %5, double** %12
%19 = load i32, i32* %8
%20 = load i32, i32* %10
%21 = add nsw i32 %19, %20
store i32 %21, i32* %13
%22 = load i32, i32* %9
%23 = load i32, i32* %11
%24 = add nsw i32 %22, %23
store i32 %24, i32* %14
store i32 0, i32* %15
br label %25
25: ; preds = %74, %6
%26 = load i32, i32* %15
%27 = load i32, i32* %8
%28 = icmp slt i32 %26, %27
br i1 %28, label %29, label %77
29: ; preds = %25
%30 = load i32, i32* %15
%31 = load i32, i32* %14
%32 = mul nsw i32 %30, %31
store i32 %32, i32* %17
%33 = load i32, i32* %15
%34 = load i32, i32* %9
%35 = mul nsw i32 %33, %34
store i32 %35, i32* %18
store i32 0, i32* %16
br label %36
36: ; preds = %54, %29
%37 = load i32, i32* %16
%38 = load i32, i32* %9
%39 = icmp slt i32 %37, %38
br i1 %39, label %40, label %57
40: ; preds = %36
%41 = load double*, double** %7
%42 = load i32, i32* %18
%43 = load i32, i32* %16
%44 = add nsw i32 %42, %43
%45 = sext i32 %44 to i64
%46 = getelementptr inbounds double, double* %41, i64 %45
%47 = load double, double* %46
%48 = load double*, double** %12
%49 = load i32, i32* %17
%50 = load i32, i32* %16
%51 = add nsw i32 %49, %50
%52 = sext i32 %51 to i64
%53 = getelementptr inbounds double, double* %48, i64 %52
store double %47, double* %53
br label %54
54: ; preds = %40
%55 = load i32, i32* %16
%56 = add nsw i32 %55, 1
store i32 %56, i32* %16
br label %36
57: ; preds = %36
%58 = load i32, i32* %9
store i32 %58, i32* %16
br label %59
59: ; preds = %70, %57
%60 = load i32, i32* %16
%61 = load i32, i32* %14
%62 = icmp slt i32 %60, %61
br i1 %62, label %63, label %73
63: ; preds = %59
%64 = load double*, double** %12
%65 = load i32, i32* %17
%66 = load i32, i32* %16
%67 = add nsw i32 %65, %66
%68 = sext i32 %67 to i64
%69 = getelementptr inbounds double, double* %64, i64 %68
store double 0.000000e+00, double* %69
br label %70
70: ; preds = %63
%71 = load i32, i32* %16
%72 = add nsw i32 %71, 1
store i32 %72, i32* %16
br label %59
73: ; preds = %59
br label %74
74: ; preds = %73
%75 = load i32, i32* %15
%76 = add nsw i32 %75, 1
store i32 %76, i32* %15
br label %25
77: ; preds = %25
%78 = load i32, i32* %8
store i32 %78, i32* %15
br label %79
79: ; preds = %102, %77
%80 = load i32, i32* %15
%81 = load i32, i32* %13
%82 = icmp slt i32 %80, %81
br i1 %82, label %83, label %105
83: ; preds = %79
%84 = load i32, i32* %15
%85 = load i32, i32* %14
%86 = mul nsw i32 %84, %85
store i32 %86, i32* %17
store i32 0, i32* %16
br label %87
87: ; preds = %98, %83
%88 = load i32, i32* %16
%89 = load i32, i32* %14
%90 = icmp slt i32 %88, %89
br i1 %90, label %91, label %101
91: ; preds = %87
%92 = load double*, double** %12
%93 = load i32, i32* %17
%94 = load i32, i32* %16
%95 = add nsw i32 %93, %94
%96 = sext i32 %95 to i64
%97 = getelementptr inbounds double, double* %92, i64 %96
store double 0.000000e+00, double* %97
br label %98
98: ; preds = %91
%99 = load i32, i32* %16
%100 = add nsw i32 %99, 1
store i32 %100, i32* %16
br label %87
101: ; preds = %87
br label %102
102: ; preds = %101
%103 = load i32, i32* %15
%104 = add nsw i32 %103, 1
store i32 %104, i32* %15
br label %79
105: ; preds = %79
ret void
}

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