#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 */
scalar_t__ NUM_lu(double *, double *, int);
int /*<<< orphan*/ NUM_print_matrix(double *, int);
int /*<<< orphan*/ NUM_print_vec(double *, int);
scalar_t__ malloc(int);
int /*<<< orphan*/ printf(char *);
int test_lu(void) {
  int i;
  double *matrix;
  double *solution;
  int size = 8;
  printf("Testing the LU factorisation algotythm.\n");
  matrix = (double *)malloc(sizeof(double) * size * (size + 1));
  solution = (double *)malloc(sizeof(double) * size);
  matrix[0] = 4.77088e-02;
  matrix[8] = 1.17204e-01;
  matrix[16] = 1.88670e-02;
  matrix[1] = 1.17204e-01;
  matrix[9] = 4.07815e-01;
  matrix[17] = 1.25752e-02;
  matrix[2] = 1.88670e-02;
  matrix[10] = 1.25752e-02;
  matrix[18] = 4.40215e-02;
  matrix[3] = 0.00000e+00;
  matrix[11] = 0.00000e+00;
  matrix[19] = 0.00000e+00;
  matrix[4] = 0.00000e+00;
  matrix[12] = 0.00000e+00;
  matrix[20] = 0.00000e+00;
  matrix[5] = 1.00000e+00;
  matrix[13] = 0.00000e+00;
  matrix[21] = 3.00000e+00;
  matrix[6] = 2.97603e-02;
  matrix[14] = 8.67031e-02;
  matrix[22] = 2.51546e-02;
  matrix[7] = 5.15962e+01;
  matrix[15] = 1.67940e+02;
  matrix[23] = -7.46975e+01;
  matrix[24] = 0.00000e+00;
  matrix[25] = 0.00000e+00;
  matrix[26] = 0.00000e+00;
  matrix[27] = 1.28681e-02;
  matrix[28] = 0.00000e+00;
  matrix[29] = 0.00000e+00;
  matrix[30] = 6.43406e-03;
  matrix[31] = -7.23674e-01;
  matrix[32] = 0.0000e+00;
  matrix[40] = 1.00000e+00;
  matrix[48] = 2.97603e-02;
  matrix[33] = 0.0000e+00;
  matrix[41] = 0.00000e+00;
  matrix[49] = 8.67031e-02;
  matrix[34] = 0.0000e+00;
  matrix[42] = 3.00000e+00;
  matrix[50] = 2.51546e-02;
  matrix[35] = 0.0000e+00;
  matrix[43] = 0.00000e+00;
  matrix[51] = 6.43406e-03;
  matrix[36] = 0.0000e+00;
  matrix[44] = 2.00000e+00;
  matrix[52] = 0.00000e+00;
  matrix[37] = 2.0000e+00;
  matrix[45] = 0.00000e+00;
  matrix[53] = 0.00000e+00;
  matrix[38] = 0.0000e+00;
  matrix[46] = 0.00000e+00;
  matrix[54] = 1.45616e-02;
  matrix[39] = 0.0000e+00;
  matrix[47] = -9.68727e+03;
  matrix[55] = -3.06747e+01;
  matrix[56] = 5.15962e+01;
  matrix[64] = -1.12677e+01;
  matrix[57] = 1.67940e+02;
  matrix[65] = 8.29437e+00;
  matrix[58] = -7.46975e+01;
  matrix[66] = -7.39145e+01;
  matrix[59] = -7.23674e-01;
  matrix[67] = -5.99249e-01;
  matrix[60] = 0.00000e+00;
  matrix[68] = 1.58804e-02;
  matrix[61] = -9.68727e+03;
  matrix[69] = -9.62706e+03;
  matrix[62] = -3.06747e+01;
  matrix[70] = -4.63904e+01;
  matrix[63] = 4.68590e+05;
  matrix[71] = 2.37298e+05;
  NUM_print_matrix(matrix, size);
  if (NUM_lu(matrix, solution, size))
    printf("No solution: Error in the numerical method,\n");
  else
    NUM_print_vec(solution, size);
  printf("\n");
  return 0;
}


x86 assembly (clang, -O0):
<test_lu>:
push   %rbp
mov    %rsp,%rbp
sub    $0x20,%rsp
movl   $0x8,-0x1c(%rbp)
lea    0xf8a(%rip),%rdi
call   1030 <printf@plt>
movslq -0x1c(%rbp),%rax
shl    $0x3,%rax
mov    -0x1c(%rbp),%ecx
add    $0x1,%ecx
movslq %ecx,%rcx
imul   %rcx,%rax
mov    %eax,%edi
call   1060 <malloc@plt>
mov    %rax,-0x10(%rbp)
movslq -0x1c(%rbp),%rax
shl    $0x3,%rax
mov    %eax,%edi
call   1060 <malloc@plt>
mov    %rax,-0x18(%rbp)
mov    -0x10(%rbp),%rax
movsd  0xf3e(%rip),%xmm0
movsd  %xmm0,(%rax)
mov    -0x10(%rbp),%rax
movsd  0xf26(%rip),%xmm0
movsd  %xmm0,0x40(%rax)
mov    -0x10(%rbp),%rax
movsd  0xf05(%rip),%xmm0
movsd  %xmm0,0x80(%rax)
mov    -0x10(%rbp),%rax
movsd  0xf01(%rip),%xmm0
movsd  %xmm0,0x8(%rax)
mov    -0x10(%rbp),%rax
movsd  0xee8(%rip),%xmm0
movsd  %xmm0,0x48(%rax)
mov    -0x10(%rbp),%rax
movsd  0xec7(%rip),%xmm0
movsd  %xmm0,0x88(%rax)
mov    -0x10(%rbp),%rax
movsd  0xebb(%rip),%xmm0
movsd  %xmm0,0x10(%rax)
mov    -0x10(%rbp),%rax
movsd  0xea2(%rip),%xmm0
movsd  %xmm0,0x50(%rax)
mov    -0x10(%rbp),%rax
movsd  0xe89(%rip),%xmm0
movsd  %xmm0,0x90(%rax)
mov    -0x10(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,0x18(%rax)
mov    -0x10(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,0x58(%rax)
mov    -0x10(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,0x98(%rax)
mov    -0x10(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,0x20(%rax)
mov    -0x10(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,0x60(%rax)
mov    -0x10(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,0xa0(%rax)
mov    -0x10(%rbp),%rax
movsd  0xe17(%rip),%xmm0
movsd  %xmm0,0x28(%rax)
mov    -0x10(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,0x68(%rax)
mov    -0x10(%rbp),%rax
movsd  0xde2(%rip),%xmm0
movsd  %xmm0,0xa8(%rax)
mov    -0x10(%rbp),%rax
movsd  0xdde(%rip),%xmm0
movsd  %xmm0,0x30(%rax)
mov    -0x10(%rbp),%rax
movsd  0xdc5(%rip),%xmm0
movsd  %xmm0,0x70(%rax)
mov    -0x10(%rbp),%rax
movsd  0xda4(%rip),%xmm0
movsd  %xmm0,0xb0(%rax)
mov    -0x10(%rbp),%rax
movsd  0xd70(%rip),%xmm0
movsd  %xmm0,0x38(%rax)
mov    -0x10(%rbp),%rax
movsd  0xd4f(%rip),%xmm0
movsd  %xmm0,0x78(%rax)
mov    -0x10(%rbp),%rax
movsd  0xd2e(%rip),%xmm0
movsd  %xmm0,0xb8(%rax)
mov    -0x10(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,0xc0(%rax)
mov    -0x10(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,0xc8(%rax)
mov    -0x10(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,0xd0(%rax)
mov    -0x10(%rbp),%rax
movsd  0xd55(%rip),%xmm0
movsd  %xmm0,0xd8(%rax)
mov    -0x10(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,0xe0(%rax)
mov    -0x10(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,0xe8(%rax)
mov    -0x10(%rbp),%rax
movsd  0xcf3(%rip),%xmm0
movsd  %xmm0,0xf0(%rax)
mov    -0x10(%rbp),%rax
movsd  0xc97(%rip),%xmm0
movsd  %xmm0,0xf8(%rax)
mov    -0x10(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,0x100(%rax)
mov    -0x10(%rbp),%rax
movsd  0xce4(%rip),%xmm0
movsd  %xmm0,0x140(%rax)
mov    -0x10(%rbp),%rax
movsd  0xcc8(%rip),%xmm0
movsd  %xmm0,0x180(%rax)
mov    -0x10(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,0x108(%rax)
mov    -0x10(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,0x148(%rax)
mov    -0x10(%rbp),%rax
movsd  0xc8e(%rip),%xmm0
movsd  %xmm0,0x188(%rax)
mov    -0x10(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,0x110(%rax)
mov    -0x10(%rbp),%rax
movsd  0xc63(%rip),%xmm0
movsd  %xmm0,0x150(%rax)
mov    -0x10(%rbp),%rax
movsd  0xc47(%rip),%xmm0
movsd  %xmm0,0x190(%rax)
mov    -0x10(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,0x118(%rax)
mov    -0x10(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,0x158(%rax)
mov    -0x10(%rbp),%rax
movsd  0xc0d(%rip),%xmm0
movsd  %xmm0,0x198(%rax)
mov    -0x10(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,0x120(%rax)
mov    -0x10(%rbp),%rax
movsd  0xbe2(%rip),%xmm0
movsd  %xmm0,0x160(%rax)
mov    -0x10(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,0x1a0(%rax)
mov    -0x10(%rbp),%rax
movsd  0xbbf(%rip),%xmm0
movsd  %xmm0,0x128(%rax)
mov    -0x10(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,0x168(%rax)
mov    -0x10(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,0x1a8(%rax)
mov    -0x10(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,0x130(%rax)
mov    -0x10(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,0x170(%rax)
mov    -0x10(%rbp),%rax
movsd  0xb67(%rip),%xmm0
movsd  %xmm0,0x1b0(%rax)
mov    -0x10(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,0x138(%rax)
mov    -0x10(%rbp),%rax
movsd  0xaf4(%rip),%xmm0
movsd  %xmm0,0x178(%rax)
mov    -0x10(%rbp),%rax
movsd  0xad0(%rip),%xmm0
movsd  %xmm0,0x1b8(%rax)
mov    -0x10(%rbp),%rax
movsd  0xb14(%rip),%xmm0
movsd  %xmm0,0x1c0(%rax)
mov    -0x10(%rbp),%rax
movsd  0xaf8(%rip),%xmm0
movsd  %xmm0,0x200(%rax)
mov    -0x10(%rbp),%rax
movsd  0xadc(%rip),%xmm0
movsd  %xmm0,0x1c8(%rax)
mov    -0x10(%rbp),%rax
movsd  0xac0(%rip),%xmm0
movsd  %xmm0,0x208(%rax)
mov    -0x10(%rbp),%rax
movsd  0xaa4(%rip),%xmm0
movsd  %xmm0,0x1d0(%rax)
mov    -0x10(%rbp),%rax
movsd  0xa88(%rip),%xmm0
movsd  %xmm0,0x210(%rax)
mov    -0x10(%rbp),%rax
movsd  0xa6c(%rip),%xmm0
movsd  %xmm0,0x1d8(%rax)
mov    -0x10(%rbp),%rax
movsd  0xa50(%rip),%xmm0
movsd  %xmm0,0x218(%rax)
mov    -0x10(%rbp),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,0x1e0(%rax)
mov    -0x10(%rbp),%rax
movsd  0xa25(%rip),%xmm0
movsd  %xmm0,0x220(%rax)
mov    -0x10(%rbp),%rax
movsd  0xa09(%rip),%xmm0
movsd  %xmm0,0x1e8(%rax)
mov    -0x10(%rbp),%rax
movsd  0x9ed(%rip),%xmm0
movsd  %xmm0,0x228(%rax)
mov    -0x10(%rbp),%rax
movsd  0x9d1(%rip),%xmm0
movsd  %xmm0,0x1f0(%rax)
mov    -0x10(%rbp),%rax
movsd  0x9b5(%rip),%xmm0
movsd  %xmm0,0x230(%rax)
mov    -0x10(%rbp),%rax
movsd  0x999(%rip),%xmm0
movsd  %xmm0,0x1f8(%rax)
mov    -0x10(%rbp),%rax
movsd  0x97d(%rip),%xmm0
movsd  %xmm0,0x238(%rax)
mov    -0x10(%rbp),%rdi
mov    -0x1c(%rbp),%esi
call   1040 <NUM_print_matrix@plt>
mov    -0x10(%rbp),%rdi
mov    -0x18(%rbp),%rsi
mov    -0x1c(%rbp),%edx
call   1070 <NUM_lu@plt>
cmp    $0x0,%rax
je     16c2 <test_lu+0x572>
lea    0xa61(%rip),%rdi
call   1030 <printf@plt>
jmp    16ce <test_lu+0x57e>
mov    -0x18(%rbp),%rdi
mov    -0x1c(%rbp),%esi
call   1050 <NUM_print_vec@plt>
lea    0xa6f(%rip),%rdi
call   1030 <printf@plt>
xor    %eax,%eax
add    $0x20,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<test_lu>:
push   %r14
push   %rbx
push   %rax
lea    0xfbd(%rip),%rdi
call   1040 <puts@plt>
mov    $0x240,%edi
call   1070 <malloc@plt>
mov    %rax,%rbx
mov    $0x40,%edi
call   1070 <malloc@plt>
mov    %rax,%r14
movaps 0xe6f(%rip),%xmm0
movups %xmm0,(%rbx)
movabs $0x3f9351deefe4ffc9,%rax
51 93 3f
mov    %rax,0x10(%rbx)
movaps 0xe67(%rip),%xmm0
movups %xmm0,0x48(%rbx)
movaps 0xe6c(%rip),%xmm0
movups %xmm0,0x88(%rbx)
xorps  %xmm0,%xmm0
movups %xmm0,0x18(%rbx)
movups %xmm0,0x58(%rbx)
movups %xmm0,0x98(%rbx)
movaps 0xe5c(%rip),%xmm1
movups %xmm1,0x28(%rbx)
xorps  %xmm1,%xmm1
movhps 0xf1e(%rip),%xmm1
movups %xmm1,0x68(%rbx)
movaps 0xe53(%rip),%xmm1
movups %xmm1,0xa8(%rbx)
movaps 0xe55(%rip),%xmm1
movups %xmm1,0x38(%rbx)
movaps 0xe5a(%rip),%xmm1
movups %xmm1,0x78(%rbx)
movabs $0xc052aca3d70a3d71,%rax
ac 52 c0
mov    %rax,0xb8(%rbx)
movups %xmm0,0xc0(%rbx)
movq   $0x0,0xd0(%rbx)
movabs $0x3f8a5a97254b1324,%rax
5a 8a 3f
mov    %rax,0xd8(%rbx)
movups %xmm0,0xe0(%rbx)
movaps 0xe24(%rip),%xmm1
movups %xmm1,0xf0(%rbx)
movabs $0x3ff0000000000000,%rax
mov    %rax,0x140(%rbx)
movabs $0x3f9e797bece0688b,%rax
79 9e 3f
mov    %rax,0x180(%rbx)
movups %xmm0,0x100(%rbx)
movq   $0x0,0x148(%rbx)
movabs $0x4008000000000000,%rax
mov    %rax,0x150(%rbx)
movaps 0xde1(%rip),%xmm1
movups %xmm1,0x188(%rbx)
movups %xmm0,0x110(%rbx)
movq   $0x0,0x158(%rbx)
movabs $0x4000000000000000,%rax
mov    %rax,0x160(%rbx)
movsd  0xe37(%rip),%xmm1
xorps  %xmm2,%xmm2
movhps 0xe35(%rip),%xmm2
movups %xmm1,0x198(%rbx)
movups %xmm2,0x120(%rbx)
movups %xmm0,0x168(%rbx)
xorps  %xmm1,%xmm1
movhps 0xe1e(%rip),%xmm1
movups %xmm1,0x1a8(%rbx)
movups %xmm0,0x130(%rbx)
movabs $0xc0c2eba28f5c28f6,%rax
eb c2 c0
mov    %rax,0x178(%rbx)
movaps 0xd70(%rip),%xmm0
movups %xmm0,0x1b8(%rbx)
movaps 0xd72(%rip),%xmm0
movups %xmm0,0x1c8(%rbx)
movaps 0xd74(%rip),%xmm0
movups %xmm0,0x208(%rbx)
movsd  0xdd5(%rip),%xmm0
movups %xmm0,0x1d8(%rbx)
movaps 0xd67(%rip),%xmm0
movups %xmm0,0x218(%rbx)
movaps 0xd69(%rip),%xmm0
movups %xmm0,0x1e8(%rbx)
movaps 0xd6b(%rip),%xmm0
movups %xmm0,0x228(%rbx)
movaps 0xd6d(%rip),%xmm0
movups %xmm0,0x1f8(%rbx)
movabs $0x410cf79000000000,%rax
f7 0c 41
mov    %rax,0x238(%rbx)
mov    %rbx,%rdi
mov    $0x8,%esi
call   1050 <NUM_print_matrix@plt>
mov    %rbx,%rdi
mov    %r14,%rsi
mov    $0x8,%edx
call   1080 <NUM_lu@plt>
test   %rax,%rax
je     13cb <test_lu+0x26b>
lea    0xd8c(%rip),%rdi
call   1040 <puts@plt>
jmp    13d8 <test_lu+0x278>
mov    %r14,%rdi
mov    $0x8,%esi
call   1060 <NUM_print_vec@plt>
mov    $0xa,%edi
call   1030 <putchar@plt>
xor    %eax,%eax
add    $0x8,%rsp
pop    %rbx
pop    %r14
ret

LLVM IR (-O0):
@.str = [41 x i8] c"Testing the LU factorisation algotythm.\0A\00"
@.str.1 = [45 x i8] c"No solution: Error in the numerical method,\0A\00"
@.str.2 = [2 x i8] c"\0A\00"
define i32 @test_lu() {
%1 = alloca i32
%2 = alloca double*
%3 = alloca double*
%4 = alloca i32
store i32 8, i32* %4
%5 = call i32 @printf(i8* getelementptr inbounds ([41 x i8], [41 x i8]* @.str, i64 0, i64 0))
%6 = load i32, i32* %4
%7 = sext i32 %6 to i64
%8 = mul i64 8, %7
%9 = load i32, i32* %4
%10 = add nsw i32 %9, 1
%11 = sext i32 %10 to i64
%12 = mul i64 %8, %11
%13 = trunc i64 %12 to i32
%14 = call i64 @malloc(i32 %13)
%15 = inttoptr i64 %14 to double*
store double* %15, double** %2
%16 = load i32, i32* %4
%17 = sext i32 %16 to i64
%18 = mul i64 8, %17
%19 = trunc i64 %18 to i32
%20 = call i64 @malloc(i32 %19)
%21 = inttoptr i64 %20 to double*
store double* %21, double** %3
%22 = load double*, double** %2
%23 = getelementptr inbounds double, double* %22, i64 0
store double 4.770880e-02, double* %23
%24 = load double*, double** %2
%25 = getelementptr inbounds double, double* %24, i64 8
store double 1.172040e-01, double* %25
%26 = load double*, double** %2
%27 = getelementptr inbounds double, double* %26, i64 16
store double 1.886700e-02, double* %27
%28 = load double*, double** %2
%29 = getelementptr inbounds double, double* %28, i64 1
store double 1.172040e-01, double* %29
%30 = load double*, double** %2
%31 = getelementptr inbounds double, double* %30, i64 9
store double 4.078150e-01, double* %31
%32 = load double*, double** %2
%33 = getelementptr inbounds double, double* %32, i64 17
store double 1.257520e-02, double* %33
%34 = load double*, double** %2
%35 = getelementptr inbounds double, double* %34, i64 2
store double 1.886700e-02, double* %35
%36 = load double*, double** %2
%37 = getelementptr inbounds double, double* %36, i64 10
store double 1.257520e-02, double* %37
%38 = load double*, double** %2
%39 = getelementptr inbounds double, double* %38, i64 18
store double 4.402150e-02, double* %39
%40 = load double*, double** %2
%41 = getelementptr inbounds double, double* %40, i64 3
store double 0.000000e+00, double* %41
%42 = load double*, double** %2
%43 = getelementptr inbounds double, double* %42, i64 11
store double 0.000000e+00, double* %43
%44 = load double*, double** %2
%45 = getelementptr inbounds double, double* %44, i64 19
store double 0.000000e+00, double* %45
%46 = load double*, double** %2
%47 = getelementptr inbounds double, double* %46, i64 4
store double 0.000000e+00, double* %47
%48 = load double*, double** %2
%49 = getelementptr inbounds double, double* %48, i64 12
store double 0.000000e+00, double* %49
%50 = load double*, double** %2
%51 = getelementptr inbounds double, double* %50, i64 20
store double 0.000000e+00, double* %51
%52 = load double*, double** %2
%53 = getelementptr inbounds double, double* %52, i64 5
store double 1.000000e+00, double* %53
%54 = load double*, double** %2
%55 = getelementptr inbounds double, double* %54, i64 13
store double 0.000000e+00, double* %55
%56 = load double*, double** %2
%57 = getelementptr inbounds double, double* %56, i64 21
store double 3.000000e+00, double* %57
%58 = load double*, double** %2
%59 = getelementptr inbounds double, double* %58, i64 6
store double 2.976030e-02, double* %59
%60 = load double*, double** %2
%61 = getelementptr inbounds double, double* %60, i64 14
store double 8.670310e-02, double* %61
%62 = load double*, double** %2
%63 = getelementptr inbounds double, double* %62, i64 22
store double 2.515460e-02, double* %63
%64 = load double*, double** %2
%65 = getelementptr inbounds double, double* %64, i64 7
store double 5.159620e+01, double* %65
%66 = load double*, double** %2
%67 = getelementptr inbounds double, double* %66, i64 15
store double 1.679400e+02, double* %67
%68 = load double*, double** %2
%69 = getelementptr inbounds double, double* %68, i64 23
store double -7.469750e+01, double* %69
%70 = load double*, double** %2
%71 = getelementptr inbounds double, double* %70, i64 24
store double 0.000000e+00, double* %71
%72 = load double*, double** %2
%73 = getelementptr inbounds double, double* %72, i64 25
store double 0.000000e+00, double* %73
%74 = load double*, double** %2
%75 = getelementptr inbounds double, double* %74, i64 26
store double 0.000000e+00, double* %75
%76 = load double*, double** %2
%77 = getelementptr inbounds double, double* %76, i64 27
store double 1.286810e-02, double* %77
%78 = load double*, double** %2
%79 = getelementptr inbounds double, double* %78, i64 28
store double 0.000000e+00, double* %79
%80 = load double*, double** %2
%81 = getelementptr inbounds double, double* %80, i64 29
store double 0.000000e+00, double* %81
%82 = load double*, double** %2
%83 = getelementptr inbounds double, double* %82, i64 30
store double 6.434060e-03, double* %83
%84 = load double*, double** %2
%85 = getelementptr inbounds double, double* %84, i64 31
store double -7.236740e-01, double* %85
%86 = load double*, double** %2
%87 = getelementptr inbounds double, double* %86, i64 32
store double 0.000000e+00, double* %87
%88 = load double*, double** %2
%89 = getelementptr inbounds double, double* %88, i64 40
store double 1.000000e+00, double* %89
%90 = load double*, double** %2
%91 = getelementptr inbounds double, double* %90, i64 48
store double 2.976030e-02, double* %91
%92 = load double*, double** %2
%93 = getelementptr inbounds double, double* %92, i64 33
store double 0.000000e+00, double* %93
%94 = load double*, double** %2
%95 = getelementptr inbounds double, double* %94, i64 41
store double 0.000000e+00, double* %95
%96 = load double*, double** %2
%97 = getelementptr inbounds double, double* %96, i64 49
store double 8.670310e-02, double* %97
%98 = load double*, double** %2
%99 = getelementptr inbounds double, double* %98, i64 34
store double 0.000000e+00, double* %99
%100 = load double*, double** %2
%101 = getelementptr inbounds double, double* %100, i64 42
store double 3.000000e+00, double* %101
%102 = load double*, double** %2
%103 = getelementptr inbounds double, double* %102, i64 50
store double 2.515460e-02, double* %103
%104 = load double*, double** %2
%105 = getelementptr inbounds double, double* %104, i64 35
store double 0.000000e+00, double* %105
%106 = load double*, double** %2
%107 = getelementptr inbounds double, double* %106, i64 43
store double 0.000000e+00, double* %107
%108 = load double*, double** %2
%109 = getelementptr inbounds double, double* %108, i64 51
store double 6.434060e-03, double* %109
%110 = load double*, double** %2
%111 = getelementptr inbounds double, double* %110, i64 36
store double 0.000000e+00, double* %111
%112 = load double*, double** %2
%113 = getelementptr inbounds double, double* %112, i64 44
store double 2.000000e+00, double* %113
%114 = load double*, double** %2
%115 = getelementptr inbounds double, double* %114, i64 52
store double 0.000000e+00, double* %115
%116 = load double*, double** %2
%117 = getelementptr inbounds double, double* %116, i64 37
store double 2.000000e+00, double* %117
%118 = load double*, double** %2
%119 = getelementptr inbounds double, double* %118, i64 45
store double 0.000000e+00, double* %119
%120 = load double*, double** %2
%121 = getelementptr inbounds double, double* %120, i64 53
store double 0.000000e+00, double* %121
%122 = load double*, double** %2
%123 = getelementptr inbounds double, double* %122, i64 38
store double 0.000000e+00, double* %123
%124 = load double*, double** %2
%125 = getelementptr inbounds double, double* %124, i64 46
store double 0.000000e+00, double* %125
%126 = load double*, double** %2
%127 = getelementptr inbounds double, double* %126, i64 54
store double 1.456160e-02, double* %127
%128 = load double*, double** %2
%129 = getelementptr inbounds double, double* %128, i64 39
store double 0.000000e+00, double* %129
%130 = load double*, double** %2
%131 = getelementptr inbounds double, double* %130, i64 47
store double -9.687270e+03, double* %131
%132 = load double*, double** %2
%133 = getelementptr inbounds double, double* %132, i64 55
store double -3.067470e+01, double* %133
%134 = load double*, double** %2
%135 = getelementptr inbounds double, double* %134, i64 56
store double 5.159620e+01, double* %135
%136 = load double*, double** %2
%137 = getelementptr inbounds double, double* %136, i64 64
store double -1.126770e+01, double* %137
%138 = load double*, double** %2
%139 = getelementptr inbounds double, double* %138, i64 57
store double 1.679400e+02, double* %139
%140 = load double*, double** %2
%141 = getelementptr inbounds double, double* %140, i64 65
store double 8.294370e+00, double* %141
%142 = load double*, double** %2
%143 = getelementptr inbounds double, double* %142, i64 58
store double -7.469750e+01, double* %143
%144 = load double*, double** %2
%145 = getelementptr inbounds double, double* %144, i64 66
store double -7.391450e+01, double* %145
%146 = load double*, double** %2
%147 = getelementptr inbounds double, double* %146, i64 59
store double -7.236740e-01, double* %147
%148 = load double*, double** %2
%149 = getelementptr inbounds double, double* %148, i64 67
store double -5.992490e-01, double* %149
%150 = load double*, double** %2
%151 = getelementptr inbounds double, double* %150, i64 60
store double 0.000000e+00, double* %151
%152 = load double*, double** %2
%153 = getelementptr inbounds double, double* %152, i64 68
store double 1.588040e-02, double* %153
%154 = load double*, double** %2
%155 = getelementptr inbounds double, double* %154, i64 61
store double -9.687270e+03, double* %155
%156 = load double*, double** %2
%157 = getelementptr inbounds double, double* %156, i64 69
store double 0xC0C2CD87AE147AE1, double* %157
%158 = load double*, double** %2
%159 = getelementptr inbounds double, double* %158, i64 62
store double -3.067470e+01, double* %159
%160 = load double*, double** %2
%161 = getelementptr inbounds double, double* %160, i64 70
store double -4.639040e+01, double* %161
%162 = load double*, double** %2
%163 = getelementptr inbounds double, double* %162, i64 63
store double 4.685900e+05, double* %163
%164 = load double*, double** %2
%165 = getelementptr inbounds double, double* %164, i64 71
store double 2.372980e+05, double* %165
%166 = load double*, double** %2
%167 = load i32, i32* %4
%168 = call i32 @NUM_print_matrix(double* %166, i32 %167)
%169 = load double*, double** %2
%170 = load double*, double** %3
%171 = load i32, i32* %4
%172 = call i64 @NUM_lu(double* %169, double* %170, i32 %171)
%173 = icmp ne i64 %172, 0
br i1 %173, label %174, label %176
174: ; preds = %0
%175 = call i32 @printf(i8* getelementptr inbounds ([45 x i8], [45 x i8]* @.str.1, i64 0, i64 0))
br label %180
176: ; preds = %0
%177 = load double*, double** %3
%178 = load i32, i32* %4
%179 = call i32 @NUM_print_vec(double* %177, i32 %178)
br label %180
180: ; preds = %176, %174
%181 = call i32 @printf(i8* getelementptr inbounds ([2 x i8], [2 x i8]* @.str.2, i64 0, i64 0))
ret i32 0
}

LLVM IR (-O2):
@str = [40 x i8] c"Testing the LU factorisation algotythm.\00"
@str.3 = [44 x i8] c"No solution: Error in the numerical method,\00"
define i32 @test_lu() {
%1 = tail call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([40 x i8], [40 x i8]* @str, i64 0, i64 0))
%2 = tail call i64 @malloc(i32 576)
%3 = inttoptr i64 %2 to double*
%4 = tail call i64 @malloc(i32 64)
%5 = inttoptr i64 %4 to double*
%6 = inttoptr i64 %2 to <2 x double>*
store <2 x double> <double 4.770880e-02, double 1.172040e-01>, <2 x double>* %6
%7 = getelementptr inbounds double, double* %3, i64 9
%8 = getelementptr inbounds double, double* %3, i64 17
%9 = getelementptr inbounds double, double* %3, i64 2
store double 1.886700e-02, double* %9
%10 = bitcast double* %7 to <2 x double>*
store <2 x double> <double 4.078150e-01, double 1.257520e-02>, <2 x double>* %10
%11 = bitcast double* %8 to <2 x double>*
store <2 x double> <double 1.257520e-02, double 4.402150e-02>, <2 x double>* %11
%12 = getelementptr inbounds double, double* %3, i64 3
%13 = getelementptr inbounds double, double* %3, i64 11
%14 = getelementptr inbounds double, double* %3, i64 19
%15 = getelementptr inbounds double, double* %3, i64 5
%16 = bitcast double* %12 to i8*
%17 = bitcast double* %13 to i8*
%18 = bitcast double* %14 to i8*
%19 = getelementptr inbounds double, double* %3, i64 13
%20 = getelementptr inbounds double, double* %3, i64 21
%21 = bitcast double* %15 to <2 x double>*
store <2 x double> <double 1.000000e+00, double 2.976030e-02>, <2 x double>* %21
%22 = bitcast double* %19 to <2 x double>*
store <2 x double> <double 0.000000e+00, double 8.670310e-02>, <2 x double>* %22
%23 = bitcast double* %20 to <2 x double>*
store <2 x double> <double 3.000000e+00, double 2.515460e-02>, <2 x double>* %23
%24 = getelementptr inbounds double, double* %3, i64 7
%25 = bitcast double* %24 to <2 x double>*
store <2 x double> <double 5.159620e+01, double 1.172040e-01>, <2 x double>* %25
%26 = getelementptr inbounds double, double* %3, i64 15
%27 = bitcast double* %26 to <2 x double>*
store <2 x double> <double 1.679400e+02, double 1.886700e-02>, <2 x double>* %27
%28 = getelementptr inbounds double, double* %3, i64 23
store double -7.469750e+01, double* %28
%29 = getelementptr inbounds double, double* %3, i64 24
%30 = getelementptr inbounds double, double* %3, i64 27
%31 = bitcast double* %29 to i8*
store double 1.286810e-02, double* %30
%32 = getelementptr inbounds double, double* %3, i64 28
%33 = getelementptr inbounds double, double* %3, i64 30
%34 = bitcast double* %32 to i8*
%35 = bitcast double* %33 to <2 x double>*
store <2 x double> <double 6.434060e-03, double -7.236740e-01>, <2 x double>* %35
%36 = getelementptr inbounds double, double* %3, i64 32
%37 = getelementptr inbounds double, double* %3, i64 40
store double 1.000000e+00, double* %37
%38 = getelementptr inbounds double, double* %3, i64 48
store double 2.976030e-02, double* %38
%39 = bitcast double* %36 to <2 x double>*
store <2 x double> zeroinitializer, <2 x double>* %39
%40 = getelementptr inbounds double, double* %3, i64 41
store double 0.000000e+00, double* %40
%41 = getelementptr inbounds double, double* %3, i64 49
%42 = getelementptr inbounds double, double* %3, i64 34
%43 = getelementptr inbounds double, double* %3, i64 42
store double 3.000000e+00, double* %43
%44 = bitcast double* %41 to <2 x double>*
store <2 x double> <double 8.670310e-02, double 2.515460e-02>, <2 x double>* %44
%45 = bitcast double* %42 to <2 x double>*
store <2 x double> zeroinitializer, <2 x double>* %45
%46 = getelementptr inbounds double, double* %3, i64 43
store double 0.000000e+00, double* %46
%47 = getelementptr inbounds double, double* %3, i64 51
%48 = getelementptr inbounds double, double* %3, i64 36
%49 = getelementptr inbounds double, double* %3, i64 44
store double 2.000000e+00, double* %49
%50 = bitcast double* %47 to <2 x double>*
store <2 x double> <double 6.434060e-03, double 0.000000e+00>, <2 x double>* %50
%51 = bitcast double* %48 to <2 x double>*
store <2 x double> <double 0.000000e+00, double 2.000000e+00>, <2 x double>* %51
%52 = getelementptr inbounds double, double* %3, i64 45
%53 = getelementptr inbounds double, double* %3, i64 53
%54 = getelementptr inbounds double, double* %3, i64 38
%55 = bitcast double* %52 to i8*
%56 = bitcast double* %53 to <2 x double>*
store <2 x double> <double 0.000000e+00, double 1.456160e-02>, <2 x double>* %56
%57 = bitcast double* %54 to <2 x double>*
store <2 x double> zeroinitializer, <2 x double>* %57
%58 = getelementptr inbounds double, double* %3, i64 47
store double -9.687270e+03, double* %58
%59 = getelementptr inbounds double, double* %3, i64 55
%60 = bitcast double* %59 to <2 x double>*
store <2 x double> <double -3.067470e+01, double 5.159620e+01>, <2 x double>* %60
%61 = getelementptr inbounds double, double* %3, i64 57
%62 = getelementptr inbounds double, double* %3, i64 65
%63 = bitcast double* %61 to <2 x double>*
store <2 x double> <double 1.679400e+02, double -7.469750e+01>, <2 x double>* %63
%64 = bitcast double* %62 to <2 x double>*
store <2 x double> <double 8.294370e+00, double -7.391450e+01>, <2 x double>* %64
%65 = getelementptr inbounds double, double* %3, i64 59
%66 = getelementptr inbounds double, double* %3, i64 67
%67 = bitcast double* %65 to <2 x double>*
store <2 x double> <double -7.236740e-01, double 0.000000e+00>, <2 x double>* %67
%68 = bitcast double* %66 to <2 x double>*
store <2 x double> <double -5.992490e-01, double 1.588040e-02>, <2 x double>* %68
%69 = getelementptr inbounds double, double* %3, i64 61
%70 = getelementptr inbounds double, double* %3, i64 69
%71 = bitcast double* %69 to <2 x double>*
store <2 x double> <double -9.687270e+03, double -3.067470e+01>, <2 x double>* %71
%72 = bitcast double* %70 to <2 x double>*
store <2 x double> <double 0xC0C2CD87AE147AE1, double -4.639040e+01>, <2 x double>* %72
%73 = getelementptr inbounds double, double* %3, i64 63
%74 = bitcast double* %73 to <2 x double>*
store <2 x double> <double 4.685900e+05, double -1.126770e+01>, <2 x double>* %74
%75 = getelementptr inbounds double, double* %3, i64 71
store double 2.372980e+05, double* %75
%76 = tail call i32 @NUM_print_matrix(double* %3, i32 8)
%77 = tail call i64 @NUM_lu(double* %3, double* %5, i32 8)
%78 = icmp eq i64 %77, 0
br i1 %78, label %81, label %79
79: ; preds = %0
%80 = tail call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([44 x i8], [44 x i8]* @str.3, i64 0, i64 0))
br label %83
81: ; preds = %0
%82 = tail call i32 @NUM_print_vec(double* %5, i32 8)
br label %83
83: ; preds = %81, %79
%84 = tail call i32 @putchar(i32 10)
ret i32 0
}