#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
int /*<<< orphan*/ isinff(float);
int /*<<< orphan*/ isnanf(float);
void arr_float_sum_stages_pos(float *pDest, float *pSrc,
                              unsigned long int uStates,
                              unsigned long int uStages) {
  float fTol = 1E-3;
  unsigned short int i, j;
  for (i = 0; i < uStates; i++) {
    pDest[i] = 0;
    for (j = 0; j < uStages; j++) {
      float v = pSrc[j * uStates + i];
      if (!isnanf(v) && !isinff(v)) {
        if (v > fTol) {
          pDest[i] += v;
        }
      }
    }
  }
}


x86 assembly (clang, -O0):
<arr_float_sum_stages_pos>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
movss  0xec0(%rip),%xmm0
movss  %xmm0,-0x24(%rbp)
movw   $0x0,-0x26(%rbp)
movzwl -0x26(%rbp),%eax
cmp    -0x18(%rbp),%rax
jae    1218 <arr_float_sum_stages_pos+0xf8>
mov    -0x8(%rbp),%rax
movzwl -0x26(%rbp),%ecx
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax,%rcx,4)
movw   $0x0,-0x28(%rbp)
movzwl -0x28(%rbp),%eax
cmp    -0x20(%rbp),%rax
jae    1202 <arr_float_sum_stages_pos+0xe2>
mov    -0x10(%rbp),%rax
movzwl -0x28(%rbp),%ecx
imul   -0x18(%rbp),%rcx
movzwl -0x26(%rbp),%edx
add    %rdx,%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x2c(%rbp)
movss  -0x2c(%rbp),%xmm0
call   1030 <isnanf@plt>
cmp    $0x0,%eax
jne    11ec <arr_float_sum_stages_pos+0xcc>
movss  -0x2c(%rbp),%xmm0
call   1040 <isinff@plt>
cmp    $0x0,%eax
jne    11ec <arr_float_sum_stages_pos+0xcc>
movss  -0x2c(%rbp),%xmm0
ucomiss -0x24(%rbp),%xmm0
jbe    11e7 <arr_float_sum_stages_pos+0xc7>
movss  -0x2c(%rbp),%xmm0
mov    -0x8(%rbp),%rax
movzwl -0x26(%rbp),%ecx
addss  (%rax,%rcx,4),%xmm0
movss  %xmm0,(%rax,%rcx,4)
jmp    11ec <arr_float_sum_stages_pos+0xcc>
jmp    11f1 <arr_float_sum_stages_pos+0xd1>
mov    -0x28(%rbp),%ax
add    $0x1,%ax
mov    %ax,-0x28(%rbp)
jmp    116f <arr_float_sum_stages_pos+0x4f>
jmp    1207 <arr_float_sum_stages_pos+0xe7>
mov    -0x26(%rbp),%ax
add    $0x1,%ax
mov    %ax,-0x26(%rbp)
jmp    114b <arr_float_sum_stages_pos+0x2b>
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<arr_float_sum_stages_pos>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
push   %rax
test   %rdx,%rdx
je     11cc <arr_float_sum_stages_pos+0xac>
mov    %rcx,%r14
mov    %rdx,%r13
mov    %rsi,%r15
mov    %rdi,%r12
xor    %ebx,%ebx
jmp    115b <arr_float_sum_stages_pos+0x3b>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
add    $0x1,%ebx
movzwl %bx,%ebx
cmp    %r13,%rbx
jae    11cc <arr_float_sum_stages_pos+0xac>
movl   $0x0,(%r12,%rbx,4)
test   %r14,%r14
je     1150 <arr_float_sum_stages_pos+0x30>
xor    %ebp,%ebp
jmp    117b <arr_float_sum_stages_pos+0x5b>
nopl   0x0(%rax)
add    $0x1,%ebp
movzwl %bp,%ebp
cmp    %r14,%rbp
jae    1150 <arr_float_sum_stages_pos+0x30>
mov    %rbp,%rax
imul   %r13,%rax
add    %rbx,%rax
movss  (%r15,%rax,4),%xmm0
movss  %xmm0,0x4(%rsp)
call   1030 <isnanf@plt>
test   %eax,%eax
jne    1170 <arr_float_sum_stages_pos+0x50>
movss  0x4(%rsp),%xmm0
call   1040 <isinff@plt>
movss  0x4(%rsp),%xmm0
ucomiss 0xe4e(%rip),%xmm0
jbe    1170 <arr_float_sum_stages_pos+0x50>
test   %eax,%eax
jne    1170 <arr_float_sum_stages_pos+0x50>
movss  0x4(%rsp),%xmm0
addss  (%r12,%rbx,4),%xmm0
movss  %xmm0,(%r12,%rbx,4)
jmp    1170 <arr_float_sum_stages_pos+0x50>
add    $0x8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define void @arr_float_sum_stages_pos(float* nocapture %0, float* nocapture %1, i64 %2, i64 %3) {
%5 = icmp eq i64 %2, 0
br i1 %5, label %35, label %6
6: ; preds = %4
%7 = icmp eq i64 %3, 0
br label %8
8: ; preds = %6, %31
%9 = phi i64 [ 0, %6 ], [ %33, %31 ]
%10 = getelementptr inbounds float, float* %0, i64 %9
store float 0.000000e+00, float* %10
br i1 %7, label %31, label %11
11: ; preds = %8, %27
%12 = phi i64 [ %29, %27 ], [ 0, %8 ]
%13 = mul i64 %12, %2
%14 = add i64 %13, %9
%15 = getelementptr inbounds float, float* %1, i64 %14
%16 = load float, float* %15
%17 = tail call i32 @isnanf(float %16)
%18 = icmp eq i32 %17, 0
br i1 %18, label %19, label %27
19: ; preds = %11
%20 = tail call i32 @isinff(float %16)
%21 = icmp eq i32 %20, 0
%22 = fcmp ogt float %16, 0x3F50624DE0000000
%23 = and i1 %22, %21
br i1 %23, label %24, label %27
24: ; preds = %19
%25 = load float, float* %10
%26 = fadd float %16, %25
store float %26, float* %10
br label %27
27: ; preds = %24, %19, %11
%28 = add nuw nsw i64 %12, 1
%29 = and i64 %28, 65535
%30 = icmp ult i64 %29, %3
br i1 %30, label %11, label %31
31: ; preds = %27, %8
%32 = add nuw nsw i64 %9, 1
%33 = and i64 %32, 65535
%34 = icmp ult i64 %33, %2
br i1 %34, label %8, label %35
35: ; preds = %31, %4
ret void
}