unsigned int detectPeaks(float x[], unsigned int signal_size) {
  float s = 0;
  float sf[] = {0, 0};
  const float a = 0.2;
  unsigned int k;
  float y[] = {0, 0};
  char slopeWasPositive = 0;
  char slopeIsPositive = 0;
  unsigned int pulseCount = 0;
  for (k = 1; k < signal_size; k++) {
    y[1] = y[0] + (a * (x[k] - y[0]));
    s = y[1] - y[0];
    sf[1] = sf[0] + (a * (s - sf[0]));
    sf[1] += 0.0005;
    y[0] = y[1];
    sf[0] = sf[1];
    if (sf[1] > 0)
      slopeIsPositive = 1;
    else if (sf[1] < 0)
      slopeIsPositive = 0;
    if (!slopeIsPositive && slopeWasPositive)
      pulseCount += 1;
    if (slopeIsPositive)
      slopeWasPositive = 1;
    else
      slopeWasPositive = 0;
  }
  return pulseCount;
}


x86 assembly (clang, -O0):
<detectPeaks>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x10(%rbp)
lea    -0x18(%rbp),%rdi
xor    %esi,%esi
mov    $0x8,%edx
call   1030 <memset@plt>
movss  0xec1(%rip),%xmm0
movss  %xmm0,-0x1c(%rbp)
lea    -0x28(%rbp),%rdi
xor    %esi,%esi
mov    $0x8,%edx
call   1030 <memset@plt>
movb   $0x0,-0x29(%rbp)
movb   $0x0,-0x2a(%rbp)
movl   $0x0,-0x30(%rbp)
movl   $0x1,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0xc(%rbp),%eax
jae    1280 <detectPeaks+0x170>
movss  -0x28(%rbp),%xmm1
mov    -0x8(%rbp),%rax
mov    -0x20(%rbp),%ecx
movss  (%rax,%rcx,4),%xmm0
subss  -0x28(%rbp),%xmm0
movss  0xe6c(%rip),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x24(%rbp)
movss  -0x24(%rbp),%xmm0
subss  -0x28(%rbp),%xmm0
movss  %xmm0,-0x10(%rbp)
movss  -0x18(%rbp),%xmm1
movss  -0x10(%rbp),%xmm0
subss  -0x18(%rbp),%xmm0
movss  0xe39(%rip),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x14(%rbp)
movss  -0x14(%rbp),%xmm0
cvtss2sd %xmm0,%xmm0
movsd  0xe23(%rip),%xmm1
addsd  %xmm1,%xmm0
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0x14(%rbp)
movss  -0x24(%rbp),%xmm0
movss  %xmm0,-0x28(%rbp)
movss  -0x14(%rbp),%xmm0
movss  %xmm0,-0x18(%rbp)
movss  -0x14(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomiss %xmm1,%xmm0
jbe    1220 <detectPeaks+0x110>
movb   $0x1,-0x2a(%rbp)
jmp    1236 <detectPeaks+0x126>
xorps  %xmm0,%xmm0
ucomiss -0x14(%rbp),%xmm0
jbe    1231 <detectPeaks+0x121>
movb   $0x0,-0x2a(%rbp)
jmp    1236 <detectPeaks+0x126>
cmpb   $0x0,-0x2a(%rbp)
jne    1256 <detectPeaks+0x146>
movsbl -0x29(%rbp),%eax
cmp    $0x0,%eax
je     1256 <detectPeaks+0x146>
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
cmpb   $0x0,-0x2a(%rbp)
je     1269 <detectPeaks+0x159>
movb   $0x1,-0x29(%rbp)
jmp    126d <detectPeaks+0x15d>
movb   $0x0,-0x29(%rbp)
jmp    1272 <detectPeaks+0x162>
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    116a <detectPeaks+0x5a>
mov    -0x30(%rbp),%eax
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<detectPeaks>:
xor    %eax,%eax
cmp    $0x2,%esi
jb     11b5 <detectPeaks+0xb5>
mov    %esi,%r8d
xorps  %xmm0,%xmm0
xor    %eax,%eax
mov    $0x1,%edx
movss  0xee0(%rip),%xmm1
movsd  0xee0(%rip),%xmm2
xor    %r10d,%r10d
xor    %r9d,%r9d
xorps  %xmm3,%xmm3
xorps  %xmm4,%xmm4
jmp    116a <detectPeaks+0x6a>
cs nopw 0x0(%rax,%rax,1)
xor    %esi,%esi
test   %r11b,%r11b
setne  %sil
sete   %r10b
test   %r9d,%r9d
setne  %cl
and    %r10b,%cl
movzbl %cl,%ecx
add    %ecx,%eax
add    $0x1,%rdx
mov    %r11d,%r10d
mov    %esi,%r9d
cmp    %rdx,%r8
je     11b5 <detectPeaks+0xb5>
movaps %xmm3,%xmm5
movss  (%rdi,%rdx,4),%xmm3
subss  %xmm5,%xmm3
mulss  %xmm1,%xmm3
addss  %xmm5,%xmm3
movaps %xmm3,%xmm6
subss  %xmm5,%xmm6
subss  %xmm4,%xmm6
mulss  %xmm1,%xmm6
addss  %xmm4,%xmm6
xorps  %xmm4,%xmm4
cvtss2sd %xmm6,%xmm4
addsd  %xmm2,%xmm4
cvtsd2ss %xmm4,%xmm4
mov    $0x1,%r11b
ucomiss %xmm0,%xmm4
ja     1140 <detectPeaks+0x40>
ucomiss %xmm4,%xmm0
mov    %r10d,%r11d
jbe    1140 <detectPeaks+0x40>
xor    %r11d,%r11d
jmp    1140 <detectPeaks+0x40>
ret

LLVM IR (-O0):
define i32 @detectPeaks(float* %0, i32 %1) {
%3 = alloca float*
%4 = alloca i32
%5 = alloca float
%6 = alloca [2 x float]
%7 = alloca float
%8 = alloca i32
%9 = alloca [2 x float]
%10 = alloca i8
%11 = alloca i8
%12 = alloca i32
store float* %0, float** %3
store i32 %1, i32* %4
store float 0.000000e+00, float* %5
%13 = bitcast [2 x float]* %6 to i8*
store float 0x3FC99999A0000000, float* %7
%14 = bitcast [2 x float]* %9 to i8*
store i8 0, i8* %10
store i8 0, i8* %11
store i32 0, i32* %12
store i32 1, i32* %8
br label %15
15: ; preds = %82, %2
%16 = load i32, i32* %8
%17 = load i32, i32* %4
%18 = icmp ult i32 %16, %17
br i1 %18, label %19, label %85
19: ; preds = %15
%20 = getelementptr inbounds [2 x float], [2 x float]* %9, i64 0, i64 0
%21 = load float, float* %20
%22 = load float*, float** %3
%23 = load i32, i32* %8
%24 = zext i32 %23 to i64
%25 = getelementptr inbounds float, float* %22, i64 %24
%26 = load float, float* %25
%27 = getelementptr inbounds [2 x float], [2 x float]* %9, i64 0, i64 0
%28 = load float, float* %27
%29 = fsub float %26, %28
%30 = call float @llvm.fmuladd.f32(float 0x3FC99999A0000000, float %29, float %21)
%31 = getelementptr inbounds [2 x float], [2 x float]* %9, i64 0, i64 1
store float %30, float* %31
%32 = getelementptr inbounds [2 x float], [2 x float]* %9, i64 0, i64 1
%33 = load float, float* %32
%34 = getelementptr inbounds [2 x float], [2 x float]* %9, i64 0, i64 0
%35 = load float, float* %34
%36 = fsub float %33, %35
store float %36, float* %5
%37 = getelementptr inbounds [2 x float], [2 x float]* %6, i64 0, i64 0
%38 = load float, float* %37
%39 = load float, float* %5
%40 = getelementptr inbounds [2 x float], [2 x float]* %6, i64 0, i64 0
%41 = load float, float* %40
%42 = fsub float %39, %41
%43 = call float @llvm.fmuladd.f32(float 0x3FC99999A0000000, float %42, float %38)
%44 = getelementptr inbounds [2 x float], [2 x float]* %6, i64 0, i64 1
store float %43, float* %44
%45 = getelementptr inbounds [2 x float], [2 x float]* %6, i64 0, i64 1
%46 = load float, float* %45
%47 = fpext float %46 to double
%48 = fadd double %47, 5.000000e-04
%49 = fptrunc double %48 to float
store float %49, float* %45
%50 = getelementptr inbounds [2 x float], [2 x float]* %9, i64 0, i64 1
%51 = load float, float* %50
%52 = getelementptr inbounds [2 x float], [2 x float]* %9, i64 0, i64 0
store float %51, float* %52
%53 = getelementptr inbounds [2 x float], [2 x float]* %6, i64 0, i64 1
%54 = load float, float* %53
%55 = getelementptr inbounds [2 x float], [2 x float]* %6, i64 0, i64 0
store float %54, float* %55
%56 = getelementptr inbounds [2 x float], [2 x float]* %6, i64 0, i64 1
%57 = load float, float* %56
%58 = fcmp ogt float %57, 0.000000e+00
br i1 %58, label %59, label %60
59: ; preds = %19
store i8 1, i8* %11
br label %66
60: ; preds = %19
%61 = getelementptr inbounds [2 x float], [2 x float]* %6, i64 0, i64 1
%62 = load float, float* %61
%63 = fcmp olt float %62, 0.000000e+00
br i1 %63, label %64, label %65
64: ; preds = %60
store i8 0, i8* %11
br label %65
65: ; preds = %64, %60
br label %66
66: ; preds = %65, %59
%67 = load i8, i8* %11
%68 = icmp ne i8 %67, 0
br i1 %68, label %76, label %69
69: ; preds = %66
%70 = load i8, i8* %10
%71 = sext i8 %70 to i32
%72 = icmp ne i32 %71, 0
br i1 %72, label %73, label %76
73: ; preds = %69
%74 = load i32, i32* %12
%75 = add i32 %74, 1
store i32 %75, i32* %12
br label %76
76: ; preds = %73, %69, %66
%77 = load i8, i8* %11
%78 = icmp ne i8 %77, 0
br i1 %78, label %79, label %80
79: ; preds = %76
store i8 1, i8* %10
br label %81
80: ; preds = %76
store i8 0, i8* %10
br label %81
81: ; preds = %80, %79
br label %82
82: ; preds = %81
%83 = load i32, i32* %8
%84 = add i32 %83, 1
store i32 %84, i32* %8
br label %15
85: ; preds = %15
%86 = load i32, i32* %12
ret i32 %86
}

LLVM IR (-O2):
define i32 @detectPeaks(float* nocapture %0, i32 %1) {
%3 = icmp ugt i32 %1, 1
br i1 %3, label %4, label %38
4: ; preds = %2
%5 = zext i32 %1 to i64
br label %6
6: ; preds = %4, %27
%7 = phi i64 [ 1, %4 ], [ %36, %27 ]
%8 = phi i32 [ 0, %4 ], [ %33, %27 ]
%9 = phi i8 [ 0, %4 ], [ %28, %27 ]
%10 = phi i32 [ 0, %4 ], [ %35, %27 ]
%11 = phi float [ 0.000000e+00, %4 ], [ %16, %27 ]
%12 = phi float [ 0.000000e+00, %4 ], [ %22, %27 ]
%13 = getelementptr inbounds float, float* %0, i64 %7
%14 = load float, float* %13
%15 = fsub float %14, %11
%16 = tail call float @llvm.fmuladd.f32(float %15, float 0x3FC99999A0000000, float %11)
%17 = fsub float %16, %11
%18 = fsub float %17, %12
%19 = tail call float @llvm.fmuladd.f32(float %18, float 0x3FC99999A0000000, float %12)
%20 = fpext float %19 to double
%21 = fadd double %20, 5.000000e-04
%22 = fptrunc double %21 to float
%23 = fcmp ogt float %22, 0.000000e+00
br i1 %23, label %27, label %24
24: ; preds = %6
%25 = fcmp olt float %22, 0.000000e+00
br i1 %25, label %26, label %27
26: ; preds = %24
br label %27
27: ; preds = %6, %24, %26
%28 = phi i8 [ 0, %26 ], [ %9, %24 ], [ 1, %6 ]
%29 = icmp eq i8 %28, 0
%30 = icmp ne i32 %10, 0
%31 = and i1 %30, %29
%32 = zext i1 %31 to i32
%33 = add i32 %8, %32
%34 = xor i1 %29, true
%35 = zext i1 %34 to i32
%36 = add nuw nsw i64 %7, 1
%37 = icmp eq i64 %36, %5
br i1 %37, label %38, label %6
38: ; preds = %27, %2
%39 = phi i32 [ 0, %2 ], [ %33, %27 ]
ret i32 %39
}