void composite_subimage_vr(int pesize, int *pe_id, int pixn, double *n_subimage,
                           double *n_subopa, double *subimage) {
  int i, j, pe_no;
  double acc_opa, r, g, b;
  for (j = 0; j < pixn; j++) {
    pe_no = pe_id[0];
    subimage[j * 3 + 0] = n_subimage[pe_no * pixn * 3 + j * 3 + 0];
    subimage[j * 3 + 1] = n_subimage[pe_no * pixn * 3 + j * 3 + 1];
    subimage[j * 3 + 2] = n_subimage[pe_no * pixn * 3 + j * 3 + 2];
    acc_opa = n_subopa[pe_no * pixn + j];
    i = 0;
    while ((acc_opa < 0.99) && (i < pesize - 1)) {
      i++;
      pe_no = pe_id[i];
      r = n_subimage[pe_no * pixn * 3 + j * 3];
      g = n_subimage[pe_no * pixn * 3 + j * 3 + 1];
      b = n_subimage[pe_no * pixn * 3 + j * 3 + 2];
      subimage[j * 3] += r * (1.0 - acc_opa);
      subimage[j * 3 + 1] += g * (1.0 - acc_opa);
      subimage[j * 3 + 2] += b * (1.0 - acc_opa);
      acc_opa += n_subopa[pe_no * pixn + j] * (1.0 - acc_opa);
    }
  }
  return;
}


x86 assembly (clang, -O0):
<composite_subimage_vr>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9,-0x30(%rbp)
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1387 <composite_subimage_vr+0x287>
mov    -0x10(%rbp),%rax
mov    (%rax),%eax
mov    %eax,-0x3c(%rbp)
mov    -0x20(%rbp),%rax
mov    -0x3c(%rbp),%ecx
imul   -0x14(%rbp),%ecx
imul   $0x3,%ecx,%ecx
imul   $0x3,-0x38(%rbp),%edx
add    %edx,%ecx
add    $0x0,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x30(%rbp),%rax
imul   $0x3,-0x38(%rbp),%ecx
add    $0x0,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%rax
mov    -0x3c(%rbp),%ecx
imul   -0x14(%rbp),%ecx
imul   $0x3,%ecx,%ecx
imul   $0x3,-0x38(%rbp),%edx
add    %edx,%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x30(%rbp),%rax
imul   $0x3,-0x38(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%rax
mov    -0x3c(%rbp),%ecx
imul   -0x14(%rbp),%ecx
imul   $0x3,%ecx,%ecx
imul   $0x3,-0x38(%rbp),%edx
add    %edx,%ecx
add    $0x2,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x30(%rbp),%rax
imul   $0x3,-0x38(%rbp),%ecx
add    $0x2,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%rax
mov    -0x3c(%rbp),%ecx
imul   -0x14(%rbp),%ecx
add    -0x38(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x48(%rbp)
movl   $0x0,-0x34(%rbp)
xor    %eax,%eax
movsd  0xe08(%rip),%xmm0
ucomisd -0x48(%rbp),%xmm0
mov    %al,-0x61(%rbp)
jbe    1217 <composite_subimage_vr+0x117>
mov    -0x34(%rbp),%eax
mov    -0x4(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
setl   %al
mov    %al,-0x61(%rbp)
mov    -0x61(%rbp),%al
test   $0x1,%al
jne    1227 <composite_subimage_vr+0x127>
jmp    1374 <composite_subimage_vr+0x274>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x3c(%rbp)
mov    -0x20(%rbp),%rax
mov    -0x3c(%rbp),%ecx
imul   -0x14(%rbp),%ecx
imul   $0x3,%ecx,%ecx
imul   $0x3,-0x38(%rbp),%edx
add    %edx,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x20(%rbp),%rax
mov    -0x3c(%rbp),%ecx
imul   -0x14(%rbp),%ecx
imul   $0x3,%ecx,%ecx
imul   $0x3,-0x38(%rbp),%edx
add    %edx,%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x58(%rbp)
mov    -0x20(%rbp),%rax
mov    -0x3c(%rbp),%ecx
imul   -0x14(%rbp),%ecx
imul   $0x3,%ecx,%ecx
imul   $0x3,-0x38(%rbp),%edx
add    %edx,%ecx
add    $0x2,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x60(%rbp)
movsd  -0x50(%rbp),%xmm0
movsd  0xd54(%rip),%xmm2
subsd  -0x48(%rbp),%xmm2
mov    -0x30(%rbp),%rax
imul   $0x3,-0x38(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x58(%rbp),%xmm0
movsd  0xd25(%rip),%xmm2
subsd  -0x48(%rbp),%xmm2
mov    -0x30(%rbp),%rax
imul   $0x3,-0x38(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x60(%rbp),%xmm0
movsd  0xcf3(%rip),%xmm2
subsd  -0x48(%rbp),%xmm2
mov    -0x30(%rbp),%rax
imul   $0x3,-0x38(%rbp),%ecx
add    $0x2,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%rax
mov    -0x3c(%rbp),%ecx
imul   -0x14(%rbp),%ecx
add    -0x38(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  0xcb0(%rip),%xmm2
subsd  -0x48(%rbp),%xmm2
movsd  -0x48(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
jmp    11ee <composite_subimage_vr+0xee>
jmp    1379 <composite_subimage_vr+0x279>
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    1121 <composite_subimage_vr+0x21>
pop    %rbp
ret

x86 assembly (clang, -O2):
<composite_subimage_vr>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
test   %edx,%edx
jle    122c <composite_subimage_vr+0x12c>
movslq (%rsi),%rax
movslq %edx,%r11
imul   %r11,%rax
lea    (%rax,%rax,2),%r10
lea    -0x1(%rdi),%ebx
movslq %ebx,%r12
movslq %eax,%r14
mov    %edx,%r15d
xor    %r13d,%r13d
movsd  0xec9(%rip),%xmm0
movsd  0xec9(%rip),%xmm1
jmp    115d <composite_subimage_vr+0x5d>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
add    $0x1,%r13
cmp    %r15,%r13
je     122c <composite_subimage_vr+0x12c>
lea    0x0(,%r13,2),%rdx
add    %r13,%rdx
lea    (%rdx,%r10,1),%rax
movsd  (%rcx,%rax,8),%xmm2
movsd  %xmm2,(%r9,%rdx,8)
movsd  0x8(%rcx,%rax,8),%xmm5
movsd  %xmm5,0x8(%r9,%rdx,8)
movsd  0x10(%rcx,%rax,8),%xmm4
movsd  %xmm4,0x10(%r9,%rdx,8)
lea    (%r14,%r13,1),%rax
movsd  (%r8,%rax,8),%xmm3
ucomisd %xmm3,%xmm0
jbe    1150 <composite_subimage_vr+0x50>
cmp    $0x2,%edi
jl     1150 <composite_subimage_vr+0x50>
unpcklpd %xmm5,%xmm2
mov    $0x1,%eax
nop
movapd %xmm2,%xmm6
movapd %xmm3,%xmm5
movapd %xmm4,%xmm7
movslq (%rsi,%rax,4),%rbx
imul   %r11,%rbx
lea    (%rbx,%rbx,2),%ebp
add    %edx,%ebp
movslq %ebp,%rbp
movapd %xmm1,%xmm3
subsd  %xmm5,%xmm3
movupd (%rcx,%rbp,8),%xmm4
movapd %xmm3,%xmm2
unpcklpd %xmm3,%xmm2
mulpd  %xmm4,%xmm2
addpd  %xmm6,%xmm2
movsd  0x10(%rcx,%rbp,8),%xmm4
mulsd  %xmm3,%xmm4
addsd  %xmm7,%xmm4
movupd %xmm2,(%r9,%rdx,8)
movsd  %xmm4,0x10(%r9,%rdx,8)
add    %r13,%rbx
mulsd  (%r8,%rbx,8),%xmm3
addsd  %xmm5,%xmm3
ucomisd %xmm3,%xmm0
jbe    1150 <composite_subimage_vr+0x50>
lea    0x1(%rax),%rbx
cmp    %r12,%rax
mov    %rbx,%rax
jl     11b0 <composite_subimage_vr+0xb0>
jmp    1150 <composite_subimage_vr+0x50>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @composite_subimage_vr(i32 %0, i32* %1, i32 %2, double* %3, double* %4, double* %5) {
%7 = alloca i32
%8 = alloca i32*
%9 = alloca i32
%10 = alloca double*
%11 = alloca double*
%12 = alloca double*
%13 = alloca i32
%14 = alloca i32
%15 = alloca i32
%16 = alloca double
%17 = alloca double
%18 = alloca double
%19 = alloca double
store i32 %0, i32* %7
store i32* %1, i32** %8
store i32 %2, i32* %9
store double* %3, double** %10
store double* %4, double** %11
store double* %5, double** %12
store i32 0, i32* %14
br label %20
20: ; preds = %190, %6
%21 = load i32, i32* %14
%22 = load i32, i32* %9
%23 = icmp slt i32 %21, %22
br i1 %23, label %24, label %193
24: ; preds = %20
%25 = load i32*, i32** %8
%26 = getelementptr inbounds i32, i32* %25, i64 0
%27 = load i32, i32* %26
store i32 %27, i32* %15
%28 = load double*, double** %10
%29 = load i32, i32* %15
%30 = load i32, i32* %9
%31 = mul nsw i32 %29, %30
%32 = mul nsw i32 %31, 3
%33 = load i32, i32* %14
%34 = mul nsw i32 %33, 3
%35 = add nsw i32 %32, %34
%36 = add nsw i32 %35, 0
%37 = sext i32 %36 to i64
%38 = getelementptr inbounds double, double* %28, i64 %37
%39 = load double, double* %38
%40 = load double*, double** %12
%41 = load i32, i32* %14
%42 = mul nsw i32 %41, 3
%43 = add nsw i32 %42, 0
%44 = sext i32 %43 to i64
%45 = getelementptr inbounds double, double* %40, i64 %44
store double %39, double* %45
%46 = load double*, double** %10
%47 = load i32, i32* %15
%48 = load i32, i32* %9
%49 = mul nsw i32 %47, %48
%50 = mul nsw i32 %49, 3
%51 = load i32, i32* %14
%52 = mul nsw i32 %51, 3
%53 = add nsw i32 %50, %52
%54 = add nsw i32 %53, 1
%55 = sext i32 %54 to i64
%56 = getelementptr inbounds double, double* %46, i64 %55
%57 = load double, double* %56
%58 = load double*, double** %12
%59 = load i32, i32* %14
%60 = mul nsw i32 %59, 3
%61 = add nsw i32 %60, 1
%62 = sext i32 %61 to i64
%63 = getelementptr inbounds double, double* %58, i64 %62
store double %57, double* %63
%64 = load double*, double** %10
%65 = load i32, i32* %15
%66 = load i32, i32* %9
%67 = mul nsw i32 %65, %66
%68 = mul nsw i32 %67, 3
%69 = load i32, i32* %14
%70 = mul nsw i32 %69, 3
%71 = add nsw i32 %68, %70
%72 = add nsw i32 %71, 2
%73 = sext i32 %72 to i64
%74 = getelementptr inbounds double, double* %64, i64 %73
%75 = load double, double* %74
%76 = load double*, double** %12
%77 = load i32, i32* %14
%78 = mul nsw i32 %77, 3
%79 = add nsw i32 %78, 2
%80 = sext i32 %79 to i64
%81 = getelementptr inbounds double, double* %76, i64 %80
store double %75, double* %81
%82 = load double*, double** %11
%83 = load i32, i32* %15
%84 = load i32, i32* %9
%85 = mul nsw i32 %83, %84
%86 = load i32, i32* %14
%87 = add nsw i32 %85, %86
%88 = sext i32 %87 to i64
%89 = getelementptr inbounds double, double* %82, i64 %88
%90 = load double, double* %89
store double %90, double* %16
store i32 0, i32* %13
br label %91
91: ; preds = %101, %24
%92 = load double, double* %16
%93 = fcmp olt double %92, 0x3FEFAE147AE147AE
br i1 %93, label %94, label %99
94: ; preds = %91
%95 = load i32, i32* %13
%96 = load i32, i32* %7
%97 = sub nsw i32 %96, 1
%98 = icmp slt i32 %95, %97
br label %99
99: ; preds = %94, %91
%100 = phi i1 [ false, %91 ], [ %98, %94 ]
br i1 %100, label %101, label %189
101: ; preds = %99
%102 = load i32, i32* %13
%103 = add nsw i32 %102, 1
store i32 %103, i32* %13
%104 = load i32*, i32** %8
%105 = load i32, i32* %13
%106 = sext i32 %105 to i64
%107 = getelementptr inbounds i32, i32* %104, i64 %106
%108 = load i32, i32* %107
store i32 %108, i32* %15
%109 = load double*, double** %10
%110 = load i32, i32* %15
%111 = load i32, i32* %9
%112 = mul nsw i32 %110, %111
%113 = mul nsw i32 %112, 3
%114 = load i32, i32* %14
%115 = mul nsw i32 %114, 3
%116 = add nsw i32 %113, %115
%117 = sext i32 %116 to i64
%118 = getelementptr inbounds double, double* %109, i64 %117
%119 = load double, double* %118
store double %119, double* %17
%120 = load double*, double** %10
%121 = load i32, i32* %15
%122 = load i32, i32* %9
%123 = mul nsw i32 %121, %122
%124 = mul nsw i32 %123, 3
%125 = load i32, i32* %14
%126 = mul nsw i32 %125, 3
%127 = add nsw i32 %124, %126
%128 = add nsw i32 %127, 1
%129 = sext i32 %128 to i64
%130 = getelementptr inbounds double, double* %120, i64 %129
%131 = load double, double* %130
store double %131, double* %18
%132 = load double*, double** %10
%133 = load i32, i32* %15
%134 = load i32, i32* %9
%135 = mul nsw i32 %133, %134
%136 = mul nsw i32 %135, 3
%137 = load i32, i32* %14
%138 = mul nsw i32 %137, 3
%139 = add nsw i32 %136, %138
%140 = add nsw i32 %139, 2
%141 = sext i32 %140 to i64
%142 = getelementptr inbounds double, double* %132, i64 %141
%143 = load double, double* %142
store double %143, double* %19
%144 = load double, double* %17
%145 = load double, double* %16
%146 = fsub double 1.000000e+00, %145
%147 = load double*, double** %12
%148 = load i32, i32* %14
%149 = mul nsw i32 %148, 3
%150 = sext i32 %149 to i64
%151 = getelementptr inbounds double, double* %147, i64 %150
%152 = load double, double* %151
%153 = call double @llvm.fmuladd.f64(double %144, double %146, double %152)
store double %153, double* %151
%154 = load double, double* %18
%155 = load double, double* %16
%156 = fsub double 1.000000e+00, %155
%157 = load double*, double** %12
%158 = load i32, i32* %14
%159 = mul nsw i32 %158, 3
%160 = add nsw i32 %159, 1
%161 = sext i32 %160 to i64
%162 = getelementptr inbounds double, double* %157, i64 %161
%163 = load double, double* %162
%164 = call double @llvm.fmuladd.f64(double %154, double %156, double %163)
store double %164, double* %162
%165 = load double, double* %19
%166 = load double, double* %16
%167 = fsub double 1.000000e+00, %166
%168 = load double*, double** %12
%169 = load i32, i32* %14
%170 = mul nsw i32 %169, 3
%171 = add nsw i32 %170, 2
%172 = sext i32 %171 to i64
%173 = getelementptr inbounds double, double* %168, i64 %172
%174 = load double, double* %173
%175 = call double @llvm.fmuladd.f64(double %165, double %167, double %174)
store double %175, double* %173
%176 = load double*, double** %11
%177 = load i32, i32* %15
%178 = load i32, i32* %9
%179 = mul nsw i32 %177, %178
%180 = load i32, i32* %14
%181 = add nsw i32 %179, %180
%182 = sext i32 %181 to i64
%183 = getelementptr inbounds double, double* %176, i64 %182
%184 = load double, double* %183
%185 = load double, double* %16
%186 = fsub double 1.000000e+00, %185
%187 = load double, double* %16
%188 = call double @llvm.fmuladd.f64(double %184, double %186, double %187)
store double %188, double* %16
br label %91
189: ; preds = %99
br label %190
190: ; preds = %189
%191 = load i32, i32* %14
%192 = add nsw i32 %191, 1
store i32 %192, i32* %14
br label %20
193: ; preds = %20
ret void
}

LLVM IR (-O2):
define void @composite_subimage_vr(i32 %0, i32* nocapture %1, i32 %2, double* nocapture %3, double* nocapture %4, double* nocapture %5) {
%7 = icmp sgt i32 %2, 0
br i1 %7, label %8, label %86
8: ; preds = %6
%9 = load i32, i32* %1
%10 = mul nsw i32 %9, %2
%11 = mul nsw i32 %10, 3
%12 = add nsw i32 %0, -1
%13 = icmp sgt i32 %0, 1
%14 = sext i32 %12 to i64
%15 = sext i32 %11 to i64
%16 = sext i32 %10 to i64
%17 = zext i32 %2 to i64
br label %18
18: ; preds = %8, %83
%19 = phi i64 [ 0, %8 ], [ %84, %83 ]
%20 = mul nuw nsw i64 %19, 3
%21 = add nsw i64 %20, %15
%22 = getelementptr inbounds double, double* %3, i64 %21
%23 = load double, double* %22
%24 = getelementptr inbounds double, double* %5, i64 %20
store double %23, double* %24
%25 = add nsw i64 %21, 1
%26 = getelementptr inbounds double, double* %3, i64 %25
%27 = load double, double* %26
%28 = add nuw nsw i64 %20, 1
%29 = getelementptr inbounds double, double* %5, i64 %28
store double %27, double* %29
%30 = add nsw i64 %21, 2
%31 = getelementptr inbounds double, double* %3, i64 %30
%32 = load double, double* %31
%33 = add nuw nsw i64 %20, 2
%34 = getelementptr inbounds double, double* %5, i64 %33
store double %32, double* %34
%35 = add nsw i64 %19, %16
%36 = getelementptr inbounds double, double* %4, i64 %35
%37 = load double, double* %36
%38 = fcmp olt double %37, 0x3FEFAE147AE147AE
%39 = and i1 %38, %13
br i1 %39, label %40, label %83
40: ; preds = %18
%41 = trunc i64 %20 to i32
%42 = insertelement <2 x double> poison, double %23, i64 0
%43 = insertelement <2 x double> %42, double %27, i64 1
%44 = bitcast double* %24 to <2 x double>*
br label %45
45: ; preds = %40, %45
%46 = phi double [ %32, %40 ], [ %74, %45 ]
%47 = phi i64 [ 0, %40 ], [ %50, %45 ]
%48 = phi double [ %37, %40 ], [ %79, %45 ]
%49 = phi <2 x double> [ %43, %40 ], [ %72, %45 ]
%50 = add nuw nsw i64 %47, 1
%51 = getelementptr inbounds i32, i32* %1, i64 %50
%52 = load i32, i32* %51
%53 = mul nsw i32 %52, %2
%54 = mul nsw i32 %53, 3
%55 = add nsw i32 %54, %41
%56 = sext i32 %55 to i64
%57 = getelementptr inbounds double, double* %3, i64 %56
%58 = load double, double* %57
%59 = add nsw i32 %55, 1
%60 = sext i32 %59 to i64
%61 = getelementptr inbounds double, double* %3, i64 %60
%62 = load double, double* %61
%63 = add nsw i32 %55, 2
%64 = sext i32 %63 to i64
%65 = getelementptr inbounds double, double* %3, i64 %64
%66 = load double, double* %65
%67 = fsub double 1.000000e+00, %48
%68 = insertelement <2 x double> poison, double %58, i64 0
%69 = insertelement <2 x double> %68, double %62, i64 1
%70 = insertelement <2 x double> poison, double %67, i64 0
%71 = shufflevector <2 x double> %70, <2 x double> poison, <2 x i32> zeroinitializer
%72 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %69, <2 x double> %71, <2 x double> %49)
%73 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %69, <2 x double> %71, <2 x double> %49)
store <2 x double> %73, <2 x double>* %44
%74 = tail call double @llvm.fmuladd.f64(double %66, double %67, double %46)
store double %74, double* %34
%75 = sext i32 %53 to i64
%76 = add nsw i64 %19, %75
%77 = getelementptr inbounds double, double* %4, i64 %76
%78 = load double, double* %77
%79 = tail call double @llvm.fmuladd.f64(double %78, double %67, double %48)
%80 = fcmp olt double %79, 0x3FEFAE147AE147AE
%81 = icmp slt i64 %50, %14
%82 = select i1 %80, i1 %81, i1 false
br i1 %82, label %45, label %83
83: ; preds = %45, %18
%84 = add nuw nsw i64 %19, 1
%85 = icmp eq i64 %84, %17
br i1 %85, label %86, label %18
86: ; preds = %83, %6
ret void
}