#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 */
struct timeval {
  double tv_sec;
  double tv_usec;
};
typedef int /*<<< orphan*/ FILE;

/* Variables and functions */
void *calloc(int, int);
int /*<<< orphan*/ dfeast_scsrgv(char *, int *, double *, int *, int *,
                                 double *, int *, int *, int *, double *, int *,
                                 double *, double *, int *, double *, double *,
                                 int *, double *, int *);
int /*<<< orphan*/ fclose(int /*<<< orphan*/ *);
int /*<<< orphan*/ feastinit(int *);
int /*<<< orphan*/ *fopen(char *, char *);
int fscanf(int /*<<< orphan*/ *, char *, int *, int *, ...);
int /*<<< orphan*/ gettimeofday(struct timeval *, int /*<<< orphan*/ *);
int /*<<< orphan*/ printf(char *, ...);
int main() {
  int feastparam[64];
  double epsout;
  int loop;
  char UPLO = 'F';
  FILE *fp;
  char name[] = "../../system1";
  int N, nnz;
  double *sa, *sb;
  int *isa, *jsa;
  struct timeval t1, t2;
  int i, k, err;
  int M0, M, info;
  double Emin, Emax, trace;
  double *X;
  double *E, *res;
  fp = fopen(name, "r");
  err = fscanf(fp, "%d%d\n", &N, &nnz);
  sa = calloc(nnz, sizeof(double));
  sb = calloc(nnz, sizeof(double));
  isa = calloc(N + 1, sizeof(int));
  jsa = calloc(nnz, sizeof(int));
  for (i = 0; i <= N; i++) {
    *(isa + i) = 0;
  };
  *(isa) = 1;
  for (k = 0; k <= nnz - 1; k++) {
    err = fscanf(fp, "%d%d%lf%lf\n", &i, jsa + k, sa + k, sb + k);
    *(isa + i) = *(isa + i) + 1;
  };
  fclose(fp);
  for (i = 1; i <= N; i++) {
    *(isa + i) = *(isa + i) + *(isa + i - 1);
  };
  printf("sparse matrix -system1- size %.d\n", N);
  printf("nnz %d \n", nnz);
  gettimeofday(&t1, NULL);
  Emin = (double)0.18;
  Emax = (double)1.0;
  M0 = 25;
  E = calloc(M0, sizeof(double));
  res = calloc(M0, sizeof(double));
  X = calloc(N * M0, sizeof(double));
  feastinit(feastparam);
  feastparam[0] = 1;
  dfeast_scsrgv(&UPLO, &N, sa, isa, jsa, sb, isa, jsa, feastparam, &epsout,
                &loop, &Emin, &Emax, &M0, E, X, &M, res, &info);
  gettimeofday(&t2, NULL);
  printf("FEAST OUTPUT INFO %d\n", info);
  if (info == 0) {
    printf("*************************************************\n");
    printf("************** REPORT ***************************\n");
    printf("*************************************************\n");
    printf("SIMULATION TIME %f\n", (t2.tv_sec - t1.tv_sec) * 1.0 +
                                       (t2.tv_usec - t1.tv_usec) * 0.000001);
    printf("# Search interval [Emin,Emax] %.15e %.15e\n", Emin, Emax);
    printf("# mode found/subspace %d %d \n", M, M0);
    printf("# iterations %d \n", loop);
    trace = (double)0.0;
    for (i = 0; i <= M - 1; i = i + 1) {
      trace = trace + *(E + i);
    }
    printf("TRACE %.15e\n", trace);
    printf("Relative error on the Trace %.15e\n", epsout);
    printf("Eigenvalues/Residuals\n");
    for (i = 0; i <= M - 1; i = i + 1) {
      printf("   %d %.15e %.15e\n", i, *(E + i), *(res + i));
    }
  }
  return 0;
}


x86 assembly (clang, -O0):
<main>:
push   %rbp
mov    %rsp,%rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x258,%rsp
movl   $0x0,-0x2c(%rbp)
movb   $0x46,-0x13d(%rbp)
mov    0xe6f(%rip),%rax
mov    %rax,-0x156(%rbp)
mov    0xe6a(%rip),%eax
mov    %eax,-0x14e(%rbp)
mov    0xe61(%rip),%ax
mov    %ax,-0x14a(%rbp)
lea    -0x156(%rbp),%rdi
lea    0xe4e(%rip),%rsi
call   10a0 <fopen@plt>
mov    %rax,-0x148(%rbp)
mov    -0x148(%rbp),%rdi
lea    0xe36(%rip),%rsi
lea    -0x15c(%rbp),%rdx
lea    -0x160(%rbp),%rcx
mov    $0x0,%al
call   1070 <fscanf@plt>
mov    %eax,-0x1ac(%rbp)
mov    -0x160(%rbp),%edi
mov    $0x8,%esi
call   1080 <calloc@plt>
mov    %rax,-0x168(%rbp)
mov    -0x160(%rbp),%edi
mov    $0x8,%esi
call   1080 <calloc@plt>
mov    %rax,-0x170(%rbp)
mov    -0x15c(%rbp),%edi
add    $0x1,%edi
mov    $0x4,%esi
call   1080 <calloc@plt>
mov    %rax,-0x178(%rbp)
mov    -0x160(%rbp),%edi
mov    $0x4,%esi
call   1080 <calloc@plt>
mov    %rax,-0x180(%rbp)
movl   $0x0,-0x1a4(%rbp)
mov    -0x1a4(%rbp),%eax
cmp    -0x15c(%rbp),%eax
jg     12b1 <main+0x131>
mov    -0x178(%rbp),%rax
movslq -0x1a4(%rbp),%rcx
movl   $0x0,(%rax,%rcx,4)
mov    -0x1a4(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1a4(%rbp)
jmp    1276 <main+0xf6>
mov    -0x178(%rbp),%rax
movl   $0x1,(%rax)
movl   $0x0,-0x1a8(%rbp)
mov    -0x1a8(%rbp),%eax
mov    -0x160(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jg     1379 <main+0x1f9>
mov    -0x148(%rbp),%rdi
mov    -0x180(%rbp),%rcx
movslq -0x1a8(%rbp),%rax
shl    $0x2,%rax
add    %rax,%rcx
mov    -0x168(%rbp),%r8
movslq -0x1a8(%rbp),%rax
shl    $0x3,%rax
add    %rax,%r8
mov    -0x170(%rbp),%r9
movslq -0x1a8(%rbp),%rax
shl    $0x3,%rax
add    %rax,%r9
lea    0xd02(%rip),%rsi
lea    -0x1a4(%rbp),%rdx
mov    $0x0,%al
call   1070 <fscanf@plt>
mov    %eax,-0x1ac(%rbp)
mov    -0x178(%rbp),%rax
movslq -0x1a4(%rbp),%rcx
mov    (%rax,%rcx,4),%edx
add    $0x1,%edx
mov    -0x178(%rbp),%rax
movslq -0x1a4(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x1a8(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1a8(%rbp)
jmp    12c8 <main+0x148>
mov    -0x148(%rbp),%rdi
call   1030 <fclose@plt>
movl   $0x1,-0x1a4(%rbp)
mov    -0x1a4(%rbp),%eax
cmp    -0x15c(%rbp),%eax
jg     13e9 <main+0x269>
mov    -0x178(%rbp),%rax
movslq -0x1a4(%rbp),%rcx
mov    (%rax,%rcx,4),%edx
mov    -0x178(%rbp),%rax
movslq -0x1a4(%rbp),%rcx
add    -0x4(%rax,%rcx,4),%edx
mov    -0x178(%rbp),%rax
movslq -0x1a4(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x1a4(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1a4(%rbp)
jmp    138f <main+0x20f>
mov    -0x15c(%rbp),%esi
lea    0xc44(%rip),%rdi
mov    $0x0,%al
call   1040 <printf@plt>
mov    -0x160(%rbp),%esi
lea    0xc52(%rip),%rdi
mov    $0x0,%al
call   1040 <printf@plt>
lea    -0x190(%rbp),%rdi
xor    %eax,%eax
mov    %eax,%esi
call   1060 <gettimeofday@plt>
movsd  0xbdf(%rip),%xmm0
movsd  %xmm0,-0x1c0(%rbp)
ff
movsd  0xbc7(%rip),%xmm0
movsd  %xmm0,-0x1c8(%rbp)
ff
movl   $0x19,-0x1b0(%rbp)
mov    -0x1b0(%rbp),%edi
mov    $0x8,%esi
call   1080 <calloc@plt>
mov    %rax,-0x1e0(%rbp)
mov    -0x1b0(%rbp),%edi
mov    $0x8,%esi
call   1080 <calloc@plt>
mov    %rax,-0x1e8(%rbp)
mov    -0x15c(%rbp),%edi
imul   -0x1b0(%rbp),%edi
mov    $0x8,%esi
call   1080 <calloc@plt>
mov    %rax,-0x1d8(%rbp)
lea    -0x130(%rbp),%rdi
call   1050 <feastinit@plt>
movl   $0x1,-0x130(%rbp)
mov    -0x168(%rbp),%rdx
mov    -0x178(%rbp),%rcx
mov    -0x180(%rbp),%r8
mov    -0x170(%rbp),%r9
mov    -0x178(%rbp),%rax
mov    %rax,-0x218(%rbp)
mov    -0x180(%rbp),%rax
mov    %rax,-0x210(%rbp)
lea    -0x130(%rbp),%rax
mov    %rax,-0x208(%rbp)
mov    -0x1e0(%rbp),%r14
mov    -0x1d8(%rbp),%rbx
mov    -0x1e8(%rbp),%r10
lea    -0x13d(%rbp),%rdi
lea    -0x15c(%rbp),%rsi
lea    -0x138(%rbp),%rax
mov    %rax,-0x200(%rbp)
lea    -0x13c(%rbp),%rax
mov    %rax,-0x1f8(%rbp)
lea    -0x1c0(%rbp),%r13
lea    -0x1c8(%rbp),%r12
lea    -0x1b0(%rbp),%r15
lea    -0x1b4(%rbp),%r11
lea    -0x1b8(%rbp),%rax
mov    %rax,-0x1f0(%rbp)
mov    -0x218(%rbp),%rax
mov    %rax,(%rsp)
mov    -0x210(%rbp),%rax
mov    %rax,0x8(%rsp)
mov    -0x208(%rbp),%rax
mov    %rax,0x10(%rsp)
mov    -0x200(%rbp),%rax
mov    %rax,0x18(%rsp)
mov    -0x1f8(%rbp),%rax
mov    %rax,0x20(%rsp)
mov    -0x1f0(%rbp),%rax
mov    %r13,0x28(%rsp)
mov    %r12,0x30(%rsp)
mov    %r15,0x38(%rsp)
mov    %r14,0x40(%rsp)
mov    %rbx,0x48(%rsp)
mov    %r11,0x50(%rsp)
mov    %r10,0x58(%rsp)
mov    %rax,0x60(%rsp)
call   1090 <dfeast_scsrgv@plt>
lea    -0x1a0(%rbp),%rdi
xor    %eax,%eax
mov    %eax,%esi
call   1060 <gettimeofday@plt>
mov    -0x1b8(%rbp),%esi
lea    0xa7d(%rip),%rdi
mov    $0x0,%al
call   1040 <printf@plt>
cmpl   $0x0,-0x1b8(%rbp)
jne    17c1 <main+0x641>
lea    0xa78(%rip),%rdi
mov    $0x0,%al
call   1040 <printf@plt>
lea    0xa9d(%rip),%rdi
mov    $0x0,%al
call   1040 <printf@plt>
lea    0xa5c(%rip),%rdi
mov    $0x0,%al
call   1040 <printf@plt>
movsd  -0x1a0(%rbp),%xmm0
ff
subsd  -0x190(%rbp),%xmm0
ff
movsd  -0x198(%rbp),%xmm1
ff
subsd  -0x188(%rbp),%xmm1
ff
movsd  0x9c2(%rip),%xmm2
mulsd  %xmm2,%xmm1
addsd  %xmm1,%xmm0
lea    0xa84(%rip),%rdi
mov    $0x1,%al
call   1040 <printf@plt>
movsd  -0x1c0(%rbp),%xmm0
ff
movsd  -0x1c8(%rbp),%xmm1
ff
lea    0xa7a(%rip),%rdi
mov    $0x2,%al
call   1040 <printf@plt>
mov    -0x1b4(%rbp),%esi
mov    -0x1b0(%rbp),%edx
lea    0xa8b(%rip),%rdi
mov    $0x0,%al
call   1040 <printf@plt>
mov    -0x13c(%rbp),%esi
lea    0xa95(%rip),%rdi
mov    $0x0,%al
call   1040 <printf@plt>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x1d0(%rbp)
ff
movl   $0x0,-0x1a4(%rbp)
mov    -0x1a4(%rbp),%eax
mov    -0x1b4(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jg     1713 <main+0x593>
movsd  -0x1d0(%rbp),%xmm0
ff
mov    -0x1e0(%rbp),%rax
movslq -0x1a4(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x1d0(%rbp)
ff
mov    -0x1a4(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1a4(%rbp)
jmp    16c5 <main+0x545>
movsd  -0x1d0(%rbp),%xmm0
ff
lea    0xa2e(%rip),%rdi
mov    $0x1,%al
call   1040 <printf@plt>
movsd  -0x138(%rbp),%xmm0
ff
lea    0xa25(%rip),%rdi
mov    $0x1,%al
call   1040 <printf@plt>
lea    0xa3a(%rip),%rdi
mov    $0x0,%al
call   1040 <printf@plt>
movl   $0x0,-0x1a4(%rbp)
mov    -0x1a4(%rbp),%eax
mov    -0x1b4(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jg     17bc <main+0x63c>
mov    -0x1a4(%rbp),%esi
mov    -0x1e0(%rbp),%rax
movslq -0x1a4(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x1e8(%rbp),%rax
movslq -0x1a4(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
lea    0x9f6(%rip),%rdi
mov    $0x2,%al
call   1040 <printf@plt>
mov    -0x1a4(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1a4(%rbp)
jmp    1757 <main+0x5d7>
jmp    17c1 <main+0x641>
xor    %eax,%eax
add    $0x258,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

x86 assembly (clang, -O2):
<main>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x188,%rsp
movb   $0x46,0x7(%rsp)
movabs $0x316d6574737973,%rax
mov    %rax,0x56(%rsp)
movabs $0x79732f2e2e2f2e2e,%rax
2f 73 79
mov    %rax,0x50(%rsp)
lea    0xe3b(%rip),%rsi
lea    0x50(%rsp),%rdi
call   10c0 <fopen@plt>
lea    0xe2c(%rip),%rsi
xor    %r15d,%r15d
lea    0xc(%rsp),%rdx
lea    0x10(%rsp),%rcx
mov    %rax,%r12
mov    %rax,%rdi
xor    %eax,%eax
call   1090 <fscanf@plt>
mov    0x10(%rsp),%r14d
mov    %r14d,%edi
mov    $0x8,%esi
call   10a0 <calloc@plt>
mov    %rax,0x28(%rsp)
mov    %r14d,%edi
mov    $0x8,%esi
call   10a0 <calloc@plt>
mov    %rax,0x20(%rsp)
mov    0xc(%rsp),%ebp
lea    0x1(%rbp),%edi
mov    $0x4,%esi
call   10a0 <calloc@plt>
mov    %rax,%rbx
mov    %r14d,%edi
mov    $0x4,%esi
call   10a0 <calloc@plt>
mov    %rax,0x30(%rsp)
test   %ebp,%ebp
js     126f <main+0xcf>
lea    0x4(,%rbp,4),%rdx
mov    %rbx,%rdi
xor    %esi,%esi
call   1080 <memset@plt>
lea    0x1(%rbp),%r15d
mov    %r15d,(%rsp)
movl   $0x1,(%rbx)
test   %r14d,%r14d
mov    %r12,%rbp
jle    12df <main+0x13f>
mov    0x30(%rsp),%r12
mov    0x28(%rsp),%r13
mov    0x20(%rsp),%r15
xor    %r14d,%r14d
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
mov    %rbp,%rdi
lea    0xd74(%rip),%rsi
mov    %rsp,%rdx
mov    %r12,%rcx
mov    %r13,%r8
mov    %r15,%r9
xor    %eax,%eax
call   1090 <fscanf@plt>
movslq (%rsp),%rax
addl   $0x1,(%rbx,%rax,4)
add    $0x1,%r14
movslq 0x10(%rsp),%rax
add    $0x8,%r15
add    $0x8,%r13
add    $0x4,%r12
cmp    %rax,%r14
jl     12a0 <main+0x100>
mov    %rbp,%rdi
call   1040 <fclose@plt>
mov    0xc(%rsp),%esi
test   %esi,%esi
jle    1361 <main+0x1c1>
mov    (%rbx),%eax
lea    -0x1(%rsi),%rdi
mov    %esi,%ecx
and    $0x3,%ecx
mov    $0x1,%edx
cmp    $0x3,%rdi
jb     133d <main+0x19d>
mov    %esi,%edi
and    $0xfffffffc,%edi
xor    %edx,%edx
nopl   0x0(%rax)
add    0x4(%rbx,%rdx,4),%eax
mov    %eax,0x4(%rbx,%rdx,4)
add    0x8(%rbx,%rdx,4),%eax
mov    %eax,0x8(%rbx,%rdx,4)
add    0xc(%rbx,%rdx,4),%eax
mov    %eax,0xc(%rbx,%rdx,4)
add    0x10(%rbx,%rdx,4),%eax
mov    %eax,0x10(%rbx,%rdx,4)
add    $0x4,%rdx
cmp    %rdx,%rdi
jne    1310 <main+0x170>
add    $0x1,%rdx
lea    0x1(%rsi),%edi
test   %rcx,%rcx
je     1366 <main+0x1c6>
lea    (%rbx,%rdx,4),%rdx
xor    %ebp,%ebp
nopl   0x0(%rax,%rax,1)
add    (%rdx,%rbp,4),%eax
mov    %eax,(%rdx,%rbp,4)
add    $0x1,%rbp
cmp    %rbp,%rcx
jne    1350 <main+0x1b0>
jmp    1366 <main+0x1c6>
mov    $0x1,%edi
mov    %edi,(%rsp)
lea    0xcba(%rip),%rdi
xor    %eax,%eax
call   1050 <printf@plt>
mov    0x10(%rsp),%esi
lea    0xcca(%rip),%rdi
xor    %eax,%eax
call   1050 <printf@plt>
lea    0x70(%rsp),%rdi
xor    %esi,%esi
call   1070 <gettimeofday@plt>
movabs $0x3fc70a3d70a3d70a,%rax
0a c7 3f
mov    %rax,0x40(%rsp)
movabs $0x3ff0000000000000,%rax
mov    %rax,0x38(%rsp)
movl   $0x19,0x18(%rsp)
mov    $0x19,%edi
mov    $0x8,%esi
call   10a0 <calloc@plt>
mov    %rax,%r15
mov    $0x19,%edi
mov    $0x8,%esi
call   10a0 <calloc@plt>
mov    %rax,%r14
mov    0xc(%rsp),%eax
lea    (%rax,%rax,4),%eax
lea    (%rax,%rax,4),%edi
mov    $0x8,%esi
call   10a0 <calloc@plt>
mov    %rax,%rbp
lea    0x80(%rsp),%r12
mov    %r12,%rdi
call   1060 <feastinit@plt>
movl   $0x1,0x80(%rsp)
sub    $0x8,%rsp
lea    0x1c(%rsp),%rax
lea    0x10(%rsp),%r10
lea    0x20(%rsp),%r11
lea    0x40(%rsp),%r13
lea    0xf(%rsp),%rdi
lea    0x14(%rsp),%rsi
mov    0x30(%rsp),%rdx
mov    %rbx,%rcx
mov    0x38(%rsp),%r8
mov    0x28(%rsp),%r9
push   %rax
push   %r14
push   %r10
push   %rbp
push   %r15
push   %r11
push   %r13
lea    0x80(%rsp),%rax
push   %rax
lea    0x64(%rsp),%rax
push   %rax
lea    0x98(%rsp),%rax
push   %rax
push   %r12
push   %r8
push   %rbx
call   10b0 <dfeast_scsrgv@plt>
add    $0x70,%rsp
lea    0x60(%rsp),%rdi
xor    %esi,%esi
call   1070 <gettimeofday@plt>
mov    0x14(%rsp),%esi
lea    0xbc7(%rip),%rdi
xor    %eax,%eax
call   1050 <printf@plt>
cmpl   $0x0,0x14(%rsp)
jne    1648 <main+0x4a8>
lea    0xca8(%rip),%rbx
mov    %rbx,%rdi
call   1030 <puts@plt>
lea    0xc67(%rip),%rdi
call   1030 <puts@plt>
mov    %rbx,%rdi
call   1030 <puts@plt>
movapd 0x60(%rsp),%xmm1
subpd  0x70(%rsp),%xmm1
movapd %xmm1,%xmm0
unpckhpd %xmm1,%xmm0
mulsd  0xb21(%rip),%xmm0
addsd  %xmm1,%xmm0
lea    0xb81(%rip),%rdi
mov    $0x1,%al
call   1050 <printf@plt>
movsd  0x40(%rsp),%xmm0
movsd  0x38(%rsp),%xmm1
lea    0xb7b(%rip),%rdi
mov    $0x2,%al
call   1050 <printf@plt>
mov    0x8(%rsp),%esi
mov    0x18(%rsp),%edx
lea    0xb90(%rip),%rdi
xor    %eax,%eax
call   1050 <printf@plt>
mov    0x1c(%rsp),%esi
lea    0xb9c(%rip),%rdi
xor    %eax,%eax
call   1050 <printf@plt>
mov    0x8(%rsp),%eax
test   %eax,%eax
jle    1552 <main+0x3b2>
lea    -0x1(%rax),%rdx
mov    %eax,%ecx
and    $0x7,%ecx
cmp    $0x7,%rdx
jae    155a <main+0x3ba>
xorpd  %xmm0,%xmm0
xor    %edx,%edx
jmp    15b0 <main+0x410>
xorpd  %xmm0,%xmm0
xor    %eax,%eax
jmp    15ce <main+0x42e>
mov    %eax,%esi
and    $0xfffffff8,%esi
xorpd  %xmm0,%xmm0
xor    %edx,%edx
cs nopw 0x0(%rax,%rax,1)
nop
addsd  (%r15,%rdx,8),%xmm0
addsd  0x8(%r15,%rdx,8),%xmm0
addsd  0x10(%r15,%rdx,8),%xmm0
addsd  0x18(%r15,%rdx,8),%xmm0
addsd  0x20(%r15,%rdx,8),%xmm0
addsd  0x28(%r15,%rdx,8),%xmm0
addsd  0x30(%r15,%rdx,8),%xmm0
addsd  0x38(%r15,%rdx,8),%xmm0
add    $0x8,%rdx
cmp    %rdx,%rsi
jne    1570 <main+0x3d0>
test   %rcx,%rcx
je     15ce <main+0x42e>
lea    (%r15,%rdx,8),%rdx
xor    %esi,%esi
nopl   0x0(%rax,%rax,1)
addsd  (%rdx,%rsi,8),%xmm0
add    $0x1,%rsi
cmp    %rsi,%rcx
jne    15c0 <main+0x420>
mov    %eax,(%rsp)
lea    0xb02(%rip),%rdi
mov    $0x1,%al
call   1050 <printf@plt>
movsd  0x48(%rsp),%xmm0
lea    0xafb(%rip),%rdi
mov    $0x1,%al
call   1050 <printf@plt>
lea    0xb87(%rip),%rdi
call   1030 <puts@plt>
movl   $0x0,(%rsp)
cmpl   $0x0,0x8(%rsp)
jle    1648 <main+0x4a8>
xor    %eax,%eax
lea    0xaf4(%rip),%rbx
cs nopw 0x0(%rax,%rax,1)
movslq %eax,%rsi
movsd  (%r15,%rsi,8),%xmm0
movsd  (%r14,%rsi,8),%xmm1
mov    %rbx,%rdi
mov    $0x2,%al
call   1050 <printf@plt>
mov    (%rsp),%eax
add    $0x1,%eax
mov    %eax,(%rsp)
cmp    0x8(%rsp),%eax
jl     1620 <main+0x480>
xor    %eax,%eax
add    $0x188,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
%struct.timeval = type { double, double }
@__const.main.name = [14 x i8] c"../../system1\00"
@.str = [2 x i8] c"r\00"
@.str.1 = [6 x i8] c"%d%d\0A\00"
@.str.2 = [12 x i8] c"%d%d%lf%lf\0A\00"
@.str.3 = [34 x i8] c"sparse matrix -system1- size %.d\0A\00"
@.str.4 = [9 x i8] c"nnz %d \0A\00"
@.str.5 = [22 x i8] c"FEAST OUTPUT INFO %d\0A\00"
@.str.6 = [51 x i8] c"*************************************************\0A\00"
@.str.7 = [51 x i8] c"************** REPORT ***************************\0A\00"
@.str.8 = [20 x i8] c"SIMULATION TIME %f\0A\00"
@.str.9 = [43 x i8] c"# Search interval [Emin,Emax] %.15e %.15e\0A\00"
@.str.10 = [30 x i8] c"# mode found/subspace %d %d \0A\00"
@.str.11 = [18 x i8] c"# iterations %d \0A\00"
@.str.12 = [13 x i8] c"TRACE %.15e\0A\00"
@.str.13 = [35 x i8] c"Relative error on the Trace %.15e\0A\00"
@.str.14 = [23 x i8] c"Eigenvalues/Residuals\0A\00"
@.str.15 = [19 x i8] c" %d %.15e %.15e\0A\00"
define i32 @main() {
%1 = alloca i32
%2 = alloca [64 x i32]
%3 = alloca double
%4 = alloca i32
%5 = alloca i8
%6 = alloca i32*
%7 = alloca [14 x i8]
%8 = alloca i32
%9 = alloca i32
%10 = alloca double*
%11 = alloca double*
%12 = alloca i32*
%13 = alloca i32*
%14 = alloca %struct.timeval
%15 = alloca %struct.timeval
%16 = alloca i32
%17 = alloca i32
%18 = alloca i32
%19 = alloca i32
%20 = alloca i32
%21 = alloca i32
%22 = alloca double
%23 = alloca double
%24 = alloca double
%25 = alloca double*
%26 = alloca double*
%27 = alloca double*
store i32 0, i32* %1
store i8 70, i8* %5
%28 = bitcast [14 x i8]* %7 to i8*
%29 = getelementptr inbounds [14 x i8], [14 x i8]* %7, i64 0, i64 0
%30 = call i32* @fopen(i8* %29, i8* getelementptr inbounds ([2 x i8], [2 x i8]* @.str, i64 0, i64 0))
store i32* %30, i32** %6
%31 = load i32*, i32** %6
%32 = call i32 (i32*, i8*, i32*, i32*, ...) @fscanf(i32* %31, i8* getelementptr inbounds ([6 x i8], [6 x i8]* @.str.1, i64 0, i64 0), i32* %8, i32* %9)
store i32 %32, i32* %18
%33 = load i32, i32* %9
%34 = call i8* @calloc(i32 %33, i32 8)
%35 = bitcast i8* %34 to double*
store double* %35, double** %10
%36 = load i32, i32* %9
%37 = call i8* @calloc(i32 %36, i32 8)
%38 = bitcast i8* %37 to double*
store double* %38, double** %11
%39 = load i32, i32* %8
%40 = add nsw i32 %39, 1
%41 = call i8* @calloc(i32 %40, i32 4)
%42 = bitcast i8* %41 to i32*
store i32* %42, i32** %12
%43 = load i32, i32* %9
%44 = call i8* @calloc(i32 %43, i32 4)
%45 = bitcast i8* %44 to i32*
store i32* %45, i32** %13
store i32 0, i32* %16
br label %46
46: ; preds = %55, %0
%47 = load i32, i32* %16
%48 = load i32, i32* %8
%49 = icmp sle i32 %47, %48
br i1 %49, label %50, label %58
50: ; preds = %46
%51 = load i32*, i32** %12
%52 = load i32, i32* %16
%53 = sext i32 %52 to i64
%54 = getelementptr inbounds i32, i32* %51, i64 %53
store i32 0, i32* %54
br label %55
55: ; preds = %50
%56 = load i32, i32* %16
%57 = add nsw i32 %56, 1
store i32 %57, i32* %16
br label %46
58: ; preds = %46
%59 = load i32*, i32** %12
store i32 1, i32* %59
store i32 0, i32* %17
br label %60
60: ; preds = %90, %58
%61 = load i32, i32* %17
%62 = load i32, i32* %9
%63 = sub nsw i32 %62, 1
%64 = icmp sle i32 %61, %63
br i1 %64, label %65, label %93
65: ; preds = %60
%66 = load i32*, i32** %6
%67 = load i32*, i32** %13
%68 = load i32, i32* %17
%69 = sext i32 %68 to i64
%70 = getelementptr inbounds i32, i32* %67, i64 %69
%71 = load double*, double** %10
%72 = load i32, i32* %17
%73 = sext i32 %72 to i64
%74 = getelementptr inbounds double, double* %71, i64 %73
%75 = load double*, double** %11
%76 = load i32, i32* %17
%77 = sext i32 %76 to i64
%78 = getelementptr inbounds double, double* %75, i64 %77
%79 = call i32 (i32*, i8*, i32*, i32*, ...) @fscanf(i32* %66, i8* getelementptr inbounds ([12 x i8], [12 x i8]* @.str.2, i64 0, i64 0), i32* %16, i32* %70, double* %74, double* %78)
store i32 %79, i32* %18
%80 = load i32*, i32** %12
%81 = load i32, i32* %16
%82 = sext i32 %81 to i64
%83 = getelementptr inbounds i32, i32* %80, i64 %82
%84 = load i32, i32* %83
%85 = add nsw i32 %84, 1
%86 = load i32*, i32** %12
%87 = load i32, i32* %16
%88 = sext i32 %87 to i64
%89 = getelementptr inbounds i32, i32* %86, i64 %88
store i32 %85, i32* %89
br label %90
90: ; preds = %65
%91 = load i32, i32* %17
%92 = add nsw i32 %91, 1
store i32 %92, i32* %17
br label %60
93: ; preds = %60
%94 = load i32*, i32** %6
%95 = call i32 @fclose(i32* %94)
store i32 1, i32* %16
br label %96
96: ; preds = %117, %93
%97 = load i32, i32* %16
%98 = load i32, i32* %8
%99 = icmp sle i32 %97, %98
br i1 %99, label %100, label %120
100: ; preds = %96
%101 = load i32*, i32** %12
%102 = load i32, i32* %16
%103 = sext i32 %102 to i64
%104 = getelementptr inbounds i32, i32* %101, i64 %103
%105 = load i32, i32* %104
%106 = load i32*, i32** %12
%107 = load i32, i32* %16
%108 = sext i32 %107 to i64
%109 = getelementptr inbounds i32, i32* %106, i64 %108
%110 = getelementptr inbounds i32, i32* %109, i64 -1
%111 = load i32, i32* %110
%112 = add nsw i32 %105, %111
%113 = load i32*, i32** %12
%114 = load i32, i32* %16
%115 = sext i32 %114 to i64
%116 = getelementptr inbounds i32, i32* %113, i64 %115
store i32 %112, i32* %116
br label %117
117: ; preds = %100
%118 = load i32, i32* %16
%119 = add nsw i32 %118, 1
store i32 %119, i32* %16
br label %96
120: ; preds = %96
%121 = load i32, i32* %8
%122 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([34 x i8], [34 x i8]* @.str.3, i64 0, i64 0), i32 %121)
%123 = load i32, i32* %9
%124 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([9 x i8], [9 x i8]* @.str.4, i64 0, i64 0), i32 %123)
%125 = call i32 @gettimeofday(%struct.timeval* %14, i32* null)
store double 1.800000e-01, double* %22
store double 1.000000e+00, double* %23
store i32 25, i32* %19
%126 = load i32, i32* %19
%127 = call i8* @calloc(i32 %126, i32 8)
%128 = bitcast i8* %127 to double*
store double* %128, double** %26
%129 = load i32, i32* %19
%130 = call i8* @calloc(i32 %129, i32 8)
%131 = bitcast i8* %130 to double*
store double* %131, double** %27
%132 = load i32, i32* %8
%133 = load i32, i32* %19
%134 = mul nsw i32 %132, %133
%135 = call i8* @calloc(i32 %134, i32 8)
%136 = bitcast i8* %135 to double*
store double* %136, double** %25
%137 = getelementptr inbounds [64 x i32], [64 x i32]* %2, i64 0, i64 0
%138 = call i32 @feastinit(i32* %137)
%139 = getelementptr inbounds [64 x i32], [64 x i32]* %2, i64 0, i64 0
store i32 1, i32* %139
%140 = load double*, double** %10
%141 = load i32*, i32** %12
%142 = load i32*, i32** %13
%143 = load double*, double** %11
%144 = load i32*, i32** %12
%145 = load i32*, i32** %13
%146 = getelementptr inbounds [64 x i32], [64 x i32]* %2, i64 0, i64 0
%147 = load double*, double** %26
%148 = load double*, double** %25
%149 = load double*, double** %27
%150 = call i32 @dfeast_scsrgv(i8* %5, i32* %8, double* %140, i32* %141, i32* %142, double* %143, i32* %144, i32* %145, i32* %146, double* %3, i32* %4, double* %22, double* %23, i32* %19, double* %147, double* %148, i32* %20, double* %149, i32* %21)
%151 = call i32 @gettimeofday(%struct.timeval* %15, i32* null)
%152 = load i32, i32* %21
%153 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([22 x i8], [22 x i8]* @.str.5, i64 0, i64 0), i32 %152)
%154 = load i32, i32* %21
%155 = icmp eq i32 %154, 0
br i1 %155, label %156, label %225
156: ; preds = %120
%157 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([51 x i8], [51 x i8]* @.str.6, i64 0, i64 0))
%158 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([51 x i8], [51 x i8]* @.str.7, i64 0, i64 0))
%159 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([51 x i8], [51 x i8]* @.str.6, i64 0, i64 0))
%160 = getelementptr inbounds %struct.timeval, %struct.timeval* %15, i32 0, i32 0
%161 = load double, double* %160
%162 = getelementptr inbounds %struct.timeval, %struct.timeval* %14, i32 0, i32 0
%163 = load double, double* %162
%164 = fsub double %161, %163
%165 = getelementptr inbounds %struct.timeval, %struct.timeval* %15, i32 0, i32 1
%166 = load double, double* %165
%167 = getelementptr inbounds %struct.timeval, %struct.timeval* %14, i32 0, i32 1
%168 = load double, double* %167
%169 = fsub double %166, %168
%170 = fmul double %169, 0x3EB0C6F7A0B5ED8D
%171 = call double @llvm.fmuladd.f64(double %164, double 1.000000e+00, double %170)
%172 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([20 x i8], [20 x i8]* @.str.8, i64 0, i64 0), double %171)
%173 = load double, double* %22
%174 = load double, double* %23
%175 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([43 x i8], [43 x i8]* @.str.9, i64 0, i64 0), double %173, double %174)
%176 = load i32, i32* %20
%177 = load i32, i32* %19
%178 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([30 x i8], [30 x i8]* @.str.10, i64 0, i64 0), i32 %176, i32 %177)
%179 = load i32, i32* %4
%180 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([18 x i8], [18 x i8]* @.str.11, i64 0, i64 0), i32 %179)
store double 0.000000e+00, double* %24
store i32 0, i32* %16
br label %181
181: ; preds = %194, %156
%182 = load i32, i32* %16
%183 = load i32, i32* %20
%184 = sub nsw i32 %183, 1
%185 = icmp sle i32 %182, %184
br i1 %185, label %186, label %197
186: ; preds = %181
%187 = load double, double* %24
%188 = load double*, double** %26
%189 = load i32, i32* %16
%190 = sext i32 %189 to i64
%191 = getelementptr inbounds double, double* %188, i64 %190
%192 = load double, double* %191
%193 = fadd double %187, %192
store double %193, double* %24
br label %194
194: ; preds = %186
%195 = load i32, i32* %16
%196 = add nsw i32 %195, 1
store i32 %196, i32* %16
br label %181
197: ; preds = %181
%198 = load double, double* %24
%199 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([13 x i8], [13 x i8]* @.str.12, i64 0, i64 0), double %198)
%200 = load double, double* %3
%201 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([35 x i8], [35 x i8]* @.str.13, i64 0, i64 0), double %200)
%202 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([23 x i8], [23 x i8]* @.str.14, i64 0, i64 0))
store i32 0, i32* %16
br label %203
203: ; preds = %221, %197
%204 = load i32, i32* %16
%205 = load i32, i32* %20
%206 = sub nsw i32 %205, 1
%207 = icmp sle i32 %204, %206
br i1 %207, label %208, label %224
208: ; preds = %203
%209 = load i32, i32* %16
%210 = load double*, double** %26
%211 = load i32, i32* %16
%212 = sext i32 %211 to i64
%213 = getelementptr inbounds double, double* %210, i64 %212
%214 = load double, double* %213
%215 = load double*, double** %27
%216 = load i32, i32* %16
%217 = sext i32 %216 to i64
%218 = getelementptr inbounds double, double* %215, i64 %217
%219 = load double, double* %218
%220 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([19 x i8], [19 x i8]* @.str.15, i64 0, i64 0), i32 %209, double %214, double %219)
br label %221
221: ; preds = %208
%222 = load i32, i32* %16
%223 = add nsw i32 %222, 1
store i32 %223, i32* %16
br label %203
224: ; preds = %203
br label %225
225: ; preds = %224, %120
ret i32 0
}

LLVM IR (-O2):
%struct.timeval = type { double, double }
@__const.main.name = [14 x i8] c"../../system1\00"
@.str = [2 x i8] c"r\00"
@.str.1 = [6 x i8] c"%d%d\0A\00"
@.str.2 = [12 x i8] c"%d%d%lf%lf\0A\00"
@.str.3 = [34 x i8] c"sparse matrix -system1- size %.d\0A\00"
@.str.4 = [9 x i8] c"nnz %d \0A\00"
@.str.5 = [22 x i8] c"FEAST OUTPUT INFO %d\0A\00"
@.str.8 = [20 x i8] c"SIMULATION TIME %f\0A\00"
@.str.9 = [43 x i8] c"# Search interval [Emin,Emax] %.15e %.15e\0A\00"
@.str.10 = [30 x i8] c"# mode found/subspace %d %d \0A\00"
@.str.11 = [18 x i8] c"# iterations %d \0A\00"
@.str.12 = [13 x i8] c"TRACE %.15e\0A\00"
@.str.13 = [35 x i8] c"Relative error on the Trace %.15e\0A\00"
@.str.15 = [19 x i8] c" %d %.15e %.15e\0A\00"
@str.16 = [50 x i8] c"************** REPORT ***************************\00"
@str.17 = [50 x i8] c"*************************************************\00"
@str.18 = [22 x i8] c"Eigenvalues/Residuals\00"
define i32 @main() {
%1 = alloca [64 x i32]
%2 = alloca double
%3 = alloca i32
%4 = alloca i8
%5 = alloca [14 x i8]
%6 = alloca i32
%7 = alloca i32
%8 = alloca <2 x double>
%9 = bitcast <2 x double>* %8 to %struct.timeval*
%10 = alloca <2 x double>
%11 = bitcast <2 x double>* %10 to %struct.timeval*
%12 = alloca i32
%13 = alloca i32
%14 = alloca i32
%15 = alloca i32
%16 = alloca double
%17 = alloca double
%18 = bitcast [64 x i32]* %1 to i8*
%19 = bitcast double* %2 to i8*
%20 = bitcast i32* %3 to i8*
store i8 70, i8* %4
%21 = getelementptr inbounds [14 x i8], [14 x i8]* %5, i64 0, i64 0
%22 = bitcast i32* %6 to i8*
%23 = bitcast i32* %7 to i8*
%24 = bitcast <2 x double>* %8 to i8*
%25 = bitcast <2 x double>* %10 to i8*
%26 = bitcast i32* %12 to i8*
%27 = bitcast i32* %13 to i8*
%28 = bitcast i32* %14 to i8*
%29 = bitcast i32* %15 to i8*
%30 = bitcast double* %16 to i8*
%31 = bitcast double* %17 to i8*
%32 = call i32* @fopen(i8* %21, i8* getelementptr inbounds ([2 x i8], [2 x i8]* @.str, i64 0, i64 0))
%33 = call i32 (i32*, i8*, i32*, i32*, ...) @fscanf(i32* %32, i8* getelementptr inbounds ([6 x i8], [6 x i8]* @.str.1, i64 0, i64 0), i32* %6, i32* %7)
%34 = load i32, i32* %7
%35 = call i8* @calloc(i32 %34, i32 8)
%36 = bitcast i8* %35 to double*
%37 = call i8* @calloc(i32 %34, i32 8)
%38 = bitcast i8* %37 to double*
%39 = load i32, i32* %6
%40 = add i32 %39, 1
%41 = call i8* @calloc(i32 %40, i32 4)
%42 = bitcast i8* %41 to i32*
%43 = call i8* @calloc(i32 %34, i32 4)
%44 = bitcast i8* %43 to i32*
%45 = icmp slt i32 %39, 0
br i1 %45, label %51, label %46
46: ; preds = %0
%47 = zext i32 %39 to i64
%48 = shl nuw nsw i64 %47, 2
%49 = add nuw nsw i64 %48, 4
%50 = add nuw i32 %39, 1
br label %51
51: ; preds = %46, %0
%52 = phi i32 [ 0, %0 ], [ %50, %46 ]
store i32 %52, i32* %12
store i32 1, i32* %42
%53 = icmp sgt i32 %34, 0
br i1 %53, label %54, label %69
54: ; preds = %51, %54
%55 = phi i64 [ %65, %54 ], [ 0, %51 ]
%56 = getelementptr inbounds i32, i32* %44, i64 %55
%57 = getelementptr inbounds double, double* %36, i64 %55
%58 = getelementptr inbounds double, double* %38, i64 %55
%59 = call i32 (i32*, i8*, i32*, i32*, ...) @fscanf(i32* %32, i8* getelementptr inbounds ([12 x i8], [12 x i8]* @.str.2, i64 0, i64 0), i32* %12, i32* %56, double* %57, double* %58)
%60 = load i32, i32* %12
%61 = sext i32 %60 to i64
%62 = getelementptr inbounds i32, i32* %42, i64 %61
%63 = load i32, i32* %62
%64 = add nsw i32 %63, 1
store i32 %64, i32* %62
%65 = add nuw nsw i64 %55, 1
%66 = load i32, i32* %7
%67 = sext i32 %66 to i64
%68 = icmp slt i64 %65, %67
br i1 %68, label %54, label %69
69: ; preds = %54, %51
%70 = call i32 @fclose(i32* %32)
%71 = load i32, i32* %6
%72 = icmp slt i32 %71, 1
br i1 %72, label %119, label %73
73: ; preds = %69
%74 = add nuw i32 %71, 1
%75 = zext i32 %74 to i64
%76 = load i32, i32* %42
%77 = add nsw i64 %75, -1
%78 = add nsw i64 %75, -2
%79 = and i64 %77, 3
%80 = icmp ult i64 %78, 3
br i1 %80, label %105, label %81
81: ; preds = %73
%82 = and i64 %77, -4
br label %83
83: ; preds = %83, %81
%84 = phi i32 [ %76, %81 ], [ %101, %83 ]
%85 = phi i64 [ 1, %81 ], [ %102, %83 ]
%86 = phi i64 [ 0, %81 ], [ %103, %83 ]
%87 = getelementptr inbounds i32, i32* %42, i64 %85
%88 = load i32, i32* %87
%89 = add nsw i32 %84, %88
store i32 %89, i32* %87
%90 = add nuw nsw i64 %85, 1
%91 = getelementptr inbounds i32, i32* %42, i64 %90
%92 = load i32, i32* %91
%93 = add nsw i32 %89, %92
store i32 %93, i32* %91
%94 = add nuw nsw i64 %85, 2
%95 = getelementptr inbounds i32, i32* %42, i64 %94
%96 = load i32, i32* %95
%97 = add nsw i32 %93, %96
store i32 %97, i32* %95
%98 = add nuw nsw i64 %85, 3
%99 = getelementptr inbounds i32, i32* %42, i64 %98
%100 = load i32, i32* %99
%101 = add nsw i32 %97, %100
store i32 %101, i32* %99
%102 = add nuw nsw i64 %85, 4
%103 = add i64 %86, 4
%104 = icmp eq i64 %103, %82
br i1 %104, label %105, label %83
105: ; preds = %83, %73
%106 = phi i32 [ %76, %73 ], [ %101, %83 ]
%107 = phi i64 [ 1, %73 ], [ %102, %83 ]
%108 = icmp eq i64 %79, 0
br i1 %108, label %119, label %109
109: ; preds = %105, %109
%110 = phi i32 [ %115, %109 ], [ %106, %105 ]
%111 = phi i64 [ %116, %109 ], [ %107, %105 ]
%112 = phi i64 [ %117, %109 ], [ 0, %105 ]
%113 = getelementptr inbounds i32, i32* %42, i64 %111
%114 = load i32, i32* %113
%115 = add nsw i32 %110, %114
store i32 %115, i32* %113
%116 = add nuw nsw i64 %111, 1
%117 = add i64 %112, 1
%118 = icmp eq i64 %117, %79
br i1 %118, label %119, label %109
119: ; preds = %105, %109, %69
%120 = phi i32 [ 1, %69 ], [ %74, %109 ], [ %74, %105 ]
store i32 %120, i32* %12
%121 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([34 x i8], [34 x i8]* @.str.3, i64 0, i64 0), i32 %71)
%122 = load i32, i32* %7
%123 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([9 x i8], [9 x i8]* @.str.4, i64 0, i64 0), i32 %122)
%124 = call i32 @gettimeofday(%struct.timeval* %9, i32* null)
store double 1.800000e-01, double* %16
store double 1.000000e+00, double* %17
store i32 25, i32* %13
%125 = call dereferenceable_or_null(200) i8* @calloc(i32 25, i32 8)
%126 = bitcast i8* %125 to double*
%127 = call dereferenceable_or_null(200) i8* @calloc(i32 25, i32 8)
%128 = bitcast i8* %127 to double*
%129 = load i32, i32* %6
%130 = mul nsw i32 %129, 25
%131 = call i8* @calloc(i32 %130, i32 8)
%132 = bitcast i8* %131 to double*
%133 = getelementptr inbounds [64 x i32], [64 x i32]* %1, i64 0, i64 0
%134 = call i32 @feastinit(i32* %133)
store i32 1, i32* %133
%135 = call i32 @dfeast_scsrgv(i8* %4, i32* %6, double* %36, i32* %42, i32* %44, double* %38, i32* %42, i32* %44, i32* %133, double* %2, i32* %3, double* %16, double* %17, i32* %13, double* %126, double* %132, i32* %14, double* %128, i32* %15)
%136 = call i32 @gettimeofday(%struct.timeval* %11, i32* null)
%137 = load i32, i32* %15
%138 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([22 x i8], [22 x i8]* @.str.5, i64 0, i64 0), i32 %137)
%139 = load i32, i32* %15
%140 = icmp eq i32 %139, 0
br i1 %140, label %141, label %244
141: ; preds = %119
%142 = call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([50 x i8], [50 x i8]* @str.17, i64 0, i64 0))
%143 = call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([50 x i8], [50 x i8]* @str.16, i64 0, i64 0))
%144 = call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([50 x i8], [50 x i8]* @str.17, i64 0, i64 0))
%145 = load <2 x double>, <2 x double>* %10
%146 = load <2 x double>, <2 x double>* %8
%147 = fsub <2 x double> %145, %146
%148 = extractelement <2 x double> %147, i64 1
%149 = fmul double %148, 0x3EB0C6F7A0B5ED8D
%150 = extractelement <2 x double> %147, i64 0
%151 = fadd double %150, %149
%152 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([20 x i8], [20 x i8]* @.str.8, i64 0, i64 0), double %151)
%153 = load double, double* %16
%154 = load double, double* %17
%155 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([43 x i8], [43 x i8]* @.str.9, i64 0, i64 0), double %153, double %154)
%156 = load i32, i32* %14
%157 = load i32, i32* %13
%158 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([30 x i8], [30 x i8]* @.str.10, i64 0, i64 0), i32 %156, i32 %157)
%159 = load i32, i32* %3
%160 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([18 x i8], [18 x i8]* @.str.11, i64 0, i64 0), i32 %159)
%161 = load i32, i32* %14
%162 = icmp sgt i32 %161, 0
br i1 %162, label %163, label %223
163: ; preds = %141
%164 = zext i32 %161 to i64
%165 = add nsw i64 %164, -1
%166 = and i64 %164, 7
%167 = icmp ult i64 %165, 7
br i1 %167, label %208, label %168
168: ; preds = %163
%169 = and i64 %164, 4294967288
br label %170
170: ; preds = %170, %168
%171 = phi i64 [ 0, %168 ], [ %205, %170 ]
%172 = phi double [ 0.000000e+00, %168 ], [ %204, %170 ]
%173 = phi i64 [ 0, %168 ], [ %206, %170 ]
%174 = getelementptr inbounds double, double* %126, i64 %171
%175 = load double, double* %174
%176 = fadd double %172, %175
%177 = or i64 %171, 1
%178 = getelementptr inbounds double, double* %126, i64 %177
%179 = load double, double* %178
%180 = fadd double %176, %179
%181 = or i64 %171, 2
%182 = getelementptr inbounds double, double* %126, i64 %181
%183 = load double, double* %182
%184 = fadd double %180, %183
%185 = or i64 %171, 3
%186 = getelementptr inbounds double, double* %126, i64 %185
%187 = load double, double* %186
%188 = fadd double %184, %187
%189 = or i64 %171, 4
%190 = getelementptr inbounds double, double* %126, i64 %189
%191 = load double, double* %190
%192 = fadd double %188, %191
%193 = or i64 %171, 5
%194 = getelementptr inbounds double, double* %126, i64 %193
%195 = load double, double* %194
%196 = fadd double %192, %195
%197 = or i64 %171, 6
%198 = getelementptr inbounds double, double* %126, i64 %197
%199 = load double, double* %198
%200 = fadd double %196, %199
%201 = or i64 %171, 7
%202 = getelementptr inbounds double, double* %126, i64 %201
%203 = load double, double* %202
%204 = fadd double %200, %203
%205 = add nuw nsw i64 %171, 8
%206 = add i64 %173, 8
%207 = icmp eq i64 %206, %169
br i1 %207, label %208, label %170
208: ; preds = %170, %163
%209 = phi double [ undef, %163 ], [ %204, %170 ]
%210 = phi i64 [ 0, %163 ], [ %205, %170 ]
%211 = phi double [ 0.000000e+00, %163 ], [ %204, %170 ]
%212 = icmp eq i64 %166, 0
br i1 %212, label %223, label %213
213: ; preds = %208, %213
%214 = phi i64 [ %220, %213 ], [ %210, %208 ]
%215 = phi double [ %219, %213 ], [ %211, %208 ]
%216 = phi i64 [ %221, %213 ], [ 0, %208 ]
%217 = getelementptr inbounds double, double* %126, i64 %214
%218 = load double, double* %217
%219 = fadd double %215, %218
%220 = add nuw nsw i64 %214, 1
%221 = add i64 %216, 1
%222 = icmp eq i64 %221, %166
br i1 %222, label %223, label %213
223: ; preds = %208, %213, %141
%224 = phi i32 [ 0, %141 ], [ %161, %213 ], [ %161, %208 ]
%225 = phi double [ 0.000000e+00, %141 ], [ %209, %208 ], [ %219, %213 ]
store i32 %224, i32* %12
%226 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([13 x i8], [13 x i8]* @.str.12, i64 0, i64 0), double %225)
%227 = load double, double* %2
%228 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([35 x i8], [35 x i8]* @.str.13, i64 0, i64 0), double %227)
%229 = call i32 @puts(i8* dereferenceable(1) getelementptr inbounds ([22 x i8], [22 x i8]* @str.18, i64 0, i64 0))
store i32 0, i32* %12
%230 = load i32, i32* %14
%231 = icmp sgt i32 %230, 0
br i1 %231, label %232, label %244
232: ; preds = %223, %232
%233 = phi i32 [ %241, %232 ], [ 0, %223 ]
%234 = sext i32 %233 to i64
%235 = getelementptr inbounds double, double* %126, i64 %234
%236 = load double, double* %235
%237 = getelementptr inbounds double, double* %128, i64 %234
%238 = load double, double* %237
%239 = call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([19 x i8], [19 x i8]* @.str.15, i64 0, i64 0), i32 %233, double %236, double %238)
%240 = load i32, i32* %12
%241 = add nsw i32 %240, 1
store i32 %241, i32* %12
%242 = load i32, i32* %14
%243 = icmp slt i32 %241, %242
br i1 %243, label %232, label %244
244: ; preds = %232, %223, %119
ret i32 0
}