void average_values(int numvalues, int numgrids, float *grids[],
                    float *result) {
  int i, j;
  for (i = 0; i < numvalues; i++) {
    float sum = 0.0;
    int count = 0;
    for (j = 0; j < numgrids; j++) {
      float value = grids[j][i];
      if (!((value) >= 1.0e30)) {
        sum += value;
        count++;
      }
    }
    if (count > 0) {
      result[i] = sum / (float)count;
    } else {
      result[i] = 1.0e35;
    }
  }
}


x86 assembly (clang, -O0):
<average_values>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %rdx,-0x10(%rbp)
mov    %rcx,-0x18(%rbp)
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    11f9 <average_values+0xf9>
xorps  %xmm0,%xmm0
movss  %xmm0,-0x24(%rbp)
movl   $0x0,-0x28(%rbp)
movl   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x8(%rbp),%eax
jge    11a7 <average_values+0xa7>
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x1c(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x2c(%rbp)
movss  -0x2c(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xe96(%rip),%xmm1
ucomisd %xmm1,%xmm0
jae    1194 <average_values+0x94>
movss  -0x2c(%rbp),%xmm0
addss  -0x24(%rbp),%xmm0
movss  %xmm0,-0x24(%rbp)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1199 <average_values+0x99>
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    113b <average_values+0x3b>
cmpl   $0x0,-0x28(%rbp)
jle    11d1 <average_values+0xd1>
movss  -0x24(%rbp),%xmm0
cvtsi2ssl -0x28(%rbp),%xmm1
divss  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
jmp    11e6 <average_values+0xe6>
mov    -0x18(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movss  0xe1f(%rip),%xmm0
movss  %xmm0,(%rax,%rcx,4)
jmp    11eb <average_values+0xeb>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1119 <average_values+0x19>
pop    %rbp
ret

x86 assembly (clang, -O2):
<average_values>:
push   %rbx
test   %edi,%edi
jle    11fe <average_values+0xfe>
mov    %edi,%r9d
mov    %esi,%r8d
mov    %r8d,%r10d
and    $0xfffffffe,%r10d
xor    %r11d,%r11d
movss  0xedf(%rip),%xmm0
movsd  0xedf(%rip),%xmm1
movss  0xed3(%rip),%xmm2
jmp    1153 <average_values+0x53>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movss  %xmm4,(%rcx,%r11,4)
add    $0x1,%r11
cmp    %r9,%r11
je     11fe <average_values+0xfe>
movaps %xmm0,%xmm4
test   %esi,%esi
jle    1140 <average_values+0x40>
xorps  %xmm3,%xmm3
xor    %eax,%eax
xor    %edi,%edi
cmp    $0x1,%esi
jne    11c4 <average_values+0xc4>
test   $0x1,%r8b
je     118d <average_values+0x8d>
mov    (%rdx,%rax,8),%rax
movss  (%rax,%r11,4),%xmm4
xorps  %xmm5,%xmm5
cvtss2sd %xmm4,%xmm5
ucomisd %xmm1,%xmm5
jb     1186 <average_values+0x86>
movaps %xmm2,%xmm4
addss  %xmm4,%xmm3
adc    $0x0,%edi
movaps %xmm0,%xmm4
test   %edi,%edi
je     1140 <average_values+0x40>
xorps  %xmm4,%xmm4
cvtsi2ss %edi,%xmm4
divss  %xmm4,%xmm3
movaps %xmm3,%xmm4
jmp    1140 <average_values+0x40>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
addss  %xmm4,%xmm3
addss  %xmm5,%xmm3
adc    $0x0,%edi
add    $0x2,%rax
cmp    %rax,%r10
je     1166 <average_values+0x66>
mov    (%rdx,%rax,8),%rbx
movss  (%rbx,%r11,4),%xmm4
xorps  %xmm5,%xmm5
cvtss2sd %xmm4,%xmm5
ucomisd %xmm1,%xmm5
jb     11de <average_values+0xde>
movaps %xmm2,%xmm4
mov    0x8(%rdx,%rax,8),%rbx
adc    $0x0,%edi
movss  (%rbx,%r11,4),%xmm5
xorps  %xmm6,%xmm6
cvtss2sd %xmm5,%xmm6
ucomisd %xmm1,%xmm6
jb     11b0 <average_values+0xb0>
movaps %xmm2,%xmm5
jmp    11b0 <average_values+0xb0>
pop    %rbx
ret

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

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