#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 /*<<< orphan*/ trnm(double *, int);
int psinv(double *v, int n) {
  double z, *p, *q, *r, *s, *t;
  int j, k;
  for (j = 0, p = v; j < n; ++j, p += n + 1) {
    for (q = v + 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 = v + j * n, s = v + k * n, z = 0.; r < p;)
        z += *r++ * *s++;
      *q -= z;
      *q /= *p;
    }
  }
  trnm(v, n);
  for (j = 0, p = v; j < n; ++j, p += n + 1) {
    *p = 1. / *p;
    for (q = v + j, t = v; q < p; t += n + 1, q += n) {
      for (s = q, r = t, z = 0.; s < p; s += n)
        z -= *s * *r++;
      *q = z * *p;
    }
  }
  for (j = 0, p = v; j < n; ++j, p += n + 1) {
    for (q = v + j, t = p - j; q <= p; q += n) {
      for (k = j, r = p, s = q, z = 0.; k < n; ++k)
        z += *r++ * *s++;
      *t++ = (*q = z);
    }
  }
  return 0;
}


x86 assembly (clang, -O0):
<psinv>:
push   %rbp
mov    %rsp,%rbp
sub    $0x50,%rsp
mov    %rdi,-0x10(%rbp)
mov    %esi,-0x14(%rbp)
movl   $0x0,-0x4c(%rbp)
mov    -0x10(%rbp),%rax
mov    %rax,-0x28(%rbp)
mov    -0x4c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1317 <psinv+0x1f7>
mov    -0x10(%rbp),%rax
mov    -0x4c(%rbp),%ecx
imul   -0x14(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x30(%rbp)
mov    -0x30(%rbp),%rax
cmp    -0x28(%rbp),%rax
jae    11bd <psinv+0x9d>
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    1163 <psinv+0x43>
mov    -0x28(%rbp),%rax
xorps  %xmm0,%xmm0
ucomisd (%rax),%xmm0
jb     11da <psinv+0xba>
movl   $0xffffffff,-0x4(%rbp)
jmp    15a5 <psinv+0x485>
mov    -0x28(%rbp),%rax
movsd  (%rax),%xmm0
call   1040 <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 -0x14(%rbp),%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x30(%rbp)
mov    -0x50(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    12ec <psinv+0x1cc>
mov    -0x10(%rbp),%rax
mov    -0x4c(%rbp),%ecx
imul   -0x14(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x38(%rbp)
mov    -0x10(%rbp),%rax
mov    -0x50(%rbp),%ecx
imul   -0x14(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x40(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x38(%rbp),%rax
cmp    -0x28(%rbp),%rax
jae    129c <psinv+0x17c>
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  -0x20(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x20(%rbp)
jmp    1251 <psinv+0x131>
movsd  -0x20(%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    -0x14(%rbp),%ecx
mov    -0x30(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x30(%rbp)
jmp    120b <psinv+0xeb>
jmp    12f1 <psinv+0x1d1>
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
mov    -0x14(%rbp),%ecx
add    $0x1,%ecx
mov    -0x28(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x28(%rbp)
jmp    113e <psinv+0x1e>
mov    -0x10(%rbp),%rdi
mov    -0x14(%rbp),%esi
call   1030 <trnm@plt>
movl   $0x0,-0x4c(%rbp)
mov    -0x10(%rbp),%rax
mov    %rax,-0x28(%rbp)
mov    -0x4c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1475 <psinv+0x355>
mov    -0x28(%rbp),%rax
movsd  0xcb6(%rip),%xmm0
divsd  (%rax),%xmm0
mov    -0x28(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x10(%rbp),%rax
movslq -0x4c(%rbp),%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x30(%rbp)
mov    -0x10(%rbp),%rax
mov    %rax,-0x48(%rbp)
mov    -0x30(%rbp),%rax
cmp    -0x28(%rbp),%rax
jae    144a <psinv+0x32a>
mov    -0x30(%rbp),%rax
mov    %rax,-0x40(%rbp)
mov    -0x48(%rbp),%rax
mov    %rax,-0x38(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x40(%rbp),%rax
cmp    -0x28(%rbp),%rax
jae    1403 <psinv+0x2e3>
mov    -0x40(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x38(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x38(%rbp)
movsd  (%rax),%xmm2
movsd  -0x20(%rbp),%xmm1
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x14(%rbp),%ecx
mov    -0x40(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x40(%rbp)
jmp    1397 <psinv+0x277>
movsd  -0x20(%rbp),%xmm0
mov    -0x28(%rbp),%rax
mulsd  (%rax),%xmm0
mov    -0x30(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x14(%rbp),%ecx
add    $0x1,%ecx
mov    -0x48(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x48(%rbp)
mov    -0x14(%rbp),%ecx
mov    -0x30(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x30(%rbp)
jmp    1371 <psinv+0x251>
jmp    144f <psinv+0x32f>
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
mov    -0x14(%rbp),%ecx
add    $0x1,%ecx
mov    -0x28(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x28(%rbp)
jmp    1332 <psinv+0x212>
movl   $0x0,-0x4c(%rbp)
mov    -0x10(%rbp),%rax
mov    %rax,-0x28(%rbp)
mov    -0x4c(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    159e <psinv+0x47e>
mov    -0x10(%rbp),%rax
movslq -0x4c(%rbp),%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x30(%rbp)
mov    -0x28(%rbp),%rax
movslq -0x4c(%rbp),%rdx
xor    %ecx,%ecx
sub    %rdx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x48(%rbp)
mov    -0x30(%rbp),%rax
cmp    -0x28(%rbp),%rax
ja     1573 <psinv+0x453>
mov    -0x4c(%rbp),%eax
mov    %eax,-0x50(%rbp)
mov    -0x28(%rbp),%rax
mov    %rax,-0x38(%rbp)
mov    -0x30(%rbp),%rax
mov    %rax,-0x40(%rbp)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x50(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1539 <psinv+0x419>
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  -0x20(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x20(%rbp)
mov    -0x50(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x50(%rbp)
jmp    14e7 <psinv+0x3c7>
movsd  -0x20(%rbp),%xmm0
mov    -0x30(%rbp),%rax
movsd  %xmm0,(%rax)
mov    -0x48(%rbp),%rax
mov    %rax,%rcx
add    $0x8,%rcx
mov    %rcx,-0x48(%rbp)
movsd  %xmm0,(%rax)
mov    -0x14(%rbp),%ecx
mov    -0x30(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x30(%rbp)
jmp    14bb <psinv+0x39b>
jmp    1578 <psinv+0x458>
mov    -0x4c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x4c(%rbp)
mov    -0x14(%rbp),%ecx
add    $0x1,%ecx
mov    -0x28(%rbp),%rax
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x28(%rbp)
jmp    1484 <psinv+0x364>
movl   $0x0,-0x4(%rbp)
mov    -0x4(%rbp),%eax
add    $0x50,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<psinv>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x18,%rsp
test   %esi,%esi
jle    125b <psinv+0x13b>
mov    %esi,%r12d
mov    %rsi,0x10(%rsp)
lea    0x1(%rsi),%ebx
mov    $0x1,%r13d
xor    %ebp,%ebp
xorpd  %xmm2,%xmm2
mov    %rdi,%r15
mov    %rdi,0x8(%rsp)
mov    %rbx,(%rsp)
jmp    1171 <psinv+0x51>
nopl   0x0(%rax,%rax,1)
lea    (%r15,%rbx,8),%r15
add    $0x1,%r13
cmp    %r12,%rbp
je     126f <psinv+0x14f>
mov    %rbp,%rax
imul   %r12,%rax
lea    (%rdi,%rax,8),%r14
movsd  (%r15),%xmm0
cmp    %r15,%r14
jae    11aa <psinv+0x8a>
mov    %r14,%rax
nopl   0x0(%rax)
movsd  (%rax),%xmm1
mulsd  %xmm1,%xmm1
subsd  %xmm1,%xmm0
movsd  %xmm0,(%r15)
add    $0x8,%rax
cmp    %r15,%rax
jb     1190 <psinv+0x70>
ucomisd %xmm0,%xmm2
jae    1265 <psinv+0x145>
ucomisd %xmm2,%xmm0
jb     11c0 <psinv+0xa0>
sqrtsd %xmm0,%xmm0
jmp    11d2 <psinv+0xb2>
call   1040 <sqrt@plt>
xorpd  %xmm2,%xmm2
mov    (%rsp),%rbx
mov    0x8(%rsp),%rdi
movsd  %xmm0,(%r15)
add    $0x1,%rbp
cmp    %r12,%rbp
jae    1160 <psinv+0x40>
mov    %r13,%rax
mov    %r15,%rcx
jmp    1216 <psinv+0xf6>
nopl   0x0(%rax,%rax,1)
lea    (%rcx,%r12,8),%rcx
movsd  (%rcx),%xmm1
subsd  %xmm0,%xmm1
movsd  %xmm1,(%rcx)
divsd  (%r15),%xmm1
movsd  %xmm1,(%rcx)
add    $0x1,%rax
cmp    %r12,%rax
je     1160 <psinv+0x40>
xorpd  %xmm0,%xmm0
cmp    %r15,%r14
jae    11f0 <psinv+0xd0>
mov    %rax,%rdx
imul   %r12,%rdx
lea    (%rdi,%rdx,8),%rdx
xorpd  %xmm0,%xmm0
mov    %r14,%rsi
cs nopw 0x0(%rax,%rax,1)
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 <psinv+0x120>
jmp    11f0 <psinv+0xd0>
call   1030 <trnm@plt>
jmp    146b <psinv+0x34b>
mov    $0xffffffff,%eax
jmp    146d <psinv+0x34d>
mov    0x10(%rsp),%r13
mov    %r13d,%esi
mov    %rdi,%rbx
call   1030 <trnm@plt>
test   %r13d,%r13d
jle    146b <psinv+0x34b>
lea    0x0(,%r12,8),%rax
xor    %r9d,%r9d
movsd  0xd65(%rip),%xmm0
mov    %rbx,%r8
mov    %rbx,%rsi
mov    %rbx,%r10
mov    (%rsp),%rcx
jmp    12c1 <psinv+0x1a1>
nopw   0x0(%rax,%rax,1)
add    $0x1,%r9
lea    (%rsi,%rcx,8),%rsi
add    $0x8,%r8
cmp    %r12,%r9
je     131f <psinv+0x1ff>
movapd %xmm0,%xmm1
divsd  (%rsi),%xmm1
movsd  %xmm1,(%rsi)
lea    (%r10,%r9,8),%rdx
cmp    %rsi,%rdx
jae    12b0 <psinv+0x190>
mov    %r8,%rdx
mov    %r10,%rbx
nopl   0x0(%rax)
xorpd  %xmm1,%xmm1
mov    %rdx,%rdi
mov    %rbx,%rbp
nopw   0x0(%rax,%rax,1)
movsd  (%rdi),%xmm2
mulsd  0x0(%rbp),%xmm2
add    $0x8,%rbp
subsd  %xmm2,%xmm1
add    %rax,%rdi
cmp    %rsi,%rdi
jb     12f0 <psinv+0x1d0>
mulsd  (%rsi),%xmm1
movsd  %xmm1,(%rdx)
lea    (%rbx,%rcx,8),%rbx
add    %rax,%rdx
cmp    %rsi,%rdx
jb     12e0 <psinv+0x1c0>
jmp    12b0 <psinv+0x190>
test   %r13d,%r13d
jle    146b <psinv+0x34b>
xor    %r8d,%r8d
mov    %r13d,%ecx
xor    %r9d,%r9d
mov    %r10,%r15
jmp    1361 <psinv+0x241>
cs nopw 0x0(%rax,%rax,1)
add    $0x1,%r9
mov    (%rsp),%rax
lea    (%r15,%rax,8),%r15
add    $0xffffffff,%r8d
add    $0x3,%cl
cmp    %r12,%r9
mov    0x8(%rsp),%r10
je     146b <psinv+0x34b>
movzbl %cl,%ecx
and    $0x3,%ecx
lea    (%r10,%r9,8),%rbp
cmp    %r15,%rbp
ja     1340 <psinv+0x220>
mov    %r13d,%r11d
sub    %r9d,%r11d
mov    %r9d,%r14d
not    %r14d
add    %r13d,%r14d
lea    0x0(,%r9,8),%rax
mov    %r15,%rsi
sub    %rax,%rsi
mov    %r8d,%r10d
neg    %r10d
jmp    13b7 <psinv+0x297>
cs nopw 0x0(%rax,%rax,1)
nop
movsd  %xmm0,0x0(%rbp)
movsd  %xmm0,(%rsi)
add    $0x8,%rsi
lea    0x0(%rbp,%r12,8),%rbp
cmp    %r15,%rbp
ja     1340 <psinv+0x220>
test   $0x3,%r11b
je     1400 <psinv+0x2e0>
xor    %eax,%eax
xorpd  %xmm0,%xmm0
mov    %rbp,%rdi
mov    %r15,%rdx
nopl   0x0(%rax)
movsd  (%rdx),%xmm1
add    $0x8,%rdx
mulsd  (%rdi),%xmm1
add    $0x8,%rdi
addsd  %xmm1,%xmm0
add    $0x1,%eax
cmp    %eax,%ecx
jne    13d0 <psinv+0x2b0>
add    %r10d,%eax
cmp    $0x3,%r14d
jb     13a0 <psinv+0x280>
jmp    1413 <psinv+0x2f3>
cs nopw 0x0(%rax,%rax,1)
xorpd  %xmm0,%xmm0
mov    %r9d,%eax
mov    %rbp,%rdi
mov    %r15,%rdx
cmp    $0x3,%r14d
jb     13a0 <psinv+0x280>
mov    %r13d,%ebx
sub    %eax,%ebx
xor    %eax,%eax
nopw   0x0(%rax,%rax,1)
movsd  (%rdx,%rax,8),%xmm1
movsd  0x8(%rdx,%rax,8),%xmm2
mulsd  (%rdi,%rax,8),%xmm1
mulsd  0x8(%rdi,%rax,8),%xmm2
addsd  %xmm0,%xmm1
movsd  0x10(%rdx,%rax,8),%xmm3
mulsd  0x10(%rdi,%rax,8),%xmm3
addsd  %xmm1,%xmm2
movsd  0x18(%rdx,%rax,8),%xmm0
mulsd  0x18(%rdi,%rax,8),%xmm0
addsd  %xmm2,%xmm3
addsd  %xmm3,%xmm0
add    $0x4,%rax
cmp    %eax,%ebx
jne    1420 <psinv+0x300>
jmp    13a0 <psinv+0x280>
xor    %eax,%eax
add    $0x18,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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