void overwriteDenseF32(float *data, float *annodata, int *dims) {
  int i, j, k, index;
  int xdim = dims[0];
  int ydim = dims[1];
  int zdim = dims[2];
  for (i = 0; i < xdim; i++)
    for (j = 0; j < ydim; j++)
      for (k = 0; k < zdim; k++) {
        index = (i * zdim * ydim) + (j * zdim) + (k);
        if (annodata[index] != 0)
          data[index] = annodata[index];
      }
}


x86 assembly (clang, -O0):
<overwriteDenseF32>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    -0x18(%rbp),%rax
mov    (%rax),%eax
mov    %eax,-0x2c(%rbp)
mov    -0x18(%rbp),%rax
mov    0x4(%rax),%eax
mov    %eax,-0x30(%rbp)
mov    -0x18(%rbp),%rax
mov    0x8(%rax),%eax
mov    %eax,-0x34(%rbp)
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x2c(%rbp),%eax
jge    11f7 <overwriteDenseF32+0xf7>
movl   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x30(%rbp),%eax
jge    11e4 <overwriteDenseF32+0xe4>
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    11d1 <overwriteDenseF32+0xd1>
mov    -0x1c(%rbp),%eax
imul   -0x34(%rbp),%eax
imul   -0x30(%rbp),%eax
mov    -0x20(%rbp),%ecx
imul   -0x34(%rbp),%ecx
add    %ecx,%eax
add    -0x24(%rbp),%eax
mov    %eax,-0x28(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x28(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
xorps  %xmm1,%xmm1
ucomiss %xmm1,%xmm0
jne    11a4 <overwriteDenseF32+0xa4>
jp     11a4 <overwriteDenseF32+0xa4>
jmp    11be <overwriteDenseF32+0xbe>
mov    -0x10(%rbp),%rax
movslq -0x28(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x28(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
jmp    11c3 <overwriteDenseF32+0xc3>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    115a <overwriteDenseF32+0x5a>
jmp    11d6 <overwriteDenseF32+0xd6>
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    1147 <overwriteDenseF32+0x47>
jmp    11e9 <overwriteDenseF32+0xe9>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1134 <overwriteDenseF32+0x34>
pop    %rbp
ret

x86 assembly (clang, -O2):
<overwriteDenseF32>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    (%rdx),%eax
mov    %eax,-0x10(%rsp)
test   %eax,%eax
jle    134c <overwriteDenseF32+0x24c>
mov    0x4(%rdx),%ecx
mov    0x8(%rdx),%r9d
lea    -0x1(%r9),%r12
mov    %ecx,%eax
imul   %r9d,%eax
mov    %eax,-0x1c(%rsp)
lea    (%rsi,%r9,4),%rax
cmp    %rdi,%rax
seta   %al
lea    (%rdi,%r9,4),%rdx
cmp    %rsi,%rdx
seta   %r10b
and    %al,%r10b
mov    %r9d,%r11d
and    $0xfffffffc,%r11d
mov    %r9,%r13
neg    %r13
movl   $0x1,-0x20(%rsp)
xor    %eax,%eax
mov    %rax,-0x18(%rsp)
xorps  %xmm0,%xmm0
xorps  %xmm1,%xmm1
movl   $0x0,-0x24(%rsp)
jmp    11a9 <overwriteDenseF32+0xa9>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    -0x24(%rsp),%ebp
add    $0x1,%ebp
mov    -0x1c(%rsp),%eax
add    %eax,-0x20(%rsp)
mov    -0x18(%rsp),%rdx
add    %eax,%edx
mov    %rdx,-0x18(%rsp)
mov    %ebp,-0x24(%rsp)
cmp    -0x10(%rsp),%ebp
je     134c <overwriteDenseF32+0x24c>
test   %ecx,%ecx
jle    1180 <overwriteDenseF32+0x80>
mov    -0x1c(%rsp),%edx
mov    -0x24(%rsp),%eax
imul   %eax,%edx
mov    %edx,-0xc(%rsp)
imul   %ecx,%eax
mov    %rax,-0x8(%rsp)
xor    %r8d,%r8d
mov    -0x18(%rsp),%r15
mov    -0x20(%rsp),%eax
mov    %eax,%r14d
jmp    11ef <overwriteDenseF32+0xef>
cs nopw 0x0(%rax,%rax,1)
nop
add    $0x1,%r8d
add    %r9d,%r14d
add    %r9d,%r15d
cmp    %ecx,%r8d
je     1180 <overwriteDenseF32+0x80>
test   %r9d,%r9d
jle    11e0 <overwriteDenseF32+0xe0>
cmp    $0x4,%r9d
jb     1223 <overwriteDenseF32+0x123>
mov    %r9d,%eax
imul   %r8d,%eax
add    -0xc(%rsp),%eax
mov    %r12,%rdx
shr    $0x20,%rdx
setne  %dl
lea    (%rax,%r12,1),%ebp
cmp    %eax,%ebp
setl   %al
or     %dl,%al
or     %r10b,%al
je     12aa <overwriteDenseF32+0x1aa>
xor    %edx,%edx
mov    %rdx,%rbx
test   $0x1,%r9b
je     1257 <overwriteDenseF32+0x157>
mov    -0x8(%rsp),%rax
add    %r8d,%eax
imul   %r9d,%eax
add    %edx,%eax
movslq %eax,%rbp
movss  (%rsi,%rbp,4),%xmm2
ucomiss %xmm0,%xmm2
jne    124b <overwriteDenseF32+0x14b>
jnp    1250 <overwriteDenseF32+0x150>
movss  %xmm2,(%rdi,%rbp,4)
mov    %rdx,%rbx
or     $0x1,%rbx
not    %rdx
cmp    %r13,%rdx
je     11e0 <overwriteDenseF32+0xe0>
lea    (%r14,%rbx,1),%ebp
mov    %r9,%rdx
sub    %rbx,%rdx
add    %r15d,%ebx
jmp    1280 <overwriteDenseF32+0x180>
xchg   %ax,%ax
add    $0x2,%ebp
add    $0x2,%ebx
add    $0xfffffffffffffffe,%rdx
je     11e0 <overwriteDenseF32+0xe0>
movslq %ebx,%rax
movss  (%rsi,%rax,4),%xmm2
ucomiss %xmm0,%xmm2
jne    128f <overwriteDenseF32+0x18f>
jnp    1294 <overwriteDenseF32+0x194>
movss  %xmm2,(%rdi,%rax,4)
movslq %ebp,%rax
movss  (%rsi,%rax,4),%xmm2
ucomiss %xmm0,%xmm2
jne    12a3 <overwriteDenseF32+0x1a3>
jnp    1270 <overwriteDenseF32+0x170>
movss  %xmm2,(%rdi,%rax,4)
jmp    1270 <overwriteDenseF32+0x170>
mov    %r10d,%ebx
mov    %r14d,%r10d
mov    %r11,%rbp
jmp    12ca <overwriteDenseF32+0x1ca>
cs nopw 0x0(%rax,%rax,1)
nop
add    $0x4,%r10d
add    $0xfffffffffffffffc,%r11
je     1335 <overwriteDenseF32+0x235>
lea    -0x1(%r10),%eax
movslq %eax,%rdx
movups (%rsi,%rdx,4),%xmm2
ucomiss %xmm0,%xmm2
jne    12dc <overwriteDenseF32+0x1dc>
jnp    12e1 <overwriteDenseF32+0x1e1>
movss  %xmm2,(%rdi,%rdx,4)
movaps %xmm2,%xmm3
cmpneqps %xmm1,%xmm3
pextrw $0x2,%xmm3,%eax
test   $0x1,%al
je     1300 <overwriteDenseF32+0x200>
movslq %r10d,%rax
movaps %xmm2,%xmm4
shufps $0x55,%xmm2,%xmm4
movss  %xmm4,(%rdi,%rax,4)
pextrw $0x4,%xmm3,%eax
test   $0x1,%al
je     131b <overwriteDenseF32+0x21b>
lea    0x1(%r10),%eax
cltq
movaps %xmm2,%xmm4
unpckhpd %xmm2,%xmm4
movss  %xmm4,(%rdi,%rax,4)
pextrw $0x6,%xmm3,%eax
test   $0x1,%al
je     12c0 <overwriteDenseF32+0x1c0>
lea    0x2(%r10),%eax
cltq
shufps $0xff,%xmm2,%xmm2
movss  %xmm2,(%rdi,%rax,4)
jmp    12c0 <overwriteDenseF32+0x1c0>
mov    %rbp,%r11
mov    %rbp,%rdx
cmp    %r9,%rbp
mov    %ebx,%r10d
je     11e0 <overwriteDenseF32+0xe0>
jmp    1225 <overwriteDenseF32+0x125>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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