int find_pivot_row(const int m, const int n, const float *p, int col) {
  int row = -1;
  float minr;
  int in = n;
  for (int i = 1; i < m; i++, in += n) {
    float r = p[in] / p[in + col];
    if (r > 0 && (r < minr || row < 0)) {
      minr = r;
      row = i;
    }
  }
  return row;
}


x86 assembly (clang, -O0):
<find_pivot_row>:
push   %rbp
mov    %rsp,%rbp
mov    %edi,-0x4(%rbp)
mov    %esi,-0x8(%rbp)
mov    %rdx,-0x10(%rbp)
mov    %ecx,-0x14(%rbp)
movl   $0xffffffff,-0x18(%rbp)
mov    -0x8(%rbp),%eax
mov    %eax,-0x20(%rbp)
movl   $0x1,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    11af <find_pivot_row+0xaf>
mov    -0x10(%rbp),%rax
movslq -0x20(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%ecx
add    -0x14(%rbp),%ecx
movslq %ecx,%rcx
divss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x28(%rbp)
movss  -0x28(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomiss %xmm1,%xmm0
jbe    1193 <find_pivot_row+0x93>
movss  -0x28(%rbp),%xmm1
movss  -0x1c(%rbp),%xmm0
ucomiss %xmm1,%xmm0
ja     1183 <find_pivot_row+0x83>
cmpl   $0x0,-0x18(%rbp)
jge    1193 <find_pivot_row+0x93>
movss  -0x28(%rbp),%xmm0
movss  %xmm0,-0x1c(%rbp)
mov    -0x24(%rbp),%eax
mov    %eax,-0x18(%rbp)
jmp    1198 <find_pivot_row+0x98>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
mov    -0x8(%rbp),%eax
add    -0x20(%rbp),%eax
mov    %eax,-0x20(%rbp)
jmp    1125 <find_pivot_row+0x25>
mov    -0x18(%rbp),%eax
pop    %rbp
ret

x86 assembly (clang, -O2):
<find_pivot_row>:
mov    $0xffffffff,%eax
cmp    $0x2,%edi
jl     1169 <find_pivot_row+0x69>
movslq %esi,%r10
movslq %ecx,%rax
shl    $0x2,%rax
lea    (%rax,%r10,4),%r8
shl    $0x2,%r10
mov    $0x1,%esi
mov    $0xffffffff,%eax
xorps  %xmm0,%xmm0
jmp    1140 <find_pivot_row+0x40>
nopl   0x0(%rax,%rax,1)
cmovne %esi,%eax
movaps %xmm2,%xmm1
add    $0x1,%esi
add    %r10,%rdx
cmp    %esi,%edi
je     1169 <find_pivot_row+0x69>
movss  (%rdx,%r10,1),%xmm2
divss  (%rdx,%r8,1),%xmm2
ucomiss %xmm0,%xmm2
jbe    1136 <find_pivot_row+0x36>
ucomiss %xmm2,%xmm1
seta   %r9b
test   %eax,%eax
sets   %cl
or     %r9b,%cl
test   %cl,%cl
jne    1130 <find_pivot_row+0x30>
movaps %xmm1,%xmm2
jmp    1130 <find_pivot_row+0x30>
ret

LLVM IR (-O0):
define i32 @find_pivot_row(i32 %0, i32 %1, float* %2, i32 %3) {
%5 = alloca i32
%6 = alloca i32
%7 = alloca float*
%8 = alloca i32
%9 = alloca i32
%10 = alloca float
%11 = alloca i32
%12 = alloca i32
%13 = alloca float
store i32 %0, i32* %5
store i32 %1, i32* %6
store float* %2, float** %7
store i32 %3, i32* %8
store i32 -1, i32* %9
%14 = load i32, i32* %6
store i32 %14, i32* %11
store i32 1, i32* %12
br label %15
15: ; preds = %46, %4
%16 = load i32, i32* %12
%17 = load i32, i32* %5
%18 = icmp slt i32 %16, %17
br i1 %18, label %19, label %52
19: ; preds = %15
%20 = load float*, float** %7
%21 = load i32, i32* %11
%22 = sext i32 %21 to i64
%23 = getelementptr inbounds float, float* %20, i64 %22
%24 = load float, float* %23
%25 = load float*, float** %7
%26 = load i32, i32* %11
%27 = load i32, i32* %8
%28 = add nsw i32 %26, %27
%29 = sext i32 %28 to i64
%30 = getelementptr inbounds float, float* %25, i64 %29
%31 = load float, float* %30
%32 = fdiv float %24, %31
store float %32, float* %13
%33 = load float, float* %13
%34 = fcmp ogt float %33, 0.000000e+00
br i1 %34, label %35, label %45
35: ; preds = %19
%36 = load float, float* %13
%37 = load float, float* %10
%38 = fcmp olt float %36, %37
br i1 %38, label %42, label %39
39: ; preds = %35
%40 = load i32, i32* %9
%41 = icmp slt i32 %40, 0
br i1 %41, label %42, label %45
42: ; preds = %39, %35
%43 = load float, float* %13
store float %43, float* %10
%44 = load i32, i32* %12
store i32 %44, i32* %9
br label %45
45: ; preds = %42, %39, %19
br label %46
46: ; preds = %45
%47 = load i32, i32* %12
%48 = add nsw i32 %47, 1
store i32 %48, i32* %12
%49 = load i32, i32* %6
%50 = load i32, i32* %11
%51 = add nsw i32 %50, %49
store i32 %51, i32* %11
br label %15
52: ; preds = %15
%53 = load i32, i32* %9
ret i32 %53
}

LLVM IR (-O2):
define i32 @find_pivot_row(i32 %0, i32 %1, float* nocapture %2, i32 %3) {
%5 = icmp sgt i32 %0, 1
br i1 %5, label %6, label %9
6: ; preds = %4
%7 = sext i32 %1 to i64
%8 = sext i32 %3 to i64
br label %11
9: ; preds = %29, %4
%10 = phi i32 [ -1, %4 ], [ %31, %29 ]
ret i32 %10
11: ; preds = %6, %29
%12 = phi i64 [ %7, %6 ], [ %33, %29 ]
%13 = phi i32 [ -1, %6 ], [ %31, %29 ]
%14 = phi i32 [ 1, %6 ], [ %32, %29 ]
%15 = phi float [ undef, %6 ], [ %30, %29 ]
%16 = getelementptr inbounds float, float* %2, i64 %12
%17 = load float, float* %16
%18 = add nsw i64 %12, %8
%19 = getelementptr inbounds float, float* %2, i64 %18
%20 = load float, float* %19
%21 = fdiv float %17, %20
%22 = fcmp ogt float %21, 0.000000e+00
br i1 %22, label %23, label %29
23: ; preds = %11
%24 = fcmp olt float %21, %15
%25 = icmp slt i32 %13, 0
%26 = select i1 %24, i1 true, i1 %25
%27 = select i1 %26, float %21, float %15
%28 = select i1 %26, i32 %14, i32 %13
br label %29
29: ; preds = %23, %11
%30 = phi float [ %15, %11 ], [ %27, %23 ]
%31 = phi i32 [ %13, %11 ], [ %28, %23 ]
%32 = add nuw nsw i32 %14, 1
%33 = add i64 %12, %7
%34 = icmp eq i32 %32, %0
br i1 %34, label %9, label %11
}