int combine_maps(float *data, float *map, float *mask, int *sizes) {
  int i, j, k, index;
  for (i = 0; i < sizes[0]; i++) {
    for (j = 0; j < sizes[1]; j++) {
      for (k = 0; k < sizes[2]; k++) {
        index = i * sizes[2] * sizes[1] + j * sizes[2] + k;
        if (mask[index]) {
          data[index] = map[index];
        }
      }
    }
  }
  return 0;
}


x86 assembly (clang, -O0):
<combine_maps>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
mov    -0x20(%rbp),%rcx
cmp    (%rcx),%eax
jge    11f5 <combine_maps+0xf5>
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
mov    -0x20(%rbp),%rcx
cmp    0x4(%rcx),%eax
jge    11e2 <combine_maps+0xe2>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
mov    -0x20(%rbp),%rcx
cmp    0x8(%rcx),%eax
jge    11cf <combine_maps+0xcf>
mov    -0x24(%rbp),%eax
mov    -0x20(%rbp),%rcx
imul   0x8(%rcx),%eax
mov    -0x20(%rbp),%rcx
imul   0x4(%rcx),%eax
mov    -0x28(%rbp),%ecx
mov    -0x20(%rbp),%rdx
imul   0x8(%rdx),%ecx
add    %ecx,%eax
add    -0x2c(%rbp),%eax
mov    %eax,-0x30(%rbp)
mov    -0x18(%rbp),%rax
movslq -0x30(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
xorps  %xmm1,%xmm1
ucomiss %xmm1,%xmm0
jne    11a2 <combine_maps+0xa2>
jp     11a2 <combine_maps+0xa2>
jmp    11bc <combine_maps+0xbc>
mov    -0x10(%rbp),%rax
movslq -0x30(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
movslq -0x30(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
jmp    11c1 <combine_maps+0xc1>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1148 <combine_maps+0x48>
jmp    11d4 <combine_maps+0xd4>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1131 <combine_maps+0x31>
jmp    11e7 <combine_maps+0xe7>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    111b <combine_maps+0x1b>
xor    %eax,%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<combine_maps>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    %rcx,-0x8(%rsp)
mov    (%rcx),%eax
mov    %eax,-0xc(%rsp)
test   %eax,%eax
jle    1226 <combine_maps+0x126>
mov    -0x8(%rsp),%rax
mov    0x4(%rax),%r11d
xor    %r9d,%r9d
xorps  %xmm0,%xmm0
xor    %r10d,%r10d
jmp    1152 <combine_maps+0x52>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
add    $0x1,%r10d
add    %r11d,%r9d
cmp    -0xc(%rsp),%r10d
je     1226 <combine_maps+0x126>
test   %r11d,%r11d
jle    1140 <combine_maps+0x40>
mov    -0x8(%rsp),%rax
mov    0x8(%rax),%r15d
mov    %r11d,%r14d
imul   %r10d,%r14d
mov    %r15d,%eax
and    $0xfffffffe,%eax
mov    %r15d,%r12d
imul   %r9d,%r12d
xor    %r13d,%r13d
jmp    118c <combine_maps+0x8c>
nopl   0x0(%rax)
add    $0x1,%r13d
add    %r15d,%r12d
cmp    %r11d,%r13d
je     1140 <combine_maps+0x40>
test   %r15d,%r15d
jle    1180 <combine_maps+0x80>
cmp    $0x1,%r15d
jne    11d0 <combine_maps+0xd0>
xor    %ebp,%ebp
test   $0x1,%r15b
je     1180 <combine_maps+0x80>
lea    (%r14,%r13,1),%ecx
imul   %r15d,%ecx
add    %ecx,%ebp
movslq %ebp,%rcx
movss  (%rdx,%rcx,4),%xmm1
ucomiss %xmm0,%xmm1
jne    11b8 <combine_maps+0xb8>
jnp    1180 <combine_maps+0x80>
movss  (%rsi,%rcx,4),%xmm1
movss  %xmm1,(%rdi,%rcx,4)
jmp    1180 <combine_maps+0x80>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
mov    %r12d,%ebx
xor    %ebp,%ebp
jmp    11ec <combine_maps+0xec>
nopw   0x0(%rax,%rax,1)
add    $0x2,%rbp
add    $0x2,%ebx
cmp    %rbp,%rax
je     1199 <combine_maps+0x99>
movslq %ebx,%r8
movss  (%rdx,%r8,4),%xmm1
ucomiss %xmm0,%xmm1
jne    11fc <combine_maps+0xfc>
jnp    1208 <combine_maps+0x108>
movss  (%rsi,%r8,4),%xmm1
movss  %xmm1,(%rdi,%r8,4)
lea    0x1(%rbx),%ecx
movslq %ecx,%rcx
movss  (%rdx,%rcx,4),%xmm1
ucomiss %xmm0,%xmm1
jne    121a <combine_maps+0x11a>
jnp    11e0 <combine_maps+0xe0>
movss  (%rsi,%rcx,4),%xmm1
movss  %xmm1,(%rdi,%rcx,4)
jmp    11e0 <combine_maps+0xe0>
xor    %eax,%eax
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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