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


x86 assembly (clang, -O0):
<remove_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
sub    -0x14(%rbp),%eax
mov    %eax,-0x24(%rbp)
mov    -0x10(%rbp),%eax
sub    -0x18(%rbp),%eax
mov    %eax,-0x28(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    11aa <remove_zero_pad+0xaa>
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    -0x28(%rbp),%eax
jge    1197 <remove_zero_pad+0x97>
mov    -0x8(%rbp),%rax
mov    -0x30(%rbp),%ecx
add    -0x38(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x30(%rbp),%ecx
add    -0x34(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1159 <remove_zero_pad+0x59>
jmp    119c <remove_zero_pad+0x9c>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1132 <remove_zero_pad+0x32>
pop    %rbp
ret

x86 assembly (clang, -O2):
<remove_zero_pad>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    %ecx,%esi
test   %esi,%esi
jle    1375 <remove_zero_pad+0x275>
mov    %rdi,%r15
mov    %r9,%rbp
movslq %edx,%rax
sub    %r8d,%edx
mov    %edx,-0x4c(%rsp)
movslq %edx,%rbx
mov    %esi,%edx
mov    %ebx,%r12d
mov    %r12d,%ecx
and    $0xfffffffc,%ecx
mov    %rcx,-0x18(%rsp)
add    $0xfffffffffffffffc,%rcx
mov    %rcx,-0x20(%rsp)
mov    %rcx,%rsi
shr    $0x2,%rsi
add    $0x1,%rsi
mov    %rsi,%rcx
mov    %rsi,-0x28(%rsp)
and    $0xfffffffffffffffe,%rsi
mov    %rsi,-0x30(%rsp)
mov    %r12d,%edi
and    $0x3,%edi
lea    0x30(%r15),%rcx
mov    %rax,-0x8(%rsp)
lea    0x0(,%rax,8),%r9
mov    %rbp,%rax
add    $0x30,%rbp
mov    %rbx,-0x10(%rsp)
lea    0x0(,%rbx,8),%r14
xor    %r11d,%r11d
mov    %r15,%r10
mov    %rax,%r13
mov    %r15,-0x40(%rsp)
mov    %rax,-0x48(%rsp)
mov    %r14,-0x38(%rsp)
jmp    11c9 <remove_zero_pad+0xc9>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
add    $0x1,%r11
add    %r9,%rcx
add    %r14,%rbp
add    %r14,%r13
add    %r9,%r10
cmp    %rdx,%r11
je     1375 <remove_zero_pad+0x275>
cmpl   $0x0,-0x4c(%rsp)
jle    11b0 <remove_zero_pad+0xb0>
cmpl   $0x4,-0x4c(%rsp)
jae    11e0 <remove_zero_pad+0xe0>
xor    %ebx,%ebx
jmp    12e0 <remove_zero_pad+0x1e0>
xchg   %ax,%ax
mov    %r9,%r14
mov    %r11,%r8
imul   -0x10(%rsp),%r8
mov    %r11,%rbx
imul   -0x8(%rsp),%rbx
mov    %r15,%r9
mov    %rax,%r15
lea    (%rbx,%r12,1),%rax
lea    (%r15,%r8,8),%rsi
lea    (%r9,%rax,8),%rax
cmp    %rax,%rsi
jae    1232 <remove_zero_pad+0x132>
lea    (%r8,%r12,1),%rax
lea    (%r15,%rax,8),%rax
lea    (%r9,%rbx,8),%rsi
cmp    %rax,%rsi
jae    1232 <remove_zero_pad+0x132>
xor    %ebx,%ebx
mov    %r15,%rax
mov    %r9,%r15
mov    %r14,%r9
mov    -0x38(%rsp),%r14
jmp    12e0 <remove_zero_pad+0x1e0>
cmpq   $0x0,-0x20(%rsp)
je     128b <remove_zero_pad+0x18b>
mov    -0x30(%rsp),%rsi
xor    %r15d,%r15d
mov    %r14,%r9
cs nopw 0x0(%rax,%rax,1)
nop
movups -0x30(%rcx,%r15,8),%xmm0
movups -0x20(%rcx,%r15,8),%xmm1
movups %xmm0,-0x30(%rbp,%r15,8)
movups %xmm1,-0x20(%rbp,%r15,8)
movups -0x10(%rcx,%r15,8),%xmm0
movups (%rcx,%r15,8),%xmm1
movups %xmm0,-0x10(%rbp,%r15,8)
movups %xmm1,0x0(%rbp,%r15,8)
add    $0x8,%r15
add    $0xfffffffffffffffe,%rsi
jne    1250 <remove_zero_pad+0x150>
jmp    1291 <remove_zero_pad+0x191>
xor    %r15d,%r15d
mov    %r14,%r9
testb  $0x1,-0x28(%rsp)
mov    -0x38(%rsp),%r14
je     12c1 <remove_zero_pad+0x1c1>
add    %r15,%rbx
mov    -0x40(%rsp),%rax
movups (%rax,%rbx,8),%xmm0
movups 0x10(%rax,%rbx,8),%xmm1
add    %r15,%r8
mov    -0x48(%rsp),%rax
movups %xmm0,(%rax,%r8,8)
movups %xmm1,0x10(%rax,%r8,8)
mov    -0x18(%rsp),%rax
mov    %rax,%rbx
cmp    %r12,%rax
mov    -0x48(%rsp),%rax
mov    -0x40(%rsp),%r15
je     11b0 <remove_zero_pad+0xb0>
nopl   0x0(%rax)
mov    %rax,%r8
mov    %rbx,%rsi
not    %rsi
add    %r12,%rsi
test   %rdi,%rdi
je     1317 <remove_zero_pad+0x217>
mov    %rdi,%rax
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movsd  (%r10,%rbx,8),%xmm0
movsd  %xmm0,0x0(%r13,%rbx,8)
add    $0x1,%rbx
add    $0xffffffffffffffff,%rax
jne    1300 <remove_zero_pad+0x200>
cmp    $0x3,%rsi
mov    %r8,%rax
jb     11b0 <remove_zero_pad+0xb0>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movsd  (%r10,%rbx,8),%xmm0
movsd  %xmm0,0x0(%r13,%rbx,8)
movsd  0x8(%r10,%rbx,8),%xmm0
movsd  %xmm0,0x8(%r13,%rbx,8)
movsd  0x10(%r10,%rbx,8),%xmm0
movsd  %xmm0,0x10(%r13,%rbx,8)
movsd  0x18(%r10,%rbx,8),%xmm0
movsd  %xmm0,0x18(%r13,%rbx,8)
add    $0x4,%rbx
cmp    %rbx,%r12
jne    1330 <remove_zero_pad+0x230>
jmp    11b0 <remove_zero_pad+0xb0>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @remove_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 = sub nsw i32 %19, %20
store i32 %21, i32* %13
%22 = load i32, i32* %9
%23 = load i32, i32* %11
%24 = sub nsw i32 %22, %23
store i32 %24, i32* %14
store i32 0, i32* %15
br label %25
25: ; preds = %58, %6
%26 = load i32, i32* %15
%27 = load i32, i32* %13
%28 = icmp slt i32 %26, %27
br i1 %28, label %29, label %61
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* %14
%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* %16
%43 = load i32, i32* %18
%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* %16
%50 = load i32, i32* %17
%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
br label %58
58: ; preds = %57
%59 = load i32, i32* %15
%60 = add nsw i32 %59, 1
store i32 %60, i32* %15
br label %25
61: ; preds = %25
ret void
}

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