#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
int /*<<< orphan*/ dsvdcmp0(double **, int, int, double *, double **);
int /*<<< orphan*/ free(double *);
scalar_t__ malloc(int);
double sqrt(double);
double **svd_calloc_double2(int, int);
int /*<<< orphan*/ svd_errm();
int /*<<< orphan*/ svd_free_double2(double **);
void dsvdsqrtinv(double **A, int n, double *det) {
  double **V, **H, *W;
  int i, j, k;
  V = svd_calloc_double2(n, n);
  H = svd_calloc_double2(n, n);
  if (!(W = (double *)malloc(n * sizeof(double))))
    svd_errm();
  dsvdcmp0(A, n, n, W, V);
  *det = 1.;
  for (i = 0; i < n; i++)
    *det *= W[i];
  for (j = 0; j < n; j++)
    for (i = 0; i < n; i++)
      H[j][i] = A[j][i] / sqrt(W[i]);
  for (j = 0; j < n; j++)
    for (i = 0; i < n; i++) {
      for (A[i][j] = k = 0; k < n; k++)
        A[i][j] += V[i][k] * H[j][k];
    }
  svd_free_double2(V);
  svd_free_double2(H);
  free(W);
}


x86 assembly (clang, -O0):
<dsvdsqrtinv>:
push   %rbp
mov    %rsp,%rbp
sub    $0x50,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %rdx,-0x18(%rbp)
mov    -0xc(%rbp),%edi
mov    -0xc(%rbp),%esi
call   1050 <svd_calloc_double2@plt>
mov    %rax,-0x20(%rbp)
mov    -0xc(%rbp),%edi
mov    -0xc(%rbp),%esi
call   1050 <svd_calloc_double2@plt>
mov    %rax,-0x28(%rbp)
movslq -0xc(%rbp),%rax
shl    $0x3,%rax
mov    %eax,%edi
call   1070 <malloc@plt>
mov    %rax,-0x30(%rbp)
cmp    $0x0,%rax
jne    11c5 <dsvdsqrtinv+0x55>
mov    $0x0,%al
call   1080 <svd_errm@plt>
mov    -0x8(%rbp),%rdi
mov    -0xc(%rbp),%esi
mov    -0xc(%rbp),%edx
mov    -0x30(%rbp),%rcx
mov    -0x20(%rbp),%r8
call   1060 <dsvdcmp0@plt>
mov    -0x18(%rbp),%rax
movsd  0xe18(%rip),%xmm0
movsd  %xmm0,(%rax)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1226 <dsvdsqrtinv+0xb6>
mov    -0x30(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mulsd  (%rax),%xmm0
movsd  %xmm0,(%rax)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    11f3 <dsvdsqrtinv+0x83>
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    12ba <dsvdsqrtinv+0x14a>
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    12a7 <dsvdsqrtinv+0x137>
mov    -0x8(%rbp),%rax
movslq -0x38(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x48(%rbp)
mov    -0x30(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
call   1090 <sqrt@plt>
movaps %xmm0,%xmm1
movsd  -0x48(%rbp),%xmm0
divsd  %xmm1,%xmm0
mov    -0x28(%rbp),%rax
movslq -0x38(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x34(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1240 <dsvdsqrtinv+0xd0>
jmp    12ac <dsvdsqrtinv+0x13c>
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    122d <dsvdsqrtinv+0xbd>
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    138b <dsvdsqrtinv+0x21b>
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1378 <dsvdsqrtinv+0x208>
movl   $0x0,-0x3c(%rbp)
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x38(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x3c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1365 <dsvdsqrtinv+0x1f5>
mov    -0x20(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x38(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x3c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x8(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x38(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    12ff <dsvdsqrtinv+0x18f>
jmp    136a <dsvdsqrtinv+0x1fa>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    12d4 <dsvdsqrtinv+0x164>
jmp    137d <dsvdsqrtinv+0x20d>
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    12c1 <dsvdsqrtinv+0x151>
mov    -0x20(%rbp),%rdi
call   1040 <svd_free_double2@plt>
mov    -0x28(%rbp),%rdi
call   1040 <svd_free_double2@plt>
mov    -0x30(%rbp),%rdi
call   1030 <free@plt>
add    $0x50,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<dsvdsqrtinv>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x28,%rsp
mov    %rdx,%r13
mov    %esi,%ebp
mov    %rdi,%r12
mov    %esi,%edi
call   1050 <svd_calloc_double2@plt>
mov    %rax,%r15
mov    %ebp,%edi
mov    %ebp,%esi
call   1050 <svd_calloc_double2@plt>
mov    %rax,%rbx
lea    0x0(,%rbp,8),%edi
call   1070 <malloc@plt>
test   %rax,%rax
mov    %rax,(%rsp)
jne    11bc <dsvdsqrtinv+0x4c>
xor    %eax,%eax
call   1080 <svd_errm@plt>
mov    (%rsp),%rax
mov    %r12,%rdi
mov    %ebp,%esi
mov    %ebp,%edx
mov    %rax,%rcx
mov    %r15,%r8
call   1060 <dsvdcmp0@plt>
movabs $0x3ff0000000000000,%rax
mov    %rax,0x0(%r13)
test   %ebp,%ebp
mov    %rbx,%r8
jle    148e <dsvdsqrtinv+0x31e>
mov    %ebp,%r14d
lea    -0x1(%r14),%rdi
mov    %r14d,%eax
and    $0x3,%eax
cmp    $0x3,%rdi
jae    120a <dsvdsqrtinv+0x9a>
movsd  0xdfe(%rip),%xmm0
xor    %ecx,%ecx
mov    (%rsp),%rsi
jmp    1258 <dsvdsqrtinv+0xe8>
mov    %r14d,%edx
and    $0xfffffffc,%edx
movsd  0xde8(%rip),%xmm0
xor    %ecx,%ecx
mov    (%rsp),%rsi
xchg   %ax,%ax
mulsd  (%rsi,%rcx,8),%xmm0
movsd  %xmm0,0x0(%r13)
mulsd  0x8(%rsi,%rcx,8),%xmm0
movsd  %xmm0,0x0(%r13)
mulsd  0x10(%rsi,%rcx,8),%xmm0
movsd  %xmm0,0x0(%r13)
mulsd  0x18(%rsi,%rcx,8),%xmm0
movsd  %xmm0,0x0(%r13)
add    $0x4,%rcx
cmp    %rcx,%rdx
jne    1220 <dsvdsqrtinv+0xb0>
test   %rax,%rax
je     1284 <dsvdsqrtinv+0x114>
lea    (%rsi,%rcx,8),%rcx
xor    %edx,%edx
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
mulsd  (%rcx,%rdx,8),%xmm0
movsd  %xmm0,0x0(%r13)
add    $0x1,%rdx
cmp    %rdx,%rax
jne    1270 <dsvdsqrtinv+0x100>
test   %ebp,%ebp
jle    148e <dsvdsqrtinv+0x31e>
mov    %rbp,0x18(%rsp)
mov    %rdi,0x10(%rsp)
mov    %r15,0x20(%rsp)
mov    %r14d,%r13d
and    $0xfffffffe,%r13d
xor    %ebp,%ebp
xorpd  %xmm1,%xmm1
mov    (%rsp),%rcx
jmp    12e7 <dsvdsqrtinv+0x177>
xchg   %ax,%ax
movsd  %xmm2,0x8(%rsp)
call   1090 <sqrt@plt>
movsd  0x8(%rsp),%xmm2
xorpd  %xmm1,%xmm1
mov    (%rsp),%rcx
mov    %rbx,%r8
divsd  %xmm0,%xmm2
mov    (%r8,%rbp,8),%rax
movsd  %xmm2,(%rax,%r15,8)
add    $0x1,%rbp
cmp    %r14,%rbp
je     13d0 <dsvdsqrtinv+0x260>
cmpq   $0x0,0x10(%rsp)
je     13a0 <dsvdsqrtinv+0x230>
xor    %r15d,%r15d
jmp    1334 <dsvdsqrtinv+0x1c4>
nopl   0x0(%rax,%rax,1)
movsd  %xmm2,0x8(%rsp)
call   1090 <sqrt@plt>
movsd  0x8(%rsp),%xmm2
xorpd  %xmm1,%xmm1
mov    (%rsp),%rcx
mov    %rbx,%r8
divsd  %xmm0,%xmm2
mov    (%r8,%rbp,8),%rax
movsd  %xmm2,0x8(%rax,%r15,8)
add    $0x2,%r15
cmp    %r15,%r13
je     13a3 <dsvdsqrtinv+0x233>
mov    (%r12,%rbp,8),%rax
movsd  (%rax,%r15,8),%xmm2
movsd  (%rcx,%r15,8),%xmm0
ucomisd %xmm1,%xmm0
jb     1350 <dsvdsqrtinv+0x1e0>
sqrtsd %xmm0,%xmm0
jmp    136c <dsvdsqrtinv+0x1fc>
movsd  %xmm2,0x8(%rsp)
call   1090 <sqrt@plt>
movsd  0x8(%rsp),%xmm2
xorpd  %xmm1,%xmm1
mov    (%rsp),%rcx
mov    %rbx,%r8
divsd  %xmm0,%xmm2
mov    (%r8,%rbp,8),%rax
movsd  %xmm2,(%rax,%r15,8)
mov    (%r12,%rbp,8),%rax
movsd  0x8(%rax,%r15,8),%xmm2
movsd  0x8(%rcx,%r15,8),%xmm0
ucomisd %xmm1,%xmm0
jb     1300 <dsvdsqrtinv+0x190>
sqrtsd %xmm0,%xmm0
jmp    131c <dsvdsqrtinv+0x1ac>
nopl   0x0(%rax)
xor    %r15d,%r15d
test   $0x1,%r14b
je     12da <dsvdsqrtinv+0x16a>
mov    (%r12,%rbp,8),%rax
movsd  (%rax,%r15,8),%xmm2
movsd  (%rcx,%r15,8),%xmm0
ucomisd %xmm1,%xmm0
jb     12b0 <dsvdsqrtinv+0x140>
sqrtsd %xmm0,%xmm0
jmp    12cc <dsvdsqrtinv+0x15c>
cmpl   $0x0,0x18(%rsp)
mov    0x20(%rsp),%r15
mov    0x10(%rsp),%r9
jle    148e <dsvdsqrtinv+0x31e>
mov    %r14d,%eax
and    $0xfffffffe,%eax
xor    %ecx,%ecx
jmp    13fd <dsvdsqrtinv+0x28d>
nop
add    $0x1,%rcx
cmp    %r14,%rcx
je     148e <dsvdsqrtinv+0x31e>
mov    (%r8,%rcx,8),%rdx
xor    %esi,%esi
jmp    1419 <dsvdsqrtinv+0x2a9>
cs nopw 0x0(%rax,%rax,1)
nop
add    $0x1,%rsi
cmp    %r14,%rsi
je     13f0 <dsvdsqrtinv+0x280>
mov    (%r12,%rsi,8),%rdi
movq   $0x0,(%rdi,%rcx,8)
mov    (%r15,%rsi,8),%rbp
xorpd  %xmm0,%xmm0
xor    %ebx,%ebx
test   %r9,%r9
je     1472 <dsvdsqrtinv+0x302>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
movsd  0x0(%rbp,%rbx,8),%xmm1
mulsd  (%rdx,%rbx,8),%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,(%rdi,%rcx,8)
movsd  0x8(%rbp,%rbx,8),%xmm0
mulsd  0x8(%rdx,%rbx,8),%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rdi,%rcx,8)
add    $0x2,%rbx
cmp    %rbx,%rax
jne    1440 <dsvdsqrtinv+0x2d0>
test   $0x1,%r14b
je     1410 <dsvdsqrtinv+0x2a0>
movsd  0x0(%rbp,%rbx,8),%xmm1
mulsd  (%rdx,%rbx,8),%xmm1
addsd  %xmm0,%xmm1
movsd  %xmm1,(%rdi,%rcx,8)
jmp    1410 <dsvdsqrtinv+0x2a0>
mov    %r15,%rdi
mov    %r8,%rbx
call   1040 <svd_free_double2@plt>
mov    %rbx,%rdi
call   1040 <svd_free_double2@plt>
mov    (%rsp),%rdi
add    $0x28,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
jmp    1030 <free@plt>

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

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