void crossings(const double *x, const unsigned N, double *ti,
               const double threshold) {
  unsigned i, j;
  i = 0;
  for (j = 1; j < N; j++) {
    if (x[j] > threshold && x[j - 1] < threshold) {
      ti[i++] = (double)(j - 1) + (threshold - x[j - 1]) / (x[j] - x[j - 1]);
    }
  }
}


x86 assembly (clang, -O0):
<crossings>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %rdx,-0x18(%rbp)
movsd  %xmm0,-0x20(%rbp)
movl   $0x0,-0x24(%rbp)
movl   $0x1,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0xc(%rbp),%eax
jae    11d4 <crossings+0xd4>
mov    -0x8(%rbp),%rax
mov    -0x28(%rbp),%ecx
movsd  (%rax,%rcx,8),%xmm0
ucomisd -0x20(%rbp),%xmm0
jbe    11c1 <crossings+0xc1>
mov    -0x8(%rbp),%rax
mov    -0x28(%rbp),%ecx
sub    $0x1,%ecx
mov    %ecx,%ecx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x20(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    11c1 <crossings+0xc1>
mov    -0x28(%rbp),%eax
add    $0xffffffff,%eax
cvtsi2sd %rax,%xmm0
movsd  -0x20(%rbp),%xmm1
mov    -0x8(%rbp),%rax
mov    -0x28(%rbp),%ecx
sub    $0x1,%ecx
mov    %ecx,%ecx
subsd  (%rax,%rcx,8),%xmm1
mov    -0x8(%rbp),%rax
mov    -0x28(%rbp),%ecx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x8(%rbp),%rax
mov    -0x28(%rbp),%ecx
sub    $0x1,%ecx
mov    %ecx,%ecx
subsd  (%rax,%rcx,8),%xmm2
divsd  %xmm2,%xmm1
addsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
mov    -0x24(%rbp),%ecx
mov    %ecx,%edx
add    $0x1,%edx
mov    %edx,-0x24(%rbp)
mov    %ecx,%ecx
movsd  %xmm0,(%rax,%rcx,8)
jmp    11c6 <crossings+0xc6>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1122 <crossings+0x22>
pop    %rbp
ret

x86 assembly (clang, -O2):
<crossings>:
cmp    $0x2,%esi
jb     116c <crossings+0x6c>
mov    %esi,%eax
add    $0xffffffffffffffff,%rax
xor    %ecx,%ecx
xor    %r8d,%r8d
jmp    1129 <crossings+0x29>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
add    $0x1,%rcx
cmp    %rcx,%rax
je     116c <crossings+0x6c>
movsd  0x8(%rdi,%rcx,8),%xmm1
ucomisd %xmm0,%xmm1
jbe    1120 <crossings+0x20>
movsd  (%rdi,%rcx,8),%xmm2
ucomisd %xmm2,%xmm0
jbe    1120 <crossings+0x20>
mov    %ecx,%esi
xorps  %xmm3,%xmm3
cvtsi2sd %rsi,%xmm3
movapd %xmm0,%xmm4
subsd  %xmm2,%xmm4
subsd  %xmm2,%xmm1
divsd  %xmm1,%xmm4
addsd  %xmm3,%xmm4
mov    %r8d,%esi
add    $0x1,%r8d
movsd  %xmm4,(%rdx,%rsi,8)
jmp    1120 <crossings+0x20>
ret

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

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