void _bin_delta_c(double *rho, int *pixel_ind, double *pixel_weight,
                  int *radial_ind, double *radial_weight, double *out,
                  int npart, int npix, int nrad) {
#pragma omp parallel for
  for (int ipart = 0; ipart < npart; ipart++) {
    double v = rho[ipart];
    for (int ipix = 0; ipix < npix; ipix++) {
      int ip = ipart * npix + ipix;
      int pi = pixel_ind[ip];
      double pw = pixel_weight[ip];
      for (int irad = 0; irad < nrad; irad++) {
        int ir = ipart * nrad + irad;
        int ri = radial_ind[ir];
        double rw = radial_weight[ir];
        if (rw < 0)
          continue;
#pragma omp atomic
        out[ri * npix + pi] += v * pw * rw;
      }
    }
  }
}


x86 assembly (clang, -O0):
<_bin_delta_c>:
push   %rbp
mov    %rsp,%rbp
mov    0x20(%rbp),%eax
mov    0x18(%rbp),%eax
mov    0x10(%rbp),%eax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9,-0x30(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    0x10(%rbp),%eax
jge    1243 <_bin_delta_c+0x143>
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x40(%rbp)
movl   $0x0,-0x44(%rbp)
mov    -0x44(%rbp),%eax
cmp    0x18(%rbp),%eax
jge    1230 <_bin_delta_c+0x130>
mov    -0x34(%rbp),%eax
imul   0x18(%rbp),%eax
add    -0x44(%rbp),%eax
mov    %eax,-0x48(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x48(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x4c(%rbp)
mov    -0x18(%rbp),%rax
movslq -0x48(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x58(%rbp)
movl   $0x0,-0x5c(%rbp)
mov    -0x5c(%rbp),%eax
cmp    0x20(%rbp),%eax
jge    121d <_bin_delta_c+0x11d>
mov    -0x34(%rbp),%eax
imul   0x20(%rbp),%eax
add    -0x5c(%rbp),%eax
mov    %eax,-0x60(%rbp)
mov    -0x20(%rbp),%rax
movslq -0x60(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x64(%rbp)
mov    -0x28(%rbp),%rax
movslq -0x60(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x70(%rbp)
xorps  %xmm0,%xmm0
ucomisd -0x70(%rbp),%xmm0
jbe    11dd <_bin_delta_c+0xdd>
jmp    120f <_bin_delta_c+0x10f>
movsd  -0x40(%rbp),%xmm0
mulsd  -0x58(%rbp),%xmm0
movsd  -0x70(%rbp),%xmm2
mov    -0x30(%rbp),%rax
mov    -0x64(%rbp),%ecx
imul   0x18(%rbp),%ecx
add    -0x4c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x5c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x5c(%rbp)
jmp    1191 <_bin_delta_c+0x91>
jmp    1222 <_bin_delta_c+0x122>
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
jmp    1151 <_bin_delta_c+0x51>
jmp    1235 <_bin_delta_c+0x135>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    112c <_bin_delta_c+0x2c>
pop    %rbp
ret

x86 assembly (clang, -O2):
<_bin_delta_c>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    %rdi,-0x20(%rsp)
mov    0x38(%rsp),%edi
test   %edi,%edi
jle    1202 <_bin_delta_c+0x102>
mov    0x48(%rsp),%r11d
movslq 0x40(%rsp),%rax
movslq %r11d,%rbp
mov    %edi,%edi
mov    %rdi,-0x10(%rsp)
mov    %eax,%r12d
mov    %ebp,%ebx
lea    0x0(,%rbp,4),%rdi
mov    %rdi,-0x18(%rsp)
shl    $0x3,%rbp
mov    %rbp,-0x8(%rsp)
xor    %r15d,%r15d
xorpd  %xmm0,%xmm0
jmp    1179 <_bin_delta_c+0x79>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
add    $0x1,%r15
add    -0x18(%rsp),%rcx
add    -0x8(%rsp),%r8
cmp    -0x10(%rsp),%r15
je     1202 <_bin_delta_c+0x102>
cmpl   $0x0,0x40(%rsp)
jle    1160 <_bin_delta_c+0x60>
mov    -0x20(%rsp),%rdi
movsd  (%rdi,%r15,8),%xmm1
mov    %r15,%r10
imul   %rax,%r10
xor    %r14d,%r14d
jmp    11a9 <_bin_delta_c+0xa9>
nopw   0x0(%rax,%rax,1)
add    $0x1,%r14
cmp    %r12,%r14
je     1160 <_bin_delta_c+0x60>
test   %r11d,%r11d
jle    11a0 <_bin_delta_c+0xa0>
lea    (%r14,%r10,1),%rbp
movslq (%rsi,%rbp,4),%rdi
movsd  (%rdx,%rbp,8),%xmm2
mulsd  %xmm1,%xmm2
xor    %r13d,%r13d
jmp    11d9 <_bin_delta_c+0xd9>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
add    $0x1,%r13
cmp    %r13,%rbx
je     11a0 <_bin_delta_c+0xa0>
movsd  (%r8,%r13,8),%xmm3
ucomisd %xmm3,%xmm0
ja     11d0 <_bin_delta_c+0xd0>
movslq (%rcx,%r13,4),%rbp
imul   %rax,%rbp
add    %rdi,%rbp
mulsd  %xmm2,%xmm3
addsd  (%r9,%rbp,8),%xmm3
movsd  %xmm3,(%r9,%rbp,8)
jmp    11d0 <_bin_delta_c+0xd0>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define void @_bin_delta_c(double* nocapture %0, i32* nocapture %1, double* nocapture %2, i32* nocapture %3, double* nocapture %4, double* nocapture %5, i32 %6, i32 %7, i32 %8) {
%10 = icmp sgt i32 %6, 0
br i1 %10, label %11, label %19
11: ; preds = %9
%12 = icmp sgt i32 %7, 0
%13 = icmp sgt i32 %8, 0
%14 = sext i32 %7 to i64
%15 = sext i32 %8 to i64
%16 = zext i32 %6 to i64
%17 = zext i32 %7 to i64
%18 = zext i32 %8 to i64
br label %20
19: ; preds = %27, %9
ret void
20: ; preds = %11, %27
%21 = phi i64 [ 0, %11 ], [ %28, %27 ]
%22 = getelementptr inbounds double, double* %0, i64 %21
%23 = load double, double* %22
br i1 %12, label %24, label %27
24: ; preds = %20
%25 = mul nsw i64 %21, %14
%26 = mul nsw i64 %21, %15
br label %30
27: ; preds = %39, %20
%28 = add nuw nsw i64 %21, 1
%29 = icmp eq i64 %28, %16
br i1 %29, label %19, label %20
30: ; preds = %24, %39
%31 = phi i64 [ 0, %24 ], [ %40, %39 ]
%32 = add nsw i64 %31, %25
%33 = getelementptr inbounds i32, i32* %1, i64 %32
%34 = load i32, i32* %33
br i1 %13, label %35, label %39
35: ; preds = %30
%36 = getelementptr inbounds double, double* %2, i64 %32
%37 = load double, double* %36
%38 = fmul double %23, %37
br label %42
39: ; preds = %57, %30
%40 = add nuw nsw i64 %31, 1
%41 = icmp eq i64 %40, %17
br i1 %41, label %27, label %30
42: ; preds = %35, %57
%43 = phi i64 [ 0, %35 ], [ %58, %57 ]
%44 = add nsw i64 %43, %26
%45 = getelementptr inbounds double, double* %4, i64 %44
%46 = load double, double* %45
%47 = fcmp olt double %46, 0.000000e+00
br i1 %47, label %57, label %48
48: ; preds = %42
%49 = getelementptr inbounds i32, i32* %3, i64 %44
%50 = load i32, i32* %49
%51 = mul nsw i32 %50, %7
%52 = add nsw i32 %51, %34
%53 = sext i32 %52 to i64
%54 = getelementptr inbounds double, double* %5, i64 %53
%55 = load double, double* %54
%56 = tail call double @llvm.fmuladd.f64(double %38, double %46, double %55)
store double %56, double* %54
br label %57
57: ; preds = %42, %48
%58 = add nuw nsw i64 %43, 1
%59 = icmp eq i64 %58, %18
br i1 %59, label %39, label %42
}