#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 */
typedef int /*<<< orphan*/ g;

/* Variables and functions */
int /*<<< orphan*/ fprintf(int /*<<< orphan*/, char *, unsigned int,
                           unsigned int, double const, double);
int /*<<< orphan*/ memset(double *, int /*<<< orphan*/, unsigned int);
int sqrt(double);
int /*<<< orphan*/ stderr;
int drcholfactor(const double *a, double *g, unsigned int d) {
  unsigned int i, j, k;
  double s, s1;
  memset(g, 0, d * d * sizeof(g[0]));
  for (i = 0; i < d; i++) {
    s = a[i * d + i];
    for (j = 0; j < i; j++)
      s -= g[i * d + j] * g[i * d + j];
    if (s <= 0) {
      fprintf(stderr,
              "frcholfactor: matrix not positive definite a[%u][%u]=%g s=%g\n",
              i, i, a[i * d + i], s);
      return -1;
    }
    s = 1 / sqrt(s);
    g[i * d + i] = s;
    for (j = i + 1; j < d; j++) {
      s1 = 0;
      for (k = 0; k < i; k++)
        s1 += g[j * d + k] * g[i * d + k];
      g[j * d + i] = s * (a[j * d + i] - s1);
    }
  }
  return 0;
}


x86 assembly (clang, -O0):
<drcholfactor>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
mov    %rdi,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %edx,-0x1c(%rbp)
mov    -0x18(%rbp),%rdi
mov    -0x1c(%rbp),%eax
imul   -0x1c(%rbp),%eax
mov    %eax,%eax
shl    $0x3,%rax
mov    %eax,%edx
xor    %esi,%esi
call   1030 <memset@plt>
movl   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jae    1357 <drcholfactor+0x227>
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
add    -0x20(%rbp),%ecx
mov    %ecx,%ecx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x30(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x20(%rbp),%eax
jae    11fe <drcholfactor+0xce>
mov    -0x18(%rbp),%rax
mov    -0x20(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
add    -0x24(%rbp),%ecx
mov    %ecx,%ecx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x20(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
add    -0x24(%rbp),%ecx
mov    %ecx,%ecx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x30(%rbp),%xmm1
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x30(%rbp)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1191 <drcholfactor+0x61>
xorps  %xmm0,%xmm0
ucomisd -0x30(%rbp),%xmm0
jb     124d <drcholfactor+0x11d>
mov    0x2dcd(%rip),%rax
mov    (%rax),%edi
mov    -0x20(%rbp),%edx
mov    -0x20(%rbp),%ecx
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%esi
imul   -0x1c(%rbp),%esi
add    -0x20(%rbp),%esi
mov    %esi,%esi
movsd  (%rax,%rsi,8),%xmm0
movsd  -0x30(%rbp),%xmm1
lea    0xdc4(%rip),%rsi
call   1040 <fprintf@plt>
movl   $0xffffffff,-0x4(%rbp)
jmp    135e <drcholfactor+0x22e>
movsd  -0x30(%rbp),%xmm0
call   1050 <sqrt@plt>
mov    %eax,%ecx
mov    $0x1,%eax
cltd
idiv   %ecx
cvtsi2sd %eax,%xmm0
movsd  %xmm0,-0x30(%rbp)
movsd  -0x30(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x20(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
add    -0x20(%rbp),%ecx
mov    %ecx,%ecx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jae    1344 <drcholfactor+0x214>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x38(%rbp)
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x20(%rbp),%eax
jae    12fe <drcholfactor+0x1ce>
mov    -0x18(%rbp),%rax
mov    -0x24(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
add    -0x28(%rbp),%ecx
mov    %ecx,%ecx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x20(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
add    -0x28(%rbp),%ecx
mov    %ecx,%ecx
movsd  (%rax,%rcx,8),%xmm2
movsd  -0x38(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    12a8 <drcholfactor+0x178>
movsd  -0x30(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
add    -0x20(%rbp),%ecx
mov    %ecx,%ecx
movsd  (%rax,%rcx,8),%xmm1
subsd  -0x38(%rbp),%xmm1
mulsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
mov    -0x24(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
add    -0x20(%rbp),%ecx
mov    %ecx,%ecx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    128d <drcholfactor+0x15d>
jmp    1349 <drcholfactor+0x219>
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    1164 <drcholfactor+0x34>
movl   $0x0,-0x4(%rbp)
mov    -0x4(%rbp),%eax
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<drcholfactor>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x28,%rsp
mov    %edx,%r15d
mov    %rsi,%rbx
mov    %rdi,%r14
imul   %edx,%edx
shl    $0x3,%edx
xor    %ebp,%ebp
mov    %rsi,%rdi
xor    %esi,%esi
call   1030 <memset@plt>
test   %r15d,%r15d
je     13bd <drcholfactor+0x28d>
mov    %r15d,%r11d
mov    $0x1,%eax
mov    %rax,0x10(%rsp)
mov    $0x1,%r12d
mov    %r15d,0x8(%rsp)
xor    %edi,%edi
xor    %r8d,%r8d
mov    %r14,0x18(%rsp)
jmp    11ae <drcholfactor+0x7e>
cs nopw 0x0(%rax,%rax,1)
addq   $0x1,0x10(%rsp)
mov    0xc(%rsp),%edi
add    %r15d,%edi
add    %r15d,0x8(%rsp)
add    %r15d,%r12d
cmp    %r11,%r8
je     1399 <drcholfactor+0x269>
mov    %r8,%r13
mov    %r13d,%eax
imul   %r15d,%eax
mov    %rax,0x20(%rsp)
add    %eax,%r8d
movsd  (%r14,%r8,8),%xmm0
movapd %xmm0,%xmm1
test   %r13,%r13
je     1278 <drcholfactor+0x148>
lea    -0x1(%r13),%rax
cmp    $0x3,%rax
jae    11f0 <drcholfactor+0xc0>
xor    %eax,%eax
movapd %xmm0,%xmm1
jmp    124b <drcholfactor+0x11b>
cs nopw 0x0(%rax,%rax,1)
nop
mov    %r13,%rcx
and    $0xfffffffffffffffc,%rcx
mov    %edi,%edx
xor    %eax,%eax
movapd %xmm0,%xmm1
nop
mov    %edx,%esi
movsd  (%rbx,%rsi,8),%xmm2
mulsd  %xmm2,%xmm2
subsd  %xmm2,%xmm1
lea    0x1(%rdx),%esi
movsd  (%rbx,%rsi,8),%xmm2
mulsd  %xmm2,%xmm2
subsd  %xmm2,%xmm1
lea    0x2(%rdx),%esi
movsd  (%rbx,%rsi,8),%xmm2
mulsd  %xmm2,%xmm2
subsd  %xmm2,%xmm1
lea    0x3(%rdx),%esi
movsd  (%rbx,%rsi,8),%xmm2
mulsd  %xmm2,%xmm2
subsd  %xmm2,%xmm1
add    $0x4,%rax
add    $0x4,%edx
cmp    %rax,%rcx
jne    1200 <drcholfactor+0xd0>
test   $0x3,%r13b
je     1278 <drcholfactor+0x148>
mov    %r13d,%ecx
and    $0x3,%ecx
add    %edi,%eax
nopl   0x0(%rax)
mov    %eax,%edx
movsd  (%rbx,%rdx,8),%xmm2
mulsd  %xmm2,%xmm2
subsd  %xmm2,%xmm1
add    $0x1,%eax
add    $0xffffffffffffffff,%rcx
jne    1260 <drcholfactor+0x130>
xorpd  %xmm2,%xmm2
ucomisd %xmm1,%xmm2
jae    139d <drcholfactor+0x26d>
mov    %edi,0xc(%rsp)
movapd %xmm1,%xmm0
mov    %r11,%rbp
mov    %r12d,%r14d
mov    %r8,%r12
call   1050 <sqrt@plt>
mov    %rbp,%r11
lea    0x1(%rax),%ecx
cmp    $0x3,%ecx
mov    $0x0,%ecx
cmovae %ecx,%eax
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
movsd  %xmm0,(%rbx,%r12,8)
mov    %r14d,%r12d
mov    0x18(%rsp),%r14
lea    0x1(%r13),%r8
cmp    %rbp,%r8
jae    1190 <drcholfactor+0x60>
mov    %r13,%rcx
and    $0xfffffffffffffffe,%rcx
mov    0x8(%rsp),%r9d
mov    0x10(%rsp),%r10
jmp    1316 <drcholfactor+0x1e6>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
add    %r13d,%edi
movsd  (%r14,%rdi,8),%xmm2
subsd  %xmm1,%xmm2
mulsd  %xmm0,%xmm2
movsd  %xmm2,(%rbx,%rdi,8)
add    $0x1,%r10
add    %r15d,%r9d
cmp    %r15d,%r10d
je     1190 <drcholfactor+0x60>
mov    %r10d,%edi
imul   %r15d,%edi
xorpd  %xmm1,%xmm1
test   %r13,%r13
je     12f0 <drcholfactor+0x1c0>
cmp    $0x1,%r13
jne    1330 <drcholfactor+0x200>
xor    %eax,%eax
jmp    1375 <drcholfactor+0x245>
mov    %r12d,%ebp
mov    %r9d,%edx
xor    %eax,%eax
nopl   0x0(%rax,%rax,1)
mov    %edx,%esi
movsd  (%rbx,%rsi,8),%xmm2
lea    -0x1(%rbp),%esi
mulsd  (%rbx,%rsi,8),%xmm2
addsd  %xmm1,%xmm2
lea    0x1(%rdx),%esi
movsd  (%rbx,%rsi,8),%xmm1
mov    %ebp,%esi
mulsd  (%rbx,%rsi,8),%xmm1
addsd  %xmm2,%xmm1
add    $0x2,%rax
add    $0x2,%edx
add    $0x2,%ebp
cmp    %rax,%rcx
jne    1340 <drcholfactor+0x210>
test   $0x1,%r13b
je     12f0 <drcholfactor+0x1c0>
lea    (%rdi,%rax,1),%edx
movsd  (%rbx,%rdx,8),%xmm2
add    0x20(%rsp),%eax
mulsd  (%rbx,%rax,8),%xmm2
addsd  %xmm2,%xmm1
jmp    12f0 <drcholfactor+0x1c0>
xor    %ebp,%ebp
jmp    13bd <drcholfactor+0x28d>
mov    0x2c3c(%rip),%rax
mov    (%rax),%edi
lea    0xc5b(%rip),%rsi
mov    %r13d,%edx
mov    %r13d,%ecx
call   1040 <fprintf@plt>
mov    $0xffffffff,%ebp
mov    %ebp,%eax
add    $0x28,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

LLVM IR (-O2):
@stderr = global i32 0
@.str = [62 x i8] c"frcholfactor: matrix not positive definite a[%u][%u]=%g s=%g\0A\00"
define i32 @drcholfactor(double* nocapture %0, double* %1, i32 %2) {
%4 = shl i32 %2, 3
%5 = mul i32 %4, %2
%6 = tail call i32 @memset(double* %1, i32 0, i32 %5)
%7 = icmp eq i32 %2, 0
br i1 %7, label %167, label %8
8: ; preds = %3
%9 = zext i32 %2 to i64
br label %13
10: ; preds = %155, %93
%11 = add nuw nsw i64 %15, 1
%12 = icmp eq i64 %100, %9
br i1 %12, label %167, label %13
13: ; preds = %8, %10
%14 = phi i64 [ 0, %8 ], [ %100, %10 ]
%15 = phi i64 [ 1, %8 ], [ %11, %10 ]
%16 = trunc i64 %14 to i32
%17 = mul i32 %16, %2
%18 = add i32 %17, %16
%19 = zext i32 %18 to i64
%20 = getelementptr inbounds double, double* %0, i64 %19
%21 = load double, double* %20
%22 = icmp eq i64 %14, 0
br i1 %22, label %86, label %23
23: ; preds = %13
%24 = add nsw i64 %14, -1
%25 = and i64 %14, 3
%26 = icmp ult i64 %24, 3
br i1 %26, label %67, label %27
27: ; preds = %23
%28 = and i64 %14, 9223372036854775804
br label %29
29: ; preds = %29, %27
%30 = phi i64 [ 0, %27 ], [ %64, %29 ]
%31 = phi double [ %21, %27 ], [ %63, %29 ]
%32 = phi i64 [ 0, %27 ], [ %65, %29 ]
%33 = trunc i64 %30 to i32
%34 = add i32 %17, %33
%35 = zext i32 %34 to i64
%36 = getelementptr inbounds double, double* %1, i64 %35
%37 = load double, double* %36
%38 = fneg double %37
%39 = tail call double @llvm.fmuladd.f64(double %38, double %37, double %31)
%40 = trunc i64 %30 to i32
%41 = or i32 %40, 1
%42 = add i32 %17, %41
%43 = zext i32 %42 to i64
%44 = getelementptr inbounds double, double* %1, i64 %43
%45 = load double, double* %44
%46 = fneg double %45
%47 = tail call double @llvm.fmuladd.f64(double %46, double %45, double %39)
%48 = trunc i64 %30 to i32
%49 = or i32 %48, 2
%50 = add i32 %17, %49
%51 = zext i32 %50 to i64
%52 = getelementptr inbounds double, double* %1, i64 %51
%53 = load double, double* %52
%54 = fneg double %53
%55 = tail call double @llvm.fmuladd.f64(double %54, double %53, double %47)
%56 = trunc i64 %30 to i32
%57 = or i32 %56, 3
%58 = add i32 %17, %57
%59 = zext i32 %58 to i64
%60 = getelementptr inbounds double, double* %1, i64 %59
%61 = load double, double* %60
%62 = fneg double %61
%63 = tail call double @llvm.fmuladd.f64(double %62, double %61, double %55)
%64 = add nuw nsw i64 %30, 4
%65 = add i64 %32, 4
%66 = icmp eq i64 %65, %28
br i1 %66, label %67, label %29
67: ; preds = %29, %23
%68 = phi double [ undef, %23 ], [ %63, %29 ]
%69 = phi i64 [ 0, %23 ], [ %64, %29 ]
%70 = phi double [ %21, %23 ], [ %63, %29 ]
%71 = icmp eq i64 %25, 0
br i1 %71, label %86, label %72
72: ; preds = %67, %72
%73 = phi i64 [ %83, %72 ], [ %69, %67 ]
%74 = phi double [ %82, %72 ], [ %70, %67 ]
%75 = phi i64 [ %84, %72 ], [ 0, %67 ]
%76 = trunc i64 %73 to i32
%77 = add i32 %17, %76
%78 = zext i32 %77 to i64
%79 = getelementptr inbounds double, double* %1, i64 %78
%80 = load double, double* %79
%81 = fneg double %80
%82 = tail call double @llvm.fmuladd.f64(double %81, double %80, double %74)
%83 = add nuw nsw i64 %73, 1
%84 = add i64 %75, 1
%85 = icmp eq i64 %84, %25
br i1 %85, label %86, label %72
86: ; preds = %67, %72, %13
%87 = phi double [ %21, %13 ], [ %68, %67 ], [ %82, %72 ]
%88 = fcmp ugt double %87, 0.000000e+00
br i1 %88, label %93, label %89
89: ; preds = %86
%90 = trunc i64 %14 to i32
%91 = load i32, i32* @stderr
%92 = tail call i32 @fprintf(i32 %91, i8* getelementptr inbounds ([62 x i8], [62 x i8]* @.str, i64 0, i64 0), i32 %90, i32 %90, double %21, double %87)
br label %167
93: ; preds = %86
%94 = tail call i32 @sqrt(double %87)
%95 = add i32 %94, 1
%96 = icmp ult i32 %95, 3
%97 = select i1 %96, i32 %94, i32 0
%98 = sitofp i32 %97 to double
%99 = getelementptr inbounds double, double* %1, i64 %19
store double %98, double* %99
%100 = add nuw nsw i64 %14, 1
%101 = icmp ult i64 %100, %9
br i1 %101, label %102, label %10
102: ; preds = %93
%103 = and i64 %14, 1
%104 = icmp eq i64 %14, 1
%105 = and i64 %14, 9223372036854775806
%106 = icmp eq i64 %103, 0
br label %107
107: ; preds = %102, %155
%108 = phi i64 [ %164, %155 ], [ %15, %102 ]
%109 = trunc i64 %108 to i32
%110 = mul i32 %109, %2
br i1 %22, label %155, label %111
111: ; preds = %107
br i1 %104, label %140, label %112
112: ; preds = %111, %112
%113 = phi i64 [ %137, %112 ], [ 0, %111 ]
%114 = phi double [ %136, %112 ], [ 0.000000e+00, %111 ]
%115 = phi i64 [ %138, %112 ], [ 0, %111 ]
%116 = trunc i64 %113 to i32
%117 = add i32 %110, %116
%118 = zext i32 %117 to i64
%119 = getelementptr inbounds double, double* %1, i64 %118
%120 = load double, double* %119
%121 = add i32 %17, %116
%122 = zext i32 %121 to i64
%123 = getelementptr inbounds double, double* %1, i64 %122
%124 = load double, double* %123
%125 = tail call double @llvm.fmuladd.f64(double %120, double %124, double %114)
%126 = trunc i64 %113 to i32
%127 = or i32 %126, 1
%128 = add i32 %110, %127
%129 = zext i32 %128 to i64
%130 = getelementptr inbounds double, double* %1, i64 %129
%131 = load double, double* %130
%132 = add i32 %17, %127
%133 = zext i32 %132 to i64
%134 = getelementptr inbounds double, double* %1, i64 %133
%135 = load double, double* %134
%136 = tail call double @llvm.fmuladd.f64(double %131, double %135, double %125)
%137 = add nuw nsw i64 %113, 2
%138 = add i64 %115, 2
%139 = icmp eq i64 %138, %105
br i1 %139, label %140, label %112
140: ; preds = %112, %111
%141 = phi double [ undef, %111 ], [ %136, %112 ]
%142 = phi i64 [ 0, %111 ], [ %137, %112 ]
%143 = phi double [ 0.000000e+00, %111 ], [ %136, %112 ]
br i1 %106, label %155, label %144
144: ; preds = %140
%145 = trunc i64 %142 to i32
%146 = add i32 %110, %145
%147 = zext i32 %146 to i64
%148 = getelementptr inbounds double, double* %1, i64 %147
%149 = load double, double* %148
%150 = add i32 %17, %145
%151 = zext i32 %150 to i64
%152 = getelementptr inbounds double, double* %1, i64 %151
%153 = load double, double* %152
%154 = tail call double @llvm.fmuladd.f64(double %149, double %153, double %143)
br label %155
155: ; preds = %144, %140, %107
%156 = phi double [ 0.000000e+00, %107 ], [ %141, %140 ], [ %154, %144 ]
%157 = add i32 %110, %16
%158 = zext i32 %157 to i64
%159 = getelementptr inbounds double, double* %0, i64 %158
%160 = load double, double* %159
%161 = fsub double %160, %156
%162 = fmul double %161, %98
%163 = getelementptr inbounds double, double* %1, i64 %158
store double %162, double* %163
%164 = add nuw nsw i64 %108, 1
%165 = trunc i64 %164 to i32
%166 = icmp eq i32 %165, %2
br i1 %166, label %10, label %107
167: ; preds = %10, %3, %89
%168 = phi i32 [ -1, %89 ], [ 0, %3 ], [ 0, %10 ]
ret i32 %168
}