#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);
double find_tauint_flat_spot(double *xs, int N) {
  double sumi = 0.0, sumi2 = 0.0, sumj = 0.0, sumj2 = 0.0, sumij = 0.0;
  double meani, meanj, stdi, stdj;
  double x, xi, xim, xj, xjp;
  int i, j, k, t, winsz, winszm1, denom;
  double autocorr_t, tauint;
  sumi = 0.0;
  sumi2 = 0.0;
  for (k = 0; k < N; k++) {
    x = xs[k];
    sumi += x;
    sumi2 += x * x;
  }
  sumj = sumi;
  sumj2 = sumi2;
  tauint = 1.0;
  for (t = 1; t < N - 1; t++) {
    winsz = N - t;
    winszm1 = winsz - 1;
    denom = winszm1;
    if (winszm1 == 0)
      denom = 1;
    i = t;
    j = N - t - 1;
    xi = xs[i];
    xj = xs[j];
    xim = xs[i - 1];
    xjp = xs[j + 1];
    sumi -= xim;
    sumi2 -= xim * xim;
    sumj -= xjp;
    sumj2 -= xjp * xjp;
    sumij = 0.0;
    for (k = 0; k < winsz; k++)
      sumij += xs[k] * xs[t + k];
    meani = sumi / winsz;
    meanj = sumj / winsz;
    stdi = sqrt((sumi2 - (sumi * sumi / winsz)) / denom);
    stdj = sqrt((sumj2 - (sumj * sumj / winsz)) / denom);
    autocorr_t = 0.0;
    if ((stdi != 0.0) && (stdj != 0.0))
      autocorr_t = (sumij / winsz - (meani * meanj)) / (stdi * stdj);
    tauint += 2.0 * autocorr_t;
    if (autocorr_t < 0.0)
      return tauint;
  }
  return tauint;
}


x86 assembly (clang, -O0):
<find_tauint_flat_spot>:
push   %rbp
mov    %rsp,%rbp
sub    $0xc0,%rsp
mov    %rdi,-0x10(%rbp)
mov    %esi,-0x14(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x20(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x28(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x30(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x38(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x40(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x20(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x28(%rbp)
movl   $0x0,-0x94(%rbp)
mov    -0x94(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    11c7 <find_tauint_flat_spot+0xb7>
mov    -0x10(%rbp),%rax
movslq -0x94(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  -0x68(%rbp),%xmm0
addsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x20(%rbp)
movsd  -0x68(%rbp),%xmm0
movsd  -0x68(%rbp),%xmm2
movsd  -0x28(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0x94(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x94(%rbp)
jmp    1164 <find_tauint_flat_spot+0x54>
movsd  -0x20(%rbp),%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x28(%rbp),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  0xe1d(%rip),%xmm0
movsd  %xmm0,-0xb8(%rbp)
ff
movl   $0x1,-0x98(%rbp)
mov    -0x98(%rbp),%eax
mov    -0x14(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    1557 <find_tauint_flat_spot+0x447>
mov    -0x14(%rbp),%eax
sub    -0x98(%rbp),%eax
mov    %eax,-0x9c(%rbp)
mov    -0x9c(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0xa0(%rbp)
mov    -0xa0(%rbp),%eax
mov    %eax,-0xa4(%rbp)
cmpl   $0x0,-0xa0(%rbp)
jne    124a <find_tauint_flat_spot+0x13a>
movl   $0x1,-0xa4(%rbp)
mov    -0x98(%rbp),%eax
mov    %eax,-0x8c(%rbp)
mov    -0x14(%rbp),%eax
sub    -0x98(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x90(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x8c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x70(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x90(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x80(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x8c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x78(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x90(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
movsd  -0x78(%rbp),%xmm1
movsd  -0x20(%rbp),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x20(%rbp)
movsd  -0x78(%rbp),%xmm0
movsd  -0x78(%rbp),%xmm2
movsd  -0x28(%rbp),%xmm1
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x28(%rbp)
movsd  -0x88(%rbp),%xmm1
ff
movsd  -0x30(%rbp),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x88(%rbp),%xmm0
ff
movsd  -0x88(%rbp),%xmm2
ff
movsd  -0x38(%rbp),%xmm1
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x40(%rbp)
movl   $0x0,-0x94(%rbp)
mov    -0x94(%rbp),%eax
cmp    -0x9c(%rbp),%eax
jge    13d0 <find_tauint_flat_spot+0x2c0>
mov    -0x10(%rbp),%rax
movslq -0x94(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x98(%rbp),%ecx
add    -0x94(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x40(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
mov    -0x94(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x94(%rbp)
jmp    1370 <find_tauint_flat_spot+0x260>
movsd  -0x20(%rbp),%xmm0
cvtsi2sdl -0x9c(%rbp),%xmm1
ff
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x30(%rbp),%xmm0
cvtsi2sdl -0x9c(%rbp),%xmm1
ff
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x28(%rbp),%xmm0
movsd  -0x20(%rbp),%xmm1
mulsd  -0x20(%rbp),%xmm1
cvtsi2sdl -0x9c(%rbp),%xmm2
ff
divsd  %xmm2,%xmm1
subsd  %xmm1,%xmm0
cvtsi2sdl -0xa4(%rbp),%xmm1
ff
divsd  %xmm1,%xmm0
call   1030 <sqrt@plt>
movsd  %xmm0,-0x58(%rbp)
movsd  -0x38(%rbp),%xmm0
movsd  -0x30(%rbp),%xmm1
mulsd  -0x30(%rbp),%xmm1
cvtsi2sdl -0x9c(%rbp),%xmm2
ff
divsd  %xmm2,%xmm1
subsd  %xmm1,%xmm0
cvtsi2sdl -0xa4(%rbp),%xmm1
ff
divsd  %xmm1,%xmm0
call   1030 <sqrt@plt>
movsd  %xmm0,-0x60(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0xb0(%rbp)
ff
movsd  -0x58(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    148e <find_tauint_flat_spot+0x37e>
jp     148e <find_tauint_flat_spot+0x37e>
jmp    14fb <find_tauint_flat_spot+0x3eb>
movsd  -0x60(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    14ab <find_tauint_flat_spot+0x39b>
jp     14ab <find_tauint_flat_spot+0x39b>
jmp    14fb <find_tauint_flat_spot+0x3eb>
movsd  -0x40(%rbp),%xmm1
cvtsi2sdl -0x9c(%rbp),%xmm0
ff
divsd  %xmm0,%xmm1
movsd  -0x48(%rbp),%xmm0
movsd  -0x50(%rbp),%xmm2
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  -0x58(%rbp),%xmm1
mulsd  -0x60(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0xb0(%rbp)
ff
movsd  -0xb0(%rbp),%xmm0
ff
movsd  -0xb8(%rbp),%xmm1
ff
addsd  %xmm0,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0xb8(%rbp)
ff
xorps  %xmm0,%xmm0
ucomisd -0xb0(%rbp),%xmm0
ff
jbe    153e <find_tauint_flat_spot+0x42e>
movsd  -0xb8(%rbp),%xmm0
ff
movsd  %xmm0,-0x8(%rbp)
jmp    1564 <find_tauint_flat_spot+0x454>
jmp    1543 <find_tauint_flat_spot+0x433>
mov    -0x98(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x98(%rbp)
jmp    11f5 <find_tauint_flat_spot+0xe5>
movsd  -0xb8(%rbp),%xmm0
ff
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
add    $0xc0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<find_tauint_flat_spot>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x78,%rsp
movsd  0xed9(%rip),%xmm8
test   %esi,%esi
jle    14cf <find_tauint_flat_spot+0x3bf>
mov    %rdi,%rbx
mov    %esi,%edi
lea    -0x1(%rdi),%rcx
mov    %edi,%eax
and    $0x3,%eax
cmp    $0x3,%rcx
jae    114f <find_tauint_flat_spot+0x3f>
xorpd  %xmm4,%xmm4
xor    %ecx,%ecx
xorpd  %xmm5,%xmm5
jmp    11b0 <find_tauint_flat_spot+0xa0>
mov    %edi,%edx
and    $0xfffffffc,%edx
xorpd  %xmm4,%xmm4
xor    %ecx,%ecx
xorpd  %xmm5,%xmm5
xchg   %ax,%ax
movsd  (%rbx,%rcx,8),%xmm0
movsd  0x8(%rbx,%rcx,8),%xmm1
addsd  %xmm0,%xmm4
mulsd  %xmm0,%xmm0
addsd  %xmm5,%xmm0
addsd  %xmm1,%xmm4
mulsd  %xmm1,%xmm1
addsd  %xmm0,%xmm1
movsd  0x10(%rbx,%rcx,8),%xmm0
addsd  %xmm0,%xmm4
mulsd  %xmm0,%xmm0
addsd  %xmm1,%xmm0
movsd  0x18(%rbx,%rcx,8),%xmm5
addsd  %xmm5,%xmm4
mulsd  %xmm5,%xmm5
addsd  %xmm0,%xmm5
add    $0x4,%rcx
cmp    %rcx,%rdx
jne    1160 <find_tauint_flat_spot+0x50>
test   %rax,%rax
je     11da <find_tauint_flat_spot+0xca>
lea    (%rbx,%rcx,8),%rcx
xor    %edx,%edx
nopl   0x0(%rax,%rax,1)
movsd  (%rcx,%rdx,8),%xmm0
addsd  %xmm0,%xmm4
mulsd  %xmm0,%xmm0
addsd  %xmm0,%xmm5
add    $0x1,%rdx
cmp    %rdx,%rax
jne    11c0 <find_tauint_flat_spot+0xb0>
cmp    $0x3,%esi
jl     14cf <find_tauint_flat_spot+0x3bf>
mov    %esi,%eax
add    $0xffffffff,%eax
mov    %eax,%ecx
mov    %rcx,0x58(%rsp)
unpcklpd %xmm5,%xmm5
lea    0x20(%rbx),%rbp
lea    0x8(%rbx),%r13
movsd  0xdfc(%rip),%xmm8
mov    $0x1,%r14d
xor    %r15d,%r15d
xorpd  %xmm10,%xmm10
mov    %eax,0x14(%rsp)
mov    %eax,%r12d
movapd %xmm4,%xmm9
mov    %rdi,0x18(%rsp)
jmp    1254 <find_tauint_flat_spot+0x144>
cs nopw 0x0(%rax,%rax,1)
nop
addsd  %xmm10,%xmm8
add    $0x1,%r14
add    $0x1,%r15d
add    $0x8,%rbp
add    $0xffffffff,%r12d
add    $0x8,%r13
cmp    0x58(%rsp),%r14
je     14cf <find_tauint_flat_spot+0x3bf>
mov    %rdi,%rcx
sub    %r14,%rcx
mov    %ecx,%r8d
add    $0xffffffff,%r8d
mov    $0x1,%eax
cmove  %eax,%r8d
xorpd  %xmm3,%xmm3
test   %rcx,%rcx
jle    1337 <find_tauint_flat_spot+0x227>
mov    0x14(%rsp),%eax
mov    %eax,%edx
sub    %r15d,%edx
lea    -0x1(%rdx),%rax
cmp    $0x3,%rax
jae    12a0 <find_tauint_flat_spot+0x190>
xorpd  %xmm3,%xmm3
xor    %esi,%esi
jmp    12f8 <find_tauint_flat_spot+0x1e8>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
mov    %r12d,%edi
and    $0xfffffffc,%edi
xorpd  %xmm3,%xmm3
xor    %esi,%esi
nopl   0x0(%rax)
movsd  (%rbx,%rsi,8),%xmm0
movsd  0x8(%rbx,%rsi,8),%xmm1
mulsd  -0x18(%rbp,%rsi,8),%xmm0
mulsd  -0x10(%rbp,%rsi,8),%xmm1
addsd  %xmm3,%xmm0
movsd  0x10(%rbx,%rsi,8),%xmm2
mulsd  -0x8(%rbp,%rsi,8),%xmm2
addsd  %xmm0,%xmm1
movsd  0x18(%rbx,%rsi,8),%xmm3
mulsd  0x0(%rbp,%rsi,8),%xmm3
addsd  %xmm1,%xmm2
addsd  %xmm2,%xmm3
add    $0x4,%rsi
cmp    %rsi,%rdi
jne    12b0 <find_tauint_flat_spot+0x1a0>
test   $0x3,%dl
je     1337 <find_tauint_flat_spot+0x227>
mov    %r12d,%edx
and    $0x3,%edx
lea    0x0(,%rsi,8),%rdi
add    %r13,%rdi
lea    (%rbx,%rsi,8),%rsi
xor    %eax,%eax
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movsd  (%rsi,%rax,8),%xmm0
mulsd  (%rdi,%rax,8),%xmm0
addsd  %xmm0,%xmm3
add    $0x1,%rax
cmp    %rax,%rdx
jne    1320 <find_tauint_flat_spot+0x210>
movsd  -0x8(%rbx,%r14,8),%xmm0
movsd  (%rbx,%rcx,8),%xmm2
subsd  %xmm0,%xmm4
unpcklpd %xmm2,%xmm0
mulpd  %xmm0,%xmm0
subpd  %xmm0,%xmm5
xorps  %xmm7,%xmm7
cvtsi2sd %ecx,%xmm7
movapd %xmm4,%xmm1
mulsd  %xmm4,%xmm1
divsd  %xmm7,%xmm1
movapd %xmm5,%xmm0
xorps  %xmm6,%xmm6
cvtsi2sd %r8d,%xmm6
subsd  %xmm1,%xmm0
divsd  %xmm6,%xmm0
ucomisd %xmm10,%xmm0
movsd  %xmm8,0x38(%rsp)
movsd  %xmm4,0x30(%rsp)
movapd %xmm5,0x60(%rsp)
movsd  %xmm3,0x28(%rsp)
movsd  %xmm7,0x20(%rsp)
jb     13b0 <find_tauint_flat_spot+0x2a0>
xorps  %xmm11,%xmm11
sqrtsd %xmm0,%xmm11
jmp    13f1 <find_tauint_flat_spot+0x2e1>
nopl   0x0(%rax,%rax,1)
movsd  %xmm9,0x8(%rsp)
movapd %xmm2,0x40(%rsp)
movsd  %xmm6,0x50(%rsp)
call   1030 <sqrt@plt>
movsd  0x50(%rsp),%xmm6
movapd 0x40(%rsp),%xmm2
movsd  0x20(%rsp),%xmm7
movsd  0x8(%rsp),%xmm9
xorpd  %xmm10,%xmm10
movapd 0x60(%rsp),%xmm5
movapd %xmm0,%xmm11
subsd  %xmm2,%xmm9
movapd %xmm9,%xmm1
mulsd  %xmm9,%xmm1
divsd  %xmm7,%xmm1
movapd %xmm5,%xmm0
unpckhpd %xmm5,%xmm0
subsd  %xmm1,%xmm0
divsd  %xmm6,%xmm0
ucomisd %xmm10,%xmm0
jb     1430 <find_tauint_flat_spot+0x320>
sqrtsd %xmm0,%xmm0
jmp    1462 <find_tauint_flat_spot+0x352>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movsd  %xmm9,0x8(%rsp)
movsd  %xmm11,0x40(%rsp)
call   1030 <sqrt@plt>
movsd  0x40(%rsp),%xmm11
movsd  0x20(%rsp),%xmm7
movsd  0x8(%rsp),%xmm9
xorpd  %xmm10,%xmm10
movapd 0x60(%rsp),%xmm5
mov    0x18(%rsp),%rdi
movsd  0x38(%rsp),%xmm8
movsd  0x30(%rsp),%xmm4
movsd  0x28(%rsp),%xmm3
ucomisd %xmm10,%xmm11
jne    1487 <find_tauint_flat_spot+0x377>
jnp    1230 <find_tauint_flat_spot+0x120>
ucomisd %xmm10,%xmm0
jne    1494 <find_tauint_flat_spot+0x384>
jnp    1230 <find_tauint_flat_spot+0x120>
movapd %xmm9,%xmm1
divsd  %xmm7,%xmm1
divsd  %xmm7,%xmm3
movapd %xmm4,%xmm2
divsd  %xmm7,%xmm2
mulsd  %xmm1,%xmm2
subsd  %xmm2,%xmm3
mulsd  %xmm0,%xmm11
divsd  %xmm11,%xmm3
ucomisd %xmm3,%xmm10
addsd  %xmm3,%xmm3
addsd  %xmm3,%xmm8
jbe    1235 <find_tauint_flat_spot+0x125>
movapd %xmm8,%xmm0
add    $0x78,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
define double @find_tauint_flat_spot(double* nocapture %0, i32 %1) {
%3 = icmp sgt i32 %1, 0
br i1 %3, label %4, label %188
4: ; preds = %2
%5 = zext i32 %1 to i64
%6 = add nsw i64 %5, -1
%7 = and i64 %5, 3
%8 = icmp ult i64 %6, 3
br i1 %8, label %11, label %9
9: ; preds = %4
%10 = and i64 %5, 4294967292
br label %40
11: ; preds = %40, %4
%12 = phi double [ undef, %4 ], [ %62, %40 ]
%13 = phi double [ undef, %4 ], [ %63, %40 ]
%14 = phi i64 [ 0, %4 ], [ %64, %40 ]
%15 = phi double [ 0.000000e+00, %4 ], [ %62, %40 ]
%16 = phi double [ 0.000000e+00, %4 ], [ %63, %40 ]
%17 = icmp eq i64 %7, 0
br i1 %17, label %30, label %18
18: ; preds = %11, %18
%19 = phi i64 [ %27, %18 ], [ %14, %11 ]
%20 = phi double [ %25, %18 ], [ %15, %11 ]
%21 = phi double [ %26, %18 ], [ %16, %11 ]
%22 = phi i64 [ %28, %18 ], [ 0, %11 ]
%23 = getelementptr inbounds double, double* %0, i64 %19
%24 = load double, double* %23
%25 = fadd double %20, %24
%26 = tail call double @llvm.fmuladd.f64(double %24, double %24, double %21)
%27 = add nuw nsw i64 %19, 1
%28 = add i64 %22, 1
%29 = icmp eq i64 %28, %7
br i1 %29, label %30, label %18
30: ; preds = %18, %11
%31 = phi double [ %12, %11 ], [ %25, %18 ]
%32 = phi double [ %13, %11 ], [ %26, %18 ]
%33 = add i32 %1, -1
%34 = icmp sgt i32 %1, 2
br i1 %34, label %35, label %188
35: ; preds = %30
%36 = zext i32 %1 to i64
%37 = zext i32 %33 to i64
%38 = insertelement <2 x double> poison, double %32, i64 0
%39 = shufflevector <2 x double> %38, <2 x double> poison, <2 x i32> zeroinitializer
br label %67
40: ; preds = %40, %9
%41 = phi i64 [ 0, %9 ], [ %64, %40 ]
%42 = phi double [ 0.000000e+00, %9 ], [ %62, %40 ]
%43 = phi double [ 0.000000e+00, %9 ], [ %63, %40 ]
%44 = phi i64 [ 0, %9 ], [ %65, %40 ]
%45 = getelementptr inbounds double, double* %0, i64 %41
%46 = load double, double* %45
%47 = fadd double %42, %46
%48 = tail call double @llvm.fmuladd.f64(double %46, double %46, double %43)
%49 = or i64 %41, 1
%50 = getelementptr inbounds double, double* %0, i64 %49
%51 = load double, double* %50
%52 = fadd double %47, %51
%53 = tail call double @llvm.fmuladd.f64(double %51, double %51, double %48)
%54 = or i64 %41, 2
%55 = getelementptr inbounds double, double* %0, i64 %54
%56 = load double, double* %55
%57 = fadd double %52, %56
%58 = tail call double @llvm.fmuladd.f64(double %56, double %56, double %53)
%59 = or i64 %41, 3
%60 = getelementptr inbounds double, double* %0, i64 %59
%61 = load double, double* %60
%62 = fadd double %57, %61
%63 = tail call double @llvm.fmuladd.f64(double %61, double %61, double %58)
%64 = add nuw nsw i64 %41, 4
%65 = add i64 %44, 4
%66 = icmp eq i64 %65, %10
br i1 %66, label %11, label %40
67: ; preds = %35, %183
%68 = phi i32 [ 0, %35 ], [ %187, %183 ]
%69 = phi i64 [ 1, %35 ], [ %185, %183 ]
%70 = phi double [ 1.000000e+00, %35 ], [ %184, %183 ]
%71 = phi double [ %31, %35 ], [ %86, %183 ]
%72 = phi double [ %31, %35 ], [ %87, %183 ]
%73 = phi <2 x double> [ %39, %35 ], [ %91, %183 ]
%74 = sub i32 %33, %68
%75 = zext i32 %74 to i64
%76 = sub nsw i64 %36, %69
%77 = trunc i64 %76 to i32
%78 = add i32 %77, -1
%79 = icmp eq i32 %78, 0
%80 = select i1 %79, i32 1, i32 %78
%81 = add nsw i64 %69, -1
%82 = getelementptr inbounds double, double* %0, i64 %81
%83 = load double, double* %82
%84 = getelementptr inbounds double, double* %0, i64 %76
%85 = load double, double* %84
%86 = fsub double %71, %83
%87 = fsub double %72, %85
%88 = insertelement <2 x double> poison, double %83, i64 0
%89 = insertelement <2 x double> %88, double %85, i64 1
%90 = fneg <2 x double> %89
%91 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %90, <2 x double> %89, <2 x double> %73)
%92 = icmp sgt i64 %76, 0
br i1 %92, label %93, label %151
93: ; preds = %67
%94 = add nsw i64 %75, -1
%95 = and i64 %75, 3
%96 = icmp ult i64 %94, 3
br i1 %96, label %133, label %97
97: ; preds = %93
%98 = and i64 %75, 4294967292
br label %99
99: ; preds = %99, %97
%100 = phi i64 [ 0, %97 ], [ %130, %99 ]
%101 = phi double [ 0.000000e+00, %97 ], [ %129, %99 ]
%102 = phi i64 [ 0, %97 ], [ %131, %99 ]
%103 = getelementptr inbounds double, double* %0, i64 %100
%104 = load double, double* %103
%105 = add nuw nsw i64 %100, %69
%106 = getelementptr inbounds double, double* %0, i64 %105
%107 = load double, double* %106
%108 = tail call double @llvm.fmuladd.f64(double %104, double %107, double %101)
%109 = or i64 %100, 1
%110 = getelementptr inbounds double, double* %0, i64 %109
%111 = load double, double* %110
%112 = add nuw nsw i64 %109, %69
%113 = getelementptr inbounds double, double* %0, i64 %112
%114 = load double, double* %113
%115 = tail call double @llvm.fmuladd.f64(double %111, double %114, double %108)
%116 = or i64 %100, 2
%117 = getelementptr inbounds double, double* %0, i64 %116
%118 = load double, double* %117
%119 = add nuw nsw i64 %116, %69
%120 = getelementptr inbounds double, double* %0, i64 %119
%121 = load double, double* %120
%122 = tail call double @llvm.fmuladd.f64(double %118, double %121, double %115)
%123 = or i64 %100, 3
%124 = getelementptr inbounds double, double* %0, i64 %123
%125 = load double, double* %124
%126 = add nuw nsw i64 %123, %69
%127 = getelementptr inbounds double, double* %0, i64 %126
%128 = load double, double* %127
%129 = tail call double @llvm.fmuladd.f64(double %125, double %128, double %122)
%130 = add nuw nsw i64 %100, 4
%131 = add i64 %102, 4
%132 = icmp eq i64 %131, %98
br i1 %132, label %133, label %99
133: ; preds = %99, %93
%134 = phi double [ undef, %93 ], [ %129, %99 ]
%135 = phi i64 [ 0, %93 ], [ %130, %99 ]
%136 = phi double [ 0.000000e+00, %93 ], [ %129, %99 ]
%137 = icmp eq i64 %95, 0
br i1 %137, label %151, label %138
138: ; preds = %133, %138
%139 = phi i64 [ %148, %138 ], [ %135, %133 ]
%140 = phi double [ %147, %138 ], [ %136, %133 ]
%141 = phi i64 [ %149, %138 ], [ 0, %133 ]
%142 = getelementptr inbounds double, double* %0, i64 %139
%143 = load double, double* %142
%144 = add nuw nsw i64 %139, %69
%145 = getelementptr inbounds double, double* %0, i64 %144
%146 = load double, double* %145
%147 = tail call double @llvm.fmuladd.f64(double %143, double %146, double %140)
%148 = add nuw nsw i64 %139, 1
%149 = add i64 %141, 1
%150 = icmp eq i64 %149, %95
br i1 %150, label %151, label %138
151: ; preds = %133, %138, %67
%152 = phi double [ 0.000000e+00, %67 ], [ %134, %133 ], [ %147, %138 ]
%153 = trunc i64 %76 to i32
%154 = sitofp i32 %153 to double
%155 = fmul double %86, %86
%156 = fdiv double %155, %154
%157 = extractelement <2 x double> %91, i64 0
%158 = fsub double %157, %156
%159 = sitofp i32 %80 to double
%160 = fdiv double %158, %159
%161 = tail call double @sqrt(double %160)
%162 = fmul double %87, %87
%163 = fdiv double %162, %154
%164 = extractelement <2 x double> %91, i64 1
%165 = fsub double %164, %163
%166 = fdiv double %165, %159
%167 = tail call double @sqrt(double %166)
%168 = fcmp une double %161, 0.000000e+00
%169 = fcmp une double %167, 0.000000e+00
%170 = select i1 %168, i1 %169, i1 false
br i1 %170, label %173, label %171
171: ; preds = %151
%172 = fadd double %70, 0.000000e+00
br label %183
173: ; preds = %151
%174 = fdiv double %87, %154
%175 = fdiv double %152, %154
%176 = fneg double %86
%177 = fdiv double %176, %154
%178 = tail call double @llvm.fmuladd.f64(double %177, double %174, double %175)
%179 = fmul double %161, %167
%180 = fdiv double %178, %179
%181 = tail call double @llvm.fmuladd.f64(double %180, double 2.000000e+00, double %70)
%182 = fcmp olt double %180, 0.000000e+00
br i1 %182, label %188, label %183
183: ; preds = %171, %173
%184 = phi double [ %172, %171 ], [ %181, %173 ]
%185 = add nuw nsw i64 %69, 1
%186 = icmp eq i64 %185, %37
%187 = add i32 %68, 1
br i1 %186, label %188, label %67
188: ; preds = %173, %183, %2, %30
%189 = phi double [ 1.000000e+00, %30 ], [ 1.000000e+00, %2 ], [ %184, %183 ], [ %181, %173 ]
ret double %189
}