void calc_fuzz(double *g_dst, double *fuzz, int npart, int npix) {
  int i, j, k;
  double gsum, tmp;
  for (k = 0; k < npix; ++k) {
    for (i = 0; i < npart; ++i) {
      gsum = 0.0;
      tmp = *(g_dst + k * npart + i);
      for (j = 0; j < npart; ++j) {
        gsum += tmp / (*(g_dst + k * npart + j));
      }
      *(fuzz + k * npart + i) = 1.0 / gsum;
    }
  }
}


x86 assembly (clang, -O0):
<calc_fuzz>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x18(%rbp),%eax
jge    11fd <calc_fuzz+0xfd>
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    11ea <calc_fuzz+0xea>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x24(%rbp),%ecx
imul   -0x14(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
movslq -0x1c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x38(%rbp)
movl   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    11b1 <calc_fuzz+0xb1>
movsd  -0x38(%rbp),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x24(%rbp),%ecx
imul   -0x14(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
movslq -0x20(%rbp),%rcx
divsd  (%rax,%rcx,8),%xmm0
addsd  -0x30(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    116a <calc_fuzz+0x6a>
movsd  0xe47(%rip),%xmm0
divsd  -0x30(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
imul   -0x14(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
movslq -0x1c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    112c <calc_fuzz+0x2c>
jmp    11ef <calc_fuzz+0xef>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1119 <calc_fuzz+0x19>
pop    %rbp
ret

x86 assembly (clang, -O2):
<calc_fuzz>:
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
test   %ecx,%ecx
jle    11d6 <calc_fuzz+0xd6>
movslq %edx,%r8
mov    %ecx,%r9d
mov    %r8d,%r10d
mov    %r10d,%eax
and    $0xfffffffe,%eax
lea    0x8(%rdi),%rcx
lea    0x0(,%r8,8),%r11
xor    %r14d,%r14d
movsd  0xec9(%rip),%xmm0
jmp    1150 <calc_fuzz+0x50>
nopl   0x0(%rax)
add    $0x1,%r14
add    %r11,%rcx
cmp    %r9,%r14
je     11d6 <calc_fuzz+0xd6>
test   %edx,%edx
jle    1140 <calc_fuzz+0x40>
mov    %r14,%rbx
imul   %r8,%rbx
lea    (%rdi,%rbx,8),%r12
lea    (%rsi,%rbx,8),%r15
xor    %r13d,%r13d
jmp    1187 <calc_fuzz+0x87>
nopl   0x0(%rax,%rax,1)
movapd %xmm0,%xmm1
divsd  %xmm2,%xmm1
movsd  %xmm1,(%r15,%r13,8)
add    $0x1,%r13
cmp    %r10,%r13
je     1140 <calc_fuzz+0x40>
movsd  (%r12,%r13,8),%xmm1
xorpd  %xmm2,%xmm2
xor    %ebx,%ebx
cmp    $0x1,%edx
je     11c4 <calc_fuzz+0xc4>
nopl   0x0(%rax,%rax,1)
movapd %xmm1,%xmm3
divsd  -0x8(%rcx,%rbx,8),%xmm3
addsd  %xmm2,%xmm3
movapd %xmm1,%xmm2
divsd  (%rcx,%rbx,8),%xmm2
addsd  %xmm3,%xmm2
add    $0x2,%rbx
cmp    %rbx,%rax
jne    11a0 <calc_fuzz+0xa0>
test   $0x1,%r10b
je     1170 <calc_fuzz+0x70>
divsd  (%r12,%rbx,8),%xmm1
addsd  %xmm1,%xmm2
jmp    1170 <calc_fuzz+0x70>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
ret

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

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