#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 */
int /*<<< orphan*/ exit(int);
scalar_t__ fabs(double);
int /*<<< orphan*/ fprintf(int /*<<< orphan*/, char *, int);
int /*<<< orphan*/ free(double *);
scalar_t__ malloc(int);
int /*<<< orphan*/ memcpy(double *, double *, int);
double nthd(int, double *, double);
int /*<<< orphan*/ stderr;
double getTukey_d(int ndat, double *dat, double *med, double *mad) {
  int niter = 10;
  int n, i;
  double t, c, s;
  double u;
  double sdenom, snumer;
  double *a;
  if ((a = (double *)malloc(sizeof(double) * ndat)) == NULL) {
    fprintf(stderr, "Cannot allocate a in getTukey_d, ndat=%d\n", ndat);
    exit(-1);
  }
  memcpy(a, dat, sizeof(double) * ndat);
  *med = nthd(ndat, a, (double)(ndat - 1) * 0.5);
  for (i = 0; i < ndat; i++) {
    a[i] = (double)fabs((double)(a[i] - (*med)));
  }
  *mad = nthd(ndat, a, ((double)(ndat - 1) * 0.5));
  free(a);
  s = (*mad) / 0.6745;
  t = (*med);
  c = 6.0;
  if (*mad < 0.01)
    return t;
  for (n = 0; n < niter; n++) {
    snumer = 0.;
    sdenom = 0.;
    for (i = 0; i < ndat; i++) {
      u = (dat[i] - t) / (c * s);
      if (u < 1.0 && u > -1.0) {
        snumer += dat[i] * (1.0 - u * u) * (1.0 - u * u);
        sdenom += (1.0 - u * u) * (1.0 - u * u);
      }
    }
    if (sdenom != 0) {
      t = (double)(snumer / sdenom);
    } else {
      c *= 1.1;
    }
  }
  return t;
}


x86 assembly (clang, -O0):
<getTukey_d>:
push   %rbp
mov    %rsp,%rbp
sub    $0x70,%rsp
mov    %edi,-0xc(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %rdx,-0x20(%rbp)
mov    %rcx,-0x28(%rbp)
movl   $0xa,-0x2c(%rbp)
movslq -0xc(%rbp),%rax
shl    $0x3,%rax
mov    %eax,%edi
call   1070 <malloc@plt>
mov    %rax,-0x70(%rbp)
cmp    $0x0,%rax
jne    11cd <getTukey_d+0x5d>
mov    0x2e2e(%rip),%rax
mov    (%rax),%edi
mov    -0xc(%rbp),%edx
lea    0xe7a(%rip),%rsi
call   1050 <fprintf@plt>
mov    $0xffffffff,%edi
call   1080 <exit@plt>
mov    -0x70(%rbp),%rdi
mov    -0x18(%rbp),%rsi
movslq -0xc(%rbp),%rax
shl    $0x3,%rax
mov    %eax,%edx
call   1060 <memcpy@plt>
mov    -0xc(%rbp),%edi
mov    -0x70(%rbp),%rsi
mov    -0xc(%rbp),%eax
sub    $0x1,%eax
cvtsi2sd %eax,%xmm0
movsd  0xe03(%rip),%xmm1
mulsd  %xmm1,%xmm0
call   1090 <nthd@plt>
mov    -0x20(%rbp),%rax
movsd  %xmm0,(%rax)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    125b <getTukey_d+0xeb>
mov    -0x70(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
subsd  (%rax),%xmm0
call   1040 <fabs@plt>
cvtsi2sd %rax,%xmm0
mov    -0x70(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1215 <getTukey_d+0xa5>
mov    -0xc(%rbp),%edi
mov    -0x70(%rbp),%rsi
mov    -0xc(%rbp),%eax
sub    $0x1,%eax
cvtsi2sd %eax,%xmm0
movsd  0xd8c(%rip),%xmm1
mulsd  %xmm1,%xmm0
call   1090 <nthd@plt>
mov    -0x28(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x70(%rbp),%rdi
call   1030 <free@plt>
mov    -0x28(%rbp),%rax
movsd  (%rax),%xmm0
movsd  0xd7a(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
mov    -0x20(%rbp),%rax
movsd  (%rax),%xmm0
movsd  %xmm0,-0x40(%rbp)
movsd  0xd54(%rip),%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x28(%rbp),%rax
movsd  0xd3b(%rip),%xmm0
ucomisd (%rax),%xmm0
jbe    12e6 <getTukey_d+0x176>
movsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    14d7 <getTukey_d+0x367>
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0x2c(%rbp),%eax
jge    14cd <getTukey_d+0x35d>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x68(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x60(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1477 <getTukey_d+0x307>
mov    -0x18(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  -0x40(%rbp),%xmm0
movsd  -0x48(%rbp),%xmm1
mulsd  -0x50(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
movsd  0xcdf(%rip),%xmm0
ucomisd -0x58(%rbp),%xmm0
jbe    1464 <getTukey_d+0x2f4>
movsd  -0x58(%rbp),%xmm0
movsd  0xccf(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    1464 <getTukey_d+0x2f4>
mov    -0x18(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  -0x58(%rbp),%xmm1
movsd  -0x58(%rbp),%xmm2
movq   %xmm1,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm1
mulsd  %xmm2,%xmm1
movsd  0xc83(%rip),%xmm2
addsd  %xmm2,%xmm1
mulsd  %xmm1,%xmm0
movsd  -0x58(%rbp),%xmm2
movsd  -0x58(%rbp),%xmm1
movq   %xmm2,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm2
mulsd  %xmm1,%xmm2
movsd  0xc4e(%rip),%xmm1
addsd  %xmm1,%xmm2
movsd  -0x68(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  -0x58(%rbp),%xmm0
movsd  -0x58(%rbp),%xmm1
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm1,%xmm0
movsd  0xc0b(%rip),%xmm1
addsd  %xmm1,%xmm0
movsd  -0x58(%rbp),%xmm2
movsd  -0x58(%rbp),%xmm1
movq   %xmm2,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm2
mulsd  %xmm1,%xmm2
movsd  0xbda(%rip),%xmm1
addsd  %xmm1,%xmm2
movsd  -0x60(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x60(%rbp)
jmp    1469 <getTukey_d+0x2f9>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1310 <getTukey_d+0x1a0>
movsd  -0x60(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1494 <getTukey_d+0x324>
jp     1494 <getTukey_d+0x324>
jmp    14a8 <getTukey_d+0x338>
movsd  -0x68(%rbp),%xmm0
divsd  -0x60(%rbp),%xmm0
movsd  %xmm0,-0x40(%rbp)
jmp    14ba <getTukey_d+0x34a>
movsd  0xb70(%rip),%xmm0
mulsd  -0x48(%rbp),%xmm0
movsd  %xmm0,-0x48(%rbp)
jmp    14bf <getTukey_d+0x34f>
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    12ed <getTukey_d+0x17d>
movsd  -0x40(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
add    $0x70,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<getTukey_d>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x18,%rsp
mov    %rcx,%r12
mov    %rdx,%r15
mov    %rsi,%rbx
mov    %edi,%r14d
lea    0x0(,%r14,8),%ebp
mov    %ebp,%edi
call   1070 <malloc@plt>
test   %rax,%rax
je     13de <getTukey_d+0x26e>
mov    %rax,%r13
mov    %rax,%rdi
mov    %rbx,%rsi
mov    %ebp,%edx
call   1060 <memcpy@plt>
lea    -0x1(%r14),%eax
cvtsi2sd %eax,%xmm0
mulsd  0xe3e(%rip),%xmm0
mov    %r14d,%edi
mov    %r13,%rsi
movsd  %xmm0,0x10(%rsp)
call   1090 <nthd@plt>
movsd  %xmm0,(%r15)
test   %r14d,%r14d
jle    12bd <getTukey_d+0x14d>
movsd  0x0(%r13),%xmm1
subsd  %xmm0,%xmm1
movapd %xmm1,%xmm0
call   1040 <fabs@plt>
xorps  %xmm0,%xmm0
cvtsi2sd %rax,%xmm0
movsd  %xmm0,0x0(%r13)
cmp    $0x1,%r14d
je     12bd <getTukey_d+0x14d>
mov    %r12,0x8(%rsp)
mov    %r14d,%eax
add    $0xffffffffffffffff,%rax
mov    %rax,(%rsp)
mov    $0x1,%ebp
cmp    $0x2,%r14d
je     1292 <getTukey_d+0x122>
mov    (%rsp),%r12
and    $0xfffffffffffffffe,%r12
neg    %r12
mov    $0x1,%ebp
nopw   0x0(%rax,%rax,1)
movsd  0x0(%r13,%rbp,8),%xmm0
subsd  (%r15),%xmm0
call   1040 <fabs@plt>
xorps  %xmm0,%xmm0
cvtsi2sd %rax,%xmm0
movsd  %xmm0,0x0(%r13,%rbp,8)
movsd  0x8(%r13,%rbp,8),%xmm0
subsd  (%r15),%xmm0
call   1040 <fabs@plt>
xorps  %xmm0,%xmm0
cvtsi2sd %rax,%xmm0
movsd  %xmm0,0x8(%r13,%rbp,8)
lea    (%r12,%rbp,1),%rax
add    $0x2,%rax
add    $0x2,%rbp
cmp    $0x1,%rax
jne    1240 <getTukey_d+0xd0>
testb  $0x1,(%rsp)
mov    0x8(%rsp),%r12
je     12bd <getTukey_d+0x14d>
movsd  0x0(%r13,%rbp,8),%xmm0
subsd  (%r15),%xmm0
call   1040 <fabs@plt>
xorps  %xmm0,%xmm0
cvtsi2sd %rax,%xmm0
movsd  %xmm0,0x0(%r13,%rbp,8)
mov    %r14d,%edi
mov    %r13,%rsi
movsd  0x10(%rsp),%xmm0
call   1090 <nthd@plt>
movsd  %xmm0,(%r12)
mov    %r13,%rdi
call   1030 <free@plt>
movsd  (%r12),%xmm10
movsd  (%r15),%xmm0
movsd  0xd21(%rip),%xmm1
ucomisd %xmm10,%xmm1
jbe    1305 <getTukey_d+0x195>
add    $0x18,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
divsd  0xcfa(%rip),%xmm10
mov    %r14d,%eax
movsd  0xcfe(%rip),%xmm11
xor    %ecx,%ecx
movsd  0xcfc(%rip),%xmm2
movsd  0xcfc(%rip),%xmm4
xorpd  %xmm9,%xmm9
movsd  0xcf6(%rip),%xmm8
jmp    134d <getTukey_d+0x1dd>
nopl   0x0(%rax)
mulsd  %xmm8,%xmm11
add    $0x1,%ecx
cmp    $0xa,%ecx
je     12f6 <getTukey_d+0x186>
test   %r14d,%r14d
jle    1340 <getTukey_d+0x1d0>
movapd %xmm10,%xmm1
mulsd  %xmm11,%xmm1
xorpd  %xmm7,%xmm7
xor    %edx,%edx
jmp    1379 <getTukey_d+0x209>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
add    $0x1,%rdx
cmp    %rdx,%rax
je     13c0 <getTukey_d+0x250>
movsd  (%rbx,%rdx,8),%xmm6
movapd %xmm6,%xmm5
subsd  %xmm0,%xmm5
divsd  %xmm1,%xmm5
ucomisd %xmm5,%xmm2
jbe    1370 <getTukey_d+0x200>
ucomisd %xmm4,%xmm5
jbe    1370 <getTukey_d+0x200>
mulsd  %xmm5,%xmm5
movapd %xmm2,%xmm3
subsd  %xmm5,%xmm3
mulsd  %xmm3,%xmm6
movapd %xmm3,%xmm5
unpcklpd %xmm6,%xmm5
unpcklpd %xmm3,%xmm3
mulpd  %xmm5,%xmm3
addpd  %xmm3,%xmm7
jmp    1370 <getTukey_d+0x200>
nopl   0x0(%rax)
ucomisd %xmm9,%xmm7
jne    13cd <getTukey_d+0x25d>
jnp    1340 <getTukey_d+0x1d0>
movapd %xmm7,%xmm0
unpckhpd %xmm7,%xmm0
divsd  %xmm7,%xmm0
jmp    1345 <getTukey_d+0x1d5>
mov    0x2bfb(%rip),%rax
mov    (%rax),%edi
lea    0xc4a(%rip),%rsi
mov    %r14d,%edx
call   1050 <fprintf@plt>
mov    $0xffffffff,%edi
call   1080 <exit@plt>

LLVM IR (-O0):
@stderr = global i32 0
@.str = [42 x i8] c"Cannot allocate a in getTukey_d, ndat=%d\0A\00"
define double @getTukey_d(i32 %0, double* %1, double* %2, double* %3) {
%5 = alloca double
%6 = alloca i32
%7 = alloca double*
%8 = alloca double*
%9 = alloca double*
%10 = alloca i32
%11 = alloca i32
%12 = alloca i32
%13 = alloca double
%14 = alloca double
%15 = alloca double
%16 = alloca double
%17 = alloca double
%18 = alloca double
%19 = alloca double*
store i32 %0, i32* %6
store double* %1, double** %7
store double* %2, double** %8
store double* %3, double** %9
store i32 10, i32* %10
%20 = load i32, i32* %6
%21 = sext i32 %20 to i64
%22 = mul i64 8, %21
%23 = trunc i64 %22 to i32
%24 = call i64 @malloc(i32 %23)
%25 = inttoptr i64 %24 to double*
store double* %25, double** %19
%26 = icmp eq double* %25, null
br i1 %26, label %27, label %32
27: ; preds = %4
%28 = load i32, i32* @stderr
%29 = load i32, i32* %6
%30 = call i32 @fprintf(i32 %28, i8* getelementptr inbounds ([42 x i8], [42 x i8]* @.str, i64 0, i64 0), i32 %29)
%31 = call i32 @exit(i32 -1)
unreachable
32: ; preds = %4
%33 = load double*, double** %19
%34 = load double*, double** %7
%35 = load i32, i32* %6
%36 = sext i32 %35 to i64
%37 = mul i64 8, %36
%38 = trunc i64 %37 to i32
%39 = call i32 @memcpy(double* %33, double* %34, i32 %38)
%40 = load i32, i32* %6
%41 = load double*, double** %19
%42 = load i32, i32* %6
%43 = sub nsw i32 %42, 1
%44 = sitofp i32 %43 to double
%45 = fmul double %44, 5.000000e-01
%46 = call double @nthd(i32 %40, double* %41, double %45)
%47 = load double*, double** %8
store double %46, double* %47
store i32 0, i32* %12
br label %48
48: ; preds = %67, %32
%49 = load i32, i32* %12
%50 = load i32, i32* %6
%51 = icmp slt i32 %49, %50
br i1 %51, label %52, label %70
52: ; preds = %48
%53 = load double*, double** %19
%54 = load i32, i32* %12
%55 = sext i32 %54 to i64
%56 = getelementptr inbounds double, double* %53, i64 %55
%57 = load double, double* %56
%58 = load double*, double** %8
%59 = load double, double* %58
%60 = fsub double %57, %59
%61 = call i64 @fabs(double %60)
%62 = sitofp i64 %61 to double
%63 = load double*, double** %19
%64 = load i32, i32* %12
%65 = sext i32 %64 to i64
%66 = getelementptr inbounds double, double* %63, i64 %65
store double %62, double* %66
br label %67
67: ; preds = %52
%68 = load i32, i32* %12
%69 = add nsw i32 %68, 1
store i32 %69, i32* %12
br label %48
70: ; preds = %48
%71 = load i32, i32* %6
%72 = load double*, double** %19
%73 = load i32, i32* %6
%74 = sub nsw i32 %73, 1
%75 = sitofp i32 %74 to double
%76 = fmul double %75, 5.000000e-01
%77 = call double @nthd(i32 %71, double* %72, double %76)
%78 = load double*, double** %9
store double %77, double* %78
%79 = load double*, double** %19
%80 = call i32 @free(double* %79)
%81 = load double*, double** %9
%82 = load double, double* %81
%83 = fdiv double %82, 6.745000e-01
store double %83, double* %15
%84 = load double*, double** %8
%85 = load double, double* %84
store double %85, double* %13
store double 6.000000e+00, double* %14
%86 = load double*, double** %9
%87 = load double, double* %86
%88 = fcmp olt double %87, 1.000000e-02
br i1 %88, label %89, label %91
89: ; preds = %70
%90 = load double, double* %13
store double %90, double* %5
br label %165
91: ; preds = %70
store i32 0, i32* %11
br label %92
92: ; preds = %160, %91
%93 = load i32, i32* %11
%94 = load i32, i32* %10
%95 = icmp slt i32 %93, %94
br i1 %95, label %96, label %163
96: ; preds = %92
store double 0.000000e+00, double* %18
store double 0.000000e+00, double* %17
store i32 0, i32* %12
br label %97
97: ; preds = %146, %96
%98 = load i32, i32* %12
%99 = load i32, i32* %6
%100 = icmp slt i32 %98, %99
br i1 %100, label %101, label %149
101: ; preds = %97
%102 = load double*, double** %7
%103 = load i32, i32* %12
%104 = sext i32 %103 to i64
%105 = getelementptr inbounds double, double* %102, i64 %104
%106 = load double, double* %105
%107 = load double, double* %13
%108 = fsub double %106, %107
%109 = load double, double* %14
%110 = load double, double* %15
%111 = fmul double %109, %110
%112 = fdiv double %108, %111
store double %112, double* %16
%113 = load double, double* %16
%114 = fcmp olt double %113, 1.000000e+00
br i1 %114, label %115, label %145
115: ; preds = %101
%116 = load double, double* %16
%117 = fcmp ogt double %116, -1.000000e+00
br i1 %117, label %118, label %145
118: ; preds = %115
%119 = load double*, double** %7
%120 = load i32, i32* %12
%121 = sext i32 %120 to i64
%122 = getelementptr inbounds double, double* %119, i64 %121
%123 = load double, double* %122
%124 = load double, double* %16
%125 = load double, double* %16
%126 = fneg double %124
%127 = call double @llvm.fmuladd.f64(double %126, double %125, double 1.000000e+00)
%128 = fmul double %123, %127
%129 = load double, double* %16
%130 = load double, double* %16
%131 = fneg double %129
%132 = call double @llvm.fmuladd.f64(double %131, double %130, double 1.000000e+00)
%133 = load double, double* %18
%134 = call double @llvm.fmuladd.f64(double %128, double %132, double %133)
store double %134, double* %18
%135 = load double, double* %16
%136 = load double, double* %16
%137 = fneg double %135
%138 = call double @llvm.fmuladd.f64(double %137, double %136, double 1.000000e+00)
%139 = load double, double* %16
%140 = load double, double* %16
%141 = fneg double %139
%142 = call double @llvm.fmuladd.f64(double %141, double %140, double 1.000000e+00)
%143 = load double, double* %17
%144 = call double @llvm.fmuladd.f64(double %138, double %142, double %143)
store double %144, double* %17
br label %145
145: ; preds = %118, %115, %101
br label %146
146: ; preds = %145
%147 = load i32, i32* %12
%148 = add nsw i32 %147, 1
store i32 %148, i32* %12
br label %97
149: ; preds = %97
%150 = load double, double* %17
%151 = fcmp une double %150, 0.000000e+00
br i1 %151, label %152, label %156
152: ; preds = %149
%153 = load double, double* %18
%154 = load double, double* %17
%155 = fdiv double %153, %154
store double %155, double* %13
br label %159
156: ; preds = %149
%157 = load double, double* %14
%158 = fmul double %157, 1.100000e+00
store double %158, double* %14
br label %159
159: ; preds = %156, %152
br label %160
160: ; preds = %159
%161 = load i32, i32* %11
%162 = add nsw i32 %161, 1
store i32 %162, i32* %11
br label %92
163: ; preds = %92
%164 = load double, double* %13
store double %164, double* %5
br label %165
165: ; preds = %163, %89
%166 = load double, double* %5
ret double %166
}

LLVM IR (-O2):
@stderr = global i32 0
@.str = [42 x i8] c"Cannot allocate a in getTukey_d, ndat=%d\0A\00"
define double @getTukey_d(i32 %0, double* %1, double* nocapture %2, double* nocapture %3) {
%5 = shl i32 %0, 3
%6 = tail call i64 @malloc(i32 %5)
%7 = inttoptr i64 %6 to double*
%8 = icmp eq i64 %6, 0
br i1 %8, label %9, label %13
9: ; preds = %4
%10 = load i32, i32* @stderr
%11 = tail call i32 @fprintf(i32 %10, i8* getelementptr inbounds ([42 x i8], [42 x i8]* @.str, i64 0, i64 0), i32 %0)
%12 = tail call i32 @exit(i32 -1)
unreachable
13: ; preds = %4
%14 = tail call i32 @memcpy(double* %7, double* %1, i32 %5)
%15 = add nsw i32 %0, -1
%16 = sitofp i32 %15 to double
%17 = fmul double %16, 5.000000e-01
%18 = tail call double @nthd(i32 %0, double* %7, double %17)
store double %18, double* %2
%19 = icmp sgt i32 %0, 0
br i1 %19, label %20, label %62
20: ; preds = %13
%21 = load double, double* %7
%22 = fsub double %21, %18
%23 = tail call i64 @fabs(double %22)
%24 = sitofp i64 %23 to double
store double %24, double* %7
%25 = icmp eq i32 %0, 1
br i1 %25, label %62, label %26
26: ; preds = %20
%27 = zext i32 %0 to i64
%28 = add nsw i64 %27, -1
%29 = and i64 %28, 1
%30 = icmp eq i32 %0, 2
br i1 %30, label %52, label %31
31: ; preds = %26
%32 = and i64 %28, -2
br label %33
33: ; preds = %33, %31
%34 = phi i64 [ 1, %31 ], [ %49, %33 ]
%35 = phi i64 [ 0, %31 ], [ %50, %33 ]
%36 = load double, double* %2
%37 = getelementptr inbounds double, double* %7, i64 %34
%38 = load double, double* %37
%39 = fsub double %38, %36
%40 = tail call i64 @fabs(double %39)
%41 = sitofp i64 %40 to double
store double %41, double* %37
%42 = add nuw nsw i64 %34, 1
%43 = load double, double* %2
%44 = getelementptr inbounds double, double* %7, i64 %42
%45 = load double, double* %44
%46 = fsub double %45, %43
%47 = tail call i64 @fabs(double %46)
%48 = sitofp i64 %47 to double
store double %48, double* %44
%49 = add nuw nsw i64 %34, 2
%50 = add i64 %35, 2
%51 = icmp eq i64 %50, %32
br i1 %51, label %52, label %33
52: ; preds = %33, %26
%53 = phi i64 [ 1, %26 ], [ %49, %33 ]
%54 = icmp eq i64 %29, 0
br i1 %54, label %62, label %55
55: ; preds = %52
%56 = load double, double* %2
%57 = getelementptr inbounds double, double* %7, i64 %53
%58 = load double, double* %57
%59 = fsub double %58, %56
%60 = tail call i64 @fabs(double %59)
%61 = sitofp i64 %60 to double
store double %61, double* %57
br label %62
62: ; preds = %55, %52, %20, %13
%63 = tail call double @nthd(i32 %0, double* %7, double %17)
store double %63, double* %3
%64 = tail call i32 @free(double* %7)
%65 = load double, double* %3
%66 = fdiv double %65, 6.745000e-01
%67 = load double, double* %2
%68 = fcmp olt double %65, 1.000000e-02
br i1 %68, label %112, label %69
69: ; preds = %62
%70 = zext i32 %0 to i64
br label %71
71: ; preds = %69, %107
%72 = phi i32 [ 0, %69 ], [ %110, %107 ]
%73 = phi double [ 6.000000e+00, %69 ], [ %109, %107 ]
%74 = phi double [ %67, %69 ], [ %108, %107 ]
br i1 %19, label %75, label %105
75: ; preds = %71
%76 = fmul double %66, %73
br label %77
77: ; preds = %75, %95
%78 = phi i64 [ 0, %75 ], [ %97, %95 ]
%79 = phi <2 x double> [ zeroinitializer, %75 ], [ %96, %95 ]
%80 = getelementptr inbounds double, double* %1, i64 %78
%81 = load double, double* %80
%82 = fsub double %81, %74
%83 = fdiv double %82, %76
%84 = fcmp olt double %83, 1.000000e+00
%85 = fcmp ogt double %83, -1.000000e+00
%86 = and i1 %84, %85
br i1 %86, label %87, label %95
87: ; preds = %77
%88 = fneg double %83
%89 = tail call double @llvm.fmuladd.f64(double %88, double %83, double 1.000000e+00)
%90 = fmul double %81, %89
%91 = insertelement <2 x double> poison, double %89, i64 0
%92 = insertelement <2 x double> %91, double %90, i64 1
%93 = shufflevector <2 x double> %91, <2 x double> poison, <2 x i32> zeroinitializer
%94 = call <2 x double> @llvm.fmuladd.v2f64(<2 x double> %92, <2 x double> %93, <2 x double> %79)
br label %95
95: ; preds = %77, %87
%96 = phi <2 x double> [ %94, %87 ], [ %79, %77 ]
%97 = add nuw nsw i64 %78, 1
%98 = icmp eq i64 %97, %70
br i1 %98, label %99, label %77
99: ; preds = %95
%100 = extractelement <2 x double> %96, i64 0
%101 = fcmp une double %100, 0.000000e+00
br i1 %101, label %102, label %105
102: ; preds = %99
%103 = extractelement <2 x double> %96, i64 1
%104 = fdiv double %103, %100
br label %107
105: ; preds = %71, %99
%106 = fmul double %73, 1.100000e+00
br label %107
107: ; preds = %102, %105
%108 = phi double [ %104, %102 ], [ %74, %105 ]
%109 = phi double [ %73, %102 ], [ %106, %105 ]
%110 = add nuw nsw i32 %72, 1
%111 = icmp eq i32 %110, 10
br i1 %111, label %112, label %71
112: ; preds = %107, %62
%113 = phi double [ %67, %62 ], [ %108, %107 ]
ret double %113
}