void dp_all_indexes(double *p, int np, double *tv, int ntv, int *idxv) {
  int pi;
  int i;
  pi = 0;
  for (i = 0; i < ntv; ++i) {
    while (pi < np - 2 && tv[i] >= p[pi + 1])
      ++pi;
    idxv[i] = pi;
  }
}


x86 assembly (clang, -O0):
<dp_all_indexes>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    %r8,-0x28(%rbp)
movl   $0x0,-0x2c(%rbp)
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    11aa <dp_all_indexes+0xaa>
jmp    1135 <dp_all_indexes+0x35>
mov    -0x2c(%rbp),%ecx
mov    -0xc(%rbp),%edx
sub    $0x2,%edx
xor    %eax,%eax
cmp    %edx,%ecx
mov    %al,-0x31(%rbp)
jge    1170 <dp_all_indexes+0x70>
mov    -0x18(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x2c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
ucomisd (%rax,%rcx,8),%xmm0
setae  %al
mov    %al,-0x31(%rbp)
mov    -0x31(%rbp),%al
test   $0x1,%al
jne    1180 <dp_all_indexes+0x80>
jmp    118e <dp_all_indexes+0x8e>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1135 <dp_all_indexes+0x35>
mov    -0x2c(%rbp),%edx
mov    -0x28(%rbp),%rax
movslq -0x30(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1124 <dp_all_indexes+0x24>
pop    %rbp
ret

x86 assembly (clang, -O2):
<dp_all_indexes>:
push   %rbx
test   %ecx,%ecx
jle    115c <dp_all_indexes+0x5c>
add    $0xfffffffe,%esi
movslq %esi,%r10
mov    %ecx,%r9d
xor    %r11d,%r11d
xor    %eax,%eax
jmp    112f <dp_all_indexes+0x2f>
cs nopw 0x0(%rax,%rax,1)
nop
mov    %ebx,%eax
mov    %eax,(%r8,%r11,4)
add    $0x1,%r11
cmp    %r9,%r11
je     115c <dp_all_indexes+0x5c>
movslq %eax,%rcx
cmp    %esi,%ecx
mov    %esi,%ebx
cmovg  %ecx,%ebx
lea    -0x1(%rcx),%eax
nopl   0x0(%rax)
cmp    %r10,%rcx
jge    1120 <dp_all_indexes+0x20>
movsd  (%rdx,%r11,8),%xmm0
add    $0x1,%eax
ucomisd 0x8(%rdi,%rcx,8),%xmm0
lea    0x1(%rcx),%rcx
jae    1140 <dp_all_indexes+0x40>
jmp    1122 <dp_all_indexes+0x22>
pop    %rbx
ret

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

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