#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*/ cov_mat(double *, double *, int, int);
int /*<<< orphan*/ denormF(double *, double *, double *);
int /*<<< orphan*/ free(double *);
int /*<<< orphan*/ lap_eig(double *, double *, int);
int /*<<< orphan*/ lin_fm(double const *, double *, int const *, int);
int /*<<< orphan*/ lin_fmN(double const *, double *, int const *, int, double *,
                           double *);
scalar_t__ malloc(int);
int /*<<< orphan*/ normu(double const *, int const *, int, double *, double *);
int /*<<< orphan*/ scalmul(double *, double const, int, int);
int /*<<< orphan*/ singulF(double *);
int /*<<< orphan*/ svduv(double *, double *, double *, int, double *, int);
int /*<<< orphan*/ trnm(double *, int);
void u2fw(const double *u, const int *inl, const double *w, int len, double *F,
          double *buffer) {
  double A1[3], A2[3];
  double *Z, V[9 * 9], U[8 * 8], D[9], *p;
  int i, j;
  if (buffer == NULL)
    Z = (double *)malloc(sizeof(double) * 9 * len);
  else
    Z = buffer;
  if (len > 8) {
    normu(u, inl, len, A1, A2);
    lin_fmN(u, Z, inl, len, A1, A2);
    for (i = 0; i < len; i++) {
      j = inl[i];
      scalmul(Z + 9 * i, w[j], 9, 1);
    }
    cov_mat(V, Z, len, 9);
    lap_eig(V, D, 9);
    trnm(V, 9);
  } else {
    lin_fm(u, Z, inl, len);
    for (i = 0; i < len; i++) {
      j = inl[i];
      scalmul(Z + i, w[j], 9, 9);
    }
    svduv(D, Z, V, 9, U, 8);
  }
  if (len > 8) {
    j = 0;
    for (i = 1; i < 9; i++)
      if (D[i] < D[j])
        j = i;
    p = V + j;
  } else
    p = V + 8;
  for (i = 0; i < 9; i++) {
    F[i] = *p;
    p += 9;
  }
  singulF(F);
  if (len > 8)
    denormF(F, A1, A2);
  if (buffer == NULL)
    free(Z);
}


x86 assembly (clang, -O0):
<u2fw>:
push   %rbp
mov    %rsp,%rbp
sub    $0x560,%rsp
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9,-0x30(%rbp)
cmpq   $0x0,-0x30(%rbp)
jne    1205 <u2fw+0x45>
movslq -0x1c(%rbp),%rax
imul   $0x48,%rax,%rax
mov    %eax,%edi
call   10a0 <malloc@plt>
mov    %rax,-0x78(%rbp)
jmp    120d <u2fw+0x4d>
mov    -0x30(%rbp),%rax
mov    %rax,-0x78(%rbp)
cmpl   $0x8,-0x1c(%rbp)
jle    1305 <u2fw+0x145>
mov    -0x8(%rbp),%rdi
mov    -0x10(%rbp),%rsi
mov    -0x1c(%rbp),%edx
lea    -0x50(%rbp),%rcx
lea    -0x70(%rbp),%r8
call   10c0 <normu@plt>
mov    -0x8(%rbp),%rdi
mov    -0x78(%rbp),%rsi
mov    -0x10(%rbp),%rdx
mov    -0x1c(%rbp),%ecx
lea    -0x50(%rbp),%r8
lea    -0x70(%rbp),%r9
call   1050 <lin_fmN@plt>
movl   $0x0,-0x55c(%rbp)
mov    -0x55c(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    12bf <u2fw+0xff>
mov    -0x10(%rbp),%rax
movslq -0x55c(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x560(%rbp)
mov    -0x78(%rbp),%rdi
imul   $0x9,-0x55c(%rbp),%eax
cltq
shl    $0x3,%rax
add    %rax,%rdi
mov    -0x18(%rbp),%rax
movslq -0x560(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    $0x9,%esi
mov    $0x1,%edx
call   1080 <scalmul@plt>
mov    -0x55c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x55c(%rbp)
jmp    1255 <u2fw+0x95>
lea    -0x300(%rbp),%rdi
mov    -0x78(%rbp),%rsi
mov    -0x1c(%rbp),%edx
mov    $0x9,%ecx
call   10e0 <cov_mat@plt>
lea    -0x300(%rbp),%rdi
lea    -0x550(%rbp),%rsi
mov    $0x9,%edx
call   1070 <lap_eig@plt>
lea    -0x300(%rbp),%rdi
mov    $0x9,%esi
call   1060 <trnm@plt>
jmp    13b1 <u2fw+0x1f1>
mov    -0x8(%rbp),%rdi
mov    -0x78(%rbp),%rsi
mov    -0x10(%rbp),%rdx
mov    -0x1c(%rbp),%ecx
call   10d0 <lin_fm@plt>
movl   $0x0,-0x55c(%rbp)
mov    -0x55c(%rbp),%eax
cmp    -0x1c(%rbp),%eax
jge    1388 <u2fw+0x1c8>
mov    -0x10(%rbp),%rax
movslq -0x55c(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x560(%rbp)
mov    -0x78(%rbp),%rdi
movslq -0x55c(%rbp),%rax
shl    $0x3,%rax
add    %rax,%rdi
mov    -0x18(%rbp),%rax
movslq -0x560(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    $0x9,%edx
mov    %edx,%esi
call   1080 <scalmul@plt>
mov    -0x55c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x55c(%rbp)
jmp    1323 <u2fw+0x163>
lea    -0x550(%rbp),%rdi
mov    -0x78(%rbp),%rsi
lea    -0x300(%rbp),%rdx
lea    -0x500(%rbp),%r8
mov    $0x9,%ecx
mov    $0x8,%r9d
call   1090 <svduv@plt>
cmpl   $0x8,-0x1c(%rbp)
jle    144c <u2fw+0x28c>
movl   $0x0,-0x560(%rbp)
movl   $0x1,-0x55c(%rbp)
cmpl   $0x9,-0x55c(%rbp)
jge    142b <u2fw+0x26b>
movslq -0x55c(%rbp),%rax
movsd  -0x550(%rbp,%rax,8),%xmm1
ff ff
movslq -0x560(%rbp),%rax
movsd  -0x550(%rbp,%rax,8),%xmm0
ff ff
ucomisd %xmm1,%xmm0
jbe    1412 <u2fw+0x252>
mov    -0x55c(%rbp),%eax
mov    %eax,-0x560(%rbp)
jmp    1417 <u2fw+0x257>
mov    -0x55c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x55c(%rbp)
jmp    13cf <u2fw+0x20f>
lea    -0x300(%rbp),%rax
movslq -0x560(%rbp),%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x558(%rbp)
jmp    145e <u2fw+0x29e>
lea    -0x300(%rbp),%rax
add    $0x40,%rax
mov    %rax,-0x558(%rbp)
movl   $0x0,-0x55c(%rbp)
cmpl   $0x9,-0x55c(%rbp)
jge    14b6 <u2fw+0x2f6>
mov    -0x558(%rbp),%rax
movsd  (%rax),%xmm0
mov    -0x28(%rbp),%rax
movslq -0x55c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x558(%rbp),%rax
add    $0x48,%rax
mov    %rax,-0x558(%rbp)
mov    -0x55c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x55c(%rbp)
jmp    1468 <u2fw+0x2a8>
mov    -0x28(%rbp),%rdi
call   10b0 <singulF@plt>
cmpl   $0x8,-0x1c(%rbp)
jle    14da <u2fw+0x31a>
mov    -0x28(%rbp),%rdi
lea    -0x50(%rbp),%rsi
lea    -0x70(%rbp),%rdx
call   1040 <denormF@plt>
cmpq   $0x0,-0x30(%rbp)
jne    14ee <u2fw+0x32e>
mov    -0x78(%rbp),%rdi
call   1030 <free@plt>
add    $0x560,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<u2fw>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x548,%rsp
mov    %r9,%rbx
mov    %r8,%r13
mov    %ecx,%r14d
mov    %rdx,%r15
mov    %rsi,%r12
mov    %r9,%rbp
test   %r9,%r9
jne    1201 <u2fw+0x41>
lea    0x0(,%r14,8),%eax
mov    %rdi,%rbp
lea    (%rax,%rax,8),%edi
call   10a0 <malloc@plt>
mov    %rbp,%rdi
mov    %rax,%rbp
cmp    $0x9,%r14d
mov    %rbx,0x18(%rsp)
mov    %rbp,0x10(%rsp)
mov    %r14,0x8(%rsp)
jl     12d6 <u2fw+0x116>
lea    0x90(%rsp),%rcx
lea    0x70(%rsp),%r8
mov    %rdi,%rbx
mov    %r12,%rsi
mov    %r14d,%edx
call   10c0 <normu@plt>
mov    %rbx,%rdi
mov    %rbp,%rsi
mov    %r12,%rdx
mov    %r14d,%ecx
lea    0x90(%rsp),%r8
lea    0x70(%rsp),%r9
call   1050 <lin_fmN@plt>
mov    %r14d,%r14d
mov    %rbp,%rbx
xor    %ebp,%ebp
nopl   0x0(%rax,%rax,1)
movslq (%r12,%rbp,4),%rax
movsd  (%r15,%rax,8),%xmm0
mov    %rbx,%rdi
mov    $0x9,%esi
mov    $0x1,%edx
call   1080 <scalmul@plt>
add    $0x1,%rbp
add    $0x48,%rbx
cmp    %rbp,%r14
jne    1260 <u2fw+0xa0>
lea    0xb0(%rsp),%rbx
mov    %rbx,%rdi
mov    0x10(%rsp),%r14
mov    %r14,%rsi
mov    0x8(%rsp),%rbp
mov    %ebp,%edx
mov    $0x9,%ecx
call   10e0 <cov_mat@plt>
lea    0x20(%rsp),%rsi
mov    %rbx,%rdi
mov    $0x9,%edx
call   1070 <lap_eig@plt>
mov    %rbx,%rdi
mov    $0x9,%esi
call   1060 <trnm@plt>
mov    0x18(%rsp),%rbx
jmp    1360 <u2fw+0x1a0>
mov    %rbp,%rsi
mov    %r12,%rdx
mov    %r14d,%ecx
call   10d0 <lin_fm@plt>
test   %r14d,%r14d
jle    1329 <u2fw+0x169>
mov    0x8(%rsp),%r14d
xor    %ebp,%ebp
mov    0x10(%rsp),%rbx
cs nopw 0x0(%rax,%rax,1)
nop
movslq (%r12,%rbp,4),%rax
movsd  (%r15,%rax,8),%xmm0
mov    %rbx,%rdi
mov    $0x9,%esi
mov    $0x9,%edx
call   1080 <scalmul@plt>
add    $0x1,%rbp
add    $0x8,%rbx
cmp    %rbp,%r14
jne    1300 <u2fw+0x140>
lea    0x20(%rsp),%rdi
lea    0xb0(%rsp),%rdx
lea    0x340(%rsp),%r8
mov    0x10(%rsp),%r14
mov    %r14,%rsi
mov    $0x9,%ecx
mov    $0x8,%r9d
call   1090 <svduv@plt>
mov    0x18(%rsp),%rbx
mov    0x8(%rsp),%rbp
mov    $0x8,%eax
cmp    $0x9,%ebp
jl     140c <u2fw+0x24c>
movsd  0x20(%rsp),%xmm0
xor    %eax,%eax
ucomisd 0x28(%rsp),%xmm0
seta   %al
movsd  0x20(%rsp,%rax,8),%xmm0
ucomisd 0x30(%rsp),%xmm0
mov    $0x2,%ecx
cmovbe %eax,%ecx
movsd  0x20(%rsp,%rcx,8),%xmm0
ucomisd 0x38(%rsp),%xmm0
mov    $0x3,%eax
cmovbe %ecx,%eax
movsd  0x20(%rsp,%rax,8),%xmm0
ucomisd 0x40(%rsp),%xmm0
mov    $0x4,%ecx
cmovbe %eax,%ecx
movsd  0x20(%rsp,%rcx,8),%xmm0
ucomisd 0x48(%rsp),%xmm0
mov    $0x5,%eax
cmovbe %ecx,%eax
movsd  0x20(%rsp,%rax,8),%xmm0
ucomisd 0x50(%rsp),%xmm0
mov    $0x6,%ecx
cmovbe %eax,%ecx
movsd  0x20(%rsp,%rcx,8),%xmm0
ucomisd 0x58(%rsp),%xmm0
mov    $0x7,%edx
cmovbe %ecx,%edx
movsd  0x20(%rsp,%rdx,8),%xmm0
ucomisd 0x60(%rsp),%xmm0
mov    $0x8,%eax
cmovbe %rdx,%rax
movsd  0xb0(%rsp,%rax,8),%xmm0
movsd  %xmm0,0x0(%r13)
movsd  0xf8(%rsp,%rax,8),%xmm0
movsd  %xmm0,0x8(%r13)
movsd  0x140(%rsp,%rax,8),%xmm0
movsd  %xmm0,0x10(%r13)
movsd  0x188(%rsp,%rax,8),%xmm0
movsd  %xmm0,0x18(%r13)
movsd  0x1d0(%rsp,%rax,8),%xmm0
movsd  %xmm0,0x20(%r13)
movsd  0x218(%rsp,%rax,8),%xmm0
movsd  %xmm0,0x28(%r13)
movsd  0x260(%rsp,%rax,8),%xmm0
movsd  %xmm0,0x30(%r13)
movsd  0x2a8(%rsp,%rax,8),%xmm0
movsd  %xmm0,0x38(%r13)
movsd  0x2f0(%rsp,%rax,8),%xmm0
movsd  %xmm0,0x40(%r13)
mov    %r13,%rdi
call   10b0 <singulF@plt>
cmp    $0x9,%ebp
jl     14b5 <u2fw+0x2f5>
lea    0x90(%rsp),%rsi
lea    0x70(%rsp),%rdx
mov    %r13,%rdi
call   1040 <denormF@plt>
test   %rbx,%rbx
jne    14c2 <u2fw+0x302>
mov    %r14,%rdi
call   1030 <free@plt>
add    $0x548,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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