#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 */
float fabsf(float);
float sqrtf(float);
void _sigfil(float *ima, short *bpm, float nsigma, long stride, long npix) {
  long i, k, n;
  float avg1, avg2, val, std;
  for (i = 0; i < npix; i++) {
    avg1 = avg2 = 0.0f;
    n = 0;
    k = i - 1;
    if ((k >= 0) && (k < npix)) {
      avg1 += ima[k];
      avg2 += ima[k] * ima[k];
      n++;
    }
    k = i + 1;
    if ((k >= 0) && (k < npix)) {
      avg1 += ima[k];
      avg2 += ima[k] * ima[k];
      n++;
    }
    k = i - stride;
    if ((k >= 0) && (k < npix)) {
      avg1 += ima[k];
      avg2 += ima[k] * ima[k];
      n++;
    }
    k = i - stride - 1;
    if ((k >= 0) && (k < npix)) {
      avg1 += ima[k];
      avg2 += ima[k] * ima[k];
      n++;
    }
    k = i - stride + 1;
    if ((k >= 0) && (k < npix)) {
      avg1 += ima[k];
      avg2 += ima[k] * ima[k];
      n++;
    }
    k = i + stride;
    if ((k >= 0) && (k < npix)) {
      avg1 += ima[k];
      avg2 += ima[k] * ima[k];
      n++;
    }
    k = i + stride - 1;
    if ((k >= 0) && (k < npix)) {
      avg1 += ima[k];
      avg2 += ima[k] * ima[k];
      n++;
    }
    k = i + stride + 1;
    if ((k >= 0) && (k < npix)) {
      avg1 += ima[k];
      avg2 += ima[k] * ima[k];
      n++;
    }
    avg1 = avg1 / (float)n;
    avg2 = avg2 / (float)n;
    std = avg2 - avg1 * avg1;
    if (std < 0.)
      std = 0.0f;
    else
      std = sqrtf(std);
    val = fabsf(ima[i] - avg1);
    if (val > (nsigma * std))
      bpm[i] = 1;
    else
      bpm[i] = 0;
  }
}


x86 assembly (clang, -O0):
<_sigfil>:
push   %rbp
mov    %rsp,%rbp
sub    $0x50,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
movss  %xmm0,-0x14(%rbp)
mov    %rdx,-0x20(%rbp)
mov    %rcx,-0x28(%rbp)
movq   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%rax
cmp    -0x28(%rbp),%rax
jge    15ff <_sigfil+0x4ef>
xorps  %xmm0,%xmm0
movss  %xmm0,-0x48(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x44(%rbp)
movq   $0x0,-0x40(%rbp)
mov    -0x30(%rbp),%rax
sub    $0x1,%rax
mov    %rax,-0x38(%rbp)
cmpq   $0x0,-0x38(%rbp)
jl     11cf <_sigfil+0xbf>
mov    -0x38(%rbp),%rax
cmp    -0x28(%rbp),%rax
jge    11cf <_sigfil+0xbf>
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
addss  -0x44(%rbp),%xmm0
movss  %xmm0,-0x44(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
movss  -0x48(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x48(%rbp)
mov    -0x40(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x40(%rbp)
mov    -0x30(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
cmpq   $0x0,-0x38(%rbp)
jl     1243 <_sigfil+0x133>
mov    -0x38(%rbp),%rax
cmp    -0x28(%rbp),%rax
jge    1243 <_sigfil+0x133>
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
addss  -0x44(%rbp),%xmm0
movss  %xmm0,-0x44(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
movss  -0x48(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x48(%rbp)
mov    -0x40(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x40(%rbp)
mov    -0x30(%rbp),%rax
sub    -0x20(%rbp),%rax
mov    %rax,-0x38(%rbp)
cmpq   $0x0,-0x38(%rbp)
jl     12b7 <_sigfil+0x1a7>
mov    -0x38(%rbp),%rax
cmp    -0x28(%rbp),%rax
jge    12b7 <_sigfil+0x1a7>
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
addss  -0x44(%rbp),%xmm0
movss  %xmm0,-0x44(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
movss  -0x48(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x48(%rbp)
mov    -0x40(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x40(%rbp)
mov    -0x30(%rbp),%rax
sub    -0x20(%rbp),%rax
sub    $0x1,%rax
mov    %rax,-0x38(%rbp)
cmpq   $0x0,-0x38(%rbp)
jl     132f <_sigfil+0x21f>
mov    -0x38(%rbp),%rax
cmp    -0x28(%rbp),%rax
jge    132f <_sigfil+0x21f>
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
addss  -0x44(%rbp),%xmm0
movss  %xmm0,-0x44(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
movss  -0x48(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x48(%rbp)
mov    -0x40(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x40(%rbp)
mov    -0x30(%rbp),%rax
sub    -0x20(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
cmpq   $0x0,-0x38(%rbp)
jl     13a7 <_sigfil+0x297>
mov    -0x38(%rbp),%rax
cmp    -0x28(%rbp),%rax
jge    13a7 <_sigfil+0x297>
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
addss  -0x44(%rbp),%xmm0
movss  %xmm0,-0x44(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
movss  -0x48(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x48(%rbp)
mov    -0x40(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x40(%rbp)
mov    -0x30(%rbp),%rax
add    -0x20(%rbp),%rax
mov    %rax,-0x38(%rbp)
cmpq   $0x0,-0x38(%rbp)
jl     141b <_sigfil+0x30b>
mov    -0x38(%rbp),%rax
cmp    -0x28(%rbp),%rax
jge    141b <_sigfil+0x30b>
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
addss  -0x44(%rbp),%xmm0
movss  %xmm0,-0x44(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
movss  -0x48(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x48(%rbp)
mov    -0x40(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x40(%rbp)
mov    -0x30(%rbp),%rax
add    -0x20(%rbp),%rax
sub    $0x1,%rax
mov    %rax,-0x38(%rbp)
cmpq   $0x0,-0x38(%rbp)
jl     1493 <_sigfil+0x383>
mov    -0x38(%rbp),%rax
cmp    -0x28(%rbp),%rax
jge    1493 <_sigfil+0x383>
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
addss  -0x44(%rbp),%xmm0
movss  %xmm0,-0x44(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
movss  -0x48(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x48(%rbp)
mov    -0x40(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x40(%rbp)
mov    -0x30(%rbp),%rax
add    -0x20(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x38(%rbp)
cmpq   $0x0,-0x38(%rbp)
jl     150b <_sigfil+0x3fb>
mov    -0x38(%rbp),%rax
cmp    -0x28(%rbp),%rax
jge    150b <_sigfil+0x3fb>
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
addss  -0x44(%rbp),%xmm0
movss  %xmm0,-0x44(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x8(%rbp),%rax
mov    -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
movss  -0x48(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x48(%rbp)
mov    -0x40(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x40(%rbp)
movss  -0x44(%rbp),%xmm0
cvtsi2ssq -0x40(%rbp),%xmm1
divss  %xmm1,%xmm0
movss  %xmm0,-0x44(%rbp)
movss  -0x48(%rbp),%xmm0
cvtsi2ssq -0x40(%rbp),%xmm1
divss  %xmm1,%xmm0
movss  %xmm0,-0x48(%rbp)
movss  -0x48(%rbp),%xmm1
movss  -0x44(%rbp),%xmm0
movss  -0x44(%rbp),%xmm2
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x50(%rbp)
movss  -0x50(%rbp),%xmm0
cvtss2sd %xmm0,%xmm1
xorps  %xmm0,%xmm0
ucomisd %xmm1,%xmm0
jbe    157f <_sigfil+0x46f>
xorps  %xmm0,%xmm0
movss  %xmm0,-0x50(%rbp)
jmp    158e <_sigfil+0x47e>
movss  -0x50(%rbp),%xmm0
call   1030 <sqrtf@plt>
movss  %xmm0,-0x50(%rbp)
mov    -0x8(%rbp),%rax
mov    -0x30(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
subss  -0x44(%rbp),%xmm0
movaps 0xa59(%rip),%xmm1
pand   %xmm1,%xmm0
movss  %xmm0,-0x4c(%rbp)
movss  -0x4c(%rbp),%xmm0
movss  -0x14(%rbp),%xmm1
mulss  -0x50(%rbp),%xmm1
ucomiss %xmm1,%xmm0
jbe    15db <_sigfil+0x4cb>
mov    -0x10(%rbp),%rax
mov    -0x30(%rbp),%rcx
movw   $0x1,(%rax,%rcx,2)
jmp    15e9 <_sigfil+0x4d9>
mov    -0x10(%rbp),%rax
mov    -0x30(%rbp),%rcx
movw   $0x0,(%rax,%rcx,2)
jmp    15ee <_sigfil+0x4de>
mov    -0x30(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x30(%rbp)
jmp    1135 <_sigfil+0x25>
add    $0x50,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<_sigfil>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x28,%rsp
test   %rcx,%rcx
jle    1357 <_sigfil+0x247>
mov    %rdx,%r13
movaps %xmm0,%xmm2
mov    %rcx,%r15
neg    %r15
mov    %rdx,%r9
neg    %r9
lea    0x0(,%rdx,4),%rax
mov    %rdi,%r12
sub    %rax,%r12
lea    (%rdi,%rdx,4),%r14
mov    $0x1,%ebp
movaps 0xea9(%rip),%xmm3
xorps  %xmm4,%xmm4
movss  %xmm0,0x8(%rsp)
mov    %rcx,0x20(%rsp)
mov    %r9,0x18(%rsp)
mov    %r14,0x10(%rsp)
jmp    11f3 <_sigfil+0xe3>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
mov    %rsi,%r14
mov    %rdi,%rbx
movss  %xmm5,0xc(%rsp)
call   1030 <sqrtf@plt>
movss  0xc(%rsp),%xmm5
xorps  %xmm4,%xmm4
movaps 0xe5f(%rip),%xmm3
mov    0x20(%rsp),%rcx
mov    %rbx,%rdi
mov    0x18(%rsp),%r9
mov    %r14,%rsi
mov    0x10(%rsp),%r14
movss  0x8(%rsp),%xmm2
movaps %xmm0,%xmm1
movss  -0x4(%rdi,%rbp,4),%xmm0
subss  %xmm5,%xmm0
andps  %xmm3,%xmm0
mulss  %xmm2,%xmm1
xor    %eax,%eax
ucomiss %xmm1,%xmm0
seta   %al
mov    %ax,-0x2(%rsi,%rbp,2)
lea    (%r15,%rbp,1),%rax
add    $0x1,%rax
add    $0x1,%rbp
cmp    $0x1,%rax
je     1357 <_sigfil+0x247>
cmp    $0x1,%rbp
jne    1210 <_sigfil+0x100>
xor    %eax,%eax
xorps  %xmm5,%xmm5
xorps  %xmm0,%xmm0
cmp    %rcx,%rbp
jl     122f <_sigfil+0x11f>
jmp    1244 <_sigfil+0x134>
nopl   0x0(%rax,%rax,1)
movss  -0x8(%rdi,%rbp,4),%xmm0
movaps %xmm0,%xmm5
addss  %xmm4,%xmm5
mulss  %xmm0,%xmm0
addss  %xmm4,%xmm0
mov    $0x1,%eax
cmp    %rcx,%rbp
jge    1244 <_sigfil+0x134>
movss  (%rdi,%rbp,4),%xmm1
addss  %xmm1,%xmm5
mulss  %xmm1,%xmm1
addss  %xmm1,%xmm0
add    $0x1,%rax
lea    (%r9,%rbp,1),%r10
mov    %r10,%rbx
add    $0xffffffffffffffff,%rbx
js     126d <_sigfil+0x15d>
cmp    %rcx,%rbx
jge    126d <_sigfil+0x15d>
movss  -0x4(%r12,%rbp,4),%xmm1
addss  %xmm1,%xmm5
mulss  %xmm1,%xmm1
addss  %xmm1,%xmm0
add    $0x1,%rax
test   %rbx,%rbx
jle    128e <_sigfil+0x17e>
cmp    %rcx,%rbx
jg     128e <_sigfil+0x17e>
movss  -0x8(%r12,%rbp,4),%xmm1
addss  %xmm1,%xmm5
mulss  %xmm1,%xmm1
addss  %xmm1,%xmm0
add    $0x1,%rax
cmp    $0xffffffffffffffff,%rbx
jl     12af <_sigfil+0x19f>
cmp    %rcx,%r10
jge    12af <_sigfil+0x19f>
movss  (%r12,%rbp,4),%xmm1
addss  %xmm1,%xmm5
mulss  %xmm1,%xmm1
addss  %xmm1,%xmm0
add    $0x1,%rax
mov    %r13,%r10
add    %rbp,%r10
mov    %r10,%rbx
add    $0xffffffffffffffff,%rbx
js     12da <_sigfil+0x1ca>
cmp    %rcx,%rbx
jge    12da <_sigfil+0x1ca>
movss  -0x4(%r14,%rbp,4),%xmm1
addss  %xmm1,%xmm5
mulss  %xmm1,%xmm1
addss  %xmm1,%xmm0
add    $0x1,%rax
test   %rbx,%rbx
jle    12fb <_sigfil+0x1eb>
cmp    %rcx,%rbx
jg     12fb <_sigfil+0x1eb>
movss  -0x8(%r14,%rbp,4),%xmm1
addss  %xmm1,%xmm5
mulss  %xmm1,%xmm1
addss  %xmm1,%xmm0
add    $0x1,%rax
cmp    $0xffffffffffffffff,%rbx
jl     131c <_sigfil+0x20c>
cmp    %rcx,%r10
jge    131c <_sigfil+0x20c>
movss  (%r14,%rbp,4),%xmm1
addss  %xmm1,%xmm5
mulss  %xmm1,%xmm1
addss  %xmm1,%xmm0
add    $0x1,%rax
xorps  %xmm1,%xmm1
cvtsi2ss %eax,%xmm1
divss  %xmm1,%xmm5
divss  %xmm1,%xmm0
movaps %xmm5,%xmm1
mulss  %xmm5,%xmm1
subss  %xmm1,%xmm0
xorps  %xmm1,%xmm1
ucomiss %xmm0,%xmm1
ja     11bf <_sigfil+0xaf>
ucomiss %xmm4,%xmm0
jb     1180 <_sigfil+0x70>
xorps  %xmm1,%xmm1
sqrtss %xmm0,%xmm1
jmp    11bf <_sigfil+0xaf>
add    $0x28,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @_sigfil(float* %0, i16* %1, float %2, i64 %3, i64 %4) {
%6 = alloca float*
%7 = alloca i16*
%8 = alloca float
%9 = alloca i64
%10 = alloca i64
%11 = alloca i64
%12 = alloca i64
%13 = alloca i64
%14 = alloca float
%15 = alloca float
%16 = alloca float
%17 = alloca float
store float* %0, float** %6
store i16* %1, i16** %7
store float %2, float* %8
store i64 %3, i64* %9
store i64 %4, i64* %10
store i64 0, i64* %11
br label %18
18: ; preds = %299, %5
%19 = load i64, i64* %11
%20 = load i64, i64* %10
%21 = icmp slt i64 %19, %20
br i1 %21, label %22, label %302
22: ; preds = %18
store float 0.000000e+00, float* %15
store float 0.000000e+00, float* %14
store i64 0, i64* %13
%23 = load i64, i64* %11
%24 = sub nsw i64 %23, 1
store i64 %24, i64* %12
%25 = load i64, i64* %12
%26 = icmp sge i64 %25, 0
br i1 %26, label %27, label %50
27: ; preds = %22
%28 = load i64, i64* %12
%29 = load i64, i64* %10
%30 = icmp slt i64 %28, %29
br i1 %30, label %31, label %50
31: ; preds = %27
%32 = load float*, float** %6
%33 = load i64, i64* %12
%34 = getelementptr inbounds float, float* %32, i64 %33
%35 = load float, float* %34
%36 = load float, float* %14
%37 = fadd float %36, %35
store float %37, float* %14
%38 = load float*, float** %6
%39 = load i64, i64* %12
%40 = getelementptr inbounds float, float* %38, i64 %39
%41 = load float, float* %40
%42 = load float*, float** %6
%43 = load i64, i64* %12
%44 = getelementptr inbounds float, float* %42, i64 %43
%45 = load float, float* %44
%46 = load float, float* %15
%47 = call float @llvm.fmuladd.f32(float %41, float %45, float %46)
store float %47, float* %15
%48 = load i64, i64* %13
%49 = add nsw i64 %48, 1
store i64 %49, i64* %13
br label %50
50: ; preds = %31, %27, %22
%51 = load i64, i64* %11
%52 = add nsw i64 %51, 1
store i64 %52, i64* %12
%53 = load i64, i64* %12
%54 = icmp sge i64 %53, 0
br i1 %54, label %55, label %78
55: ; preds = %50
%56 = load i64, i64* %12
%57 = load i64, i64* %10
%58 = icmp slt i64 %56, %57
br i1 %58, label %59, label %78
59: ; preds = %55
%60 = load float*, float** %6
%61 = load i64, i64* %12
%62 = getelementptr inbounds float, float* %60, i64 %61
%63 = load float, float* %62
%64 = load float, float* %14
%65 = fadd float %64, %63
store float %65, float* %14
%66 = load float*, float** %6
%67 = load i64, i64* %12
%68 = getelementptr inbounds float, float* %66, i64 %67
%69 = load float, float* %68
%70 = load float*, float** %6
%71 = load i64, i64* %12
%72 = getelementptr inbounds float, float* %70, i64 %71
%73 = load float, float* %72
%74 = load float, float* %15
%75 = call float @llvm.fmuladd.f32(float %69, float %73, float %74)
store float %75, float* %15
%76 = load i64, i64* %13
%77 = add nsw i64 %76, 1
store i64 %77, i64* %13
br label %78
78: ; preds = %59, %55, %50
%79 = load i64, i64* %11
%80 = load i64, i64* %9
%81 = sub nsw i64 %79, %80
store i64 %81, i64* %12
%82 = load i64, i64* %12
%83 = icmp sge i64 %82, 0
br i1 %83, label %84, label %107
84: ; preds = %78
%85 = load i64, i64* %12
%86 = load i64, i64* %10
%87 = icmp slt i64 %85, %86
br i1 %87, label %88, label %107
88: ; preds = %84
%89 = load float*, float** %6
%90 = load i64, i64* %12
%91 = getelementptr inbounds float, float* %89, i64 %90
%92 = load float, float* %91
%93 = load float, float* %14
%94 = fadd float %93, %92
store float %94, float* %14
%95 = load float*, float** %6
%96 = load i64, i64* %12
%97 = getelementptr inbounds float, float* %95, i64 %96
%98 = load float, float* %97
%99 = load float*, float** %6
%100 = load i64, i64* %12
%101 = getelementptr inbounds float, float* %99, i64 %100
%102 = load float, float* %101
%103 = load float, float* %15
%104 = call float @llvm.fmuladd.f32(float %98, float %102, float %103)
store float %104, float* %15
%105 = load i64, i64* %13
%106 = add nsw i64 %105, 1
store i64 %106, i64* %13
br label %107
107: ; preds = %88, %84, %78
%108 = load i64, i64* %11
%109 = load i64, i64* %9
%110 = sub nsw i64 %108, %109
%111 = sub nsw i64 %110, 1
store i64 %111, i64* %12
%112 = load i64, i64* %12
%113 = icmp sge i64 %112, 0
br i1 %113, label %114, label %137
114: ; preds = %107
%115 = load i64, i64* %12
%116 = load i64, i64* %10
%117 = icmp slt i64 %115, %116
br i1 %117, label %118, label %137
118: ; preds = %114
%119 = load float*, float** %6
%120 = load i64, i64* %12
%121 = getelementptr inbounds float, float* %119, i64 %120
%122 = load float, float* %121
%123 = load float, float* %14
%124 = fadd float %123, %122
store float %124, float* %14
%125 = load float*, float** %6
%126 = load i64, i64* %12
%127 = getelementptr inbounds float, float* %125, i64 %126
%128 = load float, float* %127
%129 = load float*, float** %6
%130 = load i64, i64* %12
%131 = getelementptr inbounds float, float* %129, i64 %130
%132 = load float, float* %131
%133 = load float, float* %15
%134 = call float @llvm.fmuladd.f32(float %128, float %132, float %133)
store float %134, float* %15
%135 = load i64, i64* %13
%136 = add nsw i64 %135, 1
store i64 %136, i64* %13
br label %137
137: ; preds = %118, %114, %107
%138 = load i64, i64* %11
%139 = load i64, i64* %9
%140 = sub nsw i64 %138, %139
%141 = add nsw i64 %140, 1
store i64 %141, i64* %12
%142 = load i64, i64* %12
%143 = icmp sge i64 %142, 0
br i1 %143, label %144, label %167
144: ; preds = %137
%145 = load i64, i64* %12
%146 = load i64, i64* %10
%147 = icmp slt i64 %145, %146
br i1 %147, label %148, label %167
148: ; preds = %144
%149 = load float*, float** %6
%150 = load i64, i64* %12
%151 = getelementptr inbounds float, float* %149, i64 %150
%152 = load float, float* %151
%153 = load float, float* %14
%154 = fadd float %153, %152
store float %154, float* %14
%155 = load float*, float** %6
%156 = load i64, i64* %12
%157 = getelementptr inbounds float, float* %155, i64 %156
%158 = load float, float* %157
%159 = load float*, float** %6
%160 = load i64, i64* %12
%161 = getelementptr inbounds float, float* %159, i64 %160
%162 = load float, float* %161
%163 = load float, float* %15
%164 = call float @llvm.fmuladd.f32(float %158, float %162, float %163)
store float %164, float* %15
%165 = load i64, i64* %13
%166 = add nsw i64 %165, 1
store i64 %166, i64* %13
br label %167
167: ; preds = %148, %144, %137
%168 = load i64, i64* %11
%169 = load i64, i64* %9
%170 = add nsw i64 %168, %169
store i64 %170, i64* %12
%171 = load i64, i64* %12
%172 = icmp sge i64 %171, 0
br i1 %172, label %173, label %196
173: ; preds = %167
%174 = load i64, i64* %12
%175 = load i64, i64* %10
%176 = icmp slt i64 %174, %175
br i1 %176, label %177, label %196
177: ; preds = %173
%178 = load float*, float** %6
%179 = load i64, i64* %12
%180 = getelementptr inbounds float, float* %178, i64 %179
%181 = load float, float* %180
%182 = load float, float* %14
%183 = fadd float %182, %181
store float %183, float* %14
%184 = load float*, float** %6
%185 = load i64, i64* %12
%186 = getelementptr inbounds float, float* %184, i64 %185
%187 = load float, float* %186
%188 = load float*, float** %6
%189 = load i64, i64* %12
%190 = getelementptr inbounds float, float* %188, i64 %189
%191 = load float, float* %190
%192 = load float, float* %15
%193 = call float @llvm.fmuladd.f32(float %187, float %191, float %192)
store float %193, float* %15
%194 = load i64, i64* %13
%195 = add nsw i64 %194, 1
store i64 %195, i64* %13
br label %196
196: ; preds = %177, %173, %167
%197 = load i64, i64* %11
%198 = load i64, i64* %9
%199 = add nsw i64 %197, %198
%200 = sub nsw i64 %199, 1
store i64 %200, i64* %12
%201 = load i64, i64* %12
%202 = icmp sge i64 %201, 0
br i1 %202, label %203, label %226
203: ; preds = %196
%204 = load i64, i64* %12
%205 = load i64, i64* %10
%206 = icmp slt i64 %204, %205
br i1 %206, label %207, label %226
207: ; preds = %203
%208 = load float*, float** %6
%209 = load i64, i64* %12
%210 = getelementptr inbounds float, float* %208, i64 %209
%211 = load float, float* %210
%212 = load float, float* %14
%213 = fadd float %212, %211
store float %213, float* %14
%214 = load float*, float** %6
%215 = load i64, i64* %12
%216 = getelementptr inbounds float, float* %214, i64 %215
%217 = load float, float* %216
%218 = load float*, float** %6
%219 = load i64, i64* %12
%220 = getelementptr inbounds float, float* %218, i64 %219
%221 = load float, float* %220
%222 = load float, float* %15
%223 = call float @llvm.fmuladd.f32(float %217, float %221, float %222)
store float %223, float* %15
%224 = load i64, i64* %13
%225 = add nsw i64 %224, 1
store i64 %225, i64* %13
br label %226
226: ; preds = %207, %203, %196
%227 = load i64, i64* %11
%228 = load i64, i64* %9
%229 = add nsw i64 %227, %228
%230 = add nsw i64 %229, 1
store i64 %230, i64* %12
%231 = load i64, i64* %12
%232 = icmp sge i64 %231, 0
br i1 %232, label %233, label %256
233: ; preds = %226
%234 = load i64, i64* %12
%235 = load i64, i64* %10
%236 = icmp slt i64 %234, %235
br i1 %236, label %237, label %256
237: ; preds = %233
%238 = load float*, float** %6
%239 = load i64, i64* %12
%240 = getelementptr inbounds float, float* %238, i64 %239
%241 = load float, float* %240
%242 = load float, float* %14
%243 = fadd float %242, %241
store float %243, float* %14
%244 = load float*, float** %6
%245 = load i64, i64* %12
%246 = getelementptr inbounds float, float* %244, i64 %245
%247 = load float, float* %246
%248 = load float*, float** %6
%249 = load i64, i64* %12
%250 = getelementptr inbounds float, float* %248, i64 %249
%251 = load float, float* %250
%252 = load float, float* %15
%253 = call float @llvm.fmuladd.f32(float %247, float %251, float %252)
store float %253, float* %15
%254 = load i64, i64* %13
%255 = add nsw i64 %254, 1
store i64 %255, i64* %13
br label %256
256: ; preds = %237, %233, %226
%257 = load float, float* %14
%258 = load i64, i64* %13
%259 = sitofp i64 %258 to float
%260 = fdiv float %257, %259
store float %260, float* %14
%261 = load float, float* %15
%262 = load i64, i64* %13
%263 = sitofp i64 %262 to float
%264 = fdiv float %261, %263
store float %264, float* %15
%265 = load float, float* %15
%266 = load float, float* %14
%267 = load float, float* %14
%268 = fneg float %266
%269 = call float @llvm.fmuladd.f32(float %268, float %267, float %265)
store float %269, float* %17
%270 = load float, float* %17
%271 = fpext float %270 to double
%272 = fcmp olt double %271, 0.000000e+00
br i1 %272, label %273, label %274
273: ; preds = %256
store float 0.000000e+00, float* %17
br label %277
274: ; preds = %256
%275 = load float, float* %17
%276 = call float @sqrtf(float %275)
store float %276, float* %17
br label %277
277: ; preds = %274, %273
%278 = load float*, float** %6
%279 = load i64, i64* %11
%280 = getelementptr inbounds float, float* %278, i64 %279
%281 = load float, float* %280
%282 = load float, float* %14
%283 = fsub float %281, %282
%284 = call float @llvm.fabs.f32(float %283)
store float %284, float* %16
%285 = load float, float* %16
%286 = load float, float* %8
%287 = load float, float* %17
%288 = fmul float %286, %287
%289 = fcmp ogt float %285, %288
br i1 %289, label %290, label %294
290: ; preds = %277
%291 = load i16*, i16** %7
%292 = load i64, i64* %11
%293 = getelementptr inbounds i16, i16* %291, i64 %292
store i16 1, i16* %293
br label %298
294: ; preds = %277
%295 = load i16*, i16** %7
%296 = load i64, i64* %11
%297 = getelementptr inbounds i16, i16* %295, i64 %296
store i16 0, i16* %297
br label %298
298: ; preds = %294, %290
br label %299
299: ; preds = %298
%300 = load i64, i64* %11
%301 = add nsw i64 %300, 1
store i64 %301, i64* %11
br label %18
302: ; preds = %18
ret void
}

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