void andc_imatrix(int **m1, long nrl, long nrh, long ncl, long nch, int c,
                  int **m2) {
  long i, j;
  int x, y, *Xi, *Yi;
  for (i = nrl; i <= nrh; i++) {
    Xi = m1[i];
    Yi = m2[i];
    for (j = ncl; j <= nch; j++) {
      x = Xi[j];
      y = x & c;
      Yi[j] = y;
    }
  }
}


x86 assembly (clang, -O0):
<andc_imatrix>:
push   %rbp
mov    %rsp,%rbp
mov    0x10(%rbp),%rax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9d,-0x2c(%rbp)
mov    -0x10(%rbp),%rax
mov    %rax,-0x38(%rbp)
mov    -0x38(%rbp),%rax
cmp    -0x18(%rbp),%rax
jg     11b8 <andc_imatrix+0xb8>
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    %rax,-0x50(%rbp)
mov    0x10(%rbp),%rax
mov    -0x38(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    %rax,-0x58(%rbp)
mov    -0x20(%rbp),%rax
mov    %rax,-0x40(%rbp)
mov    -0x40(%rbp),%rax
cmp    -0x28(%rbp),%rax
jg     11a2 <andc_imatrix+0xa2>
mov    -0x50(%rbp),%rax
mov    -0x40(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x44(%rbp)
mov    -0x44(%rbp),%eax
and    -0x2c(%rbp),%eax
mov    %eax,-0x48(%rbp)
mov    -0x48(%rbp),%edx
mov    -0x58(%rbp),%rax
mov    -0x40(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x40(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x40(%rbp)
jmp    115e <andc_imatrix+0x5e>
jmp    11a7 <andc_imatrix+0xa7>
mov    -0x38(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
jmp    1128 <andc_imatrix+0x28>
pop    %rbp
ret

x86 assembly (clang, -O2):
<andc_imatrix>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
cmp    %rdx,%rsi
jle    111a <andc_imatrix+0x1a>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
mov    %rcx,%r15
mov    0x38(%rsp),%r14
lea    0x1(%r8),%r11
mov    %r11,%rax
sub    %rcx,%rax
lea    -0x8(%rax),%rcx
mov    %rcx,-0x8(%rsp)
shr    $0x3,%rcx
add    $0x1,%rcx
mov    %rax,-0x38(%rsp)
and    $0xfffffffffffffff8,%rax
mov    %rax,-0x18(%rsp)
add    %r15,%rax
mov    %rax,-0x20(%rsp)
movd   %r9d,%xmm0
pshufd $0x0,%xmm0,%xmm0
mov    %rcx,%rax
mov    %rcx,-0x10(%rsp)
and    $0xfffffffffffffffe,%rcx
mov    %rcx,-0x28(%rsp)
lea    0xc(%r15),%rax
mov    %rax,-0x30(%rsp)
jmp    118c <andc_imatrix+0x8c>
nopl   0x0(%rax)
lea    0x1(%rsi),%rax
cmp    %rdx,%rsi
mov    %rax,%rsi
je     110f <andc_imatrix+0xf>
cmp    %r8,%r15
jg     1180 <andc_imatrix+0x80>
mov    (%rdi,%rsi,8),%r12
mov    (%r14,%rsi,8),%r13
mov    %r15,%rcx
cmpq   $0x8,-0x38(%rsp)
jb     12b0 <andc_imatrix+0x1b0>
lea    0x0(,%r15,4),%rax
add    %r13,%rax
lea    (%r12,%r11,4),%rcx
cmp    %rcx,%rax
jae    11d7 <andc_imatrix+0xd7>
lea    0x0(,%r11,4),%rax
add    %r13,%rax
lea    (%r12,%r15,4),%rbp
mov    %r15,%rcx
cmp    %rax,%rbp
jb     12b0 <andc_imatrix+0x1b0>
mov    %r14,%rax
cmpq   $0x8,-0x8(%rsp)
jae    11e6 <andc_imatrix+0xe6>
xor    %ecx,%ecx
jmp    1260 <andc_imatrix+0x160>
mov    -0x30(%rsp),%rcx
lea    0x0(,%rcx,4),%r14
add    %r13,%r14
lea    (%r12,%rcx,4),%r10
mov    -0x28(%rsp),%rbx
xor    %ecx,%ecx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movdqu -0x30(%r10,%rcx,4),%xmm1
movdqu -0x20(%r10,%rcx,4),%xmm2
pand   %xmm0,%xmm1
pand   %xmm0,%xmm2
movdqu %xmm1,-0x30(%r14,%rcx,4)
movdqu %xmm2,-0x20(%r14,%rcx,4)
movdqu -0x10(%r10,%rcx,4),%xmm1
movdqu (%r10,%rcx,4),%xmm2
pand   %xmm0,%xmm1
pand   %xmm0,%xmm2
movdqu %xmm1,-0x10(%r14,%rcx,4)
movdqu %xmm2,(%r14,%rcx,4)
add    $0x10,%rcx
add    $0xfffffffffffffffe,%rbx
jne    1210 <andc_imatrix+0x110>
testb  $0x1,-0x10(%rsp)
mov    %rax,%r14
je     1290 <andc_imatrix+0x190>
add    %r15,%rcx
movdqu (%r12,%rcx,4),%xmm1
movdqu 0x10(%r12,%rcx,4),%xmm2
pand   %xmm0,%xmm1
pand   %xmm0,%xmm2
movdqu %xmm1,0x0(%r13,%rcx,4)
movdqu %xmm2,0x10(%r13,%rcx,4)
mov    -0x20(%rsp),%rcx
mov    -0x18(%rsp),%rax
cmp    %rax,-0x38(%rsp)
je     1180 <andc_imatrix+0x80>
cs nopw 0x0(%rax,%rax,1)
nop
mov    %r11d,%eax
sub    %ecx,%eax
mov    %r8,%rbx
sub    %rcx,%rbx
and    $0x3,%rax
je     12e6 <andc_imatrix+0x1e6>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    (%r12,%rcx,4),%ebp
and    %r9d,%ebp
mov    %ebp,0x0(%r13,%rcx,4)
add    $0x1,%rcx
add    $0xffffffffffffffff,%rax
jne    12d0 <andc_imatrix+0x1d0>
cmp    $0x3,%rbx
jb     1180 <andc_imatrix+0x80>
mov    %r8,%r10
sub    %rcx,%r10
lea    0xc(,%rcx,4),%rax
add    %r13,%rax
lea    (%r12,%rcx,4),%rcx
add    $0xc,%rcx
mov    $0xffffffffffffffff,%rbp
mov    -0x8(%rcx,%rbp,4),%ebx
and    %r9d,%ebx
mov    %ebx,-0x8(%rax,%rbp,4)
mov    -0x4(%rcx,%rbp,4),%ebx
and    %r9d,%ebx
mov    %ebx,-0x4(%rax,%rbp,4)
mov    (%rcx,%rbp,4),%ebx
and    %r9d,%ebx
mov    %ebx,(%rax,%rbp,4)
mov    0x4(%rcx,%rbp,4),%ebx
and    %r9d,%ebx
mov    %ebx,0x4(%rax,%rbp,4)
add    $0x4,%rbp
cmp    %rbp,%r10
jne    1310 <andc_imatrix+0x210>
jmp    1180 <andc_imatrix+0x80>

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

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