#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, float const, float);
int /*<<< orphan*/ memset(float *, int /*<<< orphan*/, unsigned int);
int sqrt(float);
int /*<<< orphan*/ stderr;
int frcholfactor(const float *a, float *g, unsigned int d) {
  unsigned int i, j, k;
  float 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):
<frcholfactor>:
push   %rbp
mov    %rsp,%rbp
sub    $0x30,%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    $0x2,%rax
mov    %eax,%edx
xor    %esi,%esi
call   1030 <memset@plt>
movl   $0x0,-0x20(%rbp)
mov    -0x20(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jae    134c <frcholfactor+0x21c>
mov    -0x10(%rbp),%rax
mov    -0x20(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
add    -0x20(%rbp),%ecx
mov    %ecx,%ecx
movss  (%rax,%rcx,4),%xmm0
movss  %xmm0,-0x2c(%rbp)
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x20(%rbp),%eax
jae    11f4 <frcholfactor+0xc4>
mov    -0x18(%rbp),%rax
mov    -0x20(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
add    -0x24(%rbp),%ecx
mov    %ecx,%ecx
movss  (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x20(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
add    -0x24(%rbp),%ecx
mov    %ecx,%ecx
movss  (%rax,%rcx,4),%xmm2
movss  -0x2c(%rbp),%xmm1
movd   %xmm0,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm0
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x2c(%rbp)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1191 <frcholfactor+0x61>
xorps  %xmm0,%xmm0
ucomiss -0x2c(%rbp),%xmm0
jb     1242 <frcholfactor+0x112>
mov    0x2dd8(%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
movss  (%rax,%rsi,4),%xmm0
movss  -0x2c(%rbp),%xmm1
lea    0xdcf(%rip),%rsi
call   1040 <fprintf@plt>
movl   $0xffffffff,-0x4(%rbp)
jmp    1353 <frcholfactor+0x223>
movss  -0x2c(%rbp),%xmm0
call   1050 <sqrt@plt>
mov    %eax,%ecx
mov    $0x1,%eax
cltd
idiv   %ecx
cvtsi2ss %eax,%xmm0
movss  %xmm0,-0x2c(%rbp)
movss  -0x2c(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x20(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
add    -0x20(%rbp),%ecx
mov    %ecx,%ecx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jae    1339 <frcholfactor+0x209>
xorps  %xmm0,%xmm0
movss  %xmm0,-0x30(%rbp)
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x20(%rbp),%eax
jae    12f3 <frcholfactor+0x1c3>
mov    -0x18(%rbp),%rax
mov    -0x24(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
add    -0x28(%rbp),%ecx
mov    %ecx,%ecx
movss  (%rax,%rcx,4),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x20(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
add    -0x28(%rbp),%ecx
mov    %ecx,%ecx
movss  (%rax,%rcx,4),%xmm2
movss  -0x30(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x30(%rbp)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    129d <frcholfactor+0x16d>
movss  -0x2c(%rbp),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x24(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
add    -0x20(%rbp),%ecx
mov    %ecx,%ecx
movss  (%rax,%rcx,4),%xmm1
subss  -0x30(%rbp),%xmm1
mulss  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
mov    -0x24(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
add    -0x20(%rbp),%ecx
mov    %ecx,%ecx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1282 <frcholfactor+0x152>
jmp    133e <frcholfactor+0x20e>
mov    -0x20(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x20(%rbp)
jmp    1164 <frcholfactor+0x34>
movl   $0x0,-0x4(%rbp)
mov    -0x4(%rbp),%eax
add    $0x30,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<frcholfactor>:
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    $0x2,%edx
xor    %ebp,%ebp
mov    %rsi,%rdi
xor    %esi,%esi
call   1030 <memset@plt>
test   %r15d,%r15d
je     13ad <frcholfactor+0x27d>
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 <frcholfactor+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     1389 <frcholfactor+0x259>
mov    %r8,%r13
mov    %r13d,%eax
imul   %r15d,%eax
mov    %rax,0x20(%rsp)
add    %eax,%r8d
movss  (%r14,%r8,4),%xmm0
movaps %xmm0,%xmm1
test   %r13,%r13
je     1278 <frcholfactor+0x148>
lea    -0x1(%r13),%rax
cmp    $0x3,%rax
jae    11f0 <frcholfactor+0xc0>
xor    %eax,%eax
movaps %xmm0,%xmm1
jmp    124b <frcholfactor+0x11b>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
mov    %r13,%rcx
and    $0xfffffffffffffffc,%rcx
mov    %edi,%edx
xor    %eax,%eax
movaps %xmm0,%xmm1
xchg   %ax,%ax
mov    %edx,%esi
movss  (%rbx,%rsi,4),%xmm2
mulss  %xmm2,%xmm2
subss  %xmm2,%xmm1
lea    0x1(%rdx),%esi
movss  (%rbx,%rsi,4),%xmm2
mulss  %xmm2,%xmm2
subss  %xmm2,%xmm1
lea    0x2(%rdx),%esi
movss  (%rbx,%rsi,4),%xmm2
mulss  %xmm2,%xmm2
subss  %xmm2,%xmm1
lea    0x3(%rdx),%esi
movss  (%rbx,%rsi,4),%xmm2
mulss  %xmm2,%xmm2
subss  %xmm2,%xmm1
add    $0x4,%rax
add    $0x4,%edx
cmp    %rax,%rcx
jne    1200 <frcholfactor+0xd0>
test   $0x3,%r13b
je     1278 <frcholfactor+0x148>
mov    %r13d,%ecx
and    $0x3,%ecx
add    %edi,%eax
nopl   0x0(%rax)
mov    %eax,%edx
movss  (%rbx,%rdx,4),%xmm2
mulss  %xmm2,%xmm2
subss  %xmm2,%xmm1
add    $0x1,%eax
add    $0xffffffffffffffff,%rcx
jne    1260 <frcholfactor+0x130>
xorps  %xmm2,%xmm2
ucomiss %xmm1,%xmm2
jae    138d <frcholfactor+0x25d>
mov    %edi,0xc(%rsp)
movaps %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
cvtsi2ss %eax,%xmm0
movss  %xmm0,(%rbx,%r12,4)
mov    %r14d,%r12d
mov    0x18(%rsp),%r14
lea    0x1(%r13),%r8
cmp    %rbp,%r8
jae    1190 <frcholfactor+0x60>
mov    %r13,%rcx
and    $0xfffffffffffffffe,%rcx
mov    0x8(%rsp),%r9d
mov    0x10(%rsp),%r10
jmp    1306 <frcholfactor+0x1d6>
nop
add    %r13d,%edi
movss  (%r14,%rdi,4),%xmm2
subss  %xmm1,%xmm2
mulss  %xmm0,%xmm2
movss  %xmm2,(%rbx,%rdi,4)
add    $0x1,%r10
add    %r15d,%r9d
cmp    %r15d,%r10d
je     1190 <frcholfactor+0x60>
mov    %r10d,%edi
imul   %r15d,%edi
xorps  %xmm1,%xmm1
test   %r13,%r13
je     12e0 <frcholfactor+0x1b0>
cmp    $0x1,%r13
jne    1320 <frcholfactor+0x1f0>
xor    %eax,%eax
jmp    1365 <frcholfactor+0x235>
nop
mov    %r12d,%ebp
mov    %r9d,%edx
xor    %eax,%eax
nopl   0x0(%rax,%rax,1)
mov    %edx,%esi
movss  (%rbx,%rsi,4),%xmm2
lea    -0x1(%rbp),%esi
mulss  (%rbx,%rsi,4),%xmm2
addss  %xmm1,%xmm2
lea    0x1(%rdx),%esi
movss  (%rbx,%rsi,4),%xmm1
mov    %ebp,%esi
mulss  (%rbx,%rsi,4),%xmm1
addss  %xmm2,%xmm1
add    $0x2,%rax
add    $0x2,%edx
add    $0x2,%ebp
cmp    %rax,%rcx
jne    1330 <frcholfactor+0x200>
test   $0x1,%r13b
je     12e0 <frcholfactor+0x1b0>
lea    (%rdi,%rax,1),%edx
movss  (%rbx,%rdx,4),%xmm2
add    0x20(%rsp),%eax
mulss  (%rbx,%rax,4),%xmm2
addss  %xmm2,%xmm1
jmp    12e0 <frcholfactor+0x1b0>
xor    %ebp,%ebp
jmp    13ad <frcholfactor+0x27d>
mov    0x2c4c(%rip),%rax
mov    (%rax),%edi
lea    0xc67(%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 @frcholfactor(float* %0, float* %1, i32 %2) {
%4 = alloca i32
%5 = alloca float*
%6 = alloca float*
%7 = alloca i32
%8 = alloca i32
%9 = alloca i32
%10 = alloca i32
%11 = alloca float
%12 = alloca float
store float* %0, float** %5
store float* %1, float** %6
store i32 %2, i32* %7
%13 = load float*, float** %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, 4
%19 = trunc i64 %18 to i32
%20 = call i32 @memset(float* %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 float*, float** %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 float, float* %26, i64 %32
%34 = load float, float* %33
store float %34, float* %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 float*, float** %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 float, float* %40, i64 %46
%48 = load float, float* %47
%49 = load float*, float** %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 float, float* %49, i64 %55
%57 = load float, float* %56
%58 = load float, float* %11
%59 = fneg float %48
%60 = call float @llvm.fmuladd.f32(float %59, float %57, float %58)
store float %60, float* %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 float, float* %11
%66 = fcmp ole float %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 float*, float** %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 float, float* %71, i64 %77
%79 = load float, float* %78
%80 = load float, float* %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, float %79, float %80)
store i32 -1, i32* %4
br label %161
82: ; preds = %64
%83 = load float, float* %11
%84 = call i32 @sqrt(float %83)
%85 = sdiv i32 1, %84
%86 = sitofp i32 %85 to float
store float %86, float* %11
%87 = load float, float* %11
%88 = load float*, float** %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 float, float* %88, i64 %94
store float %87, float* %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 float 0.000000e+00, float* %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 float*, float** %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 float, float* %108, i64 %114
%116 = load float, float* %115
%117 = load float*, float** %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 float, float* %117, i64 %123
%125 = load float, float* %124
%126 = load float, float* %12
%127 = call float @llvm.fmuladd.f32(float %116, float %125, float %126)
store float %127, float* %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 float, float* %11
%133 = load float*, float** %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 float, float* %133, i64 %139
%141 = load float, float* %140
%142 = load float, float* %12
%143 = fsub float %141, %142
%144 = fmul float %132, %143
%145 = load float*, float** %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 float, float* %145, i64 %151
store float %144, float* %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 @frcholfactor(float* nocapture %0, float* %1, i32 %2) {
%4 = shl i32 %2, 2
%5 = mul i32 %4, %2
%6 = tail call i32 @memset(float* %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 float, float* %0, i64 %19
%21 = load float, float* %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 float [ %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 float, float* %1, i64 %35
%37 = load float, float* %36
%38 = fneg float %37
%39 = tail call float @llvm.fmuladd.f32(float %38, float %37, float %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 float, float* %1, i64 %43
%45 = load float, float* %44
%46 = fneg float %45
%47 = tail call float @llvm.fmuladd.f32(float %46, float %45, float %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 float, float* %1, i64 %51
%53 = load float, float* %52
%54 = fneg float %53
%55 = tail call float @llvm.fmuladd.f32(float %54, float %53, float %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 float, float* %1, i64 %59
%61 = load float, float* %60
%62 = fneg float %61
%63 = tail call float @llvm.fmuladd.f32(float %62, float %61, float %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 float [ undef, %23 ], [ %63, %29 ]
%69 = phi i64 [ 0, %23 ], [ %64, %29 ]
%70 = phi float [ %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 float [ %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 float, float* %1, i64 %78
%80 = load float, float* %79
%81 = fneg float %80
%82 = tail call float @llvm.fmuladd.f32(float %81, float %80, float %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 float [ %21, %13 ], [ %68, %67 ], [ %82, %72 ]
%88 = fcmp ugt float %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, float %21, float %87)
br label %167
93: ; preds = %86
%94 = tail call i32 @sqrt(float %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 float
%99 = getelementptr inbounds float, float* %1, i64 %19
store float %98, float* %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 float [ %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 float, float* %1, i64 %118
%120 = load float, float* %119
%121 = add i32 %17, %116
%122 = zext i32 %121 to i64
%123 = getelementptr inbounds float, float* %1, i64 %122
%124 = load float, float* %123
%125 = tail call float @llvm.fmuladd.f32(float %120, float %124, float %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 float, float* %1, i64 %129
%131 = load float, float* %130
%132 = add i32 %17, %127
%133 = zext i32 %132 to i64
%134 = getelementptr inbounds float, float* %1, i64 %133
%135 = load float, float* %134
%136 = tail call float @llvm.fmuladd.f32(float %131, float %135, float %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 float [ undef, %111 ], [ %136, %112 ]
%142 = phi i64 [ 0, %111 ], [ %137, %112 ]
%143 = phi float [ 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 float, float* %1, i64 %147
%149 = load float, float* %148
%150 = add i32 %17, %145
%151 = zext i32 %150 to i64
%152 = getelementptr inbounds float, float* %1, i64 %151
%153 = load float, float* %152
%154 = tail call float @llvm.fmuladd.f32(float %149, float %153, float %143)
br label %155
155: ; preds = %144, %140, %107
%156 = phi float [ 0.000000e+00, %107 ], [ %141, %140 ], [ %154, %144 ]
%157 = add i32 %110, %16
%158 = zext i32 %157 to i64
%159 = getelementptr inbounds float, float* %0, i64 %158
%160 = load float, float* %159
%161 = fsub float %160, %156
%162 = fmul float %161, %98
%163 = getelementptr inbounds float, float* %1, i64 %158
store float %162, float* %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
}