#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
double sqrt(double);
void lsdf(double *data, int size, int fwhm, double f, double A, double M,
          double ratio) {
  int channel, i, j;
  double L, R, S;
  int width;
  double dhelp;
  width = (int)(f * fwhm);
  for (channel = width; channel < (size - width); channel++) {
    i = width;
    while (i > 0) {
      L = 0;
      R = 0;
      for (j = channel - i; j < channel; j++) {
        L += data[j];
      }
      for (j = channel + 1; j < channel + i; j++) {
        R += data[j];
      }
      S = data[channel] + L + R;
      if (S < M) {
        data[channel] = S / (2 * i + 1);
        break;
      }
      dhelp = (R + 1) / (L + 1);
      if ((dhelp < ratio) && (dhelp > (1 / ratio))) {
        if (S < (A * sqrt(data[channel]))) {
          data[channel] = S / (2 * i + 1);
          break;
        }
      }
      i = i - 1;
    }
  }
}


x86 assembly (clang, -O0):
<lsdf>:
push   %rbp
mov    %rsp,%rbp
sub    $0x80,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
movsd  %xmm0,-0x18(%rbp)
movsd  %xmm1,-0x20(%rbp)
movsd  %xmm2,-0x28(%rbp)
movsd  %xmm3,-0x30(%rbp)
movsd  -0x18(%rbp),%xmm0
cvtsi2sdl -0x10(%rbp),%xmm1
mulsd  %xmm1,%xmm0
cvttsd2si %xmm0,%eax
mov    %eax,-0x5c(%rbp)
mov    -0x5c(%rbp),%eax
mov    %eax,-0x34(%rbp)
mov    -0x34(%rbp),%eax
mov    -0xc(%rbp),%ecx
sub    -0x5c(%rbp),%ecx
cmp    %ecx,%eax
jge    1338 <lsdf+0x228>
mov    -0x5c(%rbp),%eax
mov    %eax,-0x38(%rbp)
cmpl   $0x0,-0x38(%rbp)
jle    1325 <lsdf+0x215>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x48(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x34(%rbp),%eax
sub    -0x38(%rbp),%eax
mov    %eax,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x34(%rbp),%eax
jge    11bf <lsdf+0xaf>
mov    -0x8(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    118e <lsdf+0x7e>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
mov    -0x34(%rbp),%ecx
add    -0x38(%rbp),%ecx
cmp    %ecx,%eax
jge    11fe <lsdf+0xee>
mov    -0x8(%rbp),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    11c8 <lsdf+0xb8>
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
addsd  -0x48(%rbp),%xmm0
addsd  -0x50(%rbp),%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  -0x58(%rbp),%xmm1
movsd  -0x28(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    1256 <lsdf+0x146>
movsd  -0x58(%rbp),%xmm0
mov    -0x38(%rbp),%eax
shl    $0x1,%eax
add    $0x1,%eax
cvtsi2sd %eax,%xmm1
divsd  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    1325 <lsdf+0x215>
movsd  0xda2(%rip),%xmm0
addsd  -0x50(%rbp),%xmm0
movsd  0xd95(%rip),%xmm1
addsd  -0x48(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  -0x68(%rbp),%xmm1
movsd  -0x30(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    1317 <lsdf+0x207>
movsd  -0x68(%rbp),%xmm0
movsd  0xd66(%rip),%xmm1
divsd  -0x30(%rbp),%xmm1
ucomisd %xmm1,%xmm0
jbe    1317 <lsdf+0x207>
movsd  -0x58(%rbp),%xmm0
movsd  %xmm0,-0x78(%rbp)
movsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
call   1030 <sqrt@plt>
movsd  -0x78(%rbp),%xmm1
movaps %xmm0,%xmm2
movsd  -0x70(%rbp),%xmm0
mulsd  %xmm2,%xmm0
ucomisd %xmm1,%xmm0
jbe    1312 <lsdf+0x202>
movsd  -0x58(%rbp),%xmm0
mov    -0x38(%rbp),%eax
shl    $0x1,%eax
add    $0x1,%eax
cvtsi2sd %eax,%xmm1
divsd  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    1325 <lsdf+0x215>
jmp    1317 <lsdf+0x207>
mov    -0x38(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    116b <lsdf+0x5b>
jmp    132a <lsdf+0x21a>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1154 <lsdf+0x44>
add    $0x80,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<lsdf>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x48,%rsp
cvtsi2sd %edx,%xmm4
mulsd  %xmm0,%xmm4
cvttsd2si %xmm4,%ecx
mov    %esi,0xc(%rsp)
mov    %esi,%eax
sub    %ecx,%eax
mov    %rcx,0x10(%rsp)
cmp    %ecx,%eax
jle    12bf <lsdf+0x1af>
mov    %rdi,%rbx
movsd  0xeb6(%rip),%xmm6
movapd %xmm6,%xmm9
divsd  %xmm3,%xmm9
mov    0x10(%rsp),%rax
movslq %eax,%r15
add    %eax,%eax
sub    %eax,0xc(%rsp)
xor    %r12d,%r12d
xorpd  %xmm8,%xmm8
mov    %r15,0x18(%rsp)
movsd  %xmm3,0x38(%rsp)
movsd  %xmm2,0x30(%rsp)
movsd  %xmm1,0x28(%rsp)
movsd  %xmm9,0x20(%rsp)
jmp    11b7 <lsdf+0xa7>
nopw   0x0(%rax,%rax,1)
lea    (%r14,%r14,1),%eax
add    $0x1,%eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
divsd  %xmm0,%xmm7
movsd  %xmm7,(%rbx,%r13,8)
add    $0x1,%r12d
cmp    0xc(%rsp),%r12d
je     12bf <lsdf+0x1af>
mov    %r15,%r13
add    $0x1,%r15
cmpl   $0x0,0x10(%rsp)
jle    11a8 <lsdf+0x98>
mov    0x18(%rsp),%r14
mov    %r12d,%ebp
jmp    1221 <lsdf+0x111>
nop
movsd  %xmm7,0x40(%rsp)
call   1030 <sqrt@plt>
movsd  0x40(%rsp),%xmm7
xorpd  %xmm8,%xmm8
movsd  0x20(%rsp),%xmm9
movsd  0xe0b(%rip),%xmm6
movsd  0x28(%rsp),%xmm1
movsd  0x30(%rsp),%xmm2
movsd  0x38(%rsp),%xmm3
mulsd  %xmm1,%xmm0
ucomisd %xmm7,%xmm0
ja     1190 <lsdf+0x80>
add    $0xffffffffffffffff,%r14
add    $0x1,%ebp
cmp    %r13d,%ebp
je     11a8 <lsdf+0x98>
movslq %ebp,%rax
xorpd  %xmm5,%xmm5
nopl   0x0(%rax,%rax,1)
addsd  (%rbx,%rax,8),%xmm5
add    $0x1,%rax
cmp    %r13,%rax
jl     1230 <lsdf+0x120>
xorpd  %xmm4,%xmm4
cmp    $0x2,%r14
jl     126e <lsdf+0x15e>
lea    (%r14,%r13,1),%rax
xorpd  %xmm4,%xmm4
mov    %r15,%rcx
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
addsd  (%rbx,%rcx,8),%xmm4
add    $0x1,%rcx
cmp    %rax,%rcx
jl     1260 <lsdf+0x150>
movsd  (%rbx,%r13,8),%xmm0
movapd %xmm5,%xmm7
addsd  %xmm0,%xmm7
addsd  %xmm4,%xmm7
ucomisd %xmm7,%xmm2
ja     1190 <lsdf+0x80>
addsd  %xmm6,%xmm4
addsd  %xmm6,%xmm5
divsd  %xmm5,%xmm4
ucomisd %xmm4,%xmm3
jbe    1215 <lsdf+0x105>
ucomisd %xmm9,%xmm4
jbe    1215 <lsdf+0x105>
ucomisd %xmm8,%xmm0
jb     11d0 <lsdf+0xc0>
sqrtsd %xmm0,%xmm0
jmp    1207 <lsdf+0xf7>
add    $0x48,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define void @lsdf(double* nocapture %0, i32 %1, i32 %2, double %3, double %4, double %5, double %6) {
%8 = sitofp i32 %2 to double
%9 = fmul double %8, %3
%10 = fptosi double %9 to i32
%11 = sub nsw i32 %1, %10
%12 = icmp sgt i32 %11, %10
br i1 %12, label %13, label %81
13: ; preds = %7
%14 = icmp sgt i32 %10, 0
%15 = fdiv double 1.000000e+00, %6
%16 = sext i32 %10 to i64
%17 = mul i32 %10, -2
%18 = add i32 %17, %1
br label %19
19: ; preds = %13, %77
%20 = phi i64 [ %16, %13 ], [ %78, %77 ]
%21 = phi i64 [ %16, %13 ], [ %23, %77 ]
%22 = phi i32 [ 0, %13 ], [ %79, %77 ]
%23 = add nsw i64 %21, 1
br i1 %14, label %24, label %77
24: ; preds = %19
%25 = getelementptr inbounds double, double* %0, i64 %20
%26 = trunc i64 %20 to i32
br label %27
27: ; preds = %67, %24
%28 = phi i64 [ %16, %24 ], [ %68, %67 ]
%29 = phi i32 [ %22, %24 ], [ %69, %67 ]
%30 = sext i32 %29 to i64
br label %34
31: ; preds = %34
%32 = add nsw i64 %28, %20
%33 = icmp sgt i64 %28, 1
br i1 %33, label %42, label %50
34: ; preds = %27, %34
%35 = phi i64 [ %30, %27 ], [ %40, %34 ]
%36 = phi double [ 0.000000e+00, %27 ], [ %39, %34 ]
%37 = getelementptr inbounds double, double* %0, i64 %35
%38 = load double, double* %37
%39 = fadd double %36, %38
%40 = add nsw i64 %35, 1
%41 = icmp slt i64 %40, %20
br i1 %41, label %34, label %31
42: ; preds = %31, %42
%43 = phi i64 [ %48, %42 ], [ %23, %31 ]
%44 = phi double [ %47, %42 ], [ 0.000000e+00, %31 ]
%45 = getelementptr inbounds double, double* %0, i64 %43
%46 = load double, double* %45
%47 = fadd double %44, %46
%48 = add nsw i64 %43, 1
%49 = icmp slt i64 %48, %32
br i1 %49, label %42, label %50
50: ; preds = %42, %31
%51 = phi double [ 0.000000e+00, %31 ], [ %47, %42 ]
%52 = load double, double* %25
%53 = fadd double %39, %52
%54 = fadd double %51, %53
%55 = fcmp olt double %54, %5
br i1 %55, label %71, label %56
56: ; preds = %50
%57 = fadd double %51, 1.000000e+00
%58 = fadd double %39, 1.000000e+00
%59 = fdiv double %57, %58
%60 = fcmp olt double %59, %6
%61 = fcmp ogt double %59, %15
%62 = and i1 %60, %61
br i1 %62, label %63, label %67
63: ; preds = %56
%64 = tail call double @sqrt(double %52)
%65 = fmul double %64, %4
%66 = fcmp olt double %54, %65
br i1 %66, label %71, label %67
67: ; preds = %63, %56
%68 = add nsw i64 %28, -1
%69 = add i32 %29, 1
%70 = icmp eq i32 %69, %26
br i1 %70, label %77, label %27
71: ; preds = %63, %50
%72 = trunc i64 %28 to i32
%73 = shl nsw i32 %72, 1
%74 = or i32 %73, 1
%75 = sitofp i32 %74 to double
%76 = fdiv double %54, %75
store double %76, double* %25
br label %77
77: ; preds = %67, %71, %19
%78 = add nsw i64 %20, 1
%79 = add i32 %22, 1
%80 = icmp eq i32 %79, %18
br i1 %80, label %81, label %19
81: ; preds = %77, %7
ret void
}