void arr_float_last_nonzero(unsigned short int *puIndex, float *pFloat,
                            unsigned long int num, unsigned long int uStride) {
  float fTol = 1E-6;
  unsigned short int uIndex = ((unsigned short int)(0 - 1));
  unsigned short int i;
  for (i = 0; i < num; i++) {
    float v = pFloat[i * uStride];
    if (v > fTol) {
      uIndex = i;
    }
  }
  *puIndex = uIndex;
}


x86 assembly (clang, -O0):
<arr_float_last_nonzero>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
movss  0xee4(%rip),%xmm0
movss  %xmm0,-0x24(%rbp)
movw   $0xffff,-0x26(%rbp)
movw   $0x0,-0x28(%rbp)
movzwl -0x28(%rbp),%eax
cmp    -0x18(%rbp),%rax
jae    117f <arr_float_last_nonzero+0x7f>
mov    -0x10(%rbp),%rax
movzwl -0x28(%rbp),%ecx
imul   -0x20(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x2c(%rbp)
movss  -0x2c(%rbp),%xmm0
ucomiss -0x24(%rbp),%xmm0
jbe    1169 <arr_float_last_nonzero+0x69>
mov    -0x28(%rbp),%ax
mov    %ax,-0x26(%rbp)
jmp    116e <arr_float_last_nonzero+0x6e>
mov    -0x28(%rbp),%ax
add    $0x1,%ax
mov    %ax,-0x28(%rbp)
jmp    112d <arr_float_last_nonzero+0x2d>
mov    -0x26(%rbp),%cx
mov    -0x8(%rbp),%rax
mov    %cx,(%rax)
pop    %rbp
ret

x86 assembly (clang, -O2):
<arr_float_last_nonzero>:
test   %rdx,%rdx
je     1145 <arr_float_last_nonzero+0x45>
shl    $0x2,%rcx
mov    $0xffff,%r9w
xor    %r8d,%r8d
movss  0xee7(%rip),%xmm0
nopl   0x0(%rax)
movss  (%rsi),%xmm1
ucomiss %xmm0,%xmm1
mov    %r8d,%eax
cmovbe %r9d,%eax
add    $0x1,%r8
movzwl %r8w,%r10d
add    %rcx,%rsi
mov    %eax,%r9d
cmp    %rdx,%r10
jb     1120 <arr_float_last_nonzero+0x20>
mov    %ax,(%rdi)
ret
mov    $0xffff,%ax
mov    %ax,(%rdi)
ret

LLVM IR (-O0):
define void @arr_float_last_nonzero(i16* %0, float* %1, i64 %2, i64 %3) {
%5 = alloca i16*
%6 = alloca float*
%7 = alloca i64
%8 = alloca i64
%9 = alloca float
%10 = alloca i16
%11 = alloca i16
%12 = alloca float
store i16* %0, i16** %5
store float* %1, float** %6
store i64 %2, i64* %7
store i64 %3, i64* %8
store float 0x3EB0C6F7A0000000, float* %9
store i16 -1, i16* %10
store i16 0, i16* %11
br label %13
13: ; preds = %32, %4
%14 = load i16, i16* %11
%15 = zext i16 %14 to i64
%16 = load i64, i64* %7
%17 = icmp ult i64 %15, %16
br i1 %17, label %18, label %35
18: ; preds = %13
%19 = load float*, float** %6
%20 = load i16, i16* %11
%21 = zext i16 %20 to i64
%22 = load i64, i64* %8
%23 = mul i64 %21, %22
%24 = getelementptr inbounds float, float* %19, i64 %23
%25 = load float, float* %24
store float %25, float* %12
%26 = load float, float* %12
%27 = load float, float* %9
%28 = fcmp ogt float %26, %27
br i1 %28, label %29, label %31
29: ; preds = %18
%30 = load i16, i16* %11
store i16 %30, i16* %10
br label %31
31: ; preds = %29, %18
br label %32
32: ; preds = %31
%33 = load i16, i16* %11
%34 = add i16 %33, 1
store i16 %34, i16* %11
br label %13
35: ; preds = %13
%36 = load i16, i16* %10
%37 = load i16*, i16** %5
store i16 %36, i16* %37
ret void
}

LLVM IR (-O2):
define void @arr_float_last_nonzero(i16* nocapture writeonly %0, float* nocapture %1, i64 %2, i64 %3) {
%5 = icmp eq i64 %2, 0
br i1 %5, label %18, label %6
6: ; preds = %4, %6
%7 = phi i64 [ %15, %6 ], [ 0, %4 ]
%8 = phi i16 [ %14, %6 ], [ -1, %4 ]
%9 = mul i64 %7, %3
%10 = getelementptr inbounds float, float* %1, i64 %9
%11 = load float, float* %10
%12 = fcmp ogt float %11, 0x3EB0C6F7A0000000
%13 = trunc i64 %7 to i16
%14 = select i1 %12, i16 %13, i16 %8
%15 = add i64 %7, 1
%16 = and i64 %15, 65535
%17 = icmp ult i64 %16, %2
br i1 %17, label %6, label %18
18: ; preds = %6, %4
%19 = phi i16 [ -1, %4 ], [ %14, %6 ]
store i16 %19, i16* %0
ret void
}