#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 */
typedef int /*<<< orphan*/ FILE;

/* Variables and functions */
int /*<<< orphan*/ fclose(int /*<<< orphan*/ *);
int /*<<< orphan*/ feof(int /*<<< orphan*/ *);
int /*<<< orphan*/ fgets(char *, int, int /*<<< orphan*/ *);
int /*<<< orphan*/ *fopen(char *, char *);
int /*<<< orphan*/ free_ivector(int *, int, int);
int /*<<< orphan*/ free_matrix(float **, int, int, int, int);
int /*<<< orphan*/ free_vector(float *, int, int);
int /*<<< orphan*/ fscanf(int /*<<< orphan*/ *, char *, ...);
int /*<<< orphan*/ getchar();
int *ivector(int, int);
int /*<<< orphan*/ lubksb(float **, int, int *, float *);
int /*<<< orphan*/ ludcmp(float **, int, int *, float *);
float **matrix(int, int, int, int);
int /*<<< orphan*/ nrerror(char *);
int /*<<< orphan*/ printf(char *, ...);
float *vector(int, int);
int main(void) {
  int j, k, l, m, n, *indx;
  float p, *x, **a, **b, **c;
  char dummy[80];
  FILE *fp;
  indx = ivector(1, 20);
  x = vector(1, 20);
  a = matrix(1, 20, 1, 20);
  b = matrix(1, 20, 1, 20);
  c = matrix(1, 20, 1, 20);
  if ((fp = fopen("matrx1.dat", "r")) == NULL)
    nrerror("Data file matrx1.dat not found\n");
  while (!feof(fp)) {
    fgets(dummy, 80, fp);
    fgets(dummy, 80, fp);
    fscanf(fp, "%d %d ", &n, &m);
    fgets(dummy, 80, fp);
    for (k = 1; k <= n; k++)
      for (l = 1; l <= n; l++)
        fscanf(fp, "%f ", &a[k][l]);
    fgets(dummy, 80, fp);
    for (l = 1; l <= m; l++)
      for (k = 1; k <= n; k++)
        fscanf(fp, "%f ", &b[k][l]);
    for (l = 1; l <= n; l++)
      for (k = 1; k <= n; k++)
        c[k][l] = a[k][l];
    ludcmp(c, n, indx, &p);
    for (k = 1; k <= m; k++) {
      for (l = 1; l <= n; l++)
        x[l] = b[l][k];
      lubksb(c, n, indx, x);
      printf("right-hand side vector:\n");
      for (l = 1; l <= n; l++)
        printf("%12.6f", b[l][k]);
      printf("\n%s%s\n", "result of matrix applied", " to sol'n vector");
      for (l = 1; l <= n; l++) {
        b[l][k] = 0.0;
        for (j = 1; j <= n; j++)
          b[l][k] += (a[l][j] * x[j]);
      }
      for (l = 1; l <= n; l++)
        printf("%12.6f", b[l][k]);
      printf("\n*********************************\n");
    }
    printf("press RETURN for next problem:\n");
    (void)getchar();
  }
  fclose(fp);
  free_matrix(c, 1, 20, 1, 20);
  free_matrix(b, 1, 20, 1, 20);
  free_matrix(a, 1, 20, 1, 20);
  free_vector(x, 1, 20);
  free_ivector(indx, 1, 20);
  return 0;
}


x86 assembly (clang, -O0):
<main>:
push   %rbp
mov    %rsp,%rbp
sub    $0xb0,%rsp
movl   $0x0,-0x4(%rbp)
mov    $0x1,%edi
mov    $0x14,%esi
call   1110 <ivector@plt>
mov    %rax,-0x20(%rbp)
mov    $0x1,%edi
mov    $0x14,%esi
call   10f0 <vector@plt>
mov    %rax,-0x30(%rbp)
mov    $0x1,%edx
mov    $0x14,%ecx
mov    %edx,%edi
mov    %ecx,%esi
call   10e0 <matrix@plt>
mov    %rax,-0x38(%rbp)
mov    $0x1,%edx
mov    $0x14,%ecx
mov    %edx,%edi
mov    %ecx,%esi
call   10e0 <matrix@plt>
mov    %rax,-0x40(%rbp)
mov    $0x1,%edx
mov    $0x14,%ecx
mov    %edx,%edi
mov    %ecx,%esi
call   10e0 <matrix@plt>
mov    %rax,-0x48(%rbp)
lea    0xd7c(%rip),%rdi
lea    0xdfa(%rip),%rsi
call   1120 <fopen@plt>
mov    %rax,-0xa8(%rbp)
cmp    $0x0,%rax
jne    12ad <main+0xad>
lea    0xd63(%rip),%rdi
call   1050 <nrerror@plt>
jmp    12b2 <main+0xb2>
mov    -0xa8(%rbp),%rdi
call   10c0 <feof@plt>
cmp    $0x0,%eax
setne  %al
xor    $0xff,%al
test   $0x1,%al
jne    12d3 <main+0xd3>
jmp    16b8 <main+0x4b8>
lea    -0xa0(%rbp),%rdi
mov    -0xa8(%rbp),%rdx
mov    $0x50,%esi
call   10a0 <fgets@plt>
lea    -0xa0(%rbp),%rdi
mov    -0xa8(%rbp),%rdx
mov    $0x50,%esi
call   10a0 <fgets@plt>
mov    -0xa8(%rbp),%rdi
lea    0xd1a(%rip),%rsi
lea    -0x18(%rbp),%rdx
lea    -0x14(%rbp),%rcx
mov    $0x0,%al
call   1090 <fscanf@plt>
lea    -0xa0(%rbp),%rdi
mov    -0xa8(%rbp),%rdx
mov    $0x50,%esi
call   10a0 <fgets@plt>
movl   $0x1,-0xc(%rbp)
mov    -0xc(%rbp),%eax
cmp    -0x18(%rbp),%eax
jg     13ab <main+0x1ab>
movl   $0x1,-0x10(%rbp)
mov    -0x10(%rbp),%eax
cmp    -0x18(%rbp),%eax
jg     1398 <main+0x198>
mov    -0xa8(%rbp),%rdi
mov    -0x38(%rbp),%rax
movslq -0xc(%rbp),%rcx
mov    (%rax,%rcx,8),%rdx
movslq -0x10(%rbp),%rax
shl    $0x2,%rax
add    %rax,%rdx
lea    0xcaf(%rip),%rsi
mov    $0x0,%al
call   1090 <fscanf@plt>
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    1352 <main+0x152>
jmp    139d <main+0x19d>
mov    -0xc(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xc(%rbp)
jmp    133f <main+0x13f>
lea    -0xa0(%rbp),%rdi
mov    -0xa8(%rbp),%rdx
mov    $0x50,%esi
call   10a0 <fgets@plt>
movl   $0x1,-0x10(%rbp)
mov    -0x10(%rbp),%eax
cmp    -0x14(%rbp),%eax
jg     1436 <main+0x236>
movl   $0x1,-0xc(%rbp)
mov    -0xc(%rbp),%eax
cmp    -0x18(%rbp),%eax
jg     1423 <main+0x223>
mov    -0xa8(%rbp),%rdi
mov    -0x40(%rbp),%rax
movslq -0xc(%rbp),%rcx
mov    (%rax,%rcx,8),%rdx
movslq -0x10(%rbp),%rax
shl    $0x2,%rax
add    %rax,%rdx
lea    0xc24(%rip),%rsi
mov    $0x0,%al
call   1090 <fscanf@plt>
mov    -0xc(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xc(%rbp)
jmp    13dd <main+0x1dd>
jmp    1428 <main+0x228>
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    13ca <main+0x1ca>
movl   $0x1,-0x10(%rbp)
mov    -0x10(%rbp),%eax
cmp    -0x18(%rbp),%eax
jg     14a7 <main+0x2a7>
movl   $0x1,-0xc(%rbp)
mov    -0xc(%rbp),%eax
cmp    -0x18(%rbp),%eax
jg     1494 <main+0x294>
mov    -0x38(%rbp),%rax
movslq -0xc(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x10(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x48(%rbp),%rax
movslq -0xc(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x10(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0xc(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xc(%rbp)
jmp    1450 <main+0x250>
jmp    1499 <main+0x299>
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    143d <main+0x23d>
mov    -0x48(%rbp),%rdi
mov    -0x18(%rbp),%esi
mov    -0x20(%rbp),%rdx
lea    -0x24(%rbp),%rcx
call   1030 <ludcmp@plt>
movl   $0x1,-0xc(%rbp)
mov    -0xc(%rbp),%eax
cmp    -0x14(%rbp),%eax
jg     169e <main+0x49e>
movl   $0x1,-0x10(%rbp)
mov    -0x10(%rbp),%eax
cmp    -0x18(%rbp),%eax
jg     1511 <main+0x311>
mov    -0x40(%rbp),%rax
movslq -0x10(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0xc(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x30(%rbp),%rax
movslq -0x10(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    14d5 <main+0x2d5>
mov    -0x48(%rbp),%rdi
mov    -0x18(%rbp),%esi
mov    -0x20(%rbp),%rdx
mov    -0x30(%rbp),%rcx
call   10d0 <lubksb@plt>
lea    0xb0a(%rip),%rdi
mov    $0x0,%al
call   1070 <printf@plt>
movl   $0x1,-0x10(%rbp)
mov    -0x10(%rbp),%eax
cmp    -0x18(%rbp),%eax
jg     157b <main+0x37b>
mov    -0x40(%rbp),%rax
movslq -0x10(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0xc(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
lea    0xae9(%rip),%rdi
mov    $0x1,%al
call   1070 <printf@plt>
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    153a <main+0x33a>
lea    0xad4(%rip),%rdi
lea    0xad4(%rip),%rsi
lea    0xae6(%rip),%rdx
mov    $0x0,%al
call   1070 <printf@plt>
movl   $0x1,-0x10(%rbp)
mov    -0x10(%rbp),%eax
cmp    -0x18(%rbp),%eax
jg     163a <main+0x43a>
mov    -0x40(%rbp),%rax
movslq -0x10(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0xc(%rbp),%rcx
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax,%rcx,4)
movl   $0x1,-0x8(%rbp)
mov    -0x8(%rbp),%eax
cmp    -0x18(%rbp),%eax
jg     1627 <main+0x427>
mov    -0x38(%rbp),%rax
movslq -0x10(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x8(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x30(%rbp),%rax
movslq -0x8(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
mov    -0x40(%rbp),%rax
movslq -0x10(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0xc(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x8(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x8(%rbp)
jmp    15c9 <main+0x3c9>
jmp    162c <main+0x42c>
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    159e <main+0x39e>
movl   $0x1,-0x10(%rbp)
mov    -0x10(%rbp),%eax
cmp    -0x18(%rbp),%eax
jg     1682 <main+0x482>
mov    -0x40(%rbp),%rax
movslq -0x10(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0xc(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
lea    0x9e2(%rip),%rdi
mov    $0x1,%al
call   1070 <printf@plt>
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    1641 <main+0x441>
lea    0x9fe(%rip),%rdi
mov    $0x0,%al
call   1070 <printf@plt>
mov    -0xc(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xc(%rbp)
jmp    14c2 <main+0x2c2>
lea    0xa06(%rip),%rdi
mov    $0x0,%al
call   1070 <printf@plt>
mov    $0x0,%al
call   10b0 <getchar@plt>
jmp    12b2 <main+0xb2>
mov    -0xa8(%rbp),%rdi
call   1060 <fclose@plt>
mov    -0x48(%rbp),%rdi
mov    $0x1,%ecx
mov    $0x14,%r8d
mov    %ecx,%esi
mov    %r8d,%edx
call   1040 <free_matrix@plt>
mov    -0x40(%rbp),%rdi
mov    $0x1,%ecx
mov    $0x14,%r8d
mov    %ecx,%esi
mov    %r8d,%edx
call   1040 <free_matrix@plt>
mov    -0x38(%rbp),%rdi
mov    $0x1,%ecx
mov    $0x14,%r8d
mov    %ecx,%esi
mov    %r8d,%edx
call   1040 <free_matrix@plt>
mov    -0x30(%rbp),%rdi
mov    $0x1,%esi
mov    $0x14,%edx
call   1080 <free_vector@plt>
mov    -0x20(%rbp),%rdi
mov    $0x1,%esi
mov    $0x14,%edx
call   1100 <free_ivector@plt>
xor    %eax,%eax
add    $0xb0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<main>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x98,%rsp
mov    $0x1,%edi
mov    $0x14,%esi
call   1120 <ivector@plt>
mov    %rax,0x18(%rsp)
mov    $0x1,%edi
mov    $0x14,%esi
call   1100 <vector@plt>
mov    %rax,%rbx
mov    $0x1,%edi
mov    $0x14,%esi
mov    $0x1,%edx
mov    $0x14,%ecx
call   10f0 <matrix@plt>
mov    %rax,%r15
mov    $0x1,%edi
mov    $0x14,%esi
mov    $0x1,%edx
mov    $0x14,%ecx
call   10f0 <matrix@plt>
mov    %rax,%r12
mov    $0x1,%edi
mov    $0x14,%esi
mov    $0x1,%edx
mov    $0x14,%ecx
call   10f0 <matrix@plt>
mov    %rax,%r13
lea    0xd5e(%rip),%rdi
lea    0xdc3(%rip),%rsi
call   1130 <fopen@plt>
mov    %rax,%rbp
test   %rax,%rax
jne    12c2 <main+0xb2>
lea    0xd4e(%rip),%rdi
call   1060 <nrerror@plt>
mov    %rbp,%rdi
call   10d0 <feof@plt>
test   %eax,%eax
je     136b <main+0x15b>
mov    %rbp,%rdi
call   1070 <fclose@plt>
mov    %r13,%rdi
mov    $0x1,%esi
mov    $0x14,%edx
mov    $0x1,%ecx
mov    $0x14,%r8d
call   1050 <free_matrix@plt>
mov    %r12,%rdi
mov    $0x1,%esi
mov    $0x14,%edx
mov    $0x1,%ecx
mov    $0x14,%r8d
call   1050 <free_matrix@plt>
mov    %r15,%rdi
mov    $0x1,%esi
mov    $0x14,%edx
mov    $0x1,%ecx
mov    $0x14,%r8d
call   1050 <free_matrix@plt>
mov    %rbx,%rdi
mov    $0x1,%esi
mov    $0x14,%edx
call   1090 <free_vector@plt>
mov    0x18(%rsp),%rdi
mov    $0x1,%esi
mov    $0x14,%edx
call   1110 <free_ivector@plt>
xor    %eax,%eax
add    $0x98,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
mov    %r13,0x20(%rsp)
mov    %rbp,0x10(%rsp)
mov    %r15,0x30(%rsp)
mov    %r12,0x28(%rsp)
jmp    13b8 <main+0x1a8>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
lea    0xcd7(%rip),%rdi
call   1040 <puts@plt>
xor    %eax,%eax
call   10c0 <getchar@plt>
mov    0x10(%rsp),%rbp
mov    %rbp,%rdi
call   10d0 <feof@plt>
test   %eax,%eax
jne    12d2 <main+0xc2>
lea    0x40(%rsp),%r14
mov    %r14,%rdi
mov    $0x50,%esi
mov    %rbp,%rdx
call   10b0 <fgets@plt>
mov    %r14,%rdi
mov    $0x50,%esi
mov    %rbp,%rdx
call   10b0 <fgets@plt>
mov    %rbp,%rdi
lea    0xc44(%rip),%rsi
lea    0x8(%rsp),%rdx
lea    0xc(%rsp),%rcx
xor    %eax,%eax
call   10a0 <fscanf@plt>
mov    %r14,%rdi
mov    $0x50,%esi
mov    %rbp,%rdx
call   10b0 <fgets@plt>
mov    0x8(%rsp),%eax
test   %eax,%eax
jle    1490 <main+0x280>
mov    $0x1,%r14d
jmp    142f <main+0x21f>
nopl   0x0(%rax)
movslq %eax,%rcx
lea    0x1(%r14),%rdx
cmp    %rcx,%r14
mov    %rdx,%r14
jge    1490 <main+0x280>
test   %eax,%eax
jle    1420 <main+0x210>
mov    $0x4,%ebp
xor    %r13d,%r13d
mov    %r15,%r12
mov    0x10(%rsp),%r15
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
mov    (%r12,%r14,8),%rdx
add    %rbp,%rdx
mov    %r15,%rdi
lea    0xbd1(%rip),%rsi
xor    %eax,%eax
call   10a0 <fscanf@plt>
movslq 0x8(%rsp),%rax
add    $0x1,%r13
add    $0x4,%rbp
cmp    %rax,%r13
jl     1450 <main+0x240>
mov    0x20(%rsp),%r13
mov    %r12,%r15
mov    0x28(%rsp),%r12
jmp    1420 <main+0x210>
nopl   0x0(%rax)
lea    0x40(%rsp),%rdi
mov    $0x50,%esi
mov    0x10(%rsp),%rdx
call   10b0 <fgets@plt>
mov    0xc(%rsp),%eax
mov    0x8(%rsp),%esi
test   %eax,%eax
jle    1530 <main+0x320>
mov    $0x1,%r14d
jmp    14cf <main+0x2bf>
nopl   0x0(%rax)
lea    0x1(%r14),%rcx
movslq %eax,%rdx
cmp    %rdx,%r14
mov    %rcx,%r14
jge    1530 <main+0x320>
test   %esi,%esi
jle    14c0 <main+0x2b0>
xor    %ebp,%ebp
mov    0x10(%rsp),%r15
lea    0xb51(%rip),%r13
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
lea    0x0(,%r14,4),%rdx
add    0x8(%r12,%rbp,8),%rdx
mov    %r15,%rdi
mov    %r13,%rsi
xor    %eax,%eax
call   10a0 <fscanf@plt>
movslq 0x8(%rsp),%rsi
add    $0x1,%rbp
cmp    %rsi,%rbp
jl     14f0 <main+0x2e0>
mov    0xc(%rsp),%eax
mov    0x30(%rsp),%r15
mov    0x20(%rsp),%r13
jmp    14c0 <main+0x2b0>
nopl   0x0(%rax,%rax,1)
test   %esi,%esi
jle    1630 <main+0x420>
mov    %esi,%r10d
lea    0x1(%r10),%r8
lea    -0x1(%r10),%r9
mov    %r10d,%r11d
and    $0x3,%r11d
and    $0xfffffffc,%r10d
mov    $0x1,%ecx
jmp    156d <main+0x35d>
cs nopw 0x0(%rax,%rax,1)
nop
add    $0x1,%rcx
cmp    %r8,%rcx
je     1630 <main+0x420>
mov    $0x1,%ebp
cmp    $0x3,%r9
jb     15dd <main+0x3cd>
xor    %ebp,%ebp
nopw   0x0(%rax,%rax,1)
mov    0x8(%r15,%rbp,8),%rax
movss  (%rax,%rcx,4),%xmm0
mov    0x8(%r13,%rbp,8),%rax
movss  %xmm0,(%rax,%rcx,4)
mov    0x10(%r15,%rbp,8),%rax
movss  (%rax,%rcx,4),%xmm0
mov    0x10(%r13,%rbp,8),%rax
movss  %xmm0,(%rax,%rcx,4)
mov    0x18(%r15,%rbp,8),%rax
movss  (%rax,%rcx,4),%xmm0
mov    0x18(%r13,%rbp,8),%rax
movss  %xmm0,(%rax,%rcx,4)
mov    0x20(%r15,%rbp,8),%rax
movss  (%rax,%rcx,4),%xmm0
mov    0x20(%r13,%rbp,8),%rax
movss  %xmm0,(%rax,%rcx,4)
add    $0x4,%rbp
cmp    %rbp,%r10
jne    1580 <main+0x370>
add    $0x1,%rbp
test   %r11,%r11
je     1560 <main+0x350>
lea    0x0(,%rbp,8),%rdx
add    %r13,%rdx
lea    (%r15,%rbp,8),%rbp
xor    %eax,%eax
nopw   0x0(%rax,%rax,1)
mov    0x0(%rbp,%rax,8),%rdi
movss  (%rdi,%rcx,4),%xmm0
mov    (%rdx,%rax,8),%rdi
movss  %xmm0,(%rdi,%rcx,4)
add    $0x1,%rax
cmp    %rax,%r11
jne    1600 <main+0x3f0>
jmp    1560 <main+0x350>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    %r13,%rdi
mov    0x18(%rsp),%rdx
lea    0x3c(%rsp),%rcx
call   1030 <ludcmp@plt>
cmpl   $0x0,0xc(%rsp)
jle    1390 <main+0x180>
mov    $0x1,%r14d
jmp    1686 <main+0x476>
cs nopw 0x0(%rax,%rax,1)
nop
lea    0xa3e(%rip),%rdi
call   1040 <puts@plt>
lea    0x1(%r14),%rax
movslq 0xc(%rsp),%rcx
cmp    %rcx,%r14
mov    %rax,%r14
mov    0x20(%rsp),%r13
jge    1390 <main+0x180>
mov    0x8(%rsp),%esi
test   %esi,%esi
jle    1729 <main+0x519>
lea    -0x1(%rsi),%rdx
mov    %esi,%eax
and    $0x3,%eax
mov    $0x1,%ecx
cmp    $0x3,%rdx
jb     1701 <main+0x4f1>
mov    %esi,%edx
and    $0xfffffffc,%edx
xor    %ecx,%ecx
nopl   (%rax)
mov    0x8(%r12,%rcx,8),%rdi
movss  (%rdi,%r14,4),%xmm0
movss  %xmm0,0x4(%rbx,%rcx,4)
mov    0x10(%r12,%rcx,8),%rdi
movss  (%rdi,%r14,4),%xmm0
movss  %xmm0,0x8(%rbx,%rcx,4)
mov    0x18(%r12,%rcx,8),%rdi
movss  (%rdi,%r14,4),%xmm0
movss  %xmm0,0xc(%rbx,%rcx,4)
mov    0x20(%r12,%rcx,8),%rdi
movss  (%rdi,%r14,4),%xmm0
movss  %xmm0,0x10(%rbx,%rcx,4)
add    $0x4,%rcx
cmp    %rcx,%rdx
jne    16b0 <main+0x4a0>
add    $0x1,%rcx
test   %rax,%rax
je     1729 <main+0x519>
lea    (%rbx,%rcx,4),%rdx
lea    (%r12,%rcx,8),%rcx
xor    %edi,%edi
mov    (%rcx,%rdi,8),%rbp
movss  0x0(%rbp,%r14,4),%xmm0
movss  %xmm0,(%rdx,%rdi,4)
add    $0x1,%rdi
cmp    %rdi,%rax
jne    1710 <main+0x500>
mov    %r13,%rdi
mov    0x18(%rsp),%rdx
mov    %rbx,%rcx
call   10e0 <lubksb@plt>
lea    0x94d(%rip),%rdi
call   1040 <puts@plt>
cmpl   $0x0,0x8(%rsp)
lea    0x8e5(%rip),%r13
jle    1787 <main+0x577>
xor    %ebp,%ebp
cs nopw 0x0(%rax,%rax,1)
nop
mov    0x8(%r12,%rbp,8),%rax
movss  (%rax,%r14,4),%xmm0
cvtss2sd %xmm0,%xmm0
mov    %r13,%rdi
mov    $0x1,%al
call   1080 <printf@plt>
movslq 0x8(%rsp),%rax
add    $0x1,%rbp
cmp    %rax,%rbp
jl     1760 <main+0x550>
lea    0x8af(%rip),%rdi
lea    0x8af(%rip),%rsi
lea    0x8c1(%rip),%rdx
xor    %eax,%eax
call   1080 <printf@plt>
mov    0x8(%rsp),%r8d
test   %r8d,%r8d
jle    1660 <main+0x450>
lea    0x1(%r8),%r9
mov    %r8d,%edx
and    $0xfffffffe,%edx
neg    %rdx
mov    $0x1,%esi
jmp    17dd <main+0x5cd>
cs nopw 0x0(%rax,%rax,1)
nop
add    $0x1,%rsi
cmp    %r9,%rsi
je     1880 <main+0x670>
mov    (%r12,%rsi,8),%rdi
movl   $0x0,(%rdi,%r14,4)
mov    (%r15,%rsi,8),%rbp
xorps  %xmm0,%xmm0
mov    $0x1,%eax
cmp    $0x2,%r9d
je     184d <main+0x63d>
mov    $0x1,%eax
xorps  %xmm0,%xmm0
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movss  0x0(%rbp,%rax,4),%xmm1
mulss  (%rbx,%rax,4),%xmm1
addss  %xmm0,%xmm1
movss  %xmm1,(%rdi,%r14,4)
movss  0x4(%rbp,%rax,4),%xmm0
mulss  0x4(%rbx,%rax,4),%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,(%rdi,%r14,4)
lea    (%rdx,%rax,1),%rcx
add    $0x2,%rcx
add    $0x2,%rax
cmp    $0x1,%rcx
jne    1810 <main+0x600>
test   $0x1,%r8b
je     17d0 <main+0x5c0>
movss  0x0(%rbp,%rax,4),%xmm1
mulss  (%rbx,%rax,4),%xmm1
addss  %xmm0,%xmm1
movss  %xmm1,(%rdi,%r14,4)
jmp    17d0 <main+0x5c0>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
test   %r8d,%r8d
lea    0x7ac(%rip),%r13
jle    1660 <main+0x450>
xor    %ebp,%ebp
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
mov    0x8(%r12,%rbp,8),%rax
movss  (%rax,%r14,4),%xmm0
cvtss2sd %xmm0,%xmm0
mov    %r13,%rdi
mov    $0x1,%al
call   1080 <printf@plt>
movslq 0x8(%rsp),%rax
add    $0x1,%rbp
cmp    %rax,%rbp
jl     18a0 <main+0x690>
jmp    1660 <main+0x450>

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

LLVM IR (-O2):
@.str = [11 x i8] c"matrx1.dat\00"
@.str.1 = [2 x i8] c"r\00"
@.str.2 = [32 x i8] c"Data file matrx1.dat not found\0A\00"
@.str.3 = [7 x i8] c"%d %d \00"
@.str.4 = [4 x i8] c"%f \00"
@.str.6 = [7 x i8] c"%12.6f\00"
@.str.7 = [7 x i8] c"\0A%s%s\0A\00"
@.str.8 = [25 x i8] c"result of matrix applied\00"
@.str.9 = [17 x i8] c" to sol'n vector\00"
@str = [31 x i8] c"press RETURN for next problem:\00"
@str.12 = [24 x i8] c"right-hand side vector:\00"
@str.13 = [35 x i8] c"\0A*********************************\00"
define i32 @main() {
%1 = alloca i32
%2 = alloca i32
%3 = alloca float
%4 = alloca [80 x i8]
%5 = bitcast i32* %1 to i8*
%6 = bitcast i32* %2 to i8*
%7 = bitcast float* %3 to i8*
%8 = getelementptr inbounds [80 x i8], [80 x i8]* %4, i64 0, i64 0
%9 = tail call i32* @ivector(i32 1, i32 20)
%10 = tail call float* @vector(i32 1, i32 20)
%11 = tail call float** @matrix(i32 1, i32 20, i32 1, i32 20)
%12 = tail call float** @matrix(i32 1, i32 20, i32 1, i32 20)
%13 = tail call float** @matrix(i32 1, i32 20, i32 1, i32 20)
%14 = tail call i32* @fopen(i8* getelementptr inbounds ([11 x i8], [11 x i8]* @.str, i64 0, i64 0), i8* getelementptr inbounds ([2 x i8], [2 x i8]* @.str.1, i64 0, i64 0))
%15 = icmp eq i32* %14, null
br i1 %15, label %16, label %18
16: ; preds = %0
%17 = tail call i32 @nrerror(i8* getelementptr inbounds ([32 x i8], [32 x i8]* @.str.2, i64 0, i64 0))
br label %18
18: ; preds = %16, %0
%19 = call i32 @feof(i32* %14)
%20 = icmp eq i32 %19, 0
br i1 %20, label %21, label %295
21: ; preds = %18, %290
%22 = call i32 @fgets(i8* %8, i32 80, i32* %14)
%23 = call i32 @fgets(i8* %8, i32 80, i32* %14)
%24 = call i32 (i32*, i8*, ...) @fscanf(i32* %14, i8* getelementptr inbounds ([7 x i8], [7 x i8]* @.str.3, i64 0, i64 0), i32* %2, i32* %1)
%25 = call i32 @fgets(i8* %8, i32 80, i32* %14)
%26 = load i32, i32* %2
%27 = icmp slt i32 %26, 1
br i1 %27, label %48, label %28
28: ; preds = %21, %43
%29 = phi i32 [ %44, %43 ], [ %26, %21 ]
%30 = phi i64 [ %46, %43 ], [ 1, %21 ]
%31 = icmp slt i32 %29, 1
br i1 %31, label %43, label %32
32: ; preds = %28
%33 = getelementptr inbounds float*, float** %11, i64 %30
br label %34
34: ; preds = %32, %34
%35 = phi i64 [ 1, %32 ], [ %39, %34 ]
%36 = load float*, float** %33
%37 = getelementptr inbounds float, float* %36, i64 %35
%38 = call i32 (i32*, i8*, ...) @fscanf(i32* %14, i8* getelementptr inbounds ([4 x i8], [4 x i8]* @.str.4, i64 0, i64 0), float* %37)
%39 = add nuw nsw i64 %35, 1
%40 = load i32, i32* %2
%41 = sext i32 %40 to i64
%42 = icmp slt i64 %35, %41
br i1 %42, label %34, label %43
43: ; preds = %34, %28
%44 = phi i32 [ %29, %28 ], [ %40, %34 ]
%45 = sext i32 %44 to i64
%46 = add nuw nsw i64 %30, 1
%47 = icmp slt i64 %30, %45
br i1 %47, label %28, label %48
48: ; preds = %43, %21
%49 = call i32 @fgets(i8* %8, i32 80, i32* %14)
%50 = load i32, i32* %1
%51 = icmp slt i32 %50, 1
%52 = load i32, i32* %2
br i1 %51, label %53, label %65
53: ; preds = %82, %48
%54 = phi i32 [ %52, %48 ], [ %84, %82 ]
%55 = icmp slt i32 %54, 1
br i1 %55, label %145, label %56
56: ; preds = %53
%57 = add nuw i32 %54, 1
%58 = zext i32 %57 to i64
%59 = add nsw i64 %58, -1
%60 = add nsw i64 %58, -2
%61 = and i64 %59, 3
%62 = icmp ult i64 %60, 3
%63 = and i64 %59, -4
%64 = icmp eq i64 %61, 0
br label %88
65: ; preds = %48, %82
%66 = phi i32 [ %83, %82 ], [ %50, %48 ]
%67 = phi i32 [ %84, %82 ], [ %52, %48 ]
%68 = phi i64 [ %85, %82 ], [ 1, %48 ]
%69 = icmp slt i32 %67, 1
br i1 %69, label %82, label %70
70: ; preds = %65, %70
%71 = phi i64 [ %76, %70 ], [ 1, %65 ]
%72 = getelementptr inbounds float*, float** %12, i64 %71
%73 = load float*, float** %72
%74 = getelementptr inbounds float, float* %73, i64 %68
%75 = call i32 (i32*, i8*, ...) @fscanf(i32* %14, i8* getelementptr inbounds ([4 x i8], [4 x i8]* @.str.4, i64 0, i64 0), float* %74)
%76 = add nuw nsw i64 %71, 1
%77 = load i32, i32* %2
%78 = sext i32 %77 to i64
%79 = icmp slt i64 %71, %78
br i1 %79, label %70, label %80
80: ; preds = %70
%81 = load i32, i32* %1
br label %82
82: ; preds = %80, %65
%83 = phi i32 [ %81, %80 ], [ %66, %65 ]
%84 = phi i32 [ %77, %80 ], [ %67, %65 ]
%85 = add nuw nsw i64 %68, 1
%86 = sext i32 %83 to i64
%87 = icmp slt i64 %68, %86
br i1 %87, label %65, label %53
88: ; preds = %142, %56
%89 = phi i64 [ 1, %56 ], [ %143, %142 ]
br i1 %62, label %127, label %90
90: ; preds = %88, %90
%91 = phi i64 [ %124, %90 ], [ 1, %88 ]
%92 = phi i64 [ %125, %90 ], [ 0, %88 ]
%93 = getelementptr inbounds float*, float** %11, i64 %91
%94 = load float*, float** %93
%95 = getelementptr inbounds float, float* %94, i64 %89
%96 = load float, float* %95
%97 = getelementptr inbounds float*, float** %13, i64 %91
%98 = load float*, float** %97
%99 = getelementptr inbounds float, float* %98, i64 %89
store float %96, float* %99
%100 = add nuw nsw i64 %91, 1
%101 = getelementptr inbounds float*, float** %11, i64 %100
%102 = load float*, float** %101
%103 = getelementptr inbounds float, float* %102, i64 %89
%104 = load float, float* %103
%105 = getelementptr inbounds float*, float** %13, i64 %100
%106 = load float*, float** %105
%107 = getelementptr inbounds float, float* %106, i64 %89
store float %104, float* %107
%108 = add nuw nsw i64 %91, 2
%109 = getelementptr inbounds float*, float** %11, i64 %108
%110 = load float*, float** %109
%111 = getelementptr inbounds float, float* %110, i64 %89
%112 = load float, float* %111
%113 = getelementptr inbounds float*, float** %13, i64 %108
%114 = load float*, float** %113
%115 = getelementptr inbounds float, float* %114, i64 %89
store float %112, float* %115
%116 = add nuw nsw i64 %91, 3
%117 = getelementptr inbounds float*, float** %11, i64 %116
%118 = load float*, float** %117
%119 = getelementptr inbounds float, float* %118, i64 %89
%120 = load float, float* %119
%121 = getelementptr inbounds float*, float** %13, i64 %116
%122 = load float*, float** %121
%123 = getelementptr inbounds float, float* %122, i64 %89
store float %120, float* %123
%124 = add nuw nsw i64 %91, 4
%125 = add i64 %92, 4
%126 = icmp eq i64 %125, %63
br i1 %126, label %127, label %90
127: ; preds = %90, %88
%128 = phi i64 [ 1, %88 ], [ %124, %90 ]
br i1 %64, label %142, label %129
129: ; preds = %127, %129
%130 = phi i64 [ %139, %129 ], [ %128, %127 ]
%131 = phi i64 [ %140, %129 ], [ 0, %127 ]
%132 = getelementptr inbounds float*, float** %11, i64 %130
%133 = load float*, float** %132
%134 = getelementptr inbounds float, float* %133, i64 %89
%135 = load float, float* %134
%136 = getelementptr inbounds float*, float** %13, i64 %130
%137 = load float*, float** %136
%138 = getelementptr inbounds float, float* %137, i64 %89
store float %135, float* %138
%139 = add nuw nsw i64 %130, 1
%140 = add i64 %131, 1
%141 = icmp eq i64 %140, %61
br i1 %141, label %142, label %129
142: ; preds = %129, %127
%143 = add nuw nsw i64 %89, 1
%144 = icmp eq i64 %143, %58
br i1 %144, label %145, label %88
145: ; preds = %142, %53
%146 = call i32 @ludcmp(float** %13, i32 %54, i32* %9, float* %3)
%147 = load i32, i32* %1
%148 = icmp slt i32 %147, 1
br i1 %148, label %290, label %149
149: ; preds = %145, %284
%150 = phi i64 [ %286, %284 ], [ 1, %145 ]
%151 = load i32, i32* %2
%152 = icmp slt i32 %151, 1
br i1 %152, label %205, label %153
153: ; preds = %149
%154 = add nuw i32 %151, 1
%155 = zext i32 %154 to i64
%156 = add nsw i64 %155, -1
%157 = add nsw i64 %155, -2
%158 = and i64 %156, 3
%159 = icmp ult i64 %157, 3
br i1 %159, label %191, label %160
160: ; preds = %153
%161 = and i64 %156, -4
br label %162
162: ; preds = %162, %160
%163 = phi i64 [ 1, %160 ], [ %188, %162 ]
%164 = phi i64 [ 0, %160 ], [ %189, %162 ]
%165 = getelementptr inbounds float*, float** %12, i64 %163
%166 = load float*, float** %165
%167 = getelementptr inbounds float, float* %166, i64 %150
%168 = load float, float* %167
%169 = getelementptr inbounds float, float* %10, i64 %163
store float %168, float* %169
%170 = add nuw nsw i64 %163, 1
%171 = getelementptr inbounds float*, float** %12, i64 %170
%172 = load float*, float** %171
%173 = getelementptr inbounds float, float* %172, i64 %150
%174 = load float, float* %173
%175 = getelementptr inbounds float, float* %10, i64 %170
store float %174, float* %175
%176 = add nuw nsw i64 %163, 2
%177 = getelementptr inbounds float*, float** %12, i64 %176
%178 = load float*, float** %177
%179 = getelementptr inbounds float, float* %178, i64 %150
%180 = load float, float* %179
%181 = getelementptr inbounds float, float* %10, i64 %176
store float %180, float* %181
%182 = add nuw nsw i64 %163, 3
%183 = getelementptr inbounds float*, float** %12, i64 %182
%184 = load float*, float** %183
%185 = getelementptr inbounds float, float* %184, i64 %150
%186 = load float, float* %185
%187 = getelementptr inbounds float, float* %10, i64 %182
store float %186, float* %187
%188 = add nuw nsw i64 %163, 4
%189 = add i64 %164, 4
%190 = icmp eq i64 %189, %161
br i1 %190, label %191, label %162
191: ; preds = %162, %153
%192 = phi i64 [ 1, %153 ], [ %188, %162 ]
%193 = icmp eq i64 %158, 0
br i1 %193, label %205, label %194
194: ; preds = %191, %194
%195 = phi i64 [ %202, %194 ], [ %192, %191 ]
%196 = phi i64 [ %203, %194 ], [ 0, %191 ]
%197 = getelementptr inbounds float*, float** %12, i64 %195
%198 = load float*, float** %197
%199 = getelementptr inbounds float, float* %198, i64 %150
%200 = load float, float* %199
%201 = getelementptr inbounds float, float* %10, i64 %195
store float %200, float* %201
%202 = add nuw nsw i64 %195, 1
%203 = add i64 %196, 1
%204 = icmp eq i64 %203, %158
br i1 %204, label %205, label %194
205: ; preds = %191, %194, %149
%206 = call i32 @lubksb(float** %13, i32 %151, i32* %9, float* %10)
%207 = call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([24 x i8], [24 x i8]* @str.12, i64 0, i64 0))
%208 = load i32, i32* %2
%209 = icmp slt i32 %208, 1
br i1 %209, label %222, label %210
210: ; preds = %205, %210
%211 = phi i64 [ %218, %210 ], [ 1, %205 ]
%212 = getelementptr inbounds float*, float** %12, i64 %211
%213 = load float*, float** %212
%214 = getelementptr inbounds float, float* %213, i64 %150
%215 = load float, float* %214
%216 = fpext float %215 to double
%217 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([7 x i8], [7 x i8]* @.str.6, i64 0, i64 0), double %216)
%218 = add nuw nsw i64 %211, 1
%219 = load i32, i32* %2
%220 = sext i32 %219 to i64
%221 = icmp slt i64 %211, %220
br i1 %221, label %210, label %222
222: ; preds = %210, %205
%223 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([7 x i8], [7 x i8]* @.str.7, i64 0, i64 0), i8* getelementptr inbounds ([25 x i8], [25 x i8]* @.str.8, i64 0, i64 0), i8* getelementptr inbounds ([17 x i8], [17 x i8]* @.str.9, i64 0, i64 0))
%224 = load i32, i32* %2
%225 = icmp slt i32 %224, 1
br i1 %225, label %284, label %226
226: ; preds = %222
%227 = add nuw i32 %224, 1
%228 = zext i32 %227 to i64
%229 = add nsw i64 %228, -1
%230 = and i64 %229, 1
%231 = icmp eq i32 %227, 2
%232 = and i64 %229, -2
%233 = icmp eq i64 %230, 0
br label %235
234: ; preds = %269
br i1 %225, label %284, label %272
235: ; preds = %269, %226
%236 = phi i64 [ 1, %226 ], [ %270, %269 ]
%237 = getelementptr inbounds float*, float** %12, i64 %236
%238 = load float*, float** %237
%239 = getelementptr inbounds float, float* %238, i64 %150
store float 0.000000e+00, float* %239
%240 = getelementptr inbounds float*, float** %11, i64 %236
%241 = load float*, float** %240
br i1 %231, label %260, label %242
242: ; preds = %235, %242
%243 = phi float [ %256, %242 ], [ 0.000000e+00, %235 ]
%244 = phi i64 [ %257, %242 ], [ 1, %235 ]
%245 = phi i64 [ %258, %242 ], [ 0, %235 ]
%246 = getelementptr inbounds float, float* %241, i64 %244
%247 = load float, float* %246
%248 = getelementptr inbounds float, float* %10, i64 %244
%249 = load float, float* %248
%250 = call float @llvm.fmuladd.f32(float %247, float %249, float %243)
store float %250, float* %239
%251 = add nuw nsw i64 %244, 1
%252 = getelementptr inbounds float, float* %241, i64 %251
%253 = load float, float* %252
%254 = getelementptr inbounds float, float* %10, i64 %251
%255 = load float, float* %254
%256 = call float @llvm.fmuladd.f32(float %253, float %255, float %250)
store float %256, float* %239
%257 = add nuw nsw i64 %244, 2
%258 = add i64 %245, 2
%259 = icmp eq i64 %258, %232
br i1 %259, label %260, label %242
260: ; preds = %242, %235
%261 = phi float [ 0.000000e+00, %235 ], [ %256, %242 ]
%262 = phi i64 [ 1, %235 ], [ %257, %242 ]
br i1 %233, label %269, label %263
263: ; preds = %260
%264 = getelementptr inbounds float, float* %241, i64 %262
%265 = load float, float* %264
%266 = getelementptr inbounds float, float* %10, i64 %262
%267 = load float, float* %266
%268 = call float @llvm.fmuladd.f32(float %265, float %267, float %261)
store float %268, float* %239
br label %269
269: ; preds = %260, %263
%270 = add nuw nsw i64 %236, 1
%271 = icmp eq i64 %270, %228
br i1 %271, label %234, label %235
272: ; preds = %234, %272
%273 = phi i64 [ %280, %272 ], [ 1, %234 ]
%274 = getelementptr inbounds float*, float** %12, i64 %273
%275 = load float*, float** %274
%276 = getelementptr inbounds float, float* %275, i64 %150
%277 = load float, float* %276
%278 = fpext float %277 to double
%279 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([7 x i8], [7 x i8]* @.str.6, i64 0, i64 0), double %278)
%280 = add nuw nsw i64 %273, 1
%281 = load i32, i32* %2
%282 = sext i32 %281 to i64
%283 = icmp slt i64 %273, %282
br i1 %283, label %272, label %284
284: ; preds = %272, %222, %234
%285 = call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([35 x i8], [35 x i8]* @str.13, i64 0, i64 0))
%286 = add nuw nsw i64 %150, 1
%287 = load i32, i32* %1
%288 = sext i32 %287 to i64
%289 = icmp slt i64 %150, %288
br i1 %289, label %149, label %290
290: ; preds = %284, %145
%291 = call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([31 x i8], [31 x i8]* @str, i64 0, i64 0))
%292 = call i32 (...) @getchar()
%293 = call i32 @feof(i32* %14)
%294 = icmp eq i32 %293, 0
br i1 %294, label %21, label %295
295: ; preds = %290, %18
%296 = call i32 @fclose(i32* %14)
%297 = call i32 @free_matrix(float** %13, i32 1, i32 20, i32 1, i32 20)
%298 = call i32 @free_matrix(float** %12, i32 1, i32 20, i32 1, i32 20)
%299 = call i32 @free_matrix(float** %11, i32 1, i32 20, i32 1, i32 20)
%300 = call i32 @free_vector(float* %10, i32 1, i32 20)
%301 = call i32 @free_ivector(i32* %9, i32 1, i32 20)
ret i32 0
}