#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_matrix(float **, int, int, int, int);
int /*<<< orphan*/ free_vector(float *, int, int);
int /*<<< orphan*/ fscanf(int /*<<< orphan*/ *, char *, ...);
int /*<<< orphan*/ getchar();
float **matrix(int, int, int, int);
int /*<<< orphan*/ nrerror(char *);
int /*<<< orphan*/ printf(char *, ...);
int /*<<< orphan*/ svdcmp(float **, int, int, float *, float **);
float *vector(int, int);
int main(void) {
  int j, k, l, m, n;
  float *w, **a, **u, **v;
  char dummy[80];
  FILE *fp;
  w = vector(1, 20);
  a = matrix(1, 20, 1, 20);
  u = matrix(1, 20, 1, 20);
  v = matrix(1, 20, 1, 20);
  if ((fp = fopen("matrx3.dat", "r")) == NULL)
    nrerror("Data file matrx3.dat not found\n");
  while (!feof(fp)) {
    fgets(dummy, 80, fp);
    fgets(dummy, 80, fp);
    fscanf(fp, "%d %d ", &m, &n);
    fgets(dummy, 80, fp);
    for (k = 1; k <= m; k++)
      for (l = 1; l <= n; l++) {
        fscanf(fp, "%f ", &a[k][l]);
        u[k][l] = a[k][l];
      }
    svdcmp(u, m, n, w, v);
    printf("Decomposition matrices:\n");
    printf("Matrix u\n");
    for (k = 1; k <= m; k++) {
      for (l = 1; l <= n; l++)
        printf("%12.6f", u[k][l]);
      printf("\n");
    }
    printf("Diagonal of matrix w\n");
    for (k = 1; k <= n; k++)
      printf("%12.6f", w[k]);
    printf("\nMatrix v-transpose\n");
    for (k = 1; k <= n; k++) {
      for (l = 1; l <= n; l++)
        printf("%12.6f", v[l][k]);
      printf("\n");
    }
    printf("\nCheck product against original matrix:\n");
    printf("Original matrix:\n");
    for (k = 1; k <= m; k++) {
      for (l = 1; l <= n; l++)
        printf("%12.6f", a[k][l]);
      printf("\n");
    }
    printf("Product u*w*(v-transpose):\n");
    for (k = 1; k <= m; k++) {
      for (l = 1; l <= n; l++) {
        a[k][l] = 0.0;
        for (j = 1; j <= n; j++)
          a[k][l] += u[k][j] * w[j] * v[l][j];
      }
      for (l = 1; l <= n; l++)
        printf("%12.6f", a[k][l]);
      printf("\n");
    }
    printf("***********************************\n");
    printf("press RETURN for next problem\n");
    (void)getchar();
  }
  fclose(fp);
  free_matrix(v, 1, 20, 1, 20);
  free_matrix(u, 1, 20, 1, 20);
  free_matrix(a, 1, 20, 1, 20);
  free_vector(w, 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   10e0 <vector@plt>
mov    %rax,-0x20(%rbp)
mov    $0x1,%edx
mov    $0x14,%ecx
mov    %edx,%edi
mov    %ecx,%esi
call   10d0 <matrix@plt>
mov    %rax,-0x28(%rbp)
mov    $0x1,%edx
mov    $0x14,%ecx
mov    %edx,%edi
mov    %ecx,%esi
call   10d0 <matrix@plt>
mov    %rax,-0x30(%rbp)
mov    $0x1,%edx
mov    $0x14,%ecx
mov    %edx,%edi
mov    %ecx,%esi
call   10d0 <matrix@plt>
mov    %rax,-0x38(%rbp)
lea    0xdbf(%rip),%rdi
lea    0xdc3(%rip),%rsi
call   10f0 <fopen@plt>
mov    %rax,-0x98(%rbp)
cmp    $0x0,%rax
jne    126a <main+0x9a>
lea    0xda8(%rip),%rdi
call   1040 <nrerror@plt>
jmp    126f <main+0x9f>
mov    -0x98(%rbp),%rdi
call   10c0 <feof@plt>
cmp    $0x0,%eax
setne  %al
xor    $0xff,%al
test   $0x1,%al
jne    1290 <main+0xc0>
jmp    1707 <main+0x537>
lea    -0x90(%rbp),%rdi
mov    -0x98(%rbp),%rdx
mov    $0x50,%esi
call   1090 <fgets@plt>
lea    -0x90(%rbp),%rdi
mov    -0x98(%rbp),%rdx
mov    $0x50,%esi
call   1090 <fgets@plt>
mov    -0x98(%rbp),%rdi
lea    0xd5f(%rip),%rsi
lea    -0x14(%rbp),%rdx
lea    -0x18(%rbp),%rcx
mov    $0x0,%al
call   1080 <fscanf@plt>
lea    -0x90(%rbp),%rdi
mov    -0x98(%rbp),%rdx
mov    $0x50,%esi
call   1090 <fgets@plt>
movl   $0x1,-0xc(%rbp)
mov    -0xc(%rbp),%eax
cmp    -0x14(%rbp),%eax
jg     1392 <main+0x1c2>
movl   $0x1,-0x10(%rbp)
mov    -0x10(%rbp),%eax
cmp    -0x18(%rbp),%eax
jg     137f <main+0x1af>
mov    -0x98(%rbp),%rdi
mov    -0x28(%rbp),%rax
movslq -0xc(%rbp),%rcx
mov    (%rax,%rcx,8),%rdx
movslq -0x10(%rbp),%rax
shl    $0x2,%rax
add    %rax,%rdx
lea    0xcf4(%rip),%rsi
mov    $0x0,%al
call   1080 <fscanf@plt>
mov    -0x28(%rbp),%rax
movslq -0xc(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x10(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x30(%rbp),%rax
movslq -0xc(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x10(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    130f <main+0x13f>
jmp    1384 <main+0x1b4>
mov    -0xc(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xc(%rbp)
jmp    12fc <main+0x12c>
mov    -0x30(%rbp),%rdi
mov    -0x14(%rbp),%esi
mov    -0x18(%rbp),%edx
mov    -0x20(%rbp),%rcx
mov    -0x38(%rbp),%r8
call   10b0 <svdcmp@plt>
lea    0xc88(%rip),%rdi
mov    $0x0,%al
call   1060 <printf@plt>
lea    0xc93(%rip),%rdi
mov    $0x0,%al
call   1060 <printf@plt>
movl   $0x1,-0xc(%rbp)
mov    -0xc(%rbp),%eax
cmp    -0x14(%rbp),%eax
jg     143c <main+0x26c>
movl   $0x1,-0x10(%rbp)
mov    -0x10(%rbp),%eax
cmp    -0x18(%rbp),%eax
jg     1420 <main+0x250>
mov    -0x30(%rbp),%rax
movslq -0xc(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x10(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
lea    0xc50(%rip),%rdi
mov    $0x1,%al
call   1060 <printf@plt>
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    13df <main+0x20f>
lea    0xce0(%rip),%rdi
mov    $0x0,%al
call   1060 <printf@plt>
mov    -0xc(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xc(%rbp)
jmp    13cc <main+0x1fc>
lea    0xc1f(%rip),%rdi
mov    $0x0,%al
call   1060 <printf@plt>
movl   $0x1,-0xc(%rbp)
mov    -0xc(%rbp),%eax
cmp    -0x18(%rbp),%eax
jg     148a <main+0x2ba>
mov    -0x20(%rbp),%rax
movslq -0xc(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
lea    0xbe6(%rip),%rdi
mov    $0x1,%al
call   1060 <printf@plt>
mov    -0xc(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xc(%rbp)
jmp    1451 <main+0x281>
lea    0xbe7(%rip),%rdi
mov    $0x0,%al
call   1060 <printf@plt>
movl   $0x1,-0xc(%rbp)
mov    -0xc(%rbp),%eax
cmp    -0x18(%rbp),%eax
jg     150f <main+0x33f>
movl   $0x1,-0x10(%rbp)
mov    -0x10(%rbp),%eax
cmp    -0x18(%rbp),%eax
jg     14f3 <main+0x323>
mov    -0x38(%rbp),%rax
movslq -0x10(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0xc(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
lea    0xb7d(%rip),%rdi
mov    $0x1,%al
call   1060 <printf@plt>
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    14b2 <main+0x2e2>
lea    0xc0d(%rip),%rdi
mov    $0x0,%al
call   1060 <printf@plt>
mov    -0xc(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xc(%rbp)
jmp    149f <main+0x2cf>
lea    0xb77(%rip),%rdi
mov    $0x0,%al
call   1060 <printf@plt>
lea    0xb92(%rip),%rdi
mov    $0x0,%al
call   1060 <printf@plt>
movl   $0x1,-0xc(%rbp)
mov    -0xc(%rbp),%eax
cmp    -0x14(%rbp),%eax
jg     15a2 <main+0x3d2>
movl   $0x1,-0x10(%rbp)
mov    -0x10(%rbp),%eax
cmp    -0x18(%rbp),%eax
jg     1586 <main+0x3b6>
mov    -0x28(%rbp),%rax
movslq -0xc(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x10(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
lea    0xaea(%rip),%rdi
mov    $0x1,%al
call   1060 <printf@plt>
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    1545 <main+0x375>
lea    0xb7a(%rip),%rdi
mov    $0x0,%al
call   1060 <printf@plt>
mov    -0xc(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xc(%rbp)
jmp    1532 <main+0x362>
lea    0xb1f(%rip),%rdi
mov    $0x0,%al
call   1060 <printf@plt>
movl   $0x1,-0xc(%rbp)
mov    -0xc(%rbp),%eax
cmp    -0x14(%rbp),%eax
jg     16df <main+0x50f>
movl   $0x1,-0x10(%rbp)
mov    -0x10(%rbp),%eax
cmp    -0x18(%rbp),%eax
jg     167b <main+0x4ab>
mov    -0x28(%rbp),%rax
movslq -0xc(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x10(%rbp),%rcx
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax,%rcx,4)
movl   $0x1,-0x8(%rbp)
mov    -0x8(%rbp),%eax
cmp    -0x18(%rbp),%eax
jg     1668 <main+0x498>
mov    -0x30(%rbp),%rax
movslq -0xc(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x8(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x8(%rbp),%rcx
mulss  (%rax,%rcx,4),%xmm0
mov    -0x38(%rbp),%rax
movslq -0x10(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x8(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
mov    -0x28(%rbp),%rax
movslq -0xc(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x10(%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    15f5 <main+0x425>
jmp    166d <main+0x49d>
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    15ca <main+0x3fa>
movl   $0x1,-0x10(%rbp)
mov    -0x10(%rbp),%eax
cmp    -0x18(%rbp),%eax
jg     16c3 <main+0x4f3>
mov    -0x28(%rbp),%rax
movslq -0xc(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x10(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
cvtss2sd %xmm0,%xmm0
lea    0x9ad(%rip),%rdi
mov    $0x1,%al
call   1060 <printf@plt>
mov    -0x10(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x10(%rbp)
jmp    1682 <main+0x4b2>
lea    0xa3d(%rip),%rdi
mov    $0x0,%al
call   1060 <printf@plt>
mov    -0xc(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0xc(%rbp)
jmp    15b7 <main+0x3e7>
lea    0x9fe(%rip),%rdi
mov    $0x0,%al
call   1060 <printf@plt>
lea    0xa15(%rip),%rdi
mov    $0x0,%al
call   1060 <printf@plt>
mov    $0x0,%al
call   10a0 <getchar@plt>
jmp    126f <main+0x9f>
mov    -0x98(%rbp),%rdi
call   1050 <fclose@plt>
mov    -0x38(%rbp),%rdi
mov    $0x1,%ecx
mov    $0x14,%r8d
mov    %ecx,%esi
mov    %r8d,%edx
call   1030 <free_matrix@plt>
mov    -0x30(%rbp),%rdi
mov    $0x1,%ecx
mov    $0x14,%r8d
mov    %ecx,%esi
mov    %r8d,%edx
call   1030 <free_matrix@plt>
mov    -0x28(%rbp),%rdi
mov    $0x1,%ecx
mov    $0x14,%r8d
mov    %ecx,%esi
mov    %r8d,%edx
call   1030 <free_matrix@plt>
mov    -0x20(%rbp),%rdi
mov    $0x1,%esi
mov    $0x14,%edx
call   1070 <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    $0x88,%rsp
mov    $0x1,%edi
mov    $0x14,%esi
call   1100 <vector@plt>
mov    %rax,%r15
mov    $0x1,%edi
mov    $0x14,%esi
mov    $0x1,%edx
mov    $0x14,%ecx
call   10f0 <matrix@plt>
mov    %rax,%r14
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,%rbp
lea    0xd92(%rip),%rdi
lea    0xd96(%rip),%rsi
call   1110 <fopen@plt>
mov    %rax,%rbx
test   %rax,%rax
jne    128e <main+0x9e>
lea    0xd84(%rip),%rdi
call   1060 <nrerror@plt>
mov    %rbx,%rdi
call   10e0 <feof@plt>
test   %eax,%eax
je     1323 <main+0x133>
mov    %rbx,%rdi
call   1070 <fclose@plt>
mov    %rbp,%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    %r14,%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
call   1090 <free_vector@plt>
xor    %eax,%eax
add    $0x88,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
mov    %rbp,0x18(%rsp)
mov    %rbx,0x28(%rsp)
mov    %r12,0x10(%rsp)
mov    %r14,0x20(%rsp)
jmp    1379 <main+0x189>
nopl   0x0(%rax)
lea    0xd96(%rip),%rdi
call   1040 <puts@plt>
lea    0xdae(%rip),%rdi
call   1040 <puts@plt>
xor    %eax,%eax
call   10c0 <getchar@plt>
mov    0x28(%rsp),%rbx
mov    %rbx,%rdi
call   10e0 <feof@plt>
test   %eax,%eax
mov    0x18(%rsp),%rbp
jne    129e <main+0xae>
lea    0x30(%rsp),%rbp
mov    %rbp,%rdi
mov    $0x50,%esi
mov    %rbx,%rdx
call   10b0 <fgets@plt>
mov    %rbp,%rdi
mov    $0x50,%esi
mov    %rbx,%rdx
call   10b0 <fgets@plt>
mov    %rbx,%rdi
lea    0xc85(%rip),%rsi
lea    0xc(%rsp),%rdx
lea    0x8(%rsp),%rcx
xor    %eax,%eax
call   10a0 <fscanf@plt>
mov    %rbp,%rdi
mov    $0x50,%esi
mov    %rbx,%rdx
call   10b0 <fgets@plt>
mov    0xc(%rsp),%esi
mov    0x8(%rsp),%edx
test   %esi,%esi
jle    1480 <main+0x290>
mov    %rbx,%rdi
mov    $0x1,%r13d
jmp    1403 <main+0x213>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
lea    0x1(%r13),%rax
movslq %esi,%rcx
cmp    %rcx,%r13
mov    %rax,%r13
jge    1480 <main+0x290>
test   %edx,%edx
jle    13f0 <main+0x200>
mov    (%r14,%r13,8),%rax
mov    $0x4,%ebp
xor    %r12d,%r12d
mov    %r14,%rbx
cs nopw 0x0(%rax,%rax,1)
add    %rbp,%rax
mov    %rdi,%r14
lea    0xc07(%rip),%rsi
mov    %rax,%rdx
xor    %eax,%eax
call   10a0 <fscanf@plt>
mov    %r14,%rdi
mov    (%rbx,%r13,8),%rax
movss  0x4(%rax,%r12,4),%xmm0
mov    0x10(%rsp),%rcx
mov    (%rcx,%r13,8),%rcx
movss  %xmm0,0x4(%rcx,%r12,4)
movslq 0x8(%rsp),%rdx
add    $0x1,%r12
add    $0x4,%rbp
cmp    %rdx,%r12
jl     1420 <main+0x230>
mov    0xc(%rsp),%esi
mov    %rbx,%r14
mov    0x10(%rsp),%r12
jmp    13f0 <main+0x200>
nopl   0x0(%rax,%rax,1)
mov    %r12,%rdi
mov    %r15,%rcx
mov    0x18(%rsp),%r13
mov    %r13,%r8
call   10d0 <svdcmp@plt>
lea    0xba5(%rip),%rdi
call   1040 <puts@plt>
lea    0xbb1(%rip),%rdi
call   1040 <puts@plt>
cmpl   $0x0,0xc(%rsp)
lea    0xb81(%rip),%rbx
jle    1530 <main+0x340>
mov    $0x1,%r12d
jmp    14ec <main+0x2fc>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    $0xa,%edi
call   1030 <putchar@plt>
lea    0x1(%r12),%rax
movslq 0xc(%rsp),%rcx
cmp    %rcx,%r12
mov    %rax,%r12
jge    1530 <main+0x340>
cmpl   $0x0,0x8(%rsp)
mov    0x10(%rsp),%r14
jle    14d0 <main+0x2e0>
xor    %ebp,%ebp
nopw   0x0(%rax,%rax,1)
mov    (%r14,%r12,8),%rax
movss  0x4(%rax,%rbp,4),%xmm0
cvtss2sd %xmm0,%xmm0
mov    %rbx,%rdi
mov    $0x1,%al
call   1080 <printf@plt>
movslq 0x8(%rsp),%rax
add    $0x1,%rbp
cmp    %rax,%rbp
jl     1500 <main+0x310>
jmp    14d0 <main+0x2e0>
nopl   0x0(%rax,%rax,1)
lea    0xb29(%rip),%rdi
call   1040 <puts@plt>
cmpl   $0x0,0x8(%rsp)
jle    1573 <main+0x383>
xor    %ebp,%ebp
cs nopw 0x0(%rax,%rax,1)
nop
movss  0x4(%r15,%rbp,4),%xmm0
cvtss2sd %xmm0,%xmm0
mov    %rbx,%rdi
mov    $0x1,%al
call   1080 <printf@plt>
movslq 0x8(%rsp),%rax
add    $0x1,%rbp
cmp    %rax,%rbp
jl     1550 <main+0x360>
lea    0xafb(%rip),%rdi
call   1040 <puts@plt>
mov    0x8(%rsp),%eax
test   %eax,%eax
mov    0x20(%rsp),%r14
jle    1600 <main+0x410>
mov    $0x1,%r12d
jmp    15bc <main+0x3cc>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
mov    $0xa,%edi
call   1030 <putchar@plt>
lea    0x1(%r12),%rcx
movslq 0x8(%rsp),%rax
cmp    %rax,%r12
mov    %rcx,%r12
jge    1600 <main+0x410>
test   %eax,%eax
jle    15a0 <main+0x3b0>
xor    %ebp,%ebp
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
mov    0x8(%r13,%rbp,8),%rax
movss  (%rax,%r12,4),%xmm0
cvtss2sd %xmm0,%xmm0
mov    %rbx,%rdi
mov    $0x1,%al
call   1080 <printf@plt>
movslq 0x8(%rsp),%rax
add    $0x1,%rbp
cmp    %rax,%rbp
jl     15d0 <main+0x3e0>
jmp    15a0 <main+0x3b0>
nopl   0x0(%rax)
lea    0xa82(%rip),%rdi
call   1040 <puts@plt>
lea    0xa9e(%rip),%rdi
call   1040 <puts@plt>
cmpl   $0x0,0xc(%rsp)
jle    1690 <main+0x4a0>
mov    $0x1,%r12d
jmp    164c <main+0x45c>
nopw   0x0(%rax,%rax,1)
mov    $0xa,%edi
call   1030 <putchar@plt>
lea    0x1(%r12),%rax
movslq 0xc(%rsp),%rcx
cmp    %rcx,%r12
mov    %rax,%r12
jge    1690 <main+0x4a0>
cmpl   $0x0,0x8(%rsp)
jle    1630 <main+0x440>
xor    %ebp,%ebp
cs nopw 0x0(%rax,%rax,1)
nop
mov    (%r14,%r12,8),%rax
movss  0x4(%rax,%rbp,4),%xmm0
cvtss2sd %xmm0,%xmm0
mov    %rbx,%rdi
mov    $0x1,%al
call   1080 <printf@plt>
movslq 0x8(%rsp),%rax
add    $0x1,%rbp
cmp    %rax,%rbp
jl     1660 <main+0x470>
jmp    1630 <main+0x440>
nopl   0x0(%rax,%rax,1)
lea    0xa2b(%rip),%rdi
call   1040 <puts@plt>
cmpl   $0x0,0xc(%rsp)
mov    0x10(%rsp),%r12
jle    1340 <main+0x150>
mov    $0x1,%r13d
jmp    16df <main+0x4ef>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
mov    $0xa,%edi
call   1030 <putchar@plt>
lea    0x1(%r13),%rax
movslq 0xc(%rsp),%rcx
cmp    %rcx,%r13
mov    %rax,%r13
jge    1340 <main+0x150>
mov    0x8(%rsp),%r8d
test   %r8d,%r8d
jle    16c0 <main+0x4d0>
mov    (%r14,%r13,8),%rcx
lea    0x1(%r8),%r9
mov    %r12,%r10
mov    (%r12,%r13,8),%rsi
mov    %r8d,%edi
and    $0xfffffffe,%edi
neg    %rdi
mov    $0x1,%ebp
jmp    171d <main+0x52d>
nopl   0x0(%rax,%rax,1)
add    $0x1,%rbp
cmp    %r9,%rbp
je     17c0 <main+0x5d0>
movl   $0x0,(%rcx,%rbp,4)
mov    0x18(%rsp),%rax
mov    (%rax,%rbp,8),%rax
xorps  %xmm0,%xmm0
mov    $0x1,%edx
cmp    $0x2,%r9d
je     1797 <main+0x5a7>
mov    $0x1,%edx
xorps  %xmm0,%xmm0
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movss  (%rsi,%rdx,4),%xmm1
mulss  (%r15,%rdx,4),%xmm1
mulss  (%rax,%rdx,4),%xmm1
addss  %xmm0,%xmm1
movss  %xmm1,(%rcx,%rbp,4)
movss  0x4(%rsi,%rdx,4),%xmm0
mulss  0x4(%r15,%rdx,4),%xmm0
mulss  0x4(%rax,%rdx,4),%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,(%rcx,%rbp,4)
lea    (%rdi,%rdx,1),%rbx
add    $0x2,%rbx
add    $0x2,%rdx
cmp    $0x1,%rbx
jne    1750 <main+0x560>
test   $0x1,%r8b
je     1710 <main+0x520>
movss  (%rsi,%rdx,4),%xmm1
mulss  (%r15,%rdx,4),%xmm1
mulss  (%rax,%rdx,4),%xmm1
addss  %xmm0,%xmm1
movss  %xmm1,(%rcx,%rbp,4)
jmp    1710 <main+0x520>
nop
test   %r8d,%r8d
mov    %r10,%r12
lea    0x86b(%rip),%rbx
jle    16c0 <main+0x4d0>
xor    %ebp,%ebp
cs nopw 0x0(%rax,%rax,1)
nop
mov    (%r14,%r13,8),%rax
movss  0x4(%rax,%rbp,4),%xmm0
cvtss2sd %xmm0,%xmm0
mov    %rbx,%rdi
mov    $0x1,%al
call   1080 <printf@plt>
movslq 0x8(%rsp),%rax
add    $0x1,%rbp
cmp    %rax,%rbp
jl     17e0 <main+0x5f0>
jmp    16c0 <main+0x4d0>

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

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