#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 EXIT_SUCCESS;
int /*<<< orphan*/ resetMartixDoubleZero(double **, int, int);
int LUFactotisation(int n, double **A, double **U, double **L) {
  int i, j, k;
  double m;
  resetMartixDoubleZero(U, n, n);
  resetMartixDoubleZero(L, n, n);
  for (i = 0; i < n; i++) {
    U[0][i] = A[0][i];
    L[i][i] = 1.0;
  }
  for (k = 0; k < n - 1; k++) {
    for (i = k + 1; i < n; i++) {
      m = A[i][k] / A[k][k];
      L[i][k] = m;
      for (j = k + 1; j < n; j++)
        U[i][j] = A[i][j] = A[i][j] - m * A[k][j];
    }
  }
  for (k = 0; k < n; k++)
    for (i = 0; i < n; i++)
      if (k < i)
        U[i][k] = 0.0;
  return (EXIT_SUCCESS);
}


x86 assembly (clang, -O0):
<LUFactotisation>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %rcx,-0x20(%rbp)
mov    -0x18(%rbp),%rdi
mov    -0x4(%rbp),%esi
mov    -0x4(%rbp),%edx
call   1030 <resetMartixDoubleZero@plt>
mov    -0x20(%rbp),%rdi
mov    -0x4(%rbp),%esi
mov    -0x4(%rbp),%edx
call   1030 <resetMartixDoubleZero@plt>
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    11a3 <LUFactotisation+0x93>
mov    -0x10(%rbp),%rax
mov    (%rax),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    (%rax),%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x20(%rbp),%rax
movslq -0x24(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x24(%rbp),%rcx
movsd  0xe70(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    114c <LUFactotisation+0x3c>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
mov    -0x4(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    12da <LUFactotisation+0x1ca>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    12c7 <LUFactotisation+0x1b7>
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
divsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x38(%rbp),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x24(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    12b4 <LUFactotisation+0x1a4>
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x38(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x28(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm2
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x24(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x28(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x28(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1222 <LUFactotisation+0x112>
jmp    12b9 <LUFactotisation+0x1a9>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    11c4 <LUFactotisation+0xb4>
jmp    12cc <LUFactotisation+0x1bc>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    11aa <LUFactotisation+0x9a>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    134a <LUFactotisation+0x23a>
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    1337 <LUFactotisation+0x227>
mov    -0x2c(%rbp),%eax
cmp    -0x24(%rbp),%eax
jge    1324 <LUFactotisation+0x214>
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x2c(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    1329 <LUFactotisation+0x219>
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    12f4 <LUFactotisation+0x1e4>
jmp    133c <LUFactotisation+0x22c>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    12e1 <LUFactotisation+0x1d1>
mov    0x2c8f(%rip),%rax
mov    (%rax),%eax
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<LUFactotisation>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x98,%rsp
mov    %rcx,%r15
mov    %rdx,%rbx
mov    %rsi,%rbp
mov    %edi,%r14d
mov    %rdx,%rdi
mov    %r14d,%esi
mov    %r14d,%edx
call   1030 <resetMartixDoubleZero@plt>
mov    %r15,%rdi
mov    %r14d,%esi
mov    %r14d,%edx
call   1030 <resetMartixDoubleZero@plt>
test   %r14d,%r14d
jle    152e <LUFactotisation+0x41e>
mov    %rbx,%r8
mov    %r15,0x20(%rsp)
mov    0x0(%rbp),%rcx
mov    (%rbx),%rdx
mov    %r14d,%eax
mov    %rax,%r9
movabs $0x3ff0000000000000,%r10
add    $0xffffffffffffffff,%r9
je     11c3 <LUFactotisation+0xb3>
mov    %eax,%r11d
and    $0xfffffffe,%r11d
xor    %edi,%edi
mov    0x20(%rsp),%rsi
cs nopw 0x0(%rax,%rax,1)
nop
movsd  (%rcx,%rdi,8),%xmm0
movsd  %xmm0,(%rdx,%rdi,8)
mov    (%rsi,%rdi,8),%rbx
mov    %r10,(%rbx,%rdi,8)
movsd  0x8(%rcx,%rdi,8),%xmm0
movsd  %xmm0,0x8(%rdx,%rdi,8)
mov    0x8(%rsi,%rdi,8),%rbx
mov    %r10,0x8(%rbx,%rdi,8)
add    $0x2,%rdi
cmp    %rdi,%r11
jne    1190 <LUFactotisation+0x80>
jmp    11c5 <LUFactotisation+0xb5>
xor    %edi,%edi
mov    %r9,0x30(%rsp)
test   $0x1,%al
je     11e5 <LUFactotisation+0xd5>
movsd  (%rcx,%rdi,8),%xmm0
movsd  %xmm0,(%rdx,%rdi,8)
mov    0x20(%rsp),%rcx
mov    (%rcx,%rdi,8),%rcx
mov    %r10,(%rcx,%rdi,8)
mov    %r14,0x28(%rsp)
cmp    $0x2,%r14d
jl     14b2 <LUFactotisation+0x3a2>
mov    0x28(%rsp),%rcx
add    $0xffffffff,%ecx
mov    %rcx,0x50(%rsp)
mov    %rax,%rcx
neg    %rcx
mov    %rcx,0x88(%rsp)
mov    $0x1,%ecx
mov    %rcx,0x10(%rsp)
mov    $0x3,%ecx
mov    %rcx,0x18(%rsp)
xor    %ecx,%ecx
movapd 0xdd3(%rip),%xmm0
mov    0x30(%rsp),%rdx
mov    %rdx,0x38(%rsp)
mov    %r8,0x70(%rsp)
mov    %rbp,0x90(%rsp)
jmp    1272 <LUFactotisation+0x162>
cs nopw 0x0(%rax,%rax,1)
addq   $0x1,0x10(%rsp)
addq   $0xffffffffffffffff,0x38(%rsp)
addq   $0x1,0x18(%rsp)
mov    0x58(%rsp),%rcx
cmp    0x50(%rsp),%rcx
je     14b2 <LUFactotisation+0x3a2>
mov    %rcx,0x40(%rsp)
add    $0x1,%rcx
mov    %rcx,0x58(%rsp)
cmp    %rax,%rcx
jae    1250 <LUFactotisation+0x140>
mov    0x38(%rsp),%rsi
and    $0xfffffffffffffffc,%rsi
mov    0x40(%rsp),%rcx
mov    %rcx,%rdx
not    %rdx
add    %rax,%rdx
mov    0x0(%rbp,%rcx,8),%r15
mov    0x10(%rsp),%rdi
lea    (%r15,%rdi,8),%rcx
mov    %rcx,0x80(%rsp)
lea    (%r15,%rax,8),%rcx
mov    %rcx,0x78(%rsp)
mov    %rdx,0x48(%rsp)
and    $0xfffffffffffffffc,%rdx
mov    %rdx,0x68(%rsp)
lea    (%rdi,%rdx,1),%rcx
mov    %rcx,0x60(%rsp)
mov    0x18(%rsp),%rcx
lea    (%r15,%rcx,8),%r11
jmp    12f5 <LUFactotisation+0x1e5>
nopl   (%rax)
add    $0x1,%rdi
cmp    %rax,%rdi
mov    0x90(%rsp),%rbp
je     1250 <LUFactotisation+0x140>
mov    0x0(%rbp,%rdi,8),%r12
mov    0x40(%rsp),%rdx
movsd  (%r12,%rdx,8),%xmm1
divsd  (%r15,%rdx,8),%xmm1
mov    0x20(%rsp),%rcx
mov    (%rcx,%rdi,8),%rcx
movsd  %xmm1,(%rcx,%rdx,8)
xorpd  %xmm0,%xmm1
mov    (%r8,%rdi,8),%r14
mov    0x10(%rsp),%r13
cmpq   $0x4,0x48(%rsp)
jae    13c0 <LUFactotisation+0x2b0>
mov    %eax,%ecx
sub    %r13d,%ecx
mov    %r13,%r10
test   $0x1,%cl
je     135f <LUFactotisation+0x24f>
movsd  (%r15,%r13,8),%xmm2
mulsd  %xmm1,%xmm2
addsd  (%r12,%r13,8),%xmm2
movsd  %xmm2,(%r12,%r13,8)
movsd  %xmm2,(%r14,%r13,8)
lea    0x1(%r13),%r10
not    %r13
cmp    0x88(%rsp),%r13
je     12e0 <LUFactotisation+0x1d0>
movsd  (%r15,%r10,8),%xmm2
mulsd  %xmm1,%xmm2
addsd  (%r12,%r10,8),%xmm2
movsd  %xmm2,(%r12,%r10,8)
movsd  %xmm2,(%r14,%r10,8)
movsd  0x8(%r15,%r10,8),%xmm2
mulsd  %xmm1,%xmm2
addsd  0x8(%r12,%r10,8),%xmm2
movsd  %xmm2,0x8(%r12,%r10,8)
movsd  %xmm2,0x8(%r14,%r10,8)
add    $0x2,%r10
cmp    %r10,%rax
jne    1370 <LUFactotisation+0x260>
jmp    12e0 <LUFactotisation+0x1d0>
nopw   0x0(%rax,%rax,1)
mov    0x10(%rsp),%r13
lea    (%r12,%r13,8),%rcx
lea    (%r12,%rax,8),%rdx
lea    (%r14,%r13,8),%rbp
lea    (%r14,%rax,8),%rbx
cmp    %rbx,%rcx
setb   0xf(%rsp)
cmp    %rdx,%rbp
setb   %r10b
mov    0x78(%rsp),%r8
cmp    %r8,%rcx
setb   %r9b
mov    0x80(%rsp),%rcx
cmp    %rdx,%rcx
setb   0xe(%rsp)
cmp    %r8,%rbp
setb   %dl
cmp    %rbx,%rcx
setb   %bl
test   %r10b,0xf(%rsp)
mov    0x70(%rsp),%r8
jne    1332 <LUFactotisation+0x222>
mov    0x10(%rsp),%r13
and    0xe(%rsp),%r9b
jne    1332 <LUFactotisation+0x222>
mov    0x10(%rsp),%r13
and    %bl,%dl
jne    1332 <LUFactotisation+0x222>
movapd %xmm1,%xmm2
unpcklpd %xmm1,%xmm2
mov    0x18(%rsp),%rcx
lea    (%r12,%rcx,8),%rbx
lea    (%r14,%rcx,8),%rbp
xor    %edx,%edx
movupd -0x10(%rbx,%rdx,8),%xmm3
movupd (%rbx,%rdx,8),%xmm4
movupd -0x10(%r11,%rdx,8),%xmm5
movupd (%r11,%rdx,8),%xmm6
mulpd  %xmm2,%xmm5
addpd  %xmm3,%xmm5
mulpd  %xmm2,%xmm6
addpd  %xmm4,%xmm6
movupd %xmm5,-0x10(%rbx,%rdx,8)
movupd %xmm6,(%rbx,%rdx,8)
movupd %xmm5,-0x10(%rbp,%rdx,8)
movupd %xmm6,0x0(%rbp,%rdx,8)
add    $0x4,%rdx
cmp    %rdx,%rsi
jne    1450 <LUFactotisation+0x340>
mov    0x60(%rsp),%r13
mov    0x68(%rsp),%rcx
cmp    %rcx,0x48(%rsp)
je     12e0 <LUFactotisation+0x1d0>
jmp    1332 <LUFactotisation+0x222>
cmpl   $0x0,0x28(%rsp)
mov    0x30(%rsp),%rbp
jle    152e <LUFactotisation+0x41e>
mov    %eax,%ecx
and    $0xfffffffe,%ecx
xor    %edx,%edx
jmp    14d9 <LUFactotisation+0x3c9>
nopw   0x0(%rax,%rax,1)
add    $0x1,%rdx
cmp    %rax,%rdx
je     152e <LUFactotisation+0x41e>
xor    %esi,%esi
test   %rbp,%rbp
jne    1509 <LUFactotisation+0x3f9>
test   $0x1,%al
je     14d0 <LUFactotisation+0x3c0>
cmp    %rsi,%rdx
jae    14d0 <LUFactotisation+0x3c0>
mov    (%r8,%rsi,8),%rsi
movq   $0x0,(%rsi,%rdx,8)
jmp    14d0 <LUFactotisation+0x3c0>
nopw   0x0(%rax,%rax,1)
add    $0x2,%rsi
cmp    %rsi,%rcx
je     14e0 <LUFactotisation+0x3d0>
cmp    %rsi,%rdx
jae    151a <LUFactotisation+0x40a>
mov    (%r8,%rsi,8),%rdi
movq   $0x0,(%rdi,%rdx,8)
cmp    %rsi,%rdx
ja     1500 <LUFactotisation+0x3f0>
mov    0x8(%r8,%rsi,8),%rdi
movq   $0x0,(%rdi,%rdx,8)
jmp    1500 <LUFactotisation+0x3f0>
mov    0x2aab(%rip),%rax
mov    (%rax),%eax
add    $0x98,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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