double trapfilter(int *pulse, int peaking, int gap, int M, int len, int len_in,
                  double N, int bl, int pol) {
  int delayk_buffer[16384];
  int delayk_index;
  int shiftk;
  long datak_out;
  int delaym_buffer[16384];
  int delaym_index;
  int shiftm;
  long datam_in;
  long datam_out;
  double slowshaperout;
  int i = 0;
  for (i = 0; i < 16384; i++) {
    delayk_buffer[i] = 0;
    delaym_buffer[i] = 0;
  }
  delayk_index = 0;
  shiftk = peaking;
  delaym_index = 0;
  shiftm = peaking + gap;
  if (len_in < len) {
    len = len_in;
  }
  int n = 0;
  i = 0;
  double max = 0;
  double accum1 = 0;
  double accum2 = 0;
  double adder2out;
  double adder3out;
  int input;
  while (i < len) {
    input = (pulse[i] - bl) * pol;
    datak_out = delayk_buffer[delayk_index];
    delayk_buffer[delayk_index] = input;
    delayk_index++;
    if (delayk_index >= shiftk)
      delayk_index = 0;
    datam_in = input - datak_out;
    datam_out = delaym_buffer[delaym_index];
    delaym_buffer[delaym_index] = datam_in;
    delaym_index++;
    if (delaym_index >= shiftm)
      delaym_index = 0;
    adder2out = datam_in - datam_out;
    accum1 += adder2out;
    adder3out = adder2out * M + accum1 * N;
    accum2 += adder3out;
    slowshaperout = accum2 / (double)peaking;
    if (slowshaperout > max) {
      max = slowshaperout;
    }
    pulse[i] = slowshaperout;
    i++;
  }
  return max;
}


x86 assembly (clang, -O0):
<trapfilter>:
push   %rbp
mov    %rsp,%rbp
sub    $0x20020,%rsp
mov    0x18(%rbp),%eax
mov    0x10(%rbp),%eax
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
mov    %ecx,-0x14(%rbp)
mov    %r8d,-0x18(%rbp)
mov    %r9d,-0x1c(%rbp)
movsd  %xmm0,-0x28(%rbp)
movl   $0x0,-0x20064(%rbp)
movl   $0x0,-0x20064(%rbp)
cmpl   $0x4000,-0x20064(%rbp)
jge    1187 <trapfilter+0x87>
movslq -0x20064(%rbp),%rax
movl   $0x0,-0x10030(%rbp,%rax,4)
movslq -0x20064(%rbp),%rax
movl   $0x0,-0x20040(%rbp,%rax,4)
mov    -0x20064(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20064(%rbp)
jmp    113f <trapfilter+0x3f>
movl   $0x0,-0x10034(%rbp)
mov    -0xc(%rbp),%eax
mov    %eax,-0x10038(%rbp)
movl   $0x0,-0x20044(%rbp)
mov    -0xc(%rbp),%eax
add    -0x10(%rbp),%eax
mov    %eax,-0x20048(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x18(%rbp),%eax
jge    11c2 <trapfilter+0xc2>
mov    -0x1c(%rbp),%eax
mov    %eax,-0x18(%rbp)
movl   $0x0,-0x20068(%rbp)
movl   $0x0,-0x20064(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x20070(%rbp)
ff
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x20078(%rbp)
ff
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x20080(%rbp)
ff
mov    -0x20064(%rbp),%eax
cmp    -0x18(%rbp),%eax
jge    13c1 <trapfilter+0x2c1>
mov    -0x8(%rbp),%rax
movslq -0x20064(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
sub    0x10(%rbp),%eax
imul   0x18(%rbp),%eax
mov    %eax,-0x20094(%rbp)
movslq -0x10034(%rbp),%rax
movslq -0x10030(%rbp,%rax,4),%rax
ff
mov    %rax,-0x10040(%rbp)
mov    -0x20094(%rbp),%ecx
movslq -0x10034(%rbp),%rax
mov    %ecx,-0x10030(%rbp,%rax,4)
mov    -0x10034(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10034(%rbp)
mov    -0x10034(%rbp),%eax
cmp    -0x10038(%rbp),%eax
jl     1276 <trapfilter+0x176>
movl   $0x0,-0x10034(%rbp)
movslq -0x20094(%rbp),%rax
sub    -0x10040(%rbp),%rax
mov    %rax,-0x20050(%rbp)
movslq -0x20044(%rbp),%rax
movslq -0x20040(%rbp,%rax,4),%rax
ff
mov    %rax,-0x20058(%rbp)
mov    -0x20050(%rbp),%rax
mov    %eax,%ecx
movslq -0x20044(%rbp),%rax
mov    %ecx,-0x20040(%rbp,%rax,4)
mov    -0x20044(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20044(%rbp)
mov    -0x20044(%rbp),%eax
cmp    -0x20048(%rbp),%eax
jl     12e3 <trapfilter+0x1e3>
movl   $0x0,-0x20044(%rbp)
mov    -0x20050(%rbp),%rax
sub    -0x20058(%rbp),%rax
cvtsi2sd %rax,%xmm0
movsd  %xmm0,-0x20088(%rbp)
ff
movsd  -0x20088(%rbp),%xmm0
ff
addsd  -0x20078(%rbp),%xmm0
ff
movsd  %xmm0,-0x20078(%rbp)
ff
movsd  -0x20088(%rbp),%xmm0
ff
cvtsi2sdl -0x14(%rbp),%xmm2
movsd  -0x20078(%rbp),%xmm1
ff
mulsd  -0x28(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x20090(%rbp)
ff
movsd  -0x20090(%rbp),%xmm0
ff
addsd  -0x20080(%rbp),%xmm0
ff
movsd  %xmm0,-0x20080(%rbp)
ff
movsd  -0x20080(%rbp),%xmm0
ff
cvtsi2sdl -0xc(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x20060(%rbp)
ff
movsd  -0x20060(%rbp),%xmm0
ff
ucomisd -0x20070(%rbp),%xmm0
ff
jbe    1397 <trapfilter+0x297>
movsd  -0x20060(%rbp),%xmm0
ff
movsd  %xmm0,-0x20070(%rbp)
ff
cvttsd2si -0x20060(%rbp),%edx
ff
mov    -0x8(%rbp),%rax
movslq -0x20064(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x20064(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20064(%rbp)
jmp    11f7 <trapfilter+0xf7>
movsd  -0x20070(%rbp),%xmm0
ff
add    $0x20020,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<trapfilter>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x20018,%rsp
movsd  %xmm0,0x8(%rsp)
mov    %r9d,%ebp
mov    %r8d,%r15d
mov    %ecx,0x4(%rsp)
mov    %edx,%r14d
mov    %esi,%r13d
mov    %rdi,%rbx
lea    0x10010(%rsp),%rdi
xor    %r12d,%r12d
mov    $0x10000,%edx
xor    %esi,%esi
call   1030 <memset@plt>
lea    0x10(%rsp),%rdi
mov    $0x10000,%edx
xor    %esi,%esi
call   1030 <memset@plt>
cmp    %r15d,%ebp
cmovl  %ebp,%r15d
test   %r15d,%r15d
jle    123a <trapfilter+0x12a>
mov    0x20058(%rsp),%r8d
mov    0x20050(%rsp),%r9d
cvtsi2sdl 0x4(%rsp),%xmm1
cvtsi2sd %r13d,%xmm2
add    %r13d,%r14d
mov    %r15d,%edx
xorpd  %xmm3,%xmm3
xor    %esi,%esi
xor    %edi,%edi
xor    %eax,%eax
xorpd  %xmm4,%xmm4
xorpd  %xmm5,%xmm5
movsd  0x8(%rsp),%xmm7
nopl   0x0(%rax,%rax,1)
mov    (%rbx,%rsi,4),%ecx
sub    %r9d,%ecx
imul   %r8d,%ecx
movslq %edi,%rdi
movslq 0x10010(%rsp,%rdi,4),%rbp
mov    %ecx,0x10010(%rsp,%rdi,4)
add    $0x1,%edi
cmp    %r13d,%edi
cmovge %r12d,%edi
movslq %ecx,%rcx
sub    %rbp,%rcx
cltq
movslq 0x10(%rsp,%rax,4),%rbp
mov    %ecx,0x10(%rsp,%rax,4)
add    $0x1,%eax
cmp    %r14d,%eax
cmovge %r12d,%eax
sub    %rbp,%rcx
xorps  %xmm0,%xmm0
cvtsi2sd %rcx,%xmm0
addsd  %xmm0,%xmm4
movapd %xmm4,%xmm6
mulsd  %xmm7,%xmm6
mulsd  %xmm1,%xmm0
addsd  %xmm6,%xmm0
addsd  %xmm0,%xmm3
movapd %xmm3,%xmm0
divsd  %xmm2,%xmm0
cvttsd2si %xmm0,%ecx
maxsd  %xmm5,%xmm0
mov    %ecx,(%rbx,%rsi,4)
add    $0x1,%rsi
movapd %xmm0,%xmm5
cmp    %rsi,%rdx
jne    11b0 <trapfilter+0xa0>
jmp    123e <trapfilter+0x12e>
xorpd  %xmm0,%xmm0
add    $0x20018,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define double @trapfilter(i32* nocapture %0, i32 %1, i32 %2, i32 %3, i32 %4, i32 %5, double %6, i32 %7, i32 %8) {
%10 = alloca [16384 x i32]
%11 = bitcast [16384 x i32]* %10 to i8*
%12 = alloca [16384 x i32]
%13 = bitcast [16384 x i32]* %12 to i8*
%14 = add nsw i32 %2, %1
%15 = icmp slt i32 %5, %4
%16 = select i1 %15, i32 %5, i32 %4
%17 = icmp sgt i32 %16, 0
br i1 %17, label %18, label %62
18: ; preds = %9
%19 = sitofp i32 %3 to double
%20 = sitofp i32 %1 to double
%21 = zext i32 %16 to i64
br label %22
22: ; preds = %18, %22
%23 = phi i64 [ 0, %18 ], [ %60, %22 ]
%24 = phi i32 [ 0, %18 ], [ %39, %22 ]
%25 = phi i32 [ 0, %18 ], [ %49, %22 ]
%26 = phi double [ 0.000000e+00, %18 ], [ %55, %22 ]
%27 = phi double [ 0.000000e+00, %18 ], [ %52, %22 ]
%28 = phi double [ 0.000000e+00, %18 ], [ %58, %22 ]
%29 = getelementptr inbounds i32, i32* %0, i64 %23
%30 = load i32, i32* %29
%31 = sub nsw i32 %30, %7
%32 = mul nsw i32 %31, %8
%33 = sext i32 %24 to i64
%34 = getelementptr inbounds [16384 x i32], [16384 x i32]* %10, i64 0, i64 %33
%35 = load i32, i32* %34
%36 = sext i32 %35 to i64
store i32 %32, i32* %34
%37 = add nsw i32 %24, 1
%38 = icmp slt i32 %37, %1
%39 = select i1 %38, i32 %37, i32 0
%40 = sext i32 %32 to i64
%41 = sub nsw i64 %40, %36
%42 = sext i32 %25 to i64
%43 = getelementptr inbounds [16384 x i32], [16384 x i32]* %12, i64 0, i64 %42
%44 = load i32, i32* %43
%45 = sext i32 %44 to i64
%46 = trunc i64 %41 to i32
store i32 %46, i32* %43
%47 = add nsw i32 %25, 1
%48 = icmp slt i32 %47, %14
%49 = select i1 %48, i32 %47, i32 0
%50 = sub nsw i64 %41, %45
%51 = sitofp i64 %50 to double
%52 = fadd double %27, %51
%53 = fmul double %52, %6
%54 = tail call double @llvm.fmuladd.f64(double %51, double %19, double %53)
%55 = fadd double %26, %54
%56 = fdiv double %55, %20
%57 = fcmp ogt double %56, %28
%58 = select i1 %57, double %56, double %28
%59 = fptosi double %56 to i32
store i32 %59, i32* %29
%60 = add nuw nsw i64 %23, 1
%61 = icmp eq i64 %60, %21
br i1 %61, label %62, label %22
62: ; preds = %22, %9
%63 = phi double [ 0.000000e+00, %9 ], [ %58, %22 ]
ret double %63
}