#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*/ free_matrix(float **, int, int, int, int);
float **matrix(int, int, int, int);
void orthog(int n, float anu[], float alpha[], float beta[], float a[],
            float b[]) {
  int k, l;
  float **sig;
  int looptmp;
  sig = matrix(1, 2 * n + 1, 1, 2 * n + 1);
  looptmp = 2 * n;
  for (l = 3; l <= looptmp; l++)
    sig[1][l] = 0.0;
  looptmp++;
  for (l = 2; l <= looptmp; l++)
    sig[2][l] = anu[l - 1];
  a[1] = alpha[1] + anu[2] / anu[1];
  b[1] = 0.0;
  for (k = 3; k <= n + 1; k++) {
    looptmp = 2 * n - k + 3;
    for (l = k; l <= looptmp; l++) {
      sig[k][l] = sig[k - 1][l + 1] +
                  (alpha[l - 1] - a[k - 2]) * sig[k - 1][l] -
                  b[k - 2] * sig[k - 2][l] + beta[l - 1] * sig[k - 1][l - 1];
    }
    a[k - 1] = alpha[k - 1] + sig[k][k + 1] / sig[k][k] -
               sig[k - 1][k] / sig[k - 1][k - 1];
    b[k - 1] = sig[k][k] / sig[k - 1][k - 1];
  }
  free_matrix(sig, 1, 2 * n + 1, 1, 2 * n + 1);
}


x86 assembly (clang, -O0):
<orthog>:
push   %rbp
mov    %rsp,%rbp
sub    $0x50,%rsp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9,-0x30(%rbp)
mov    -0x4(%rbp),%esi
shl    $0x1,%esi
add    $0x1,%esi
mov    -0x4(%rbp),%ecx
shl    $0x1,%ecx
add    $0x1,%ecx
mov    $0x1,%edx
mov    %edx,%edi
call   1040 <matrix@plt>
mov    %rax,-0x40(%rbp)
mov    -0x4(%rbp),%eax
shl    $0x1,%eax
mov    %eax,-0x44(%rbp)
movl   $0x3,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x44(%rbp),%eax
jg     119f <orthog+0x7f>
mov    -0x40(%rbp),%rax
mov    0x8(%rax),%rax
movslq -0x38(%rbp),%rcx
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    1171 <orthog+0x51>
mov    -0x44(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x44(%rbp)
movl   $0x2,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x44(%rbp),%eax
jg     11ec <orthog+0xcc>
mov    -0x10(%rbp),%rax
mov    -0x38(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x40(%rbp),%rax
mov    0x10(%rax),%rax
movslq -0x38(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    11af <orthog+0x8f>
mov    -0x18(%rbp),%rax
movss  0x4(%rax),%xmm0
mov    -0x10(%rbp),%rax
movss  0x8(%rax),%xmm1
mov    -0x10(%rbp),%rax
divss  0x4(%rax),%xmm1
addss  %xmm1,%xmm0
mov    -0x28(%rbp),%rax
movss  %xmm0,0x4(%rax)
mov    -0x30(%rbp),%rax
xorps  %xmm0,%xmm0
movss  %xmm0,0x4(%rax)
movl   $0x3,-0x34(%rbp)
mov    -0x34(%rbp),%eax
mov    -0x4(%rbp),%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jg     1443 <orthog+0x323>
mov    -0x4(%rbp),%eax
shl    $0x1,%eax
sub    -0x34(%rbp),%eax
add    $0x3,%eax
mov    %eax,-0x44(%rbp)
mov    -0x34(%rbp),%eax
mov    %eax,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x44(%rbp),%eax
jg     135b <orthog+0x23b>
mov    -0x40(%rbp),%rax
mov    -0x34(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x38(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm1
mov    -0x18(%rbp),%rax
mov    -0x38(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x28(%rbp),%rax
mov    -0x34(%rbp),%ecx
sub    $0x2,%ecx
movslq %ecx,%rcx
subss  (%rax,%rcx,4),%xmm0
mov    -0x40(%rbp),%rax
mov    -0x34(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x30(%rbp),%rax
mov    -0x34(%rbp),%ecx
sub    $0x2,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm1
mov    -0x40(%rbp),%rax
mov    -0x34(%rbp),%ecx
sub    $0x2,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x38(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
movd   %xmm1,%eax
xor    $0x80000000,%eax
movd   %eax,%xmm1
mulss  %xmm2,%xmm1
addss  %xmm0,%xmm1
mov    -0x20(%rbp),%rax
mov    -0x38(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x40(%rbp),%rax
mov    -0x34(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x38(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x40(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x38(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    124d <orthog+0x12d>
mov    -0x18(%rbp),%rax
mov    -0x34(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x40(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x34(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm1
mov    -0x40(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x34(%rbp),%rcx
divss  (%rax,%rcx,4),%xmm1
addss  %xmm1,%xmm0
mov    -0x40(%rbp),%rax
mov    -0x34(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x34(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
mov    -0x40(%rbp),%rax
mov    -0x34(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x34(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
divss  (%rax,%rcx,4),%xmm1
subss  %xmm1,%xmm0
mov    -0x28(%rbp),%rax
mov    -0x34(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x40(%rbp),%rax
movslq -0x34(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x34(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x40(%rbp),%rax
mov    -0x34(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
mov    (%rax,%rcx,8),%rax
mov    -0x34(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
divss  (%rax,%rcx,4),%xmm0
mov    -0x30(%rbp),%rax
mov    -0x34(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1227 <orthog+0x107>
mov    -0x40(%rbp),%rdi
mov    -0x4(%rbp),%edx
shl    $0x1,%edx
add    $0x1,%edx
mov    -0x4(%rbp),%r8d
shl    $0x1,%r8d
add    $0x1,%r8d
mov    $0x1,%ecx
mov    %ecx,%esi
call   1030 <free_matrix@plt>
add    $0x50,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<orthog>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x18,%rsp
mov    %r9,%r14
mov    %r8,%r12
mov    %rcx,%r15
mov    %rdx,%r13
mov    %rsi,%rbp
mov    %edi,%ebx
lea    (%rbx,%rbx,1),%eax
mov    %rax,0x8(%rsp)
lea    0x1(%rbx,%rbx,1),%ecx
mov    $0x1,%edi
mov    %ecx,%esi
mov    $0x1,%edx
mov    %ecx,0x14(%rsp)
call   1050 <matrix@plt>
mov    %rbx,%r11
mov    %rax,(%rsp)
cmp    $0x2,%r11d
jl     11ad <orthog+0x7d>
mov    (%rsp),%rax
mov    0x8(%rax),%rdi
add    $0xc,%rdi
mov    0x8(%rsp),%rax
cmp    $0x4,%eax
mov    $0x3,%edx
cmovge %eax,%edx
add    $0xfffffffe,%edx
shl    $0x2,%rdx
xor    %esi,%esi
mov    %r11,%rbx
call   1040 <memset@plt>
mov    %rbx,%r11
mov    0x14(%rsp),%edx
cmp    $0x2,%edx
jl     1308 <orthog+0x1d8>
mov    (%rsp),%rax
mov    0x10(%rax),%rax
mov    0x8(%rsp),%rcx
add    $0x2,%ecx
lea    -0x2(%rcx),%r8
mov    $0x2,%esi
cmp    $0x8,%r8
jb     1293 <orthog+0x163>
lea    0x8(%rax),%rbx
lea    -0x4(,%rcx,4),%rdi
ff
add    %rbp,%rdi
cmp    %rdi,%rbx
jae    1202 <orthog+0xd2>
lea    (%rax,%rcx,4),%rbx
lea    0x4(%rbp),%rdi
cmp    %rbx,%rdi
jb     1293 <orthog+0x163>
mov    %r8,%rsi
and    $0xfffffffffffffff8,%rsi
lea    -0x8(%rsi),%rdi
mov    %rdi,%r9
shr    $0x3,%r9
add    $0x1,%r9
test   %rdi,%rdi
je     1462 <orthog+0x332>
mov    %r9,%rbx
and    $0xfffffffffffffffe,%rbx
xor    %edi,%edi
nopw   0x0(%rax,%rax,1)
movups 0x4(%rbp,%rdi,4),%xmm0
movups 0x14(%rbp,%rdi,4),%xmm1
movups %xmm0,0x8(%rax,%rdi,4)
movups %xmm1,0x18(%rax,%rdi,4)
movups 0x24(%rbp,%rdi,4),%xmm0
movups 0x34(%rbp,%rdi,4),%xmm1
movups %xmm0,0x28(%rax,%rdi,4)
movups %xmm1,0x38(%rax,%rdi,4)
add    $0x10,%rdi
add    $0xfffffffffffffffe,%rbx
jne    1230 <orthog+0x100>
test   $0x1,%r9b
je     128a <orthog+0x15a>
shl    $0x2,%rdi
mov    %rdi,%rbx
or     $0x4,%rbx
movups 0x0(%rbp,%rbx,1),%xmm0
movups 0x10(%rbp,%rbx,1),%xmm1
or     $0x8,%rdi
movups %xmm0,(%rax,%rdi,1)
movups %xmm1,0x10(%rax,%rdi,1)
cmp    %rsi,%r8
je     1308 <orthog+0x1d8>
or     $0x2,%rsi
mov    %ecx,%edi
sub    %esi,%edi
mov    %rsi,%rbx
not    %rbx
add    %rcx,%rbx
and    $0x3,%rdi
je     12c5 <orthog+0x195>
cs nopw 0x0(%rax,%rax,1)
movss  -0x4(%rbp,%rsi,4),%xmm0
movss  %xmm0,(%rax,%rsi,4)
add    $0x1,%rsi
add    $0xffffffffffffffff,%rdi
jne    12b0 <orthog+0x180>
cmp    $0x3,%rbx
jb     1308 <orthog+0x1d8>
nopl   0x0(%rax,%rax,1)
movss  -0x4(%rbp,%rsi,4),%xmm0
movss  %xmm0,(%rax,%rsi,4)
movss  0x0(%rbp,%rsi,4),%xmm0
movss  %xmm0,0x4(%rax,%rsi,4)
movss  0x4(%rbp,%rsi,4),%xmm0
movss  %xmm0,0x8(%rax,%rsi,4)
movss  0x8(%rbp,%rsi,4),%xmm0
movss  %xmm0,0xc(%rax,%rsi,4)
add    $0x4,%rsi
cmp    %rsi,%rcx
jne    12d0 <orthog+0x1a0>
movss  0x8(%rbp),%xmm0
divss  0x4(%rbp),%xmm0
addss  0x4(%r13),%xmm0
movss  %xmm0,0x4(%r12)
movl   $0x0,0x4(%r14)
cmp    $0x2,%r11d
jge    1351 <orthog+0x221>
mov    (%rsp),%rdi
mov    $0x1,%esi
mov    $0x1,%ecx
mov    %edx,%r8d
add    $0x18,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
jmp    1030 <free_matrix@plt>
mov    0x8(%rsp),%rax
add    $0x3,%eax
mov    %rax,0x8(%rsp)
add    $0x2,%r11d
mov    $0x3,%ecx
xorps  %xmm0,%xmm0
mov    $0x4,%r8d
mov    %edx,%r10d
jmp    13e2 <orthog+0x2b2>
cs nopw 0x0(%rax,%rax,1)
nop
mov    (%rsp),%rax
mov    -0x8(%rax,%rcx,8),%rdi
mov    (%rax,%rcx,8),%rbx
lea    -0x1(%rcx),%r9
lea    0x1(%rcx),%rax
mov    %eax,%esi
movss  (%rbx,%rsi,4),%xmm0
divss  (%rbx,%rcx,4),%xmm0
addss  0x0(%r13,%r9,4),%xmm0
movss  (%rdi,%rcx,4),%xmm1
divss  (%rdi,%r9,4),%xmm1
subss  %xmm1,%xmm0
movss  %xmm0,(%r12,%r9,4)
movss  (%rbx,%rcx,4),%xmm0
divss  (%rdi,%r9,4),%xmm0
movss  %xmm0,(%r14,%r9,4)
add    $0xffffffff,%r10d
add    $0x1,%r8
mov    %rax,%rcx
cmp    %r11,%rax
je     132d <orthog+0x1fd>
mov    %r10d,%r10d
mov    0x8(%rsp),%rax
sub    %ecx,%eax
cltq
cmp    %rax,%rcx
jg     1380 <orthog+0x250>
lea    -0x1(%rcx),%r9
mov    (%rsp),%rsi
mov    -0x10(%rsi,%rcx,8),%rax
mov    -0x8(%rsi,%rcx,8),%rdi
mov    (%rsi,%rcx,8),%rbx
mov    %r8,%rbp
nopl   0x0(%rax)
mov    %ebp,%esi
movss  -0x8(%r13,%rbp,4),%xmm1
subss  -0x8(%r12,%rcx,4),%xmm1
mulss  -0x4(%rdi,%rbp,4),%xmm1
addss  (%rdi,%rsi,4),%xmm1
mulss  -0x4(%rax,%rbp,4),%xmm0
movss  -0x8(%r15,%rbp,4),%xmm2
mulss  -0x8(%rdi,%rbp,4),%xmm2
subss  %xmm0,%xmm1
addss  %xmm1,%xmm2
movss  %xmm2,-0x4(%rbx,%rbp,4)
cmp    %rbp,%r10
je     1391 <orthog+0x261>
movss  -0x8(%r14,%rcx,4),%xmm0
add    $0x1,%rbp
jmp    1410 <orthog+0x2e0>
xor    %edi,%edi
test   $0x1,%r9b
jne    1268 <orthog+0x138>
jmp    128a <orthog+0x15a>

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

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