#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);
int solvps(double *a, double *b, int n) {
  double *p, *q, *r, *s, t;
  int j, k;
  for (j = 0, p = a; j < n; ++j, p += n + 1) {
    for (q = a + j * n; q < p; ++q)
      *p -= *q * *q;
    if (*p <= 0.)
      return -1;
    *p = sqrt(*p);
    for (k = j + 1, q = p + n; k < n; ++k, q += n) {
      for (r = a + j * n, s = a + k * n, t = 0.; r < p;)
        t += *r++ * *s++;
      *q -= t;
      *q /= *p;
    }
  }
  for (j = 0, p = a; j < n; ++j, p += n + 1) {
    for (k = 0, q = a + j * n; k < j;)
      b[j] -= b[k++] * *q++;
    b[j] /= *p;
  }
  for (j = n - 1, p = a + n * n - 1; j >= 0; --j, p -= n + 1) {
    for (k = j + 1, q = p + n; k < n; q += n)
      b[j] -= b[k++] * *q;
    b[j] /= *p;
  }
  return 0;
}


x86 assembly (clang, -O0):
<solvps>:
push   %rbp
mov    %rsp,%rbp
sub    $0x50,%rsp
mov    %rdi,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %edx,-0x1c(%rbp)
movl   $0x0,-0x4c(%rbp)
mov    -0x10(%rbp),%rax
mov    %rax,-0x28(%rbp)
mov    -0x4c(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    130b <solvps+0x1fb>
mov    -0x10(%rbp),%rax
mov    -0x4c(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x30(%rbp)
mov    -0x30(%rbp),%rax
cmp    -0x28(%rbp),%rax
jae    11b1 <solvps+0xa1>
mov    -0x30(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x30(%rbp),%rax
movsd  (%rax),%xmm2
mov    -0x28(%rbp),%rax
movsd  (%rax),%xmm1
movq   %xmm0,%rcx
movabs $0x8000000000000000,%rdx
xor    %rdx,%rcx
movq   %rcx,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax)
mov    -0x30(%rbp),%rax
add    $0x8,%rax
mov    %rax,-0x30(%rbp)
jmp    1157 <solvps+0x47>
mov    -0x28(%rbp),%rax
xorps  %xmm0,%xmm0
ucomisd (%rax),%xmm0
jb     11ce <solvps+0xbe>
movl   $0xffffffff,-0x4(%rbp)
jmp    1508 <solvps+0x3f8>
mov    -0x28(%rbp),%rax
movsd  (%rax),%xmm0
call   1030 <sqrt@plt>
mov    -0x28(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x50(%rbp)
mov    -0x28(%rbp),%rax
movslq -0x1c(%rbp),%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x30(%rbp)
mov    -0x50(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    12e0 <solvps+0x1d0>
mov    -0x10(%rbp),%rax
mov    -0x4c(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x38(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x50(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x40(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x38(%rbp),%rax
cmp    -0x28(%rbp),%rax
jae    1290 <solvps+0x180>
mov    -0x38(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x38(%rbp)
movsd  (%rax),%xmm0
mov    -0x40(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x40(%rbp)
movsd  (%rax),%xmm2
movsd  -0x48(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x48(%rbp)
jmp    1245 <solvps+0x135>
movsd  -0x48(%rbp),%xmm1
mov    -0x30(%rbp),%rax
movsd  (%rax),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,(%rax)
mov    -0x28(%rbp),%rax
movsd  (%rax),%xmm1
mov    -0x30(%rbp),%rax
movsd  (%rax),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rax)
mov    -0x50(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x50(%rbp)
mov    -0x1c(%rbp),%ecx
mov    -0x30(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x30(%rbp)
jmp    11ff <solvps+0xef>
jmp    12e5 <solvps+0x1d5>
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
mov    -0x1c(%rbp),%ecx
add    $0x1,%ecx
mov    -0x28(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x28(%rbp)
jmp    1132 <solvps+0x22>
movl   $0x0,-0x4c(%rbp)
mov    -0x10(%rbp),%rax
mov    %rax,-0x28(%rbp)
mov    -0x4c(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    13f6 <solvps+0x2e6>
movl   $0x0,-0x50(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x4c(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x30(%rbp)
mov    -0x50(%rbp),%eax
cmp    -0x4c(%rbp),%eax
jge    13b2 <solvps+0x2a2>
mov    -0x18(%rbp),%rax
mov    -0x50(%rbp),%ecx
mov    %ecx,%edx
add    $0x1,%edx
mov    %edx,-0x50(%rbp)
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x30(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x30(%rbp)
movsd  (%rax),%xmm2
mov    -0x18(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movq   %xmm0,%rdx
movabs $0x8000000000000000,%rsi
xor    %rsi,%rdx
movq   %rdx,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    1346 <solvps+0x236>
mov    -0x28(%rbp),%rax
movsd  (%rax),%xmm1
mov    -0x18(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
mov    -0x1c(%rbp),%ecx
add    $0x1,%ecx
mov    -0x28(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x28(%rbp)
jmp    131a <solvps+0x20a>
mov    -0x1c(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x4c(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x1c(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
add    $0xfffffffffffffff8,%rax
mov    %rax,-0x28(%rbp)
cmpl   $0x0,-0x4c(%rbp)
jl     1501 <solvps+0x3f1>
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x50(%rbp)
mov    -0x28(%rbp),%rax
movslq -0x1c(%rbp),%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x30(%rbp)
mov    -0x50(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    14b8 <solvps+0x3a8>
mov    -0x18(%rbp),%rax
mov    -0x50(%rbp),%ecx
mov    %ecx,%edx
add    $0x1,%edx
mov    %edx,-0x50(%rbp)
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x30(%rbp),%rax
movsd  (%rax),%xmm2
mov    -0x18(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movq   %xmm0,%rdx
movabs $0x8000000000000000,%rsi
xor    %rsi,%rdx
movq   %rdx,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x1c(%rbp),%ecx
mov    -0x30(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x30(%rbp)
jmp    1442 <solvps+0x332>
mov    -0x28(%rbp),%rax
movsd  (%rax),%xmm1
mov    -0x18(%rbp),%rax
movslq -0x4c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x4c(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x4c(%rbp)
mov    -0x1c(%rbp),%ecx
add    $0x1,%ecx
mov    -0x28(%rbp),%rax
movslq %ecx,%rdx
xor    %ecx,%ecx
sub    %rdx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x28(%rbp)
jmp    141c <solvps+0x30c>
movl   $0x0,-0x4(%rbp)
mov    -0x4(%rbp),%eax
add    $0x50,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<solvps>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x18,%rsp
test   %edx,%edx
jle    142b <solvps+0x31b>
mov    %rsi,%r12
mov    %rdi,%r13
mov    %edx,%ebx
mov    %rdx,(%rsp)
lea    0x1(%rdx),%r8d
mov    $0x1,%ebp
xor    %r9d,%r9d
xorpd  %xmm2,%xmm2
mov    %rdi,%r15
mov    %r8,0x8(%rsp)
jmp    1161 <solvps+0x51>
nopl   0x0(%rax)
lea    (%r15,%r8,8),%r15
add    $0x1,%rbp
cmp    %rbx,%r9
je     1265 <solvps+0x155>
mov    %r9,%rax
imul   %rbx,%rax
lea    0x0(,%rax,8),%r14
add    %r13,%r14
movsd  (%r15),%xmm0
cmp    %r15,%r14
jae    119a <solvps+0x8a>
mov    %r14,%rax
movsd  (%rax),%xmm1
mulsd  %xmm1,%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,(%r15)
add    $0x8,%rax
cmp    %r15,%rax
jb     1180 <solvps+0x70>
ucomisd %xmm0,%xmm2
jae    125b <solvps+0x14b>
ucomisd %xmm2,%xmm0
jb     11b0 <solvps+0xa0>
sqrtsd %xmm0,%xmm0
jmp    11c8 <solvps+0xb8>
mov    %r9,0x10(%rsp)
call   1030 <sqrt@plt>
xorpd  %xmm2,%xmm2
mov    0x10(%rsp),%r9
mov    0x8(%rsp),%r8
movsd  %xmm0,(%r15)
add    $0x1,%r9
cmp    %rbx,%r9
jae    1150 <solvps+0x40>
mov    %rbp,%rax
mov    %r15,%rcx
jmp    1216 <solvps+0x106>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
lea    (%rcx,%rbx,8),%rcx
movsd  (%rcx),%xmm1
subsd  %xmm0,%xmm1
movsd  %xmm1,(%rcx)
divsd  (%r15),%xmm1
movsd  %xmm1,(%rcx)
add    $0x1,%rax
cmp    %rbx,%rax
je     1150 <solvps+0x40>
xorpd  %xmm0,%xmm0
cmp    %r15,%r14
jae    11f0 <solvps+0xe0>
mov    %rax,%rdx
imul   %rbx,%rdx
lea    0x0(,%rdx,8),%rdx
add    %r13,%rdx
xorpd  %xmm0,%xmm0
mov    %r14,%rsi
nopl   0x0(%rax,%rax,1)
movsd  (%rsi),%xmm1
mulsd  (%rdx),%xmm1
add    $0x8,%rsi
add    $0x8,%rdx
addsd  %xmm1,%xmm0
cmp    %r15,%rsi
jb     1240 <solvps+0x130>
jmp    11f0 <solvps+0xe0>
mov    $0xffffffff,%eax
jmp    142d <solvps+0x31d>
mov    (%rsp),%r9
test   %r9d,%r9d
jle    142b <solvps+0x31b>
xor    %eax,%eax
mov    %r13,%rcx
jmp    12a1 <solvps+0x191>
nopl   0x0(%rax)
movsd  (%r12),%xmm0
divsd  (%rcx),%xmm0
movsd  %xmm0,(%r12,%rax,8)
add    $0x1,%rax
lea    (%rcx,%r8,8),%rcx
cmp    %rbx,%rax
je     1338 <solvps+0x228>
test   %rax,%rax
je     1280 <solvps+0x170>
mov    %rax,%rdx
imul   %rbx,%rdx
lea    0x0(,%rdx,8),%rdx
add    %r13,%rdx
movsd  (%r12,%rax,8),%xmm0
cmp    $0x1,%rax
jne    12d0 <solvps+0x1c0>
xor    %esi,%esi
jmp    1317 <solvps+0x207>
nopl   0x0(%rax,%rax,1)
mov    %rax,%rbp
and    $0xfffffffffffffffe,%rbp
xor    %esi,%esi
nopl   0x0(%rax)
movsd  (%r12,%rsi,8),%xmm1
mulsd  (%rdx),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,(%r12,%rax,8)
movsd  0x8(%r12,%rsi,8),%xmm1
add    $0x2,%rsi
mulsd  0x8(%rdx),%xmm1
add    $0x10,%rdx
subsd  %xmm1,%xmm0
movsd  %xmm0,(%r12,%rax,8)
cmp    %rsi,%rbp
jne    12e0 <solvps+0x1d0>
test   $0x1,%al
je     1286 <solvps+0x176>
movsd  (%r12,%rsi,8),%xmm1
mulsd  (%rdx),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,(%r12,%rax,8)
jmp    1286 <solvps+0x176>
test   %r9d,%r9d
jle    142b <solvps+0x31b>
mov    %r9d,%eax
imul   %eax,%eax
lea    -0x8(,%rax,8),%rcx
ff
add    %r13,%rcx
not    %r9d
movslq %r9d,%r8
mov    %rbx,%rdx
shl    $0x4,%rdx
xor    %esi,%esi
mov    %rbx,%rdi
jmp    138b <solvps+0x27b>
cs nopw 0x0(%rax,%rax,1)
divsd  (%rcx),%xmm0
movsd  %xmm0,(%r12,%rdi,8)
lea    (%rcx,%r8,8),%rcx
add    $0x1,%esi
cmp    $0x1,%r9
jle    142b <solvps+0x31b>
mov    %rdi,%r9
add    $0xffffffffffffffff,%rdi
movsd  -0x8(%r12,%r9,8),%xmm0
cmp    %rbx,%r9
jge    1370 <solvps+0x260>
test   $0x1,%sil
jne    13c0 <solvps+0x2b0>
mov    %r9,%rbp
mov    %rcx,%rax
cmp    $0x1,%esi
jne    13f0 <solvps+0x2e0>
jmp    1370 <solvps+0x260>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
lea    (%rcx,%rbx,8),%rax
lea    0x1(%r9),%rbp
movsd  (%r12,%r9,8),%xmm1
mulsd  (%rcx,%rbx,8),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,(%r12,%rdi,8)
cmp    $0x1,%esi
je     1370 <solvps+0x260>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movsd  (%r12,%rbp,8),%xmm1
mulsd  (%rax,%rbx,8),%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,(%r12,%rdi,8)
movsd  0x8(%r12,%rbp,8),%xmm1
add    $0x2,%rbp
mulsd  (%rax,%rdx,1),%xmm1
add    %rdx,%rax
subsd  %xmm1,%xmm0
movsd  %xmm0,(%r12,%rdi,8)
cmp    %ebp,%ebx
jne    13f0 <solvps+0x2e0>
jmp    1370 <solvps+0x260>
xor    %eax,%eax
add    $0x18,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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