#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_vector(float *, int, int);
int /*<<< orphan*/ fscanf(int /*<<< orphan*/ *, char *, ...);
int /*<<< orphan*/ getchar();
int /*<<< orphan*/ nrerror(char *);
int /*<<< orphan*/ printf(char *, ...);
int /*<<< orphan*/ tridag(float *, float *, float *, float *, float *,
                          unsigned long);
float *vector(int, int);
int main(void) {
  unsigned long k, n;
  float *diag, *superd, *subd, *rhs, *u;
  char dummy[80];
  FILE *fp;
  diag = vector(1, 20);
  superd = vector(1, 20);
  subd = vector(1, 20);
  rhs = vector(1, 20);
  u = vector(1, 20);
  if ((fp = fopen("matrx2.dat", "r")) == NULL)
    nrerror("Data file matrx2.dat not found\n");
  while (!feof(fp)) {
    fgets(dummy, 80, fp);
    fgets(dummy, 80, fp);
    fscanf(fp, "%ld ", &n);
    fgets(dummy, 80, fp);
    for (k = 1; k <= n; k++)
      fscanf(fp, "%f ", &diag[k]);
    fgets(dummy, 80, fp);
    for (k = 1; k < n; k++)
      fscanf(fp, "%f ", &superd[k]);
    fgets(dummy, 80, fp);
    for (k = 2; k <= n; k++)
      fscanf(fp, "%f ", &subd[k]);
    fgets(dummy, 80, fp);
    for (k = 1; k <= n; k++)
      fscanf(fp, "%f ", &rhs[k]);
    tridag(subd, diag, superd, rhs, u, n);
    printf("\nThe solution vector is:\n");
    for (k = 1; k <= n; k++)
      printf("%12.6f", u[k]);
    printf("\n");
    printf("\n(matrix)*(sol'n vector) should be:\n");
    for (k = 1; k <= n; k++)
      printf("%12.6f", rhs[k]);
    printf("\n");
    printf("actual result is:\n");
    for (k = 1; k <= n; k++) {
      if (k == 1)
        rhs[k] = diag[1] * u[1] + superd[1] * u[2];
      else if (k == n)
        rhs[k] = subd[n] * u[n - 1] + diag[n] * u[n];
      else
        rhs[k] = subd[k] * u[k - 1] + diag[k] * u[k] + superd[k] * u[k + 1];
    }
    for (k = 1; k <= n; k++)
      printf("%12.6f", rhs[k]);
    printf("\n");
    printf("***********************************\n");
    printf("press return for next problem:\n");
    (void)getchar();
  }
  fclose(fp);
  free_vector(u, 1, 20);
  free_vector(rhs, 1, 20);
  free_vector(subd, 1, 20);
  free_vector(superd, 1, 20);
  free_vector(diag, 1, 20);
  return 0;
}


x86 assembly (clang, -O0):
<main>:
push   %rbp
mov    %rsp,%rbp
sub    $0xa0,%rsp
movl   $0x0,-0x4(%rbp)
mov    $0x1,%edi
mov    $0x14,%esi
call   10b0 <vector@plt>
mov    %rax,-0x20(%rbp)
mov    $0x1,%edi
mov    $0x14,%esi
call   10b0 <vector@plt>
mov    %rax,-0x28(%rbp)
mov    $0x1,%edi
mov    $0x14,%esi
call   10b0 <vector@plt>
mov    %rax,-0x30(%rbp)
mov    $0x1,%edi
mov    $0x14,%esi
call   10b0 <vector@plt>
mov    %rax,-0x38(%rbp)
mov    $0x1,%edi
mov    $0x14,%esi
call   10b0 <vector@plt>
mov    %rax,-0x40(%rbp)
lea    0xdd8(%rip),%rdi
lea    0xddc(%rip),%rsi
call   10c0 <fopen@plt>
mov    %rax,-0x98(%rbp)
cmp    $0x0,%rax
jne    1251 <main+0xa1>
lea    0xdc1(%rip),%rdi
call   1030 <nrerror@plt>
jmp    1256 <main+0xa6>
mov    -0x98(%rbp),%rdi
call   10a0 <feof@plt>
cmp    $0x0,%eax
setne  %al
xor    $0xff,%al
test   $0x1,%al
jne    1277 <main+0xc7>
jmp    1701 <main+0x551>
lea    -0x90(%rbp),%rdi
mov    -0x98(%rbp),%rdx
mov    $0x50,%esi
call   1080 <fgets@plt>
lea    -0x90(%rbp),%rdi
mov    -0x98(%rbp),%rdx
mov    $0x50,%esi
call   1080 <fgets@plt>
mov    -0x98(%rbp),%rdi
lea    0xd78(%rip),%rsi
lea    -0x18(%rbp),%rdx
mov    $0x0,%al
call   1070 <fscanf@plt>
lea    -0x90(%rbp),%rdi
mov    -0x98(%rbp),%rdx
mov    $0x50,%esi
call   1080 <fgets@plt>
movq   $0x1,-0x10(%rbp)
mov    -0x10(%rbp),%rax
cmp    -0x18(%rbp),%rax
ja     1323 <main+0x173>
mov    -0x98(%rbp),%rdi
mov    -0x20(%rbp),%rdx
mov    -0x10(%rbp),%rax
shl    $0x2,%rax
add    %rax,%rdx
lea    0xd27(%rip),%rsi
mov    $0x0,%al
call   1070 <fscanf@plt>
mov    -0x10(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x10(%rbp)
jmp    12e0 <main+0x130>
lea    -0x90(%rbp),%rdi
mov    -0x98(%rbp),%rdx
mov    $0x50,%esi
call   1080 <fgets@plt>
movq   $0x1,-0x10(%rbp)
mov    -0x10(%rbp),%rax
cmp    -0x18(%rbp),%rax
jae    1386 <main+0x1d6>
mov    -0x98(%rbp),%rdi
mov    -0x28(%rbp),%rdx
mov    -0x10(%rbp),%rax
shl    $0x2,%rax
add    %rax,%rdx
lea    0xcc4(%rip),%rsi
mov    $0x0,%al
call   1070 <fscanf@plt>
mov    -0x10(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x10(%rbp)
jmp    1343 <main+0x193>
lea    -0x90(%rbp),%rdi
mov    -0x98(%rbp),%rdx
mov    $0x50,%esi
call   1080 <fgets@plt>
movq   $0x2,-0x10(%rbp)
mov    -0x10(%rbp),%rax
cmp    -0x18(%rbp),%rax
ja     13e9 <main+0x239>
mov    -0x98(%rbp),%rdi
mov    -0x30(%rbp),%rdx
mov    -0x10(%rbp),%rax
shl    $0x2,%rax
add    %rax,%rdx
lea    0xc61(%rip),%rsi
mov    $0x0,%al
call   1070 <fscanf@plt>
mov    -0x10(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x10(%rbp)
jmp    13a6 <main+0x1f6>
lea    -0x90(%rbp),%rdi
mov    -0x98(%rbp),%rdx
mov    $0x50,%esi
call   1080 <fgets@plt>
movq   $0x1,-0x10(%rbp)
mov    -0x10(%rbp),%rax
cmp    -0x18(%rbp),%rax
ja     144c <main+0x29c>
mov    -0x98(%rbp),%rdi
mov    -0x38(%rbp),%rdx
mov    -0x10(%rbp),%rax
shl    $0x2,%rax
add    %rax,%rdx
lea    0xbfe(%rip),%rsi
mov    $0x0,%al
call   1070 <fscanf@plt>
mov    -0x10(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x10(%rbp)
jmp    1409 <main+0x259>
mov    -0x30(%rbp),%rdi
mov    -0x20(%rbp),%rsi
mov    -0x28(%rbp),%rdx
mov    -0x38(%rbp),%rcx
mov    -0x40(%rbp),%r8
mov    -0x18(%rbp),%r9
call   10d0 <tridag@plt>
lea    0xbc6(%rip),%rdi
mov    $0x0,%al
call   1050 <printf@plt>
movq   $0x1,-0x10(%rbp)
mov    -0x10(%rbp),%rax
cmp    -0x18(%rbp),%rax
ja     14bd <main+0x30d>
mov    -0x40(%rbp),%rax
mov    -0x10(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
lea    0xbab(%rip),%rdi
mov    $0x1,%al
call   1050 <printf@plt>
mov    -0x10(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x10(%rbp)
jmp    147f <main+0x2cf>
lea    0xbee(%rip),%rdi
mov    $0x0,%al
call   1050 <printf@plt>
lea    0xb85(%rip),%rdi
mov    $0x0,%al
call   1050 <printf@plt>
movq   $0x1,-0x10(%rbp)
mov    -0x10(%rbp),%rax
cmp    -0x18(%rbp),%rax
ja     151f <main+0x36f>
mov    -0x38(%rbp),%rax
mov    -0x10(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
lea    0xb49(%rip),%rdi
mov    $0x1,%al
call   1050 <printf@plt>
mov    -0x10(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x10(%rbp)
jmp    14e1 <main+0x331>
lea    0xb8c(%rip),%rdi
mov    $0x0,%al
call   1050 <printf@plt>
lea    0xb48(%rip),%rdi
mov    $0x0,%al
call   1050 <printf@plt>
movq   $0x1,-0x10(%rbp)
mov    -0x10(%rbp),%rax
cmp    -0x18(%rbp),%rax
ja     1685 <main+0x4d5>
cmpq   $0x1,-0x10(%rbp)
jne    159a <main+0x3ea>
mov    -0x20(%rbp),%rax
movss  0x4(%rax),%xmm0
mov    -0x40(%rbp),%rax
movss  0x4(%rax),%xmm2
mov    -0x28(%rbp),%rax
movss  0x4(%rax),%xmm1
mov    -0x40(%rbp),%rax
mulss  0x8(%rax),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x38(%rbp),%rax
mov    -0x10(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
jmp    166f <main+0x4bf>
mov    -0x10(%rbp),%rax
cmp    -0x18(%rbp),%rax
jne    15fa <main+0x44a>
mov    -0x30(%rbp),%rax
mov    -0x18(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x40(%rbp),%rax
mov    -0x18(%rbp),%rcx
sub    $0x1,%rcx
movss  (%rax,%rcx,4),%xmm2
mov    -0x20(%rbp),%rax
mov    -0x18(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
mov    -0x40(%rbp),%rax
mov    -0x18(%rbp),%rcx
mulss  (%rax,%rcx,4),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x38(%rbp),%rax
mov    -0x10(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
jmp    166a <main+0x4ba>
mov    -0x30(%rbp),%rax
mov    -0x10(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
mov    -0x40(%rbp),%rax
mov    -0x10(%rbp),%rcx
sub    $0x1,%rcx
movss  (%rax,%rcx,4),%xmm2
mov    -0x20(%rbp),%rax
mov    -0x10(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x40(%rbp),%rax
mov    -0x10(%rbp),%rcx
mulss  (%rax,%rcx,4),%xmm0
mulss  %xmm2,%xmm1
addss  %xmm0,%xmm1
mov    -0x28(%rbp),%rax
mov    -0x10(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x40(%rbp),%rax
mov    -0x10(%rbp),%rcx
movss  0x4(%rax,%rcx,4),%xmm2
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
mov    -0x38(%rbp),%rax
mov    -0x10(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
jmp    166f <main+0x4bf>
jmp    1674 <main+0x4c4>
mov    -0x10(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x10(%rbp)
jmp    1543 <main+0x393>
movq   $0x1,-0x10(%rbp)
mov    -0x10(%rbp),%rax
cmp    -0x18(%rbp),%rax
ja     16cb <main+0x51b>
mov    -0x38(%rbp),%rax
mov    -0x10(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
lea    0x99d(%rip),%rdi
mov    $0x1,%al
call   1050 <printf@plt>
mov    -0x10(%rbp),%rax
add    $0x1,%rax
mov    %rax,-0x10(%rbp)
jmp    168d <main+0x4dd>
lea    0x9e0(%rip),%rdi
mov    $0x0,%al
call   1050 <printf@plt>
lea    0x9af(%rip),%rdi
mov    $0x0,%al
call   1050 <printf@plt>
lea    0x9c6(%rip),%rdi
mov    $0x0,%al
call   1050 <printf@plt>
mov    $0x0,%al
call   1090 <getchar@plt>
jmp    1256 <main+0xa6>
mov    -0x98(%rbp),%rdi
call   1040 <fclose@plt>
mov    -0x40(%rbp),%rdi
mov    $0x1,%esi
mov    $0x14,%edx
call   1060 <free_vector@plt>
mov    -0x38(%rbp),%rdi
mov    $0x1,%esi
mov    $0x14,%edx
call   1060 <free_vector@plt>
mov    -0x30(%rbp),%rdi
mov    $0x1,%esi
mov    $0x14,%edx
call   1060 <free_vector@plt>
mov    -0x28(%rbp),%rdi
mov    $0x1,%esi
mov    $0x14,%edx
call   1060 <free_vector@plt>
mov    -0x20(%rbp),%rdi
mov    $0x1,%esi
mov    $0x14,%edx
call   1060 <free_vector@plt>
xor    %eax,%eax
add    $0xa0,%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   10d0 <vector@plt>
mov    %rax,0x10(%rsp)
mov    $0x1,%edi
mov    $0x14,%esi
call   10d0 <vector@plt>
mov    %rax,0x8(%rsp)
mov    $0x1,%edi
mov    $0x14,%esi
call   10d0 <vector@plt>
mov    %rax,0x18(%rsp)
mov    $0x1,%edi
mov    $0x14,%esi
call   10d0 <vector@plt>
mov    %rax,%r13
mov    $0x1,%edi
mov    $0x14,%esi
call   10d0 <vector@plt>
mov    %rax,%r14
lea    0xdb8(%rip),%rdi
lea    0xdbc(%rip),%rsi
call   10e0 <fopen@plt>
mov    %rax,%r15
test   %rax,%rax
jne    1268 <main+0x98>
lea    0xdaa(%rip),%rdi
call   1050 <nrerror@plt>
mov    %r15,%rdi
call   10c0 <feof@plt>
test   %eax,%eax
je     12f0 <main+0x120>
mov    %r15,%rdi
call   1060 <fclose@plt>
mov    %r14,%rdi
mov    $0x1,%esi
mov    $0x14,%edx
call   1080 <free_vector@plt>
mov    %r13,%rdi
mov    $0x1,%esi
mov    $0x14,%edx
call   1080 <free_vector@plt>
mov    0x18(%rsp),%rdi
mov    $0x1,%esi
mov    $0x14,%edx
call   1080 <free_vector@plt>
mov    0x8(%rsp),%rdi
mov    $0x1,%esi
mov    $0x14,%edx
call   1080 <free_vector@plt>
mov    0x10(%rsp),%rdi
mov    $0x1,%esi
mov    $0x14,%edx
call   1080 <free_vector@plt>
xor    %eax,%eax
add    $0x98,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
mov    0x10(%rsp),%rax
add    $0x4,%rax
mov    %rax,0x38(%rsp)
mov    0x8(%rsp),%rax
add    $0x4,%rax
mov    %rax,0x28(%rsp)
mov    0x18(%rsp),%rax
add    $0x8,%rax
mov    %rax,0x20(%rsp)
mov    %r13,%rax
add    $0x4,%rax
mov    %rax,0x30(%rsp)
jmp    1369 <main+0x199>
nopl   0x0(%rax,%rax,1)
mov    $0xa,%edi
call   1030 <putchar@plt>
lea    0xd4b(%rip),%rdi
call   1040 <puts@plt>
lea    0xd63(%rip),%rdi
call   1040 <puts@plt>
xor    %eax,%eax
call   10b0 <getchar@plt>
mov    %r15,%rdi
call   10c0 <feof@plt>
test   %eax,%eax
jne    1274 <main+0xa4>
lea    0x40(%rsp),%r12
mov    %r12,%rdi
mov    $0x50,%esi
mov    %r15,%rdx
call   10a0 <fgets@plt>
mov    %r12,%rdi
mov    $0x50,%esi
mov    %r15,%rdx
call   10a0 <fgets@plt>
mov    %r15,%rdi
lea    0xc95(%rip),%rsi
mov    %rsp,%rdx
xor    %eax,%eax
call   1090 <fscanf@plt>
mov    %r12,%rdi
mov    $0x50,%esi
mov    %r15,%rdx
call   10a0 <fgets@plt>
cmpq   $0x0,(%rsp)
lea    0xc74(%rip),%r12
je     13ee <main+0x21e>
mov    $0x1,%ebp
mov    0x38(%rsp),%rbx
nopw   0x0(%rax,%rax,1)
mov    %r15,%rdi
mov    %r12,%rsi
mov    %rbx,%rdx
xor    %eax,%eax
call   1090 <fscanf@plt>
add    $0x1,%rbp
add    $0x4,%rbx
cmp    (%rsp),%rbp
jbe    13d0 <main+0x200>
lea    0x40(%rsp),%rdi
mov    $0x50,%esi
mov    %r15,%rdx
call   10a0 <fgets@plt>
cmpq   $0x2,(%rsp)
jb     143e <main+0x26e>
mov    $0x1,%ebp
mov    0x28(%rsp),%rbx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    %r15,%rdi
mov    %r12,%rsi
mov    %rbx,%rdx
xor    %eax,%eax
call   1090 <fscanf@plt>
add    $0x1,%rbp
add    $0x4,%rbx
cmp    (%rsp),%rbp
jb     1420 <main+0x250>
lea    0x40(%rsp),%rdi
mov    $0x50,%esi
mov    %r15,%rdx
call   10a0 <fgets@plt>
cmpq   $0x2,(%rsp)
jb     148e <main+0x2be>
mov    $0x2,%ebp
mov    0x20(%rsp),%rbx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    %r15,%rdi
mov    %r12,%rsi
mov    %rbx,%rdx
xor    %eax,%eax
call   1090 <fscanf@plt>
add    $0x1,%rbp
add    $0x4,%rbx
cmp    (%rsp),%rbp
jbe    1470 <main+0x2a0>
lea    0x40(%rsp),%rdi
mov    $0x50,%esi
mov    %r15,%rdx
call   10a0 <fgets@plt>
cmpq   $0x0,(%rsp)
je     14f0 <main+0x320>
mov    $0x1,%ebx
mov    0x30(%rsp),%r12
lea    0xb7a(%rip),%rbp
nopl   0x0(%rax,%rax,1)
mov    %r15,%rdi
mov    %rbp,%rsi
mov    %r12,%rdx
xor    %eax,%eax
call   1090 <fscanf@plt>
add    $0x1,%rbx
mov    (%rsp),%r9
add    $0x4,%r12
cmp    %r9,%rbx
jbe    14c0 <main+0x2f0>
jmp    14f3 <main+0x323>
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
xor    %r9d,%r9d
mov    0x18(%rsp),%rdi
mov    0x10(%rsp),%rsi
mov    0x8(%rsp),%rdx
mov    %r13,%rcx
mov    %r14,%r8
call   10f0 <tridag@plt>
lea    0xb29(%rip),%rdi
call   1040 <puts@plt>
cmpq   $0x0,(%rsp)
lea    0xb11(%rip),%r12
je     154e <main+0x37e>
mov    $0x1,%ebx
nopl   0x0(%rax)
movss  (%r14,%rbx,4),%xmm0
cvtss2sd %xmm0,%xmm0
mov    %r12,%rdi
mov    $0x1,%al
call   1070 <printf@plt>
add    $0x1,%rbx
cmp    (%rsp),%rbx
jbe    1530 <main+0x360>
mov    $0xa,%edi
call   1030 <putchar@plt>
lea    0xaf7(%rip),%rdi
call   1040 <puts@plt>
cmpq   $0x0,(%rsp)
je     158f <main+0x3bf>
mov    $0x1,%ebx
movss  0x0(%r13,%rbx,4),%xmm0
cvtss2sd %xmm0,%xmm0
mov    %r12,%rdi
mov    $0x1,%al
call   1070 <printf@plt>
add    $0x1,%rbx
cmp    (%rsp),%rbx
jbe    1570 <main+0x3a0>
mov    $0xa,%edi
call   1030 <putchar@plt>
lea    0xada(%rip),%rdi
call   1040 <puts@plt>
mov    (%rsp),%rax
test   %rax,%rax
je     1330 <main+0x160>
mov    0x10(%rsp),%rcx
movss  0x4(%rcx),%xmm0
mov    0x8(%rsp),%rcx
movss  0x4(%rcx),%xmm1
mulss  0x8(%r14),%xmm1
mulss  0x4(%r14),%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,0x4(%r13)
cmp    $0x1,%rax
je     1690 <main+0x4c0>
lea    -0x1(%rax),%rcx
lea    -0x2(%rax),%rdx
xor    %esi,%esi
mov    0x18(%rsp),%r8
mov    0x20(%rsp),%rbx
mov    0x8(%rsp),%rbp
mov    0x10(%rsp),%rdi
jmp    164e <main+0x47e>
cs nopw 0x0(%rax,%rax,1)
movss  (%rbx,%rsi,4),%xmm0
movss  0x8(%rdi,%rsi,4),%xmm1
mulss  0x8(%r14,%rsi,4),%xmm1
mulss  0x4(%r14,%rsi,4),%xmm0
addss  %xmm1,%xmm0
movss  0x8(%rbp,%rsi,4),%xmm1
mulss  0xc(%r14,%rsi,4),%xmm1
addss  %xmm0,%xmm1
movss  %xmm1,0x8(%r13,%rsi,4)
add    $0x1,%rsi
cmp    %rsi,%rcx
je     1680 <main+0x4b0>
cmp    %rsi,%rdx
jne    1610 <main+0x440>
movss  (%r8,%rax,4),%xmm0
movss  (%rdi,%rax,4),%xmm1
mulss  (%r14,%rax,4),%xmm1
mulss  -0x4(%r14,%rax,4),%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,0x0(%r13,%rax,4)
jmp    1645 <main+0x475>
nopl   0x0(%rax,%rax,1)
test   %rax,%rax
lea    0x9ac(%rip),%r12
je     1330 <main+0x160>
mov    $0x1,%ebx
cs nopw 0x0(%rax,%rax,1)
nop
movss  0x0(%r13,%rbx,4),%xmm0
cvtss2sd %xmm0,%xmm0
mov    %r12,%rdi
mov    $0x1,%al
call   1070 <printf@plt>
add    $0x1,%rbx
cmp    (%rsp),%rbx
jbe    16a0 <main+0x4d0>
jmp    1330 <main+0x160>

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

LLVM IR (-O2):
@.str = [11 x i8] c"matrx2.dat\00"
@.str.1 = [2 x i8] c"r\00"
@.str.2 = [32 x i8] c"Data file matrx2.dat not found\0A\00"
@.str.3 = [5 x i8] c"%ld \00"
@.str.4 = [4 x i8] c"%f \00"
@.str.6 = [7 x i8] c"%12.6f\00"
@str = [25 x i8] c"\0AThe solution vector is:\00"
@str.12 = [36 x i8] c"\0A(matrix)*(sol'n vector) should be:\00"
@str.13 = [18 x i8] c"actual result is:\00"
@str.14 = [36 x i8] c"***********************************\00"
@str.15 = [31 x i8] c"press return for next problem:\00"
define i32 @main() {
%1 = alloca i64
%2 = alloca [80 x i8]
%3 = bitcast i64* %1 to i8*
%4 = getelementptr inbounds [80 x i8], [80 x i8]* %2, i64 0, i64 0
%5 = tail call float* @vector(i32 1, i32 20)
%6 = tail call float* @vector(i32 1, i32 20)
%7 = tail call float* @vector(i32 1, i32 20)
%8 = tail call float* @vector(i32 1, i32 20)
%9 = tail call float* @vector(i32 1, i32 20)
%10 = 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))
%11 = icmp eq i32* %10, null
br i1 %11, label %12, label %14
12: ; preds = %0
%13 = tail call i32 @nrerror(i8* getelementptr inbounds ([32 x i8], [32 x i8]* @.str.2, i64 0, i64 0))
br label %14
14: ; preds = %12, %0
%15 = call i32 @feof(i32* %10)
%16 = icmp eq i32 %15, 0
br i1 %16, label %17, label %169
17: ; preds = %14
%18 = getelementptr inbounds float, float* %5, i64 1
%19 = getelementptr inbounds float, float* %9, i64 1
%20 = getelementptr inbounds float, float* %6, i64 1
%21 = getelementptr inbounds float, float* %9, i64 2
%22 = getelementptr inbounds float, float* %8, i64 1
br label %23
23: ; preds = %17, %162
%24 = call i32 @fgets(i8* %4, i32 80, i32* %10)
%25 = call i32 @fgets(i8* %4, i32 80, i32* %10)
%26 = call i32 (i32*, i8*, ...) @fscanf(i32* %10, i8* getelementptr inbounds ([5 x i8], [5 x i8]* @.str.3, i64 0, i64 0), i64* %1)
%27 = call i32 @fgets(i8* %4, i32 80, i32* %10)
%28 = load i64, i64* %1
%29 = icmp eq i64 %28, 0
br i1 %29, label %37, label %30
30: ; preds = %23, %30
%31 = phi i64 [ %34, %30 ], [ 1, %23 ]
%32 = getelementptr inbounds float, float* %5, i64 %31
%33 = call i32 (i32*, i8*, ...) @fscanf(i32* %10, i8* getelementptr inbounds ([4 x i8], [4 x i8]* @.str.4, i64 0, i64 0), float* %32)
%34 = add i64 %31, 1
%35 = load i64, i64* %1
%36 = icmp ugt i64 %34, %35
br i1 %36, label %37, label %30
37: ; preds = %30, %23
%38 = call i32 @fgets(i8* %4, i32 80, i32* %10)
%39 = load i64, i64* %1
%40 = icmp ugt i64 %39, 1
br i1 %40, label %41, label %48
41: ; preds = %37, %41
%42 = phi i64 [ %45, %41 ], [ 1, %37 ]
%43 = getelementptr inbounds float, float* %6, i64 %42
%44 = call i32 (i32*, i8*, ...) @fscanf(i32* %10, i8* getelementptr inbounds ([4 x i8], [4 x i8]* @.str.4, i64 0, i64 0), float* %43)
%45 = add nuw i64 %42, 1
%46 = load i64, i64* %1
%47 = icmp ult i64 %45, %46
br i1 %47, label %41, label %48
48: ; preds = %41, %37
%49 = call i32 @fgets(i8* %4, i32 80, i32* %10)
%50 = load i64, i64* %1
%51 = icmp ult i64 %50, 2
br i1 %51, label %59, label %52
52: ; preds = %48, %52
%53 = phi i64 [ %56, %52 ], [ 2, %48 ]
%54 = getelementptr inbounds float, float* %7, i64 %53
%55 = call i32 (i32*, i8*, ...) @fscanf(i32* %10, i8* getelementptr inbounds ([4 x i8], [4 x i8]* @.str.4, i64 0, i64 0), float* %54)
%56 = add i64 %53, 1
%57 = load i64, i64* %1
%58 = icmp ugt i64 %56, %57
br i1 %58, label %59, label %52
59: ; preds = %52, %48
%60 = call i32 @fgets(i8* %4, i32 80, i32* %10)
%61 = load i64, i64* %1
%62 = icmp eq i64 %61, 0
br i1 %62, label %70, label %63
63: ; preds = %59, %63
%64 = phi i64 [ %67, %63 ], [ 1, %59 ]
%65 = getelementptr inbounds float, float* %8, i64 %64
%66 = call i32 (i32*, i8*, ...) @fscanf(i32* %10, i8* getelementptr inbounds ([4 x i8], [4 x i8]* @.str.4, i64 0, i64 0), float* %65)
%67 = add i64 %64, 1
%68 = load i64, i64* %1
%69 = icmp ugt i64 %67, %68
br i1 %69, label %70, label %63
70: ; preds = %63, %59
%71 = phi i64 [ 0, %59 ], [ %68, %63 ]
%72 = call i32 @tridag(float* %7, float* %5, float* %6, float* %8, float* %9, i64 %71)
%73 = call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([25 x i8], [25 x i8]* @str, i64 0, i64 0))
%74 = load i64, i64* %1
%75 = icmp eq i64 %74, 0
br i1 %75, label %85, label %76
76: ; preds = %70, %76
%77 = phi i64 [ %82, %76 ], [ 1, %70 ]
%78 = getelementptr inbounds float, float* %9, i64 %77
%79 = load float, float* %78
%80 = fpext float %79 to double
%81 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([7 x i8], [7 x i8]* @.str.6, i64 0, i64 0), double %80)
%82 = add i64 %77, 1
%83 = load i64, i64* %1
%84 = icmp ugt i64 %82, %83
br i1 %84, label %85, label %76
85: ; preds = %76, %70
%86 = call i32 @putchar(i32 10)
%87 = call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([36 x i8], [36 x i8]* @str.12, i64 0, i64 0))
%88 = load i64, i64* %1
%89 = icmp eq i64 %88, 0
br i1 %89, label %99, label %90
90: ; preds = %85, %90
%91 = phi i64 [ %96, %90 ], [ 1, %85 ]
%92 = getelementptr inbounds float, float* %8, i64 %91
%93 = load float, float* %92
%94 = fpext float %93 to double
%95 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([7 x i8], [7 x i8]* @.str.6, i64 0, i64 0), double %94)
%96 = add i64 %91, 1
%97 = load i64, i64* %1
%98 = icmp ugt i64 %96, %97
br i1 %98, label %99, label %90
99: ; preds = %90, %85
%100 = call i32 @putchar(i32 10)
%101 = call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([18 x i8], [18 x i8]* @str.13, i64 0, i64 0))
%102 = load i64, i64* %1
%103 = icmp eq i64 %102, 0
br i1 %103, label %162, label %104
104: ; preds = %99
%105 = getelementptr inbounds float, float* %7, i64 %102
%106 = add i64 %102, -1
%107 = getelementptr inbounds float, float* %9, i64 %106
%108 = getelementptr inbounds float, float* %5, i64 %102
%109 = getelementptr inbounds float, float* %9, i64 %102
%110 = load float, float* %18
%111 = load float, float* %19
%112 = load float, float* %20
%113 = load float, float* %21
%114 = fmul float %112, %113
%115 = call float @llvm.fmuladd.f32(float %110, float %111, float %114)
store float %115, float* %22
%116 = icmp eq i64 %102, 1
br i1 %116, label %120, label %117
117: ; preds = %104
%118 = getelementptr inbounds float, float* %8, i64 %102
br label %121
119: ; preds = %150
br i1 %103, label %162, label %120
120: ; preds = %104, %119
br label %153
121: ; preds = %150, %117
%122 = phi i64 [ 2, %117 ], [ %151, %150 ]
%123 = icmp eq i64 %122, %102
br i1 %123, label %124, label %131
124: ; preds = %121
%125 = load float, float* %105
%126 = load float, float* %107
%127 = load float, float* %108
%128 = load float, float* %109
%129 = fmul float %127, %128
%130 = call float @llvm.fmuladd.f32(float %125, float %126, float %129)
store float %130, float* %118
br label %150
131: ; preds = %121
%132 = getelementptr inbounds float, float* %7, i64 %122
%133 = load float, float* %132
%134 = add i64 %122, -1
%135 = getelementptr inbounds float, float* %9, i64 %134
%136 = load float, float* %135
%137 = getelementptr inbounds float, float* %5, i64 %122
%138 = load float, float* %137
%139 = getelementptr inbounds float, float* %9, i64 %122
%140 = load float, float* %139
%141 = fmul float %138, %140
%142 = call float @llvm.fmuladd.f32(float %133, float %136, float %141)
%143 = getelementptr inbounds float, float* %6, i64 %122
%144 = load float, float* %143
%145 = add nuw i64 %122, 1
%146 = getelementptr inbounds float, float* %9, i64 %145
%147 = load float, float* %146
%148 = call float @llvm.fmuladd.f32(float %144, float %147, float %142)
%149 = getelementptr inbounds float, float* %8, i64 %122
store float %148, float* %149
br label %150
150: ; preds = %131, %124
%151 = add nuw i64 %122, 1
%152 = icmp eq i64 %122, %102
br i1 %152, label %119, label %121
153: ; preds = %120, %153
%154 = phi i64 [ %159, %153 ], [ 1, %120 ]
%155 = getelementptr inbounds float, float* %8, i64 %154
%156 = load float, float* %155
%157 = fpext float %156 to double
%158 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([7 x i8], [7 x i8]* @.str.6, i64 0, i64 0), double %157)
%159 = add i64 %154, 1
%160 = load i64, i64* %1
%161 = icmp ugt i64 %159, %160
br i1 %161, label %162, label %153
162: ; preds = %153, %99, %119
%163 = call i32 @putchar(i32 10)
%164 = call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([36 x i8], [36 x i8]* @str.14, i64 0, i64 0))
%165 = call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([31 x i8], [31 x i8]* @str.15, i64 0, i64 0))
%166 = call i32 (...) @getchar()
%167 = call i32 @feof(i32* %10)
%168 = icmp eq i32 %167, 0
br i1 %168, label %23, label %169
169: ; preds = %162, %14
%170 = call i32 @fclose(i32* %10)
%171 = call i32 @free_vector(float* %9, i32 1, i32 20)
%172 = call i32 @free_vector(float* %8, i32 1, i32 20)
%173 = call i32 @free_vector(float* %7, i32 1, i32 20)
%174 = call i32 @free_vector(float* %6, i32 1, i32 20)
%175 = call i32 @free_vector(float* %5, i32 1, i32 20)
ret i32 0
}