#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 */
scalar_t__ isnan(double);
double sqrt(int);
void zero_lc_averages(int Njd, double *b, double clipping) {
  int i, n;
  double ave1, ave2, std;
  long double val1, val2;
  val1 = 0.;
  val2 = 0.;
  n = 0;
  for (i = 0; i < Njd; i++) {
    if (!isnan(b[i])) {
      val1 += b[i] * b[i];
      val2 += b[i];
      n++;
    }
  }
  if (n > 0) {
    ave1 = (double)(val2 / (long double)n);
    std = sqrt((double)((val1 / (long double)n) -
                        (val2 * val2 / (long double)(n * n))));
  } else {
    ave1 = 0.;
    std = 0.;
  }
  val1 = 0.;
  val2 = 0.;
  n = 0;
  for (i = 0; i < Njd; i++) {
    if (!isnan(b[i]))
      if (((b[i] - ave1) > 0 ? (b[i] - ave1) : (-(b[i] - ave1))) <
          clipping * std) {
        val1 += b[i] * b[i];
        val2 += b[i];
        n++;
      }
  }
  if (n > 0)
    ave2 = (double)(val2 / (long double)n);
  else
    ave2 = 0.;
  for (i = 0; i < Njd; i++) {
    if (isnan(b[i]) || ((b[i] - ave1) > 0 ? (b[i] - ave1) : (-(b[i] - ave1))) >=
                           clipping * std) {
      b[i] = sqrt(-1.);
    } else {
      b[i] -= ave2;
    }
  }
}


x86 assembly (clang, -O0):
<zero_lc_averages>:
push   %rbp
mov    %rsp,%rbp
sub    $0xb0,%rsp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
movsd  %xmm0,-0x18(%rbp)
fldz
fld    %st(0)
fstpt  -0x50(%rbp)
fstpt  -0x60(%rbp)
movl   $0x0,-0x20(%rbp)
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    11cd <zero_lc_averages+0xad>
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
call   1030 <isnan@plt>
cmp    $0x0,%rax
jne    11ba <zero_lc_averages+0x9a>
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mulsd  %xmm0,%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
fldl   -0x98(%rbp)
fldt   -0x50(%rbp)
faddp  %st,%st(1)
fstpt  -0x50(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
fldl   (%rax,%rcx,8)
fldt   -0x60(%rbp)
faddp  %st,%st(1)
fstpt  -0x60(%rbp)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    11bf <zero_lc_averages+0x9f>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    114f <zero_lc_averages+0x2f>
cmpl   $0x0,-0x20(%rbp)
jle    1223 <zero_lc_averages+0x103>
fldt   -0x60(%rbp)
mov    -0x20(%rbp),%eax
mov    %eax,-0x6c(%rbp)
fildl  -0x6c(%rbp)
fdivrp %st,%st(1)
fstpl  -0x28(%rbp)
fldt   -0x50(%rbp)
mov    -0x20(%rbp),%eax
mov    %eax,-0x68(%rbp)
fildl  -0x68(%rbp)
fdivrp %st,%st(1)
fldt   -0x60(%rbp)
fmul   %st,%st(0)
imul   %eax,%eax
mov    %eax,-0x64(%rbp)
fildl  -0x64(%rbp)
fdivrp %st,%st(1)
fsubrp %st,%st(1)
fstpl  -0x78(%rbp)
movsd  -0x78(%rbp),%xmm0
cvttsd2si %xmm0,%edi
call   1040 <sqrt@plt>
movsd  %xmm0,-0x38(%rbp)
jmp    1233 <zero_lc_averages+0x113>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x28(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x38(%rbp)
fldz
fld    %st(0)
fstpt  -0x50(%rbp)
fstpt  -0x60(%rbp)
movl   $0x0,-0x20(%rbp)
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    1359 <zero_lc_averages+0x239>
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
call   1030 <isnan@plt>
cmp    $0x0,%rax
jne    1346 <zero_lc_averages+0x226>
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x28(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    12b1 <zero_lc_averages+0x191>
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0xa0(%rbp)
ff
jmp    12e2 <zero_lc_averages+0x1c2>
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x28(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  %xmm0,-0xa0(%rbp)
ff
movsd  -0xa0(%rbp),%xmm1
ff
movsd  -0x18(%rbp),%xmm0
mulsd  -0x38(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    1341 <zero_lc_averages+0x221>
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mulsd  %xmm0,%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
fldl   -0x90(%rbp)
fldt   -0x50(%rbp)
faddp  %st,%st(1)
fstpt  -0x50(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
fldl   (%rax,%rcx,8)
fldt   -0x60(%rbp)
faddp  %st,%st(1)
fstpt  -0x60(%rbp)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    1346 <zero_lc_averages+0x226>
jmp    134b <zero_lc_averages+0x22b>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    124b <zero_lc_averages+0x12b>
cmpl   $0x0,-0x20(%rbp)
jle    1389 <zero_lc_averages+0x269>
fldt   -0x60(%rbp)
mov    -0x20(%rbp),%eax
mov    %eax,-0x7c(%rbp)
fildl  -0x7c(%rbp)
fdivrp %st,%st(1)
fstpl  -0x88(%rbp)
movsd  -0x88(%rbp),%xmm0
ff
movsd  %xmm0,-0x30(%rbp)
jmp    1391 <zero_lc_averages+0x271>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x30(%rbp)
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    1495 <zero_lc_averages+0x375>
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
call   1030 <isnan@plt>
cmp    $0x0,%rax
jne    144b <zero_lc_averages+0x32b>
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x28(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jbe    13fe <zero_lc_averages+0x2de>
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0xa8(%rbp)
ff
jmp    142f <zero_lc_averages+0x30f>
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x28(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  %xmm0,-0xa8(%rbp)
ff
movsd  -0xa8(%rbp),%xmm0
ff
movsd  -0x18(%rbp),%xmm1
mulsd  -0x38(%rbp),%xmm1
ucomisd %xmm1,%xmm0
jb     1467 <zero_lc_averages+0x347>
mov    $0xffffffff,%edi
call   1040 <sqrt@plt>
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    1482 <zero_lc_averages+0x362>
movsd  -0x30(%rbp),%xmm1
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    1487 <zero_lc_averages+0x367>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    1398 <zero_lc_averages+0x278>
add    $0xb0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<zero_lc_averages>:
push   %rbp
push   %r15
push   %r14
push   %r12
push   %rbx
sub    $0x70,%rsp
movsd  %xmm0,0x38(%rsp)
test   %edi,%edi
jle    1346 <zero_lc_averages+0x226>
mov    %rsi,%rbx
mov    %edi,%r14d
mov    %edi,%r15d
xor    %ebp,%ebp
fldz
fld    %st(0)
fstpt  0x10(%rsp)
fstpt  (%rsp)
xor    %r12d,%r12d
jmp    1169 <zero_lc_averages+0x49>
cs nopw 0x0(%rax,%rax,1)
nop
add    $0x1,%rbp
cmp    %rbp,%r15
je     11ad <zero_lc_averages+0x8d>
movsd  (%rbx,%rbp,8),%xmm0
call   1030 <isnan@plt>
test   %rax,%rax
jne    1160 <zero_lc_averages+0x40>
movsd  (%rbx,%rbp,8),%xmm0
movapd %xmm0,%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,0x60(%rsp)
fldt   (%rsp)
faddl  0x60(%rsp)
fstpt  (%rsp)
movsd  %xmm0,0x68(%rsp)
fldt   0x10(%rsp)
faddl  0x68(%rsp)
fstpt  0x10(%rsp)
add    $0x1,%r12d
jmp    1160 <zero_lc_averages+0x40>
test   %r12d,%r12d
jle    11fa <zero_lc_averages+0xda>
mov    %r12d,0x30(%rsp)
fildl  0x30(%rsp)
fldt   0x10(%rsp)
fld    %st(0)
fdiv   %st(2),%st
fstpl  0x58(%rsp)
movsd  0x58(%rsp),%xmm0
movsd  %xmm0,0x20(%rsp)
fldt   (%rsp)
fdivp  %st,%st(2)
fmul   %st,%st(0)
imul   %r12d,%r12d
mov    %r12d,0x34(%rsp)
fidivl 0x34(%rsp)
fsubrp %st,%st(1)
fstpl  0x50(%rsp)
cvttsd2si 0x50(%rsp),%edi
call   1040 <sqrt@plt>
jmp    1208 <zero_lc_averages+0xe8>
xorpd  %xmm0,%xmm0
movsd  %xmm0,0x20(%rsp)
xorpd  %xmm0,%xmm0
movsd  %xmm0,0x10(%rsp)
test   %r14d,%r14d
jle    1346 <zero_lc_averages+0x226>
movsd  0x10(%rsp),%xmm0
mulsd  0x38(%rsp),%xmm0
movsd  %xmm0,0x10(%rsp)
xor    %ebp,%ebp
fldz
fstpt  (%rsp)
xor    %r12d,%r12d
jmp    1249 <zero_lc_averages+0x129>
cs nopw 0x0(%rax,%rax,1)
nop
add    $0x1,%rbp
cmp    %rbp,%r15
je     12aa <zero_lc_averages+0x18a>
movsd  (%rbx,%rbp,8),%xmm0
call   1030 <isnan@plt>
test   %rax,%rax
jne    1240 <zero_lc_averages+0x120>
movsd  (%rbx,%rbp,8),%xmm0
movapd %xmm0,%xmm1
subsd  0x20(%rsp),%xmm1
movapd %xmm1,%xmm2
xorpd  0xd8d(%rip),%xmm2
xorpd  %xmm3,%xmm3
cmpltsd %xmm1,%xmm3
andpd  %xmm3,%xmm1
andnpd %xmm2,%xmm3
orpd   %xmm1,%xmm3
movsd  0x10(%rsp),%xmm1
ucomisd %xmm3,%xmm1
jbe    1240 <zero_lc_averages+0x120>
movsd  %xmm0,0x48(%rsp)
fldt   (%rsp)
faddl  0x48(%rsp)
fstpt  (%rsp)
add    $0x1,%r12d
jmp    1240 <zero_lc_averages+0x120>
test   %r12d,%r12d
jle    12c7 <zero_lc_averages+0x1a7>
mov    %r12d,0x2c(%rsp)
fldt   (%rsp)
fidivl 0x2c(%rsp)
fstpl  0x40(%rsp)
movsd  0x40(%rsp),%xmm0
jmp    12cb <zero_lc_averages+0x1ab>
xorpd  %xmm0,%xmm0
movsd  %xmm0,(%rsp)
test   %r14d,%r14d
jle    1346 <zero_lc_averages+0x226>
xor    %ebp,%ebp
jmp    12f8 <zero_lc_averages+0x1d8>
nopl   0x0(%rax)
mov    $0xffffffff,%edi
call   1040 <sqrt@plt>
movsd  %xmm0,(%rbx,%rbp,8)
add    $0x1,%rbp
cmp    %rbp,%r15
je     1346 <zero_lc_averages+0x226>
movsd  (%rbx,%rbp,8),%xmm0
call   1030 <isnan@plt>
test   %rax,%rax
jne    12e0 <zero_lc_averages+0x1c0>
movsd  (%rbx,%rbp,8),%xmm0
movapd %xmm0,%xmm1
subsd  0x20(%rsp),%xmm1
movapd %xmm1,%xmm2
xorpd  0xcde(%rip),%xmm2
xorpd  %xmm3,%xmm3
cmpltsd %xmm1,%xmm3
andpd  %xmm3,%xmm1
andnpd %xmm2,%xmm3
orpd   %xmm1,%xmm3
ucomisd 0x10(%rsp),%xmm3
jae    12e0 <zero_lc_averages+0x1c0>
subsd  (%rsp),%xmm0
jmp    12ea <zero_lc_averages+0x1ca>
add    $0x70,%rsp
pop    %rbx
pop    %r12
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @zero_lc_averages(i32 %0, double* %1, double %2) {
%4 = alloca i32
%5 = alloca double*
%6 = alloca double
%7 = alloca i32
%8 = alloca i32
%9 = alloca double
%10 = alloca double
%11 = alloca double
%12 = alloca x86_fp80
%13 = alloca x86_fp80
store i32 %0, i32* %4
store double* %1, double** %5
store double %2, double* %6
store x86_fp80 0xK00000000000000000000, x86_fp80* %12
store x86_fp80 0xK00000000000000000000, x86_fp80* %13
store i32 0, i32* %8
store i32 0, i32* %7
br label %14
14: ; preds = %52, %3
%15 = load i32, i32* %7
%16 = load i32, i32* %4
%17 = icmp slt i32 %15, %16
br i1 %17, label %18, label %55
18: ; preds = %14
%19 = load double*, double** %5
%20 = load i32, i32* %7
%21 = sext i32 %20 to i64
%22 = getelementptr inbounds double, double* %19, i64 %21
%23 = load double, double* %22
%24 = call i64 @isnan(double %23)
%25 = icmp ne i64 %24, 0
br i1 %25, label %51, label %26
26: ; preds = %18
%27 = load double*, double** %5
%28 = load i32, i32* %7
%29 = sext i32 %28 to i64
%30 = getelementptr inbounds double, double* %27, i64 %29
%31 = load double, double* %30
%32 = load double*, double** %5
%33 = load i32, i32* %7
%34 = sext i32 %33 to i64
%35 = getelementptr inbounds double, double* %32, i64 %34
%36 = load double, double* %35
%37 = fmul double %31, %36
%38 = fpext double %37 to x86_fp80
%39 = load x86_fp80, x86_fp80* %12
%40 = fadd x86_fp80 %39, %38
store x86_fp80 %40, x86_fp80* %12
%41 = load double*, double** %5
%42 = load i32, i32* %7
%43 = sext i32 %42 to i64
%44 = getelementptr inbounds double, double* %41, i64 %43
%45 = load double, double* %44
%46 = fpext double %45 to x86_fp80
%47 = load x86_fp80, x86_fp80* %13
%48 = fadd x86_fp80 %47, %46
store x86_fp80 %48, x86_fp80* %13
%49 = load i32, i32* %8
%50 = add nsw i32 %49, 1
store i32 %50, i32* %8
br label %51
51: ; preds = %26, %18
br label %52
52: ; preds = %51
%53 = load i32, i32* %7
%54 = add nsw i32 %53, 1
store i32 %54, i32* %7
br label %14
55: ; preds = %14
%56 = load i32, i32* %8
%57 = icmp sgt i32 %56, 0
br i1 %57, label %58, label %80
58: ; preds = %55
%59 = load x86_fp80, x86_fp80* %13
%60 = load i32, i32* %8
%61 = sitofp i32 %60 to x86_fp80
%62 = fdiv x86_fp80 %59, %61
%63 = fptrunc x86_fp80 %62 to double
store double %63, double* %9
%64 = load x86_fp80, x86_fp80* %12
%65 = load i32, i32* %8
%66 = sitofp i32 %65 to x86_fp80
%67 = fdiv x86_fp80 %64, %66
%68 = load x86_fp80, x86_fp80* %13
%69 = load x86_fp80, x86_fp80* %13
%70 = fmul x86_fp80 %68, %69
%71 = load i32, i32* %8
%72 = load i32, i32* %8
%73 = mul nsw i32 %71, %72
%74 = sitofp i32 %73 to x86_fp80
%75 = fdiv x86_fp80 %70, %74
%76 = fsub x86_fp80 %67, %75
%77 = fptrunc x86_fp80 %76 to double
%78 = fptosi double %77 to i32
%79 = call double @sqrt(i32 %78)
store double %79, double* %11
br label %81
80: ; preds = %55
store double 0.000000e+00, double* %9
store double 0.000000e+00, double* %11
br label %81
81: ; preds = %80, %58
store x86_fp80 0xK00000000000000000000, x86_fp80* %12
store x86_fp80 0xK00000000000000000000, x86_fp80* %13
store i32 0, i32* %8
store i32 0, i32* %7
br label %82
82: ; preds = %153, %81
%83 = load i32, i32* %7
%84 = load i32, i32* %4
%85 = icmp slt i32 %83, %84
br i1 %85, label %86, label %156
86: ; preds = %82
%87 = load double*, double** %5
%88 = load i32, i32* %7
%89 = sext i32 %88 to i64
%90 = getelementptr inbounds double, double* %87, i64 %89
%91 = load double, double* %90
%92 = call i64 @isnan(double %91)
%93 = icmp ne i64 %92, 0
br i1 %93, label %152, label %94
94: ; preds = %86
%95 = load double*, double** %5
%96 = load i32, i32* %7
%97 = sext i32 %96 to i64
%98 = getelementptr inbounds double, double* %95, i64 %97
%99 = load double, double* %98
%100 = load double, double* %9
%101 = fsub double %99, %100
%102 = fcmp ogt double %101, 0.000000e+00
br i1 %102, label %103, label %111
103: ; preds = %94
%104 = load double*, double** %5
%105 = load i32, i32* %7
%106 = sext i32 %105 to i64
%107 = getelementptr inbounds double, double* %104, i64 %106
%108 = load double, double* %107
%109 = load double, double* %9
%110 = fsub double %108, %109
br label %120
111: ; preds = %94
%112 = load double*, double** %5
%113 = load i32, i32* %7
%114 = sext i32 %113 to i64
%115 = getelementptr inbounds double, double* %112, i64 %114
%116 = load double, double* %115
%117 = load double, double* %9
%118 = fsub double %116, %117
%119 = fneg double %118
br label %120
120: ; preds = %111, %103
%121 = phi double [ %110, %103 ], [ %119, %111 ]
%122 = load double, double* %6
%123 = load double, double* %11
%124 = fmul double %122, %123
%125 = fcmp olt double %121, %124
br i1 %125, label %126, label %151
126: ; preds = %120
%127 = load double*, double** %5
%128 = load i32, i32* %7
%129 = sext i32 %128 to i64
%130 = getelementptr inbounds double, double* %127, i64 %129
%131 = load double, double* %130
%132 = load double*, double** %5
%133 = load i32, i32* %7
%134 = sext i32 %133 to i64
%135 = getelementptr inbounds double, double* %132, i64 %134
%136 = load double, double* %135
%137 = fmul double %131, %136
%138 = fpext double %137 to x86_fp80
%139 = load x86_fp80, x86_fp80* %12
%140 = fadd x86_fp80 %139, %138
store x86_fp80 %140, x86_fp80* %12
%141 = load double*, double** %5
%142 = load i32, i32* %7
%143 = sext i32 %142 to i64
%144 = getelementptr inbounds double, double* %141, i64 %143
%145 = load double, double* %144
%146 = fpext double %145 to x86_fp80
%147 = load x86_fp80, x86_fp80* %13
%148 = fadd x86_fp80 %147, %146
store x86_fp80 %148, x86_fp80* %13
%149 = load i32, i32* %8
%150 = add nsw i32 %149, 1
store i32 %150, i32* %8
br label %151
151: ; preds = %126, %120
br label %152
152: ; preds = %151, %86
br label %153
153: ; preds = %152
%154 = load i32, i32* %7
%155 = add nsw i32 %154, 1
store i32 %155, i32* %7
br label %82
156: ; preds = %82
%157 = load i32, i32* %8
%158 = icmp sgt i32 %157, 0
br i1 %158, label %159, label %165
159: ; preds = %156
%160 = load x86_fp80, x86_fp80* %13
%161 = load i32, i32* %8
%162 = sitofp i32 %161 to x86_fp80
%163 = fdiv x86_fp80 %160, %162
%164 = fptrunc x86_fp80 %163 to double
store double %164, double* %10
br label %166
165: ; preds = %156
store double 0.000000e+00, double* %10
br label %166
166: ; preds = %165, %159
store i32 0, i32* %7
br label %167
167: ; preds = %226, %166
%168 = load i32, i32* %7
%169 = load i32, i32* %4
%170 = icmp slt i32 %168, %169
br i1 %170, label %171, label %229
171: ; preds = %167
%172 = load double*, double** %5
%173 = load i32, i32* %7
%174 = sext i32 %173 to i64
%175 = getelementptr inbounds double, double* %172, i64 %174
%176 = load double, double* %175
%177 = call i64 @isnan(double %176)
%178 = icmp ne i64 %177, 0
br i1 %178, label %211, label %179
179: ; preds = %171
%180 = load double*, double** %5
%181 = load i32, i32* %7
%182 = sext i32 %181 to i64
%183 = getelementptr inbounds double, double* %180, i64 %182
%184 = load double, double* %183
%185 = load double, double* %9
%186 = fsub double %184, %185
%187 = fcmp ogt double %186, 0.000000e+00
br i1 %187, label %188, label %196
188: ; preds = %179
%189 = load double*, double** %5
%190 = load i32, i32* %7
%191 = sext i32 %190 to i64
%192 = getelementptr inbounds double, double* %189, i64 %191
%193 = load double, double* %192
%194 = load double, double* %9
%195 = fsub double %193, %194
br label %205
196: ; preds = %179
%197 = load double*, double** %5
%198 = load i32, i32* %7
%199 = sext i32 %198 to i64
%200 = getelementptr inbounds double, double* %197, i64 %199
%201 = load double, double* %200
%202 = load double, double* %9
%203 = fsub double %201, %202
%204 = fneg double %203
br label %205
205: ; preds = %196, %188
%206 = phi double [ %195, %188 ], [ %204, %196 ]
%207 = load double, double* %6
%208 = load double, double* %11
%209 = fmul double %207, %208
%210 = fcmp oge double %206, %209
br i1 %210, label %211, label %217
211: ; preds = %205, %171
%212 = call double @sqrt(i32 -1)
%213 = load double*, double** %5
%214 = load i32, i32* %7
%215 = sext i32 %214 to i64
%216 = getelementptr inbounds double, double* %213, i64 %215
store double %212, double* %216
br label %225
217: ; preds = %205
%218 = load double, double* %10
%219 = load double*, double** %5
%220 = load i32, i32* %7
%221 = sext i32 %220 to i64
%222 = getelementptr inbounds double, double* %219, i64 %221
%223 = load double, double* %222
%224 = fsub double %223, %218
store double %224, double* %222
br label %225
225: ; preds = %217, %211
br label %226
226: ; preds = %225
%227 = load i32, i32* %7
%228 = add nsw i32 %227, 1
store i32 %228, i32* %7
br label %167
229: ; preds = %167
ret void
}

LLVM IR (-O2):
define void @zero_lc_averages(i32 %0, double* nocapture %1, double %2) {
%4 = icmp sgt i32 %0, 0
br i1 %4, label %5, label %107
5: ; preds = %3
%6 = zext i32 %0 to i64
br label %7
7: ; preds = %5, %24
%8 = phi i64 [ 0, %5 ], [ %28, %24 ]
%9 = phi x86_fp80 [ 0xK00000000000000000000, %5 ], [ %27, %24 ]
%10 = phi x86_fp80 [ 0xK00000000000000000000, %5 ], [ %26, %24 ]
%11 = phi i32 [ 0, %5 ], [ %25, %24 ]
%12 = getelementptr inbounds double, double* %1, i64 %8
%13 = load double, double* %12
%14 = tail call i64 @isnan(double %13)
%15 = icmp eq i64 %14, 0
br i1 %15, label %16, label %24
16: ; preds = %7
%17 = load double, double* %12
%18 = fmul double %17, %17
%19 = fpext double %18 to x86_fp80
%20 = fadd x86_fp80 %10, %19
%21 = fpext double %17 to x86_fp80
%22 = fadd x86_fp80 %9, %21
%23 = add nsw i32 %11, 1
br label %24
24: ; preds = %7, %16
%25 = phi i32 [ %11, %7 ], [ %23, %16 ]
%26 = phi x86_fp80 [ %10, %7 ], [ %20, %16 ]
%27 = phi x86_fp80 [ %9, %7 ], [ %22, %16 ]
%28 = add nuw nsw i64 %8, 1
%29 = icmp eq i64 %28, %6
br i1 %29, label %30, label %7
30: ; preds = %24
%31 = icmp sgt i32 %25, 0
br i1 %31, label %32, label %45
32: ; preds = %30
%33 = sitofp i32 %25 to x86_fp80
%34 = fdiv x86_fp80 %27, %33
%35 = fptrunc x86_fp80 %34 to double
%36 = fdiv x86_fp80 %26, %33
%37 = fmul x86_fp80 %27, %27
%38 = mul nsw i32 %25, %25
%39 = sitofp i32 %38 to x86_fp80
%40 = fdiv x86_fp80 %37, %39
%41 = fsub x86_fp80 %36, %40
%42 = fptrunc x86_fp80 %41 to double
%43 = fptosi double %42 to i32
%44 = tail call double @sqrt(i32 %43)
br label %45
45: ; preds = %30, %32
%46 = phi double [ %35, %32 ], [ 0.000000e+00, %30 ]
%47 = phi double [ %44, %32 ], [ 0.000000e+00, %30 ]
br i1 %4, label %48, label %107
48: ; preds = %45
%49 = fmul double %47, %2
%50 = zext i32 %0 to i64
br label %51
51: ; preds = %48, %70
%52 = phi i64 [ 0, %48 ], [ %73, %70 ]
%53 = phi x86_fp80 [ 0xK00000000000000000000, %48 ], [ %72, %70 ]
%54 = phi i32 [ 0, %48 ], [ %71, %70 ]
%55 = getelementptr inbounds double, double* %1, i64 %52
%56 = load double, double* %55
%57 = tail call i64 @isnan(double %56)
%58 = icmp eq i64 %57, 0
br i1 %58, label %59, label %70
59: ; preds = %51
%60 = load double, double* %55
%61 = fsub double %60, %46
%62 = fcmp ogt double %61, 0.000000e+00
%63 = fneg double %61
%64 = select i1 %62, double %61, double %63
%65 = fcmp olt double %64, %49
br i1 %65, label %66, label %70
66: ; preds = %59
%67 = fpext double %60 to x86_fp80
%68 = fadd x86_fp80 %53, %67
%69 = add nsw i32 %54, 1
br label %70
70: ; preds = %51, %66, %59
%71 = phi i32 [ %54, %51 ], [ %69, %66 ], [ %54, %59 ]
%72 = phi x86_fp80 [ %53, %51 ], [ %68, %66 ], [ %53, %59 ]
%73 = add nuw nsw i64 %52, 1
%74 = icmp eq i64 %73, %50
br i1 %74, label %75, label %51
75: ; preds = %70
%76 = icmp sgt i32 %71, 0
br i1 %76, label %77, label %81
77: ; preds = %75
%78 = sitofp i32 %71 to x86_fp80
%79 = fdiv x86_fp80 %72, %78
%80 = fptrunc x86_fp80 %79 to double
br label %81
81: ; preds = %75, %77
%82 = phi double [ %80, %77 ], [ 0.000000e+00, %75 ]
br i1 %4, label %83, label %107
83: ; preds = %81
%84 = fmul double %47, %2
%85 = zext i32 %0 to i64
br label %86
86: ; preds = %83, %103
%87 = phi i64 [ 0, %83 ], [ %105, %103 ]
%88 = getelementptr inbounds double, double* %1, i64 %87
%89 = load double, double* %88
%90 = tail call i64 @isnan(double %89)
%91 = icmp eq i64 %90, 0
br i1 %91, label %92, label %99
92: ; preds = %86
%93 = load double, double* %88
%94 = fsub double %93, %46
%95 = fcmp ogt double %94, 0.000000e+00
%96 = fneg double %94
%97 = select i1 %95, double %94, double %96
%98 = fcmp ult double %97, %84
br i1 %98, label %101, label %99
99: ; preds = %92, %86
%100 = tail call double @sqrt(i32 -1)
br label %103
101: ; preds = %92
%102 = fsub double %93, %82
br label %103
103: ; preds = %99, %101
%104 = phi double [ %102, %101 ], [ %100, %99 ]
store double %104, double* %88
%105 = add nuw nsw i64 %87, 1
%106 = icmp eq i64 %105, %85
br i1 %106, label %107, label %86
107: ; preds = %103, %45, %3, %81
ret void
}