double sig_duration(double *arr, int n, double tn, double tk) {
  double min = arr[0], max = arr[0];
  for (int i = 1; i < n; i++) {
    if (arr[i] < min)
      min = arr[i];
    if (arr[i] > max)
      max = arr[i];
  }
  double durat = 0;
  double dt = (tk - tn) / (n - 1);
  for (int i = 0; i < n; i++) {
    if (arr[i] >= min + 0.5 * (max - min))
      durat += dt;
  }
  return durat;
}


x86 assembly (clang, -O0):
<sig_duration>:
push   %rbp
mov    %rsp,%rbp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
movsd  %xmm0,-0x18(%rbp)
movsd  %xmm1,-0x20(%rbp)
mov    -0x8(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0x8(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x30(%rbp)
movl   $0x1,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    11ad <sig_duration+0xad>
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x28(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    1170 <sig_duration+0x70>
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
ucomisd -0x30(%rbp),%xmm0
jbe    119a <sig_duration+0x9a>
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x30(%rbp)
jmp    119f <sig_duration+0x9f>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1136 <sig_duration+0x36>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  -0x20(%rbp),%xmm0
subsd  -0x18(%rbp),%xmm0
mov    -0xc(%rbp),%eax
sub    $0x1,%eax
cvtsi2sd %eax,%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
movl   $0x0,-0x4c(%rbp)
mov    -0x4c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    123d <sig_duration+0x13d>
mov    -0x8(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x28(%rbp),%xmm2
movsd  -0x30(%rbp),%xmm1
subsd  -0x28(%rbp),%xmm1
movsd  0xdf7(%rip),%xmm3
mulsd  %xmm3,%xmm1
addsd  %xmm2,%xmm1
ucomisd %xmm1,%xmm0
jb     122a <sig_duration+0x12a>
movsd  -0x48(%rbp),%xmm0
addsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
jmp    122f <sig_duration+0x12f>
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
jmp    11d9 <sig_duration+0xd9>
movsd  -0x40(%rbp),%xmm0
pop    %rbp
ret

x86 assembly (clang, -O2):
<sig_duration>:
movsd  (%rdi),%xmm3
mov    %esi,%r8d
movapd %xmm3,%xmm2
movapd %xmm3,%xmm6
cmp    $0x2,%esi
jl     11e2 <sig_duration+0xe2>
lea    -0x1(%r8),%rcx
lea    -0x2(%r8),%rdx
mov    %ecx,%eax
and    $0x3,%eax
cmp    $0x3,%rdx
jae    113a <sig_duration+0x3a>
mov    $0x1,%edx
movapd %xmm3,%xmm5
movapd %xmm3,%xmm4
jmp    11a5 <sig_duration+0xa5>
and    $0xfffffffffffffffc,%rcx
xor    %edx,%edx
movapd %xmm3,%xmm5
movapd %xmm3,%xmm4
nopl   0x0(%rax,%rax,1)
movsd  0x8(%rdi,%rdx,8),%xmm2
movsd  0x10(%rdi,%rdx,8),%xmm6
movapd %xmm2,%xmm7
minsd  %xmm5,%xmm7
maxsd  %xmm4,%xmm2
movapd %xmm6,%xmm4
minsd  %xmm7,%xmm4
maxsd  %xmm2,%xmm6
movsd  0x18(%rdi,%rdx,8),%xmm2
movapd %xmm2,%xmm7
minsd  %xmm4,%xmm7
maxsd  %xmm6,%xmm2
movsd  0x20(%rdi,%rdx,8),%xmm4
movapd %xmm4,%xmm5
minsd  %xmm7,%xmm5
maxsd  %xmm2,%xmm4
add    $0x4,%rdx
cmp    %rdx,%rcx
jne    1150 <sig_duration+0x50>
add    $0x1,%rdx
movapd %xmm4,%xmm2
movapd %xmm5,%xmm6
test   %rax,%rax
je     11e2 <sig_duration+0xe2>
lea    (%rdi,%rdx,8),%rcx
xor    %edx,%edx
nopl   0x0(%rax,%rax,1)
movsd  (%rcx,%rdx,8),%xmm2
movapd %xmm2,%xmm6
minsd  %xmm5,%xmm6
maxsd  %xmm4,%xmm2
add    $0x1,%rdx
movapd %xmm6,%xmm5
movapd %xmm2,%xmm4
cmp    %rdx,%rax
jne    11c0 <sig_duration+0xc0>
test   %esi,%esi
jle    124c <sig_duration+0x14c>
subsd  %xmm0,%xmm1
lea    -0x1(%rsi),%eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
divsd  %xmm0,%xmm1
subsd  %xmm6,%xmm2
mulsd  0xdfc(%rip),%xmm2
addsd  %xmm6,%xmm2
xorpd  %xmm4,%xmm4
addsd  %xmm1,%xmm4
movapd %xmm2,%xmm0
cmpnlesd %xmm3,%xmm0
andnpd %xmm4,%xmm0
cmp    $0x1,%esi
je     1327 <sig_duration+0x227>
lea    -0x1(%r8),%rdx
add    $0xfffffffffffffffe,%r8
mov    %edx,%ecx
and    $0x3,%ecx
cmp    $0x3,%r8
jae    1250 <sig_duration+0x150>
mov    $0x1,%eax
test   %rcx,%rcx
jne    12f2 <sig_duration+0x1f2>
jmp    1327 <sig_duration+0x227>
xorps  %xmm0,%xmm0
ret
and    $0xfffffffffffffffc,%rdx
xor    %eax,%eax
movsd  0xdaa(%rip),%xmm3
xchg   %ax,%ax
movapd %xmm2,%xmm4
cmpnlesd 0x8(%rdi,%rax,8),%xmm4
movapd %xmm4,%xmm5
andpd  %xmm3,%xmm5
andnpd %xmm1,%xmm4
orpd   %xmm5,%xmm4
addsd  %xmm0,%xmm4
movapd %xmm2,%xmm0
cmpnlesd 0x10(%rdi,%rax,8),%xmm0
movapd %xmm0,%xmm5
andpd  %xmm3,%xmm5
andnpd %xmm1,%xmm0
orpd   %xmm5,%xmm0
addsd  %xmm4,%xmm0
movapd %xmm2,%xmm4
cmpnlesd 0x18(%rdi,%rax,8),%xmm4
movapd %xmm4,%xmm5
andpd  %xmm3,%xmm5
andnpd %xmm1,%xmm4
orpd   %xmm5,%xmm4
addsd  %xmm0,%xmm4
movapd %xmm2,%xmm0
cmpnlesd 0x20(%rdi,%rax,8),%xmm0
movapd %xmm0,%xmm5
andpd  %xmm3,%xmm5
andnpd %xmm1,%xmm0
orpd   %xmm5,%xmm0
addsd  %xmm4,%xmm0
add    $0x4,%rax
cmp    %rax,%rdx
jne    1260 <sig_duration+0x160>
add    $0x1,%rax
test   %rcx,%rcx
je     1327 <sig_duration+0x227>
lea    (%rdi,%rax,8),%rax
xor    %edx,%edx
movsd  0xd08(%rip),%xmm3
movapd %xmm2,%xmm4
cmpnlesd (%rax,%rdx,8),%xmm4
movapd %xmm4,%xmm5
andpd  %xmm3,%xmm5
andnpd %xmm1,%xmm4
orpd   %xmm5,%xmm4
addsd  %xmm4,%xmm0
add    $0x1,%rdx
cmp    %rdx,%rcx
jne    1300 <sig_duration+0x200>
ret

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

LLVM IR (-O2):
define double @sig_duration(double* nocapture %0, i32 %1, double %2, double %3) {
%5 = load double, double* %0
%6 = icmp sgt i32 %1, 1
br i1 %6, label %7, label %36
7: ; preds = %4
%8 = zext i32 %1 to i64
%9 = add nsw i64 %8, -1
%10 = add nsw i64 %8, -2
%11 = and i64 %9, 3
%12 = icmp ult i64 %10, 3
br i1 %12, label %15, label %13
13: ; preds = %7
%14 = and i64 %9, -4
br label %59
15: ; preds = %59, %7
%16 = phi double [ undef, %7 ], [ %88, %59 ]
%17 = phi double [ undef, %7 ], [ %90, %59 ]
%18 = phi i64 [ 1, %7 ], [ %91, %59 ]
%19 = phi double [ %5, %7 ], [ %88, %59 ]
%20 = phi double [ %5, %7 ], [ %90, %59 ]
%21 = icmp eq i64 %11, 0
br i1 %21, label %36, label %22
22: ; preds = %15, %22
%23 = phi i64 [ %33, %22 ], [ %18, %15 ]
%24 = phi double [ %30, %22 ], [ %19, %15 ]
%25 = phi double [ %32, %22 ], [ %20, %15 ]
%26 = phi i64 [ %34, %22 ], [ 0, %15 ]
%27 = getelementptr inbounds double, double* %0, i64 %23
%28 = load double, double* %27
%29 = fcmp olt double %28, %24
%30 = select i1 %29, double %28, double %24
%31 = fcmp ogt double %28, %25
%32 = select i1 %31, double %28, double %25
%33 = add nuw nsw i64 %23, 1
%34 = add i64 %26, 1
%35 = icmp eq i64 %34, %11
br i1 %35, label %36, label %22
36: ; preds = %15, %22, %4
%37 = phi double [ %5, %4 ], [ %17, %15 ], [ %32, %22 ]
%38 = phi double [ %5, %4 ], [ %16, %15 ], [ %30, %22 ]
%39 = fsub double %3, %2
%40 = add nsw i32 %1, -1
%41 = sitofp i32 %40 to double
%42 = fdiv double %39, %41
%43 = icmp sgt i32 %1, 0
br i1 %43, label %44, label %111
44: ; preds = %36
%45 = fsub double %37, %38
%46 = tail call double @llvm.fmuladd.f64(double %45, double 5.000000e-01, double %38)
%47 = zext i32 %1 to i64
%48 = fcmp ult double %5, %46
%49 = fadd double %42, 0.000000e+00
%50 = select i1 %48, double 0.000000e+00, double %49
%51 = icmp eq i32 %1, 1
br i1 %51, label %111, label %52
52: ; preds = %44
%53 = add nsw i64 %47, -1
%54 = add nsw i64 %47, -2
%55 = and i64 %53, 3
%56 = icmp ult i64 %54, 3
br i1 %56, label %94, label %57
57: ; preds = %52
%58 = and i64 %53, -4
br label %113
59: ; preds = %59, %13
%60 = phi i64 [ 1, %13 ], [ %91, %59 ]
%61 = phi double [ %5, %13 ], [ %88, %59 ]
%62 = phi double [ %5, %13 ], [ %90, %59 ]
%63 = phi i64 [ 0, %13 ], [ %92, %59 ]
%64 = getelementptr inbounds double, double* %0, i64 %60
%65 = load double, double* %64
%66 = fcmp olt double %65, %61
%67 = select i1 %66, double %65, double %61
%68 = fcmp ogt double %65, %62
%69 = select i1 %68, double %65, double %62
%70 = add nuw nsw i64 %60, 1
%71 = getelementptr inbounds double, double* %0, i64 %70
%72 = load double, double* %71
%73 = fcmp olt double %72, %67
%74 = select i1 %73, double %72, double %67
%75 = fcmp ogt double %72, %69
%76 = select i1 %75, double %72, double %69
%77 = add nuw nsw i64 %60, 2
%78 = getelementptr inbounds double, double* %0, i64 %77
%79 = load double, double* %78
%80 = fcmp olt double %79, %74
%81 = select i1 %80, double %79, double %74
%82 = fcmp ogt double %79, %76
%83 = select i1 %82, double %79, double %76
%84 = add nuw nsw i64 %60, 3
%85 = getelementptr inbounds double, double* %0, i64 %84
%86 = load double, double* %85
%87 = fcmp olt double %86, %81
%88 = select i1 %87, double %86, double %81
%89 = fcmp ogt double %86, %83
%90 = select i1 %89, double %86, double %83
%91 = add nuw nsw i64 %60, 4
%92 = add i64 %63, 4
%93 = icmp eq i64 %92, %14
br i1 %93, label %15, label %59
94: ; preds = %113, %52
%95 = phi double [ undef, %52 ], [ %139, %113 ]
%96 = phi i64 [ 1, %52 ], [ %140, %113 ]
%97 = phi double [ %50, %52 ], [ %139, %113 ]
%98 = icmp eq i64 %55, 0
br i1 %98, label %111, label %99
99: ; preds = %94, %99
%100 = phi i64 [ %108, %99 ], [ %96, %94 ]
%101 = phi double [ %107, %99 ], [ %97, %94 ]
%102 = phi i64 [ %109, %99 ], [ 0, %94 ]
%103 = getelementptr inbounds double, double* %0, i64 %100
%104 = load double, double* %103
%105 = fcmp ult double %104, %46
%106 = select i1 %105, double -0.000000e+00, double %42
%107 = fadd double %101, %106
%108 = add nuw nsw i64 %100, 1
%109 = add i64 %102, 1
%110 = icmp eq i64 %109, %55
br i1 %110, label %111, label %99
111: ; preds = %94, %99, %44, %36
%112 = phi double [ 0.000000e+00, %36 ], [ %50, %44 ], [ %95, %94 ], [ %107, %99 ]
ret double %112
113: ; preds = %113, %57
%114 = phi i64 [ 1, %57 ], [ %140, %113 ]
%115 = phi double [ %50, %57 ], [ %139, %113 ]
%116 = phi i64 [ 0, %57 ], [ %141, %113 ]
%117 = getelementptr inbounds double, double* %0, i64 %114
%118 = load double, double* %117
%119 = fcmp ult double %118, %46
%120 = select i1 %119, double -0.000000e+00, double %42
%121 = fadd double %115, %120
%122 = add nuw nsw i64 %114, 1
%123 = getelementptr inbounds double, double* %0, i64 %122
%124 = load double, double* %123
%125 = fcmp ult double %124, %46
%126 = select i1 %125, double -0.000000e+00, double %42
%127 = fadd double %121, %126
%128 = add nuw nsw i64 %114, 2
%129 = getelementptr inbounds double, double* %0, i64 %128
%130 = load double, double* %129
%131 = fcmp ult double %130, %46
%132 = select i1 %131, double -0.000000e+00, double %42
%133 = fadd double %127, %132
%134 = add nuw nsw i64 %114, 3
%135 = getelementptr inbounds double, double* %0, i64 %134
%136 = load double, double* %135
%137 = fcmp ult double %136, %46
%138 = select i1 %137, double -0.000000e+00, double %42
%139 = fadd double %133, %138
%140 = add nuw nsw i64 %114, 4
%141 = add i64 %116, 4
%142 = icmp eq i64 %141, %58
br i1 %142, label %94, label %113
}