void drawDots(float *img, int w, int h, int bands, int pitch, int x[], int y[],
              int n, float *rgb) {
  int i, k;
  for (i = 0; i < n; i++) {
    if (x[i] < 0 || x[i] >= w || y[i] < 0 || y[i] >= h)
      continue;
    for (k = 0; k < bands; k++)
      img[pitch * y[i] + bands * x[i] + k] = rgb[k];
  }
  return;
}


x86 assembly (clang, -O0):
<drawDots>:
push   %rbp
mov    %rsp,%rbp
mov    0x20(%rbp),%rax
mov    0x18(%rbp),%eax
mov    0x10(%rbp),%rax
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
mov    %ecx,-0x14(%rbp)
mov    %r8d,-0x18(%rbp)
mov    %r9,-0x20(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    0x18(%rbp),%eax
jge    11f8 <drawDots+0xf8>
mov    -0x20(%rbp),%rax
movslq -0x24(%rbp),%rcx
cmpl   $0x0,(%rax,%rcx,4)
jl     1183 <drawDots+0x83>
mov    -0x20(%rbp),%rax
movslq -0x24(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
cmp    -0xc(%rbp),%eax
jge    1183 <drawDots+0x83>
mov    0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
cmpl   $0x0,(%rax,%rcx,4)
jl     1183 <drawDots+0x83>
mov    0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
cmp    -0x10(%rbp),%eax
jl     1188 <drawDots+0x88>
jmp    11ea <drawDots+0xea>
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    11e5 <drawDots+0xe5>
mov    0x20(%rbp),%rax
movslq -0x28(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x18(%rbp),%ecx
mov    0x10(%rbp),%rdx
movslq -0x24(%rbp),%rsi
imul   (%rdx,%rsi,4),%ecx
mov    -0x14(%rbp),%edx
mov    -0x20(%rbp),%rsi
movslq -0x24(%rbp),%rdi
imul   (%rsi,%rdi,4),%edx
add    %edx,%ecx
add    -0x28(%rbp),%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    118f <drawDots+0x8f>
jmp    11ea <drawDots+0xea>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    112b <drawDots+0x2b>
pop    %rbp
ret

x86 assembly (clang, -O2):
<drawDots>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    %r8d,-0x54(%rsp)
mov    %rdi,-0x60(%rsp)
mov    0x40(%rsp),%eax
test   %eax,%eax
jle    1382 <drawDots+0x282>
mov    0x48(%rsp),%r8
mov    %eax,%r11d
mov    %ecx,%r13d
mov    -0x60(%rsp),%rax
lea    (%rax,%r13,4),%rdi
mov    %rdi,-0x18(%rsp)
lea    (%r8,%r13,4),%rdi
mov    %rdi,-0x8(%rsp)
mov    %r13d,%edi
and    $0xfffffff8,%edi
mov    %rdi,-0x20(%rsp)
add    $0xfffffffffffffff8,%rdi
mov    %rdi,-0x28(%rsp)
shr    $0x3,%rdi
add    $0x1,%rdi
mov    %edi,%r14d
and    $0x3,%r14d
and    $0xfffffffffffffffc,%rdi
mov    %rdi,-0x48(%rsp)
mov    %r13d,%edi
and    $0x3,%edi
lea    0x70(%rax),%rbp
mov    %rbp,-0x50(%rsp)
lea    0x10(%rax),%rbp
mov    %rbp,-0x38(%rsp)
lea    0x10(%r8),%rbp
mov    %rbp,-0x40(%rsp)
mov    %r14,-0x30(%rsp)
shl    $0x5,%r14
add    $0xc,%rax
mov    %rax,-0x10(%rsp)
xor    %r15d,%r15d
jmp    11bd <drawDots+0xbd>
cs nopw 0x0(%rax,%rax,1)
add    $0x1,%r15
cmp    %r11,%r15
je     1382 <drawDots+0x282>
mov    (%r9,%r15,4),%eax
test   %eax,%eax
js     11b0 <drawDots+0xb0>
cmp    %esi,%eax
jge    11b0 <drawDots+0xb0>
mov    0x38(%rsp),%rbp
mov    0x0(%rbp,%r15,4),%ebp
test   %ebp,%ebp
js     11b0 <drawDots+0xb0>
cmp    %edx,%ebp
jge    11b0 <drawDots+0xb0>
test   %ecx,%ecx
jle    11b0 <drawDots+0xb0>
imul   -0x54(%rsp),%ebp
imul   %ecx,%eax
add    %ebp,%eax
movslq %eax,%r12
cmp    $0x8,%ecx
jb     1217 <drawDots+0x117>
mov    -0x60(%rsp),%rax
lea    (%rax,%r12,4),%rax
cmp    -0x8(%rsp),%rax
jae    12b4 <drawDots+0x1b4>
mov    -0x18(%rsp),%rax
lea    (%rax,%r12,4),%rax
cmp    %r8,%rax
jbe    12b4 <drawDots+0x1b4>
xor    %r10d,%r10d
mov    %r10,%rax
not    %rax
add    %r13,%rax
test   %rdi,%rdi
je     1256 <drawDots+0x156>
mov    -0x60(%rsp),%rbx
lea    (%rbx,%r12,4),%rbx
mov    %rdi,%rbp
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movss  (%r8,%r10,4),%xmm0
movss  %xmm0,(%rbx,%r10,4)
add    $0x1,%r10
add    $0xffffffffffffffff,%rbp
jne    1240 <drawDots+0x140>
cmp    $0x3,%rax
jb     11b0 <drawDots+0xb0>
mov    -0x10(%rsp),%rax
lea    (%rax,%r12,4),%rax
nopl   0x0(%rax)
movss  (%r8,%r10,4),%xmm0
movss  %xmm0,-0xc(%rax,%r10,4)
movss  0x4(%r8,%r10,4),%xmm0
movss  %xmm0,-0x8(%rax,%r10,4)
movss  0x8(%r8,%r10,4),%xmm0
movss  %xmm0,-0x4(%rax,%r10,4)
movss  0xc(%r8,%r10,4),%xmm0
movss  %xmm0,(%rax,%r10,4)
add    $0x4,%r10
cmp    %r10,%r13
jne    1270 <drawDots+0x170>
jmp    11b0 <drawDots+0xb0>
cmpq   $0x18,-0x28(%rsp)
jae    12c0 <drawDots+0x1c0>
xor    %ebp,%ebp
jmp    1330 <drawDots+0x230>
mov    -0x50(%rsp),%rax
lea    (%rax,%r12,4),%rbx
mov    -0x48(%rsp),%rax
xor    %ebp,%ebp
movups (%r8,%rbp,4),%xmm0
movups 0x10(%r8,%rbp,4),%xmm1
movups %xmm0,-0x70(%rbx,%rbp,4)
movups %xmm1,-0x60(%rbx,%rbp,4)
movups 0x20(%r8,%rbp,4),%xmm0
movups 0x30(%r8,%rbp,4),%xmm1
movups %xmm0,-0x50(%rbx,%rbp,4)
movups %xmm1,-0x40(%rbx,%rbp,4)
movups 0x40(%r8,%rbp,4),%xmm0
movups 0x50(%r8,%rbp,4),%xmm1
movups %xmm0,-0x30(%rbx,%rbp,4)
movups %xmm1,-0x20(%rbx,%rbp,4)
movups 0x60(%r8,%rbp,4),%xmm0
movups 0x70(%r8,%rbp,4),%xmm1
movups %xmm0,-0x10(%rbx,%rbp,4)
movups %xmm1,(%rbx,%rbp,4)
add    $0x20,%rbp
add    $0xfffffffffffffffc,%rax
jne    12d0 <drawDots+0x1d0>
cmpq   $0x0,-0x30(%rsp)
je     136c <drawDots+0x26c>
lea    (%r12,%rbp,1),%rax
mov    -0x38(%rsp),%rbx
lea    (%rbx,%rax,4),%rax
mov    -0x40(%rsp),%rbx
lea    (%rbx,%rbp,4),%rbp
xor    %ebx,%ebx
movups -0x10(%rbp,%rbx,1),%xmm0
movups 0x0(%rbp,%rbx,1),%xmm1
movups %xmm0,-0x10(%rax,%rbx,1)
movups %xmm1,(%rax,%rbx,1)
add    $0x20,%rbx
cmp    %rbx,%r14
jne    1350 <drawDots+0x250>
mov    -0x20(%rsp),%rax
mov    %rax,%r10
cmp    %r13,%rax
je     11b0 <drawDots+0xb0>
jmp    121a <drawDots+0x11a>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @drawDots(float* %0, i32 %1, i32 %2, i32 %3, i32 %4, i32* %5, i32* %6, i32 %7, float* %8) {
%10 = alloca float*
%11 = alloca i32
%12 = alloca i32
%13 = alloca i32
%14 = alloca i32
%15 = alloca i32*
%16 = alloca i32*
%17 = alloca i32
%18 = alloca float*
%19 = alloca i32
%20 = alloca i32
store float* %0, float** %10
store i32 %1, i32* %11
store i32 %2, i32* %12
store i32 %3, i32* %13
store i32 %4, i32* %14
store i32* %5, i32** %15
store i32* %6, i32** %16
store i32 %7, i32* %17
store float* %8, float** %18
store i32 0, i32* %19
br label %21
21: ; preds = %91, %9
%22 = load i32, i32* %19
%23 = load i32, i32* %17
%24 = icmp slt i32 %22, %23
br i1 %24, label %25, label %94
25: ; preds = %21
%26 = load i32*, i32** %15
%27 = load i32, i32* %19
%28 = sext i32 %27 to i64
%29 = getelementptr inbounds i32, i32* %26, i64 %28
%30 = load i32, i32* %29
%31 = icmp slt i32 %30, 0
br i1 %31, label %55, label %32
32: ; preds = %25
%33 = load i32*, i32** %15
%34 = load i32, i32* %19
%35 = sext i32 %34 to i64
%36 = getelementptr inbounds i32, i32* %33, i64 %35
%37 = load i32, i32* %36
%38 = load i32, i32* %11
%39 = icmp sge i32 %37, %38
br i1 %39, label %55, label %40
40: ; preds = %32
%41 = load i32*, i32** %16
%42 = load i32, i32* %19
%43 = sext i32 %42 to i64
%44 = getelementptr inbounds i32, i32* %41, i64 %43
%45 = load i32, i32* %44
%46 = icmp slt i32 %45, 0
br i1 %46, label %55, label %47
47: ; preds = %40
%48 = load i32*, i32** %16
%49 = load i32, i32* %19
%50 = sext i32 %49 to i64
%51 = getelementptr inbounds i32, i32* %48, i64 %50
%52 = load i32, i32* %51
%53 = load i32, i32* %12
%54 = icmp sge i32 %52, %53
br i1 %54, label %55, label %56
55: ; preds = %47, %40, %32, %25
br label %91
56: ; preds = %47
store i32 0, i32* %20
br label %57
57: ; preds = %87, %56
%58 = load i32, i32* %20
%59 = load i32, i32* %13
%60 = icmp slt i32 %58, %59
br i1 %60, label %61, label %90
61: ; preds = %57
%62 = load float*, float** %18
%63 = load i32, i32* %20
%64 = sext i32 %63 to i64
%65 = getelementptr inbounds float, float* %62, i64 %64
%66 = load float, float* %65
%67 = load float*, float** %10
%68 = load i32, i32* %14
%69 = load i32*, i32** %16
%70 = load i32, i32* %19
%71 = sext i32 %70 to i64
%72 = getelementptr inbounds i32, i32* %69, i64 %71
%73 = load i32, i32* %72
%74 = mul nsw i32 %68, %73
%75 = load i32, i32* %13
%76 = load i32*, i32** %15
%77 = load i32, i32* %19
%78 = sext i32 %77 to i64
%79 = getelementptr inbounds i32, i32* %76, i64 %78
%80 = load i32, i32* %79
%81 = mul nsw i32 %75, %80
%82 = add nsw i32 %74, %81
%83 = load i32, i32* %20
%84 = add nsw i32 %82, %83
%85 = sext i32 %84 to i64
%86 = getelementptr inbounds float, float* %67, i64 %85
store float %66, float* %86
br label %87
87: ; preds = %61
%88 = load i32, i32* %20
%89 = add nsw i32 %88, 1
store i32 %89, i32* %20
br label %57
90: ; preds = %57
br label %91
91: ; preds = %90, %55
%92 = load i32, i32* %19
%93 = add nsw i32 %92, 1
store i32 %93, i32* %19
br label %21
94: ; preds = %21
ret void
}

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