void linear_img(short int *src0, short int *src1, short int *dst, int c0,
                int c1, int c2, int ov_flag, int npix) {
  int i, res;
  for (i = 0; i < npix; ++i) {
    res = (*src0) * c0 + (*src1) * c1 + c2;
    if (res < 0) {
      if (ov_flag)
        res = 0;
      else
        res = (((res) > 0) ? ((res)) : (-(res)));
    }
    if (res > 4096 - 1)
      res = 4096 - 1;
    *dst = res;
    dst++;
    src0++;
    src1++;
  }
}


x86 assembly (clang, -O0):
<linear_img>:
push   %rbp
mov    %rsp,%rbp
mov    0x18(%rbp),%eax
mov    0x10(%rbp),%eax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    %r8d,-0x20(%rbp)
mov    %r9d,-0x24(%rbp)
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    0x18(%rbp),%eax
jge    11ed <linear_img+0xed>
mov    -0x8(%rbp),%rax
movswl (%rax),%eax
imul   -0x1c(%rbp),%eax
mov    -0x10(%rbp),%rcx
movswl (%rcx),%ecx
imul   -0x20(%rbp),%ecx
add    %ecx,%eax
add    -0x24(%rbp),%eax
mov    %eax,-0x2c(%rbp)
cmpl   $0x0,-0x2c(%rbp)
jge    119a <linear_img+0x9a>
cmpl   $0x0,0x10(%rbp)
je     1172 <linear_img+0x72>
movl   $0x0,-0x2c(%rbp)
jmp    1195 <linear_img+0x95>
cmpl   $0x0,-0x2c(%rbp)
jle    1187 <linear_img+0x87>
mov    -0x2c(%rbp),%eax
mov    %eax,-0x30(%rbp)
jmp    118f <linear_img+0x8f>
xor    %eax,%eax
sub    -0x2c(%rbp),%eax
mov    %eax,-0x30(%rbp)
mov    -0x30(%rbp),%eax
mov    %eax,-0x2c(%rbp)
jmp    119a <linear_img+0x9a>
cmpl   $0xfff,-0x2c(%rbp)
jle    11ae <linear_img+0xae>
movl   $0xfff,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
mov    %ax,%cx
mov    -0x18(%rbp),%rax
mov    %cx,(%rax)
mov    -0x18(%rbp),%rax
add    $0x2,%rax
mov    %rax,-0x18(%rbp)
mov    -0x8(%rbp),%rax
add    $0x2,%rax
mov    %rax,-0x8(%rbp)
mov    -0x10(%rbp),%rax
add    $0x2,%rax
mov    %rax,-0x10(%rbp)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1128 <linear_img+0x28>
pop    %rbp
ret

x86 assembly (clang, -O2):
<linear_img>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    0x40(%rsp),%r14d
test   %r14d,%r14d
jle    1396 <linear_img+0x296>
mov    0x38(%rsp),%r10d
lea    -0x1(%r14),%r13d
xor    %r11d,%r11d
cmp    $0x7,%r13d
jb     1332 <linear_img+0x232>
lea    (%rdi,%r13,2),%rax
add    $0x2,%rax
cmp    %rdx,%rax
seta   %r15b
lea    (%rdx,%r13,2),%rbp
add    $0x2,%rbp
cmp    %rdi,%rbp
seta   %bl
lea    (%rsi,%r13,2),%rax
add    $0x2,%rax
cmp    %rdx,%rax
seta   %al
cmp    %rsi,%rbp
seta   %bpl
test   %bl,%r15b
jne    1332 <linear_img+0x232>
and    %bpl,%al
jne    1332 <linear_img+0x232>
add    $0x1,%r13
mov    %r13,%r12
and    $0xfffffffffffffff8,%r12
lea    (%rdi,%r12,2),%r15
lea    (%rsi,%r12,2),%rax
lea    (%rdx,%r12,2),%rbx
movd   %ecx,%xmm0
pshufd $0x0,%xmm0,%xmm8
movd   %r8d,%xmm0
pshufd $0x0,%xmm0,%xmm10
movd   %r9d,%xmm0
pshufd $0x0,%xmm0,%xmm9
xor    %ebp,%ebp
pshufd $0xf5,%xmm8,%xmm11
pshufd $0xf5,%xmm10,%xmm4
movdqa 0xe41(%rip),%xmm5
jmp    124e <linear_img+0x14e>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
pxor   %xmm2,%xmm2
pxor   %xmm1,%xmm1
pxor   %xmm3,%xmm3
pcmpgtd %xmm7,%xmm3
pcmpgtd %xmm6,%xmm0
pand   %xmm0,%xmm1
pandn  %xmm6,%xmm0
por    %xmm1,%xmm0
pand   %xmm3,%xmm2
pandn  %xmm7,%xmm3
por    %xmm2,%xmm3
movdqa %xmm5,%xmm1
pcmpgtd %xmm3,%xmm1
pand   %xmm1,%xmm3
pandn  %xmm5,%xmm1
por    %xmm3,%xmm1
movdqa %xmm5,%xmm2
pcmpgtd %xmm0,%xmm2
pand   %xmm2,%xmm0
pandn  %xmm5,%xmm2
por    %xmm0,%xmm2
pslld  $0x10,%xmm2
psrad  $0x10,%xmm2
pslld  $0x10,%xmm1
psrad  $0x10,%xmm1
packssdw %xmm2,%xmm1
movdqu %xmm1,(%rdx,%rbp,2)
add    $0x8,%rbp
cmp    %rbp,%r12
je     1391 <linear_img+0x291>
movdqu (%rdi,%rbp,2),%xmm0
punpcklwd %xmm0,%xmm1
psrad  $0x10,%xmm1
punpckhwd %xmm0,%xmm0
psrad  $0x10,%xmm0
movdqa %xmm8,%xmm2
pmuludq %xmm0,%xmm2
pshufd $0xe8,%xmm2,%xmm6
pshufd $0xf5,%xmm0,%xmm0
pmuludq %xmm11,%xmm0
pshufd $0xe8,%xmm0,%xmm0
punpckldq %xmm0,%xmm6
movdqa %xmm8,%xmm0
pmuludq %xmm1,%xmm0
pshufd $0xe8,%xmm0,%xmm7
pshufd $0xf5,%xmm1,%xmm0
pmuludq %xmm11,%xmm0
pshufd $0xe8,%xmm0,%xmm0
punpckldq %xmm0,%xmm7
movdqu (%rsi,%rbp,2),%xmm0
punpckhwd %xmm0,%xmm1
psrad  $0x10,%xmm1
punpcklwd %xmm0,%xmm0
psrad  $0x10,%xmm0
movdqa %xmm10,%xmm2
pmuludq %xmm0,%xmm2
pshufd $0xe8,%xmm2,%xmm2
pshufd $0xf5,%xmm0,%xmm0
pmuludq %xmm4,%xmm0
pshufd $0xe8,%xmm0,%xmm0
punpckldq %xmm0,%xmm2
movdqa %xmm10,%xmm0
pmuludq %xmm1,%xmm0
pshufd $0xe8,%xmm0,%xmm0
pshufd $0xf5,%xmm1,%xmm1
pmuludq %xmm4,%xmm1
pshufd $0xe8,%xmm1,%xmm1
punpckldq %xmm1,%xmm0
paddd  %xmm9,%xmm7
paddd  %xmm2,%xmm7
paddd  %xmm9,%xmm6
paddd  %xmm0,%xmm6
pxor   %xmm0,%xmm0
test   %r10d,%r10d
jne    11d0 <linear_img+0xd0>
pxor   %xmm1,%xmm1
psubd  %xmm6,%xmm1
pxor   %xmm2,%xmm2
psubd  %xmm7,%xmm2
jmp    11d8 <linear_img+0xd8>
xor    %r12d,%r12d
mov    %rdi,%r15
mov    %rsi,%rax
mov    %rdx,%rbx
sub    %r12d,%r14d
xor    %edx,%edx
jmp    135d <linear_img+0x25d>
cs nopw 0x0(%rax,%rax,1)
nop
mov    %si,(%rbx,%rdx,2)
add    $0x1,%rdx
cmp    %edx,%r14d
je     1396 <linear_img+0x296>
movswl (%r15,%rdx,2),%edi
imul   %ecx,%edi
movswl (%rax,%rdx,2),%esi
imul   %r8d,%esi
add    %r9d,%edi
add    %esi,%edi
mov    %edi,%esi
neg    %esi
test   %r10d,%r10d
cmovne %r11d,%esi
test   %edi,%edi
cmovns %edi,%esi
cmp    $0xfff,%esi
jl     1350 <linear_img+0x250>
mov    $0xfff,%esi
jmp    1350 <linear_img+0x250>
cmp    %r12,%r13
jne    133e <linear_img+0x23e>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define void @linear_img(i16* nocapture %0, i16* nocapture %1, i16* nocapture writeonly %2, i32 %3, i32 %4, i32 %5, i32 %6, i32 %7) {
%9 = icmp sgt i32 %7, 0
br i1 %9, label %10, label %99
10: ; preds = %8
%11 = icmp eq i32 %6, 0
%12 = add i32 %7, -1
%13 = zext i32 %12 to i64
%14 = add nuw nsw i64 %13, 1
%15 = icmp ult i32 %12, 7
br i1 %15, label %69, label %16
16: ; preds = %10
%17 = add i32 %7, -1
%18 = zext i32 %17 to i64
%19 = add nuw nsw i64 %18, 1
%20 = getelementptr i16, i16* %2, i64 %19
%21 = getelementptr i16, i16* %0, i64 %19
%22 = getelementptr i16, i16* %1, i64 %19
%23 = icmp ugt i16* %21, %2
%24 = icmp ugt i16* %20, %0
%25 = and i1 %23, %24
%26 = icmp ugt i16* %22, %2
%27 = icmp ugt i16* %20, %1
%28 = and i1 %26, %27
%29 = or i1 %25, %28
br i1 %29, label %69, label %30
30: ; preds = %16
%31 = and i64 %14, 8589934584
%32 = trunc i64 %31 to i32
%33 = getelementptr i16, i16* %0, i64 %31
%34 = getelementptr i16, i16* %1, i64 %31
%35 = getelementptr i16, i16* %2, i64 %31
%36 = insertelement <8 x i32> poison, i32 %3, i64 0
%37 = shufflevector <8 x i32> %36, <8 x i32> poison, <8 x i32> zeroinitializer
%38 = insertelement <8 x i32> poison, i32 %4, i64 0
%39 = shufflevector <8 x i32> %38, <8 x i32> poison, <8 x i32> zeroinitializer
%40 = insertelement <8 x i32> poison, i32 %5, i64 0
%41 = shufflevector <8 x i32> %40, <8 x i32> poison, <8 x i32> zeroinitializer
br label %42
42: ; preds = %42, %30
%43 = phi i64 [ 0, %30 ], [ %65, %42 ]
%44 = getelementptr i16, i16* %0, i64 %43
%45 = getelementptr i16, i16* %1, i64 %43
%46 = getelementptr i16, i16* %2, i64 %43
%47 = bitcast i16* %44 to <8 x i16>*
%48 = load <8 x i16>, <8 x i16>* %47
%49 = sext <8 x i16> %48 to <8 x i32>
%50 = mul nsw <8 x i32> %37, %49
%51 = bitcast i16* %45 to <8 x i16>*
%52 = load <8 x i16>, <8 x i16>* %51
%53 = sext <8 x i16> %52 to <8 x i32>
%54 = mul nsw <8 x i32> %39, %53
%55 = add <8 x i32> %50, %41
%56 = add <8 x i32> %55, %54
%57 = icmp slt <8 x i32> %56, zeroinitializer
%58 = sub nsw <8 x i32> zeroinitializer, %56
%59 = select i1 %11, <8 x i32> %58, <8 x i32> zeroinitializer
%60 = select <8 x i1> %57, <8 x i32> %59, <8 x i32> %56
%61 = icmp slt <8 x i32> %60, <i32 4095, i32 4095, i32 4095, i32 4095, i32 4095, i32 4095, i32 4095, i32 4095>
%62 = select <8 x i1> %61, <8 x i32> %60, <8 x i32> <i32 4095, i32 4095, i32 4095, i32 4095, i32 4095, i32 4095, i32 4095, i32 4095>
%63 = trunc <8 x i32> %62 to <8 x i16>
%64 = bitcast i16* %46 to <8 x i16>*
store <8 x i16> %63, <8 x i16>* %64
%65 = add nuw i64 %43, 8
%66 = icmp eq i64 %65, %31
br i1 %66, label %67, label %42
67: ; preds = %42
%68 = icmp eq i64 %14, %31
br i1 %68, label %99, label %69
69: ; preds = %16, %10, %67
%70 = phi i32 [ 0, %16 ], [ 0, %10 ], [ %32, %67 ]
%71 = phi i16* [ %0, %16 ], [ %0, %10 ], [ %33, %67 ]
%72 = phi i16* [ %1, %16 ], [ %1, %10 ], [ %34, %67 ]
%73 = phi i16* [ %2, %16 ], [ %2, %10 ], [ %35, %67 ]
br label %74
74: ; preds = %69, %74
%75 = phi i32 [ %97, %74 ], [ %70, %69 ]
%76 = phi i16* [ %95, %74 ], [ %71, %69 ]
%77 = phi i16* [ %96, %74 ], [ %72, %69 ]
%78 = phi i16* [ %94, %74 ], [ %73, %69 ]
%79 = load i16, i16* %76
%80 = sext i16 %79 to i32
%81 = mul nsw i32 %80, %3
%82 = load i16, i16* %77
%83 = sext i16 %82 to i32
%84 = mul nsw i32 %83, %4
%85 = add i32 %81, %5
%86 = add i32 %85, %84
%87 = icmp slt i32 %86, 0
%88 = sub nsw i32 0, %86
%89 = select i1 %11, i32 %88, i32 0
%90 = select i1 %87, i32 %89, i32 %86
%91 = icmp slt i32 %90, 4095
%92 = select i1 %91, i32 %90, i32 4095
%93 = trunc i32 %92 to i16
store i16 %93, i16* %78
%94 = getelementptr inbounds i16, i16* %78, i64 1
%95 = getelementptr inbounds i16, i16* %76, i64 1
%96 = getelementptr inbounds i16, i16* %77, i64 1
%97 = add nuw nsw i32 %75, 1
%98 = icmp eq i32 %97, %7
br i1 %98, label %99, label %74
99: ; preds = %74, %67, %8
ret void
}