void andc_smatrix(short **m1, long nrl, long nrh, long ncl, long nch, short c,
                  short **m2) {
  long i, j;
  short 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_smatrix>:
push   %rbp
mov    %rsp,%rbp
mov    %r9w,%ax
mov    0x10(%rbp),%r9
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %ax,-0x2a(%rbp)
mov    -0x10(%rbp),%rax
mov    %rax,-0x38(%rbp)
mov    -0x38(%rbp),%rax
cmp    -0x18(%rbp),%rax
jg     11c5 <andc_smatrix+0xc5>
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     11af <andc_smatrix+0xaf>
mov    -0x50(%rbp),%rax
mov    -0x40(%rbp),%rcx
mov    (%rax,%rcx,2),%ax
mov    %ax,-0x42(%rbp)
movswl -0x42(%rbp),%eax
movswl -0x2a(%rbp),%ecx
and    %ecx,%eax
mov    %ax,-0x44(%rbp)
mov    -0x44(%rbp),%dx
mov    -0x58(%rbp),%rax
mov    -0x40(%rbp),%rcx
mov    %dx,(%rax,%rcx,2)
mov    -0x40(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x40(%rbp)
jmp    1162 <andc_smatrix+0x62>
jmp    11b4 <andc_smatrix+0xb4>
mov    -0x38(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
jmp    112c <andc_smatrix+0x2c>
pop    %rbp
ret

x86 assembly (clang, -O2):
<andc_smatrix>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
cmp    %rdx,%rsi
jle    111a <andc_smatrix+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    -0x10(%rax),%rcx
mov    %rcx,-0x8(%rsp)
shr    $0x4,%rcx
add    $0x1,%rcx
mov    %rax,-0x38(%rsp)
and    $0xfffffffffffffff0,%rax
mov    %rax,-0x18(%rsp)
add    %r15,%rax
mov    %rax,-0x20(%rsp)
movd   %r9d,%xmm0
pshuflw $0x0,%xmm0,%xmm0
pshufd $0x0,%xmm0,%xmm0
mov    %rcx,%rax
mov    %rcx,-0x10(%rsp)
and    $0xfffffffffffffffe,%rcx
mov    %rcx,-0x28(%rsp)
lea    0x18(%r15),%rax
mov    %rax,-0x30(%rsp)
jmp    118c <andc_smatrix+0x8c>
xchg   %ax,%ax
lea    0x1(%rsi),%rax
cmp    %rdx,%rsi
mov    %rax,%rsi
je     110f <andc_smatrix+0xf>
cmp    %r8,%r15
jg     1180 <andc_smatrix+0x80>
mov    (%rdi,%rsi,8),%r12
mov    (%r14,%rsi,8),%r13
mov    %r15,%rcx
cmpq   $0x10,-0x38(%rsp)
jb     12b0 <andc_smatrix+0x1b0>
lea    0x0(,%r15,2),%rax
add    %r13,%rax
lea    (%r12,%r11,2),%rcx
cmp    %rcx,%rax
jae    11d7 <andc_smatrix+0xd7>
lea    0x0(,%r11,2),%rax
add    %r13,%rax
lea    (%r12,%r15,2),%rbp
mov    %r15,%rcx
cmp    %rax,%rbp
jb     12b0 <andc_smatrix+0x1b0>
mov    %r14,%rax
cmpq   $0x10,-0x8(%rsp)
jae    11e6 <andc_smatrix+0xe6>
xor    %ecx,%ecx
jmp    1260 <andc_smatrix+0x160>
mov    -0x30(%rsp),%rcx
lea    0x0(,%rcx,2),%r14
add    %r13,%r14
lea    (%r12,%rcx,2),%r10
mov    -0x28(%rsp),%rbx
xor    %ecx,%ecx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movdqu -0x30(%r10,%rcx,2),%xmm1
movdqu -0x20(%r10,%rcx,2),%xmm2
pand   %xmm0,%xmm1
pand   %xmm0,%xmm2
movdqu %xmm1,-0x30(%r14,%rcx,2)
movdqu %xmm2,-0x20(%r14,%rcx,2)
movdqu -0x10(%r10,%rcx,2),%xmm1
movdqu (%r10,%rcx,2),%xmm2
pand   %xmm0,%xmm1
pand   %xmm0,%xmm2
movdqu %xmm1,-0x10(%r14,%rcx,2)
movdqu %xmm2,(%r14,%rcx,2)
add    $0x20,%rcx
add    $0xfffffffffffffffe,%rbx
jne    1210 <andc_smatrix+0x110>
testb  $0x1,-0x10(%rsp)
mov    %rax,%r14
je     1290 <andc_smatrix+0x190>
add    %r15,%rcx
movdqu (%r12,%rcx,2),%xmm1
movdqu 0x10(%r12,%rcx,2),%xmm2
pand   %xmm0,%xmm1
pand   %xmm0,%xmm2
movdqu %xmm1,0x0(%r13,%rcx,2)
movdqu %xmm2,0x10(%r13,%rcx,2)
mov    -0x20(%rsp),%rcx
mov    -0x18(%rsp),%rax
cmp    %rax,-0x38(%rsp)
je     1180 <andc_smatrix+0x80>
cs nopw 0x0(%rax,%rax,1)
nop
mov    %r11d,%eax
sub    %ecx,%eax
mov    %r8,%rbx
sub    %rcx,%rbx
and    $0x3,%rax
je     12e9 <andc_smatrix+0x1e9>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movzwl (%r12,%rcx,2),%ebp
and    %r9w,%bp
mov    %bp,0x0(%r13,%rcx,2)
add    $0x1,%rcx
add    $0xffffffffffffffff,%rax
jne    12d0 <andc_smatrix+0x1d0>
cmp    $0x3,%rbx
jb     1180 <andc_smatrix+0x80>
mov    %r8,%r10
sub    %rcx,%r10
lea    0x6(,%rcx,2),%rax
add    %r13,%rax
lea    (%r12,%rcx,2),%rcx
add    $0x6,%rcx
mov    $0xffffffffffffffff,%rbp
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movzwl -0x4(%rcx,%rbp,2),%ebx
and    %r9w,%bx
mov    %bx,-0x4(%rax,%rbp,2)
movzwl -0x2(%rcx,%rbp,2),%ebx
and    %r9w,%bx
mov    %bx,-0x2(%rax,%rbp,2)
movzwl (%rcx,%rbp,2),%ebx
and    %r9w,%bx
mov    %bx,(%rax,%rbp,2)
movzwl 0x2(%rcx,%rbp,2),%ebx
and    %r9w,%bx
mov    %bx,0x2(%rax,%rbp,2)
add    $0x4,%rbp
cmp    %rbp,%r10
jne    1320 <andc_smatrix+0x220>
jmp    1180 <andc_smatrix+0x80>

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

LLVM IR (-O2):
define void @andc_smatrix(i16** nocapture %0, i64 %1, i64 %2, i64 %3, i64 %4, i16 signext %5, i16** 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, -16
%14 = lshr i64 %13, 4
%15 = add nuw nsw i64 %14, 1
%16 = add i64 %4, 1
%17 = icmp ult i64 %12, 16
%18 = and i64 %12, -16
%19 = add i64 %18, %3
%20 = insertelement <8 x i16> poison, i16 %5, i64 0
%21 = shufflevector <8 x i16> %20, <8 x i16> poison, <8 x i32> zeroinitializer
%22 = insertelement <8 x i16> poison, i16 %5, i64 0
%23 = shufflevector <8 x i16> %22, <8 x i16> poison, <8 x i32> zeroinitializer
%24 = and i64 %15, 1
%25 = icmp ult i64 %13, 16
%26 = and i64 %15, 2305843009213693950
%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 i16*, i16** %0, i64 %30
%32 = load i16*, i16** %31
%33 = getelementptr inbounds i16*, i16** %6, i64 %30
%34 = load i16*, i16** %33
br i1 %10, label %137, label %35
35: ; preds = %29
br i1 %17, label %95, label %36
36: ; preds = %35
%37 = getelementptr i16, i16* %34, i64 %3
%38 = getelementptr i16, i16* %34, i64 %11
%39 = getelementptr i16, i16* %32, i64 %3
%40 = getelementptr i16, i16* %32, i64 %11
%41 = icmp ult i16* %37, %40
%42 = icmp ult i16* %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 i16, i16* %32, i64 %48
%50 = bitcast i16* %49 to <8 x i16>*
%51 = load <8 x i16>, <8 x i16>* %50
%52 = getelementptr inbounds i16, i16* %49, i64 8
%53 = bitcast i16* %52 to <8 x i16>*
%54 = load <8 x i16>, <8 x i16>* %53
%55 = and <8 x i16> %51, %21
%56 = and <8 x i16> %54, %23
%57 = getelementptr inbounds i16, i16* %34, i64 %48
%58 = bitcast i16* %57 to <8 x i16>*
store <8 x i16> %55, <8 x i16>* %58
%59 = getelementptr inbounds i16, i16* %57, i64 8
%60 = bitcast i16* %59 to <8 x i16>*
store <8 x i16> %56, <8 x i16>* %60
%61 = or i64 %46, 16
%62 = add i64 %61, %3
%63 = getelementptr inbounds i16, i16* %32, i64 %62
%64 = bitcast i16* %63 to <8 x i16>*
%65 = load <8 x i16>, <8 x i16>* %64
%66 = getelementptr inbounds i16, i16* %63, i64 8
%67 = bitcast i16* %66 to <8 x i16>*
%68 = load <8 x i16>, <8 x i16>* %67
%69 = and <8 x i16> %65, %21
%70 = and <8 x i16> %68, %23
%71 = getelementptr inbounds i16, i16* %34, i64 %62
%72 = bitcast i16* %71 to <8 x i16>*
store <8 x i16> %69, <8 x i16>* %72
%73 = getelementptr inbounds i16, i16* %71, i64 8
%74 = bitcast i16* %73 to <8 x i16>*
store <8 x i16> %70, <8 x i16>* %74
%75 = add nuw i64 %46, 32
%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 i16, i16* %32, i64 %81
%83 = bitcast i16* %82 to <8 x i16>*
%84 = load <8 x i16>, <8 x i16>* %83
%85 = getelementptr inbounds i16, i16* %82, i64 8
%86 = bitcast i16* %85 to <8 x i16>*
%87 = load <8 x i16>, <8 x i16>* %86
%88 = and <8 x i16> %84, %21
%89 = and <8 x i16> %87, %23
%90 = getelementptr inbounds i16, i16* %34, i64 %81
%91 = bitcast i16* %90 to <8 x i16>*
store <8 x i16> %88, <8 x i16>* %91
%92 = getelementptr inbounds i16, i16* %90, i64 8
%93 = bitcast i16* %92 to <8 x i16>*
store <8 x i16> %89, <8 x i16>* %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 i16, i16* %32, i64 %102
%105 = load i16, i16* %104
%106 = and i16 %105, %5
%107 = getelementptr inbounds i16, i16* %34, i64 %102
store i16 %106, i16* %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 i16, i16* %32, i64 %115
%117 = load i16, i16* %116
%118 = and i16 %117, %5
%119 = getelementptr inbounds i16, i16* %34, i64 %115
store i16 %118, i16* %119
%120 = add i64 %115, 1
%121 = getelementptr inbounds i16, i16* %32, i64 %120
%122 = load i16, i16* %121
%123 = and i16 %122, %5
%124 = getelementptr inbounds i16, i16* %34, i64 %120
store i16 %123, i16* %124
%125 = add i64 %115, 2
%126 = getelementptr inbounds i16, i16* %32, i64 %125
%127 = load i16, i16* %126
%128 = and i16 %127, %5
%129 = getelementptr inbounds i16, i16* %34, i64 %125
store i16 %128, i16* %129
%130 = add i64 %115, 3
%131 = getelementptr inbounds i16, i16* %32, i64 %130
%132 = load i16, i16* %131
%133 = and i16 %132, %5
%134 = getelementptr inbounds i16, i16* %34, i64 %130
store i16 %133, i16* %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
}