void canny_threshold_edges(int width, int height, int *data, int *edges,
                           int *no_of_edges) {
  int x, y, i;
  *no_of_edges = 0;
  x = width - 1;
  y = height - 1;
  for (i = (width * height) - 1; i >= 0; i--) {
    if (data[i] != 0) {
      data[i] = 0;
      if (*no_of_edges < 1000000) {
        edges[*no_of_edges * 2] = x;
        edges[*no_of_edges * 2 + 1] = y;
        *no_of_edges = *no_of_edges + 1;
      }
    } else {
      data[i] = 255;
    }
    x--;
    if (x < 0) {
      x += width;
      y--;
    }
  }
}


x86 assembly (clang, -O0):
<canny_threshold_edges>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %rdx,-0x10(%rbp)
mov    %rcx,-0x18(%rbp)
mov    %r8,-0x20(%rbp)
mov    -0x20(%rbp),%rax
movl   $0x0,(%rax)
mov    -0x4(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x24(%rbp)
mov    -0x8(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x28(%rbp)
mov    -0x4(%rbp),%eax
imul   -0x8(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x2c(%rbp)
cmpl   $0x0,-0x2c(%rbp)
jl     1204 <canny_threshold_edges+0x104>
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
cmpl   $0x0,(%rax,%rcx,4)
je     11bd <canny_threshold_edges+0xbd>
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movl   $0x0,(%rax,%rcx,4)
mov    -0x20(%rbp),%rax
cmpl   $0xf4240,(%rax)
jge    11b8 <canny_threshold_edges+0xb8>
mov    -0x24(%rbp),%edx
mov    -0x18(%rbp),%rax
mov    -0x20(%rbp),%rcx
mov    (%rcx),%ecx
shl    $0x1,%ecx
movslq %ecx,%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x28(%rbp),%edx
mov    -0x18(%rbp),%rax
mov    -0x20(%rbp),%rcx
mov    (%rcx),%ecx
shl    $0x1,%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x20(%rbp),%rax
mov    (%rax),%ecx
add    $0x1,%ecx
mov    -0x20(%rbp),%rax
mov    %ecx,(%rax)
jmp    11cc <canny_threshold_edges+0xcc>
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
movl   $0xff,(%rax,%rcx,4)
mov    -0x24(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x24(%rbp)
cmpl   $0x0,-0x24(%rbp)
jge    11f1 <canny_threshold_edges+0xf1>
mov    -0x4(%rbp),%eax
add    -0x24(%rbp),%eax
mov    %eax,-0x24(%rbp)
mov    -0x28(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x28(%rbp)
jmp    11f6 <canny_threshold_edges+0xf6>
mov    -0x2c(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x2c(%rbp)
jmp    113f <canny_threshold_edges+0x3f>
pop    %rbp
ret

x86 assembly (clang, -O2):
<canny_threshold_edges>:
movl   $0x0,(%r8)
mov    %esi,%eax
imul   %edi,%eax
test   %eax,%eax
jle    117e <canny_threshold_edges+0x7e>
add    $0xffffffff,%esi
lea    -0x1(%rdi),%r9d
mov    %eax,%r11d
add    $0x1,%r11
jmp    114e <canny_threshold_edges+0x4e>
movl   $0xff,(%rdx,%r10,4)
xor    %r10d,%r10d
test   %r9d,%r9d
setle  %r10b
mov    $0x0,%eax
cmovle %edi,%eax
add    %eax,%r9d
add    $0xffffffff,%r9d
sub    %r10d,%esi
add    $0xffffffffffffffff,%r11
cmp    $0x1,%r11
jbe    117e <canny_threshold_edges+0x7e>
lea    -0x2(%r11),%r10d
cmpl   $0x0,(%rdx,%r10,4)
je     1120 <canny_threshold_edges+0x20>
movl   $0x0,(%rdx,%r10,4)
movslq (%r8),%r10
cmp    $0xf423f,%r10
jg     1128 <canny_threshold_edges+0x28>
mov    %r9d,(%rcx,%r10,8)
movslq (%r8),%rax
mov    %esi,0x4(%rcx,%rax,8)
addl   $0x1,(%r8)
jmp    1128 <canny_threshold_edges+0x28>
ret

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

LLVM IR (-O2):
define void @canny_threshold_edges(i32 %0, i32 %1, i32* nocapture %2, i32* nocapture writeonly %3, i32* nocapture %4) {
store i32 0, i32* %4
%6 = mul i32 %1, %0
%7 = icmp sgt i32 %6, 0
br i1 %7, label %8, label %44
8: ; preds = %5
%9 = add nsw i32 %1, -1
%10 = add nsw i32 %0, -1
%11 = zext i32 %6 to i64
br label %12
12: ; preds = %8, %36
%13 = phi i64 [ %11, %8 ], [ %16, %36 ]
%14 = phi i32 [ %9, %8 ], [ %42, %36 ]
%15 = phi i32 [ %10, %8 ], [ %40, %36 ]
%16 = add nsw i64 %13, -1
%17 = and i64 %16, 4294967295
%18 = getelementptr inbounds i32, i32* %2, i64 %17
%19 = load i32, i32* %18
%20 = icmp eq i32 %19, 0
br i1 %20, label %35, label %21
21: ; preds = %12
store i32 0, i32* %18
%22 = load i32, i32* %4
%23 = icmp slt i32 %22, 1000000
br i1 %23, label %24, label %36
24: ; preds = %21
%25 = shl nsw i32 %22, 1
%26 = sext i32 %25 to i64
%27 = getelementptr inbounds i32, i32* %3, i64 %26
store i32 %15, i32* %27
%28 = load i32, i32* %4
%29 = shl nsw i32 %28, 1
%30 = or i32 %29, 1
%31 = sext i32 %30 to i64
%32 = getelementptr inbounds i32, i32* %3, i64 %31
store i32 %14, i32* %32
%33 = load i32, i32* %4
%34 = add nsw i32 %33, 1
store i32 %34, i32* %4
br label %36
35: ; preds = %12
store i32 255, i32* %18
br label %36
36: ; preds = %21, %24, %35
%37 = add nsw i32 %15, -1
%38 = icmp slt i32 %15, 1
%39 = select i1 %38, i32 %0, i32 0
%40 = add nsw i32 %37, %39
%41 = sext i1 %38 to i32
%42 = add nsw i32 %14, %41
%43 = icmp ugt i64 %13, 1
br i1 %43, label %12, label %44
44: ; preds = %36, %5
ret void
}