void pl_index(int *kmax, int *kmin, double pmax, double pmin, long nlevs,
              double *plevs) {
  long k;
  long MaxLev = -1, MinLev = -1;
  for (k = 0; k < nlevs; k++)
    if (plevs[k] >= pmax) {
      MaxLev = k;
      break;
    }
  for (k = nlevs - 1; k >= 0; k--)
    if (plevs[k] < pmin) {
      MinLev = k;
      break;
    }
  *kmax = MaxLev;
  *kmin = MinLev;
}


x86 assembly (clang, -O0):
<pl_index>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
movsd  %xmm0,-0x18(%rbp)
movsd  %xmm1,-0x20(%rbp)
mov    %rdx,-0x28(%rbp)
mov    %rcx,-0x30(%rbp)
movq   $0xffffffffffffffff,-0x40(%rbp)
ff
movq   $0xffffffffffffffff,-0x48(%rbp)
ff
movq   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%rax
cmp    -0x28(%rbp),%rax
jge    117f <pl_index+0x7f>
mov    -0x30(%rbp),%rax
mov    -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
ucomisd -0x18(%rbp),%xmm0
jb     1169 <pl_index+0x69>
mov    -0x38(%rbp),%rax
mov    %rax,-0x40(%rbp)
jmp    117f <pl_index+0x7f>
jmp    116e <pl_index+0x6e>
mov    -0x38(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
jmp    1136 <pl_index+0x36>
mov    -0x28(%rbp),%rax
sub    $0x1,%rax
mov    %rax,-0x38(%rbp)
cmpq   $0x0,-0x38(%rbp)
jl     11d5 <pl_index+0xd5>
mov    -0x30(%rbp),%rax
mov    -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x20(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    11bf <pl_index+0xbf>
mov    -0x38(%rbp),%rax
mov    %rax,-0x48(%rbp)
jmp    11d5 <pl_index+0xd5>
jmp    11c4 <pl_index+0xc4>
mov    -0x38(%rbp),%rax
add    $0xffffffffffffffff,%rax
mov    %rax,-0x38(%rbp)
jmp    118b <pl_index+0x8b>
mov    -0x40(%rbp),%rax
mov    %eax,%ecx
mov    -0x8(%rbp),%rax
mov    %ecx,(%rax)
mov    -0x48(%rbp),%rax
mov    %eax,%ecx
mov    -0x10(%rbp),%rax
mov    %ecx,(%rax)
pop    %rbp
ret

x86 assembly (clang, -O2):
<pl_index>:
mov    $0xffffffffffffffff,%r8
test   %rdx,%rdx
jle    1130 <pl_index+0x30>
xor    %eax,%eax
xchg   %ax,%ax
movsd  (%rcx,%rax,8),%xmm2
ucomisd %xmm0,%xmm2
jae    1126 <pl_index+0x26>
add    $0x1,%rax
cmp    %rax,%rdx
jne    1110 <pl_index+0x10>
jmp    1130 <pl_index+0x30>
mov    %rax,%r8
nopl   0x0(%rax)
test   %rdx,%rdx
jle    1146 <pl_index+0x46>
lea    -0x1(%rdx),%rax
ucomisd -0x8(%rcx,%rdx,8),%xmm1
mov    %rax,%rdx
jbe    1130 <pl_index+0x30>
jmp    114b <pl_index+0x4b>
mov    $0xffffffff,%eax
mov    %r8d,(%rdi)
mov    %eax,(%rsi)
ret

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

LLVM IR (-O2):
define void @pl_index(i32* nocapture writeonly %0, i32* nocapture writeonly %1, double %2, double %3, i64 %4, double* nocapture %5) {
%7 = icmp sgt i64 %4, 0
br i1 %7, label %8, label %16
8: ; preds = %6, %13
%9 = phi i64 [ %14, %13 ], [ 0, %6 ]
%10 = getelementptr inbounds double, double* %5, i64 %9
%11 = load double, double* %10
%12 = fcmp ult double %11, %2
br i1 %12, label %13, label %16
13: ; preds = %8
%14 = add nuw nsw i64 %9, 1
%15 = icmp eq i64 %14, %4
br i1 %15, label %16, label %8
16: ; preds = %13, %8, %6
%17 = phi i64 [ -1, %6 ], [ %9, %8 ], [ -1, %13 ]
br label %18
18: ; preds = %22, %16
%19 = phi i64 [ %4, %16 ], [ %20, %22 ]
%20 = add nsw i64 %19, -1
%21 = icmp sgt i64 %19, 0
br i1 %21, label %22, label %28
22: ; preds = %18
%23 = getelementptr inbounds double, double* %5, i64 %20
%24 = load double, double* %23
%25 = fcmp olt double %24, %3
br i1 %25, label %26, label %18
26: ; preds = %22
%27 = trunc i64 %20 to i32
br label %28
28: ; preds = %18, %26
%29 = phi i32 [ %27, %26 ], [ -1, %18 ]
%30 = trunc i64 %17 to i32
store i32 %30, i32* %0
store i32 %29, i32* %1
ret void
}