void addandc_smatrix(short **src, long nrl, long nrh, long ncl, long nch,
                     short cte, short **dst) {
  long i, j;
  short *Xi, *Yi;
  for (i = nrl; i <= nrh; i++) {
    Xi = src[i];
    Yi = dst[i];
    for (j = ncl; j <= nch; j++) {
      if (Xi[j])
        Yi[j] = Xi[j] + cte;
    }
  }
}


x86 assembly (clang, -O0):
<addandc_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     11d0 <addandc_smatrix+0xd0>
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    %rax,-0x48(%rbp)
mov    0x10(%rbp),%rax
mov    -0x38(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    %rax,-0x50(%rbp)
mov    -0x20(%rbp),%rax
mov    %rax,-0x40(%rbp)
mov    -0x40(%rbp),%rax
cmp    -0x28(%rbp),%rax
jg     11ba <addandc_smatrix+0xba>
mov    -0x48(%rbp),%rax
mov    -0x40(%rbp),%rcx
cmpw   $0x0,(%rax,%rcx,2)
je     11a4 <addandc_smatrix+0xa4>
mov    -0x48(%rbp),%rax
mov    -0x40(%rbp),%rcx
movswl (%rax,%rcx,2),%eax
movswl -0x2a(%rbp),%ecx
add    %ecx,%eax
mov    %ax,%dx
mov    -0x50(%rbp),%rax
mov    -0x40(%rbp),%rcx
mov    %dx,(%rax,%rcx,2)
jmp    11a9 <addandc_smatrix+0xa9>
mov    -0x40(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x40(%rbp)
jmp    1162 <addandc_smatrix+0x62>
jmp    11bf <addandc_smatrix+0xbf>
mov    -0x38(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
jmp    112c <addandc_smatrix+0x2c>
pop    %rbp
ret

x86 assembly (clang, -O2):
<addandc_smatrix>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
cmp    %rdx,%rsi
jle    111a <addandc_smatrix+0x1a>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
mov    0x38(%rsp),%r15
lea    0x1(%r8),%r13
mov    %r13,%r12
sub    %rcx,%r12
mov    %r12,-0x18(%rsp)
and    $0xfffffffffffffff8,%r12
lea    (%r12,%rcx,1),%rax
mov    %rax,-0x8(%rsp)
lea    0x7(%rcx),%rax
mov    %rax,-0x10(%rsp)
pxor   %xmm0,%xmm0
pcmpeqd %xmm1,%xmm1
jmp    115c <addandc_smatrix+0x5c>
xchg   %ax,%ax
lea    0x1(%rsi),%rax
cmp    %rdx,%rsi
mov    %rax,%rsi
je     110f <addandc_smatrix+0xf>
cmp    %r8,%rcx
jg     1150 <addandc_smatrix+0x50>
mov    (%rdi,%rsi,8),%rbx
mov    (%r15,%rsi,8),%r11
mov    %rcx,%r14
cmpq   $0x8,-0x18(%rsp)
jb     12b0 <addandc_smatrix+0x1b0>
lea    (%r11,%rcx,2),%r10
lea    (%rbx,%rcx,2),%rbp
lea    (%rbx,%r13,2),%rax
cmp    %rax,%r10
jae    1199 <addandc_smatrix+0x99>
lea    (%r11,%r13,2),%rax
mov    %rcx,%r14
cmp    %rax,%rbp
jb     12b0 <addandc_smatrix+0x1b0>
mov    -0x10(%rsp),%rax
lea    (%r11,%rax,2),%r15
xor    %r14d,%r14d
jmp    11bd <addandc_smatrix+0xbd>
nopw   0x0(%rax,%rax,1)
add    $0x8,%r14
cmp    %r14,%r12
je     1295 <addandc_smatrix+0x195>
movdqu 0x0(%rbp,%r14,2),%xmm2
movdqa %xmm2,%xmm3
pcmpeqw %xmm0,%xmm3
movd   %xmm3,%r10d
not    %r10d
test   $0x1,%r10b
je     11e7 <addandc_smatrix+0xe7>
movd   %xmm2,%eax
add    %r9d,%eax
mov    %ax,-0xe(%r15,%r14,2)
shr    $0x10,%r10d
test   $0x1,%r10b
je     11ff <addandc_smatrix+0xff>
pextrw $0x1,%xmm2,%eax
add    %r9d,%eax
mov    %ax,-0xc(%r15,%r14,2)
pxor   %xmm1,%xmm3
pextrw $0x2,%xmm3,%eax
test   $0x1,%al
je     121a <addandc_smatrix+0x11a>
pextrw $0x2,%xmm2,%eax
add    %r9d,%eax
mov    %ax,-0xa(%r15,%r14,2)
pextrw $0x3,%xmm3,%eax
test   $0x1,%al
je     1231 <addandc_smatrix+0x131>
pextrw $0x3,%xmm2,%eax
add    %r9d,%eax
mov    %ax,-0x8(%r15,%r14,2)
pextrw $0x4,%xmm3,%eax
test   $0x1,%al
je     1248 <addandc_smatrix+0x148>
pextrw $0x4,%xmm2,%eax
add    %r9d,%eax
mov    %ax,-0x6(%r15,%r14,2)
pextrw $0x5,%xmm3,%eax
test   $0x1,%al
je     125f <addandc_smatrix+0x15f>
pextrw $0x5,%xmm2,%eax
add    %r9d,%eax
mov    %ax,-0x4(%r15,%r14,2)
pextrw $0x6,%xmm3,%eax
test   $0x1,%al
je     1276 <addandc_smatrix+0x176>
pextrw $0x6,%xmm2,%eax
add    %r9d,%eax
mov    %ax,-0x2(%r15,%r14,2)
pextrw $0x7,%xmm3,%eax
test   $0x1,%al
je     11b0 <addandc_smatrix+0xb0>
pextrw $0x7,%xmm2,%eax
add    %r9d,%eax
mov    %ax,(%r15,%r14,2)
jmp    11b0 <addandc_smatrix+0xb0>
mov    -0x8(%rsp),%r14
cmp    %r12,-0x18(%rsp)
mov    0x38(%rsp),%r15
je     1150 <addandc_smatrix+0x50>
nopw   0x0(%rax,%rax,1)
mov    %r13d,%eax
sub    %r14d,%eax
mov    %r14,%rbp
test   $0x1,%al
je     12d3 <addandc_smatrix+0x1d3>
movzwl (%rbx,%r14,2),%ebp
test   %bp,%bp
je     12cf <addandc_smatrix+0x1cf>
add    %r9d,%ebp
mov    %bp,(%r11,%r14,2)
lea    0x1(%r14),%rbp
cmp    %r8,%r14
jne    12ed <addandc_smatrix+0x1ed>
jmp    1150 <addandc_smatrix+0x50>
nopl   (%rax)
add    $0x2,%rbp
cmp    %rbp,%r13
je     1150 <addandc_smatrix+0x50>
movzwl (%rbx,%rbp,2),%eax
test   %ax,%ax
je     12fe <addandc_smatrix+0x1fe>
add    %r9d,%eax
mov    %ax,(%r11,%rbp,2)
movzwl 0x2(%rbx,%rbp,2),%eax
test   %ax,%ax
je     12e0 <addandc_smatrix+0x1e0>
add    %r9d,%eax
mov    %ax,0x2(%r11,%rbp,2)
jmp    12e0 <addandc_smatrix+0x1e0>

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

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