#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
float fabs(float);
int /*<<< orphan*/ free_ivector(int *, int, int);
int /*<<< orphan*/ free_vector(float *, int, int);
int *ivector(int, int);
int /*<<< orphan*/ nrerror(char *);
int /*<<< orphan*/ printf(char *, ...);
float *vector(int, int);
void solvde(int itmax, float conv, float slowc, float scalv[], int indexv[],
            int ne, int nb, int m, float **y, float ***c, float **s) {
  void bksub(int ne, int nb, int jf, int k1, int k2, float ***c);
  void difeq(int k, int k1, int k2, int jsf, int is1, int isf, int indexv[],
             int ne, float **s, float **y);
  void pinvs(int ie1, int ie2, int je1, int jsf, int jc1, int k, float ***c,
             float **s);
  void red(int iz1, int iz2, int jz1, int jz2, int jm1, int jm2, int jmf,
           int ic1, int jc1, int jcf, int kc, float ***c, float **s);
  int ic1, ic2, ic3, ic4, it, j, j1, j2, j3, j4, j5, j6, j7, j8, j9;
  int jc1, jcf, jv, k, k1, k2, km, kp, nvars, *kmax;
  float err, errj, fac, vmax, vz, *ermax;
  kmax = ivector(1, ne);
  ermax = vector(1, ne);
  k1 = 1;
  k2 = m;
  nvars = ne * m;
  j1 = 1;
  j2 = nb;
  j3 = nb + 1;
  j4 = ne;
  j5 = j4 + j1;
  j6 = j4 + j2;
  j7 = j4 + j3;
  j8 = j4 + j4;
  j9 = j8 + j1;
  ic1 = 1;
  ic2 = ne - nb;
  ic3 = ic2 + 1;
  ic4 = ne;
  jc1 = 1;
  jcf = ic3;
  for (it = 1; it <= itmax; it++) {
    k = k1;
    difeq(k, k1, k2, j9, ic3, ic4, indexv, ne, s, y);
    pinvs(ic3, ic4, j5, j9, jc1, k1, c, s);
    for (k = k1 + 1; k <= k2; k++) {
      kp = k - 1;
      difeq(k, k1, k2, j9, ic1, ic4, indexv, ne, s, y);
      red(ic1, ic4, j1, j2, j3, j4, j9, ic3, jc1, jcf, kp, c, s);
      pinvs(ic1, ic4, j3, j9, jc1, k, c, s);
    }
    k = k2 + 1;
    difeq(k, k1, k2, j9, ic1, ic2, indexv, ne, s, y);
    red(ic1, ic2, j5, j6, j7, j8, j9, ic3, jc1, jcf, k2, c, s);
    pinvs(ic1, ic2, j7, j9, jcf, k2 + 1, c, s);
    bksub(ne, nb, jcf, k1, k2, c);
    err = 0.0;
    for (j = 1; j <= ne; j++) {
      jv = indexv[j];
      errj = vmax = 0.0;
      km = 0;
      for (k = k1; k <= k2; k++) {
        vz = fabs(c[jv][1][k]);
        if (vz > vmax) {
          vmax = vz;
          km = k;
        }
        errj += vz;
      }
      err += errj / scalv[j];
      ermax[j] = c[jv][1][km] / scalv[j];
      kmax[j] = km;
    }
    err /= nvars;
    fac = (err > slowc ? slowc / err : 1.0);
    for (j = 1; j <= ne; j++) {
      jv = indexv[j];
      for (k = k1; k <= k2; k++)
        y[j][k] -= fac * c[jv][1][k];
    }
    printf("\n%8s %9s %9s\n", "Iter.", "Error", "FAC");
    printf("%6d %12.6f %11.6f\n", it, err, fac);
    if (err < conv) {
      free_vector(ermax, 1, ne);
      free_ivector(kmax, 1, ne);
      return;
    }
  }
  nrerror("Too many iterations in solvde");
}


x86 assembly (clang, -O0):
<solvde>:
push   %rbp
mov    %rsp,%rbp
push   %r15
push   %r14
push   %r12
push   %rbx
sub    $0x100,%rsp
mov    0x20(%rbp),%rax
mov    0x18(%rbp),%rax
mov    0x10(%rbp),%rax
mov    %edi,-0x24(%rbp)
movss  %xmm0,-0x28(%rbp)
movss  %xmm1,-0x2c(%rbp)
mov    %rsi,-0x38(%rbp)
mov    %rdx,-0x40(%rbp)
mov    %ecx,-0x44(%rbp)
mov    %r8d,-0x48(%rbp)
mov    %r9d,-0x4c(%rbp)
mov    -0x44(%rbp),%esi
mov    $0x1,%edi
call   10d0 <ivector@plt>
mov    %rax,-0xb8(%rbp)
mov    -0x44(%rbp),%esi
mov    $0x1,%edi
call   10b0 <vector@plt>
mov    %rax,-0xd8(%rbp)
movl   $0x1,-0x9c(%rbp)
mov    -0x4c(%rbp),%eax
mov    %eax,-0xa0(%rbp)
mov    -0x44(%rbp),%eax
imul   -0x4c(%rbp),%eax
mov    %eax,-0xac(%rbp)
movl   $0x1,-0x68(%rbp)
mov    -0x48(%rbp),%eax
mov    %eax,-0x6c(%rbp)
mov    -0x48(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x70(%rbp)
mov    -0x44(%rbp),%eax
mov    %eax,-0x74(%rbp)
mov    -0x74(%rbp),%eax
add    -0x68(%rbp),%eax
mov    %eax,-0x78(%rbp)
mov    -0x74(%rbp),%eax
add    -0x6c(%rbp),%eax
mov    %eax,-0x7c(%rbp)
mov    -0x74(%rbp),%eax
add    -0x70(%rbp),%eax
mov    %eax,-0x80(%rbp)
mov    -0x74(%rbp),%eax
add    -0x74(%rbp),%eax
mov    %eax,-0x84(%rbp)
mov    -0x84(%rbp),%eax
add    -0x68(%rbp),%eax
mov    %eax,-0x88(%rbp)
movl   $0x1,-0x50(%rbp)
mov    -0x44(%rbp),%eax
sub    -0x48(%rbp),%eax
mov    %eax,-0x54(%rbp)
mov    -0x54(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x58(%rbp)
mov    -0x44(%rbp),%eax
mov    %eax,-0x5c(%rbp)
movl   $0x1,-0x8c(%rbp)
mov    -0x58(%rbp),%eax
mov    %eax,-0x90(%rbp)
movl   $0x1,-0x60(%rbp)
mov    -0x60(%rbp),%eax
cmp    -0x24(%rbp),%eax
jg     18db <solvde+0x72b>
mov    -0x9c(%rbp),%eax
mov    %eax,-0x98(%rbp)
mov    -0x98(%rbp),%edi
mov    -0x9c(%rbp),%esi
mov    -0xa0(%rbp),%edx
mov    -0x88(%rbp),%ecx
mov    -0x58(%rbp),%r8d
mov    -0x5c(%rbp),%r9d
mov    -0x40(%rbp),%rbx
mov    -0x44(%rbp),%r11d
mov    0x20(%rbp),%r10
mov    0x10(%rbp),%rax
mov    %rbx,(%rsp)
mov    %r11d,0x8(%rsp)
mov    %r10,0x10(%rsp)
mov    %rax,0x18(%rsp)
call   10a0 <difeq@plt>
mov    -0x58(%rbp),%edi
mov    -0x5c(%rbp),%esi
mov    -0x78(%rbp),%edx
mov    -0x88(%rbp),%ecx
mov    -0x8c(%rbp),%r8d
mov    -0x9c(%rbp),%r9d
mov    0x18(%rbp),%r10
mov    0x20(%rbp),%rax
mov    %r10,(%rsp)
mov    %rax,0x8(%rsp)
call   1030 <pinvs@plt>
mov    -0x9c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x98(%rbp)
mov    -0x98(%rbp),%eax
cmp    -0xa0(%rbp),%eax
jg     1474 <solvde+0x2c4>
mov    -0x98(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0xa8(%rbp)
mov    -0x98(%rbp),%edi
mov    -0x9c(%rbp),%esi
mov    -0xa0(%rbp),%edx
mov    -0x88(%rbp),%ecx
mov    -0x50(%rbp),%r8d
mov    -0x5c(%rbp),%r9d
mov    -0x40(%rbp),%rbx
mov    -0x44(%rbp),%r11d
mov    0x20(%rbp),%r10
mov    0x10(%rbp),%rax
mov    %rbx,(%rsp)
mov    %r11d,0x8(%rsp)
mov    %r10,0x10(%rsp)
mov    %rax,0x18(%rsp)
call   10a0 <difeq@plt>
mov    -0x50(%rbp),%edi
mov    -0x5c(%rbp),%esi
mov    -0x68(%rbp),%edx
mov    -0x6c(%rbp),%ecx
mov    -0x70(%rbp),%r8d
mov    -0x74(%rbp),%r9d
mov    -0x88(%rbp),%r12d
mov    -0x58(%rbp),%r15d
mov    -0x8c(%rbp),%r14d
mov    -0x90(%rbp),%ebx
mov    -0xa8(%rbp),%r11d
mov    0x18(%rbp),%r10
mov    0x20(%rbp),%rax
mov    %r12d,(%rsp)
mov    %r15d,0x8(%rsp)
mov    %r14d,0x10(%rsp)
mov    %ebx,0x18(%rsp)
mov    %r11d,0x20(%rsp)
mov    %r10,0x28(%rsp)
mov    %rax,0x30(%rsp)
call   1050 <red@plt>
mov    -0x50(%rbp),%edi
mov    -0x5c(%rbp),%esi
mov    -0x70(%rbp),%edx
mov    -0x88(%rbp),%ecx
mov    -0x8c(%rbp),%r8d
mov    -0x98(%rbp),%r9d
mov    0x18(%rbp),%r10
mov    0x20(%rbp),%rax
mov    %r10,(%rsp)
mov    %rax,0x8(%rsp)
call   1030 <pinvs@plt>
mov    -0x98(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x98(%rbp)
jmp    1363 <solvde+0x1b3>
mov    -0xa0(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x98(%rbp)
mov    -0x98(%rbp),%edi
mov    -0x9c(%rbp),%esi
mov    -0xa0(%rbp),%edx
mov    -0x88(%rbp),%ecx
mov    -0x50(%rbp),%r8d
mov    -0x54(%rbp),%r9d
mov    -0x40(%rbp),%rbx
mov    -0x44(%rbp),%r11d
mov    0x20(%rbp),%r10
mov    0x10(%rbp),%rax
mov    %rbx,(%rsp)
mov    %r11d,0x8(%rsp)
mov    %r10,0x10(%rsp)
mov    %rax,0x18(%rsp)
call   10a0 <difeq@plt>
mov    -0x50(%rbp),%edi
mov    -0x54(%rbp),%esi
mov    -0x78(%rbp),%edx
mov    -0x7c(%rbp),%ecx
mov    -0x80(%rbp),%r8d
mov    -0x84(%rbp),%r9d
mov    -0x88(%rbp),%r12d
mov    -0x58(%rbp),%r15d
mov    -0x8c(%rbp),%r14d
mov    -0x90(%rbp),%ebx
mov    -0xa0(%rbp),%r11d
mov    0x18(%rbp),%r10
mov    0x20(%rbp),%rax
mov    %r12d,(%rsp)
mov    %r15d,0x8(%rsp)
mov    %r14d,0x10(%rsp)
mov    %ebx,0x18(%rsp)
mov    %r11d,0x20(%rsp)
mov    %r10,0x28(%rsp)
mov    %rax,0x30(%rsp)
call   1050 <red@plt>
mov    -0x50(%rbp),%edi
mov    -0x54(%rbp),%esi
mov    -0x80(%rbp),%edx
mov    -0x88(%rbp),%ecx
mov    -0x90(%rbp),%r8d
mov    -0xa0(%rbp),%r9d
add    $0x1,%r9d
mov    0x18(%rbp),%r10
mov    0x20(%rbp),%rax
mov    %r10,(%rsp)
mov    %rax,0x8(%rsp)
call   1030 <pinvs@plt>
mov    -0x44(%rbp),%edi
mov    -0x48(%rbp),%esi
mov    -0x90(%rbp),%edx
mov    -0x9c(%rbp),%ecx
mov    -0xa0(%rbp),%r8d
mov    0x18(%rbp),%r9
call   1060 <bksub@plt>
xorps  %xmm0,%xmm0
movss  %xmm0,-0xbc(%rbp)
ff
movl   $0x1,-0x64(%rbp)
mov    -0x64(%rbp),%eax
cmp    -0x44(%rbp),%eax
jg     1701 <solvde+0x551>
mov    -0x40(%rbp),%rax
movslq -0x64(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x94(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0xc8(%rbp)
ff
xorps  %xmm0,%xmm0
movss  %xmm0,-0xc0(%rbp)
ff
movl   $0x0,-0xa4(%rbp)
mov    -0x9c(%rbp),%eax
mov    %eax,-0x98(%rbp)
mov    -0x98(%rbp),%eax
cmp    -0xa0(%rbp),%eax
jg     167e <solvde+0x4ce>
mov    0x18(%rbp),%rax
movslq -0x94(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    0x8(%rax),%rax
movslq -0x98(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
call   1090 <fabs@plt>
movss  %xmm0,-0xcc(%rbp)
ff
movss  -0xcc(%rbp),%xmm0
ff
ucomiss -0xc8(%rbp),%xmm0
jbe    1652 <solvde+0x4a2>
movss  -0xcc(%rbp),%xmm0
ff
movss  %xmm0,-0xc8(%rbp)
ff
mov    -0x98(%rbp),%eax
mov    %eax,-0xa4(%rbp)
movss  -0xcc(%rbp),%xmm0
ff
addss  -0xc0(%rbp),%xmm0
ff
movss  %xmm0,-0xc0(%rbp)
ff
mov    -0x98(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x98(%rbp)
jmp    15e3 <solvde+0x433>
movss  -0xc0(%rbp),%xmm0
ff
mov    -0x38(%rbp),%rax
movslq -0x64(%rbp),%rcx
divss  (%rax,%rcx,4),%xmm0
addss  -0xbc(%rbp),%xmm0
ff
movss  %xmm0,-0xbc(%rbp)
ff
mov    0x18(%rbp),%rax
movslq -0x94(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    0x8(%rax),%rax
movslq -0xa4(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x38(%rbp),%rax
movslq -0x64(%rbp),%rcx
divss  (%rax,%rcx,4),%xmm0
mov    -0xd8(%rbp),%rax
movslq -0x64(%rbp),%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0xa4(%rbp),%edx
mov    -0xb8(%rbp),%rax
movslq -0x64(%rbp),%rcx
mov    %edx,(%rax,%rcx,4)
mov    -0x64(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x64(%rbp)
jmp    159a <solvde+0x3ea>
cvtsi2ssl -0xac(%rbp),%xmm1
ff
movss  -0xbc(%rbp),%xmm0
ff
divss  %xmm1,%xmm0
movss  %xmm0,-0xbc(%rbp)
ff
movss  -0xbc(%rbp),%xmm0
ff
ucomiss -0x2c(%rbp),%xmm0
jbe    174d <solvde+0x59d>
movss  -0x2c(%rbp),%xmm0
divss  -0xbc(%rbp),%xmm0
ff
cvtss2sd %xmm0,%xmm0
movsd  %xmm0,-0xe0(%rbp)
ff
jmp    1762 <solvde+0x5b2>
movsd  0x8ab(%rip),%xmm0
movsd  %xmm0,-0xe0(%rbp)
ff
jmp    1762 <solvde+0x5b2>
movsd  -0xe0(%rbp),%xmm0
ff
cvtsd2ss %xmm0,%xmm0
movss  %xmm0,-0xc4(%rbp)
ff
movl   $0x1,-0x64(%rbp)
mov    -0x64(%rbp),%eax
cmp    -0x44(%rbp),%eax
jg     1839 <solvde+0x689>
mov    -0x40(%rbp),%rax
movslq -0x64(%rbp),%rcx
mov    (%rax,%rcx,4),%eax
mov    %eax,-0x94(%rbp)
mov    -0x9c(%rbp),%eax
mov    %eax,-0x98(%rbp)
mov    -0x98(%rbp),%eax
cmp    -0xa0(%rbp),%eax
jg     1826 <solvde+0x676>
movss  -0xc4(%rbp),%xmm0
ff
mov    0x18(%rbp),%rax
movslq -0x94(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
mov    0x8(%rax),%rax
movslq -0x98(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm2
mov    0x10(%rbp),%rax
movslq -0x64(%rbp),%rcx
mov    (%rax,%rcx,8),%rax
movslq -0x98(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
movd   %xmm0,%edx
xor    $0x80000000,%edx
movd   %edx,%xmm0
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x98(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x98(%rbp)
jmp    17a6 <solvde+0x5f6>
jmp    182b <solvde+0x67b>
mov    -0x64(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x64(%rbp)
jmp    177d <solvde+0x5cd>
lea    0x7c8(%rip),%rdi
lea    0x7cf(%rip),%rsi
lea    0x7ce(%rip),%rdx
lea    0x7cd(%rip),%rcx
mov    $0x0,%al
call   1070 <printf@plt>
mov    -0x60(%rbp),%esi
movss  -0xbc(%rbp),%xmm0
ff
cvtss2sd %xmm0,%xmm0
movss  -0xc4(%rbp),%xmm1
ff
cvtss2sd %xmm1,%xmm1
lea    0x7a8(%rip),%rdi
mov    $0x2,%al
call   1070 <printf@plt>
movss  -0xbc(%rbp),%xmm1
ff
movss  -0x28(%rbp),%xmm0
ucomiss %xmm1,%xmm0
jbe    18c8 <solvde+0x718>
mov    -0xd8(%rbp),%rdi
mov    -0x44(%rbp),%edx
mov    $0x1,%esi
call   1080 <free_vector@plt>
mov    -0xb8(%rbp),%rdi
mov    -0x44(%rbp),%edx
mov    $0x1,%esi
call   10c0 <free_ivector@plt>
jmp    18e7 <solvde+0x737>
jmp    18cd <solvde+0x71d>
mov    -0x60(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x60(%rbp)
jmp    12c1 <solvde+0x111>
lea    0x757(%rip),%rdi
call   1040 <nrerror@plt>
add    $0x100,%rsp
pop    %rbx
pop    %r12
pop    %r14
pop    %r15
pop    %rbp
ret

x86 assembly (clang, -O2):
<solvde>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x108,%rsp
mov    %r9d,%r13d
mov    %r8,0x38(%rsp)
mov    %ecx,%r15d
mov    %rdx,0x18(%rsp)
mov    %rsi,0xa0(%rsp)
movss  %xmm1,0x20(%rsp)
movss  %xmm0,0x68(%rsp)
mov    %edi,%ebx
mov    0x150(%rsp),%rbp
mov    0x148(%rsp),%rax
mov    %rax,0x10(%rsp)
mov    0x140(%rsp),%rax
mov    %rax,0x30(%rsp)
mov    $0x1,%eax
mov    %rax,0x80(%rsp)
mov    $0x1,%edi
mov    %ecx,%esi
call   10c0 <ivector@plt>
mov    %rax,0x78(%rsp)
mov    $0x1,%edi
mov    %r15d,%esi
call   10a0 <vector@plt>
mov    %rax,0x70(%rsp)
mov    %ebx,0x64(%rsp)
test   %ebx,%ebx
jle    18fc <solvde+0x75c>
mov    %rbp,%rbx
mov    0x38(%rsp),%rsi
lea    0x1(%rsi),%eax
mov    %eax,0x6c(%rsp)
lea    0x1(%r15),%ecx
lea    (%rsi,%r15,1),%eax
mov    %eax,0x60(%rsp)
lea    0x1(%rsi,%r15,1),%eax
mov    %eax,0x24(%rsp)
lea    (%r15,%r15,1),%eax
mov    %eax,0x5c(%rsp)
lea    (%r15,%r15,1),%r8d
add    $0x1,%r8d
mov    %r15d,%edx
sub    %esi,%edx
mov    %r13d,%eax
imul   %r15d,%eax
cvtsi2ss %eax,%xmm0
movss  %xmm0,0x58(%rsp)
lea    0x1(%r13),%r14d
lea    -0x1(%r14),%rax
lea    -0x9(%r14),%rsi
mov    %rsi,0xb8(%rsp)
mov    %rsi,%rdi
shr    $0x3,%rdi
add    $0x1,%rdi
mov    %rax,%rbp
and    $0xfffffffffffffffe,%rbp
mov    %rax,%rsi
mov    %rax,0x88(%rsp)
and    $0xfffffffffffffff8,%rsi
mov    %rdi,%rax
mov    %rdi,0xe0(%rsp)
and    $0xfffffffffffffffe,%rdi
mov    %rdi,0xc8(%rsp)
mov    %r14,%rax
neg    %rax
mov    %rax,0xe8(%rsp)
mov    %rdx,0xc0(%rsp)
lea    0x1(%rdx),%eax
mov    %ecx,0x28(%rsp)
mov    %ecx,%ecx
mov    %rcx,0x90(%rsp)
mov    %rsi,0xd8(%rsp)
lea    0x1(%rsi),%rcx
mov    %rcx,0xd0(%rsp)
mov    %rbx,0x8(%rsp)
mov    %r8,0xa8(%rsp)
mov    %rax,0x40(%rsp)
mov    %rbp,0x98(%rsp)
mov    %r14,0xb0(%rsp)
mov    %r13,0x50(%rsp)
mov    %r15,0x48(%rsp)
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
mov    $0x1,%edi
mov    $0x1,%esi
mov    %r13d,%edx
mov    %r8d,%ecx
mov    %rbx,%r12
mov    %rax,%rbx
mov    %r8,%r14
mov    %ebx,%r8d
mov    %r15d,%r9d
push   0x30(%rsp)
push   %r12
push   %r15
push   0x30(%rsp)
call   1090 <difeq@plt>
add    $0x20,%rsp
mov    %ebx,%edi
mov    %r15d,%esi
mov    0x28(%rsp),%edx
mov    %r14d,%ecx
mov    $0x1,%r8d
mov    $0x1,%r9d
push   %r12
push   0x18(%rsp)
call   1030 <pinvs@plt>
add    $0x10,%rsp
cmp    $0x2,%r13d
mov    0x38(%rsp),%rbp
mov    %r14,%r15
jl     1473 <solvde+0x2d3>
mov    $0x2,%r12d
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
mov    %r12d,%edi
mov    $0x1,%esi
mov    0x50(%rsp),%rdx
mov    %r15d,%ecx
mov    $0x1,%r8d
mov    0x48(%rsp),%r13
mov    %r13d,%r9d
push   0x30(%rsp)
mov    0x10(%rsp),%rbx
push   %rbx
push   %r13
push   0x30(%rsp)
call   1090 <difeq@plt>
add    $0x20,%rsp
lea    -0x1(%r12),%eax
sub    $0x8,%rsp
mov    $0x1,%edi
mov    %r13d,%esi
mov    $0x1,%edx
mov    %ebp,%ecx
mov    0x74(%rsp),%r14d
mov    %r14d,%r8d
mov    %r13d,%r9d
push   %rbx
mov    0x20(%rsp),%rbp
push   %rbp
push   %rax
mov    0x60(%rsp),%rax
push   %rax
push   $0x1
push   %rax
push   %r15
call   1050 <red@plt>
add    $0x40,%rsp
mov    $0x1,%edi
mov    %r13d,%esi
mov    %r14d,%edx
mov    %r15d,%ecx
mov    $0x1,%r8d
mov    %r12d,%r9d
push   %rbx
push   %rbp
mov    0x48(%rsp),%rbp
call   1030 <pinvs@plt>
mov    0x60(%rsp),%r13
add    $0x10,%rsp
lea    0x1(%r12),%eax
cmp    %r13d,%r12d
mov    %eax,%r12d
jne    13c0 <solvde+0x220>
mov    0xb0(%rsp),%r14
mov    %r14d,%edi
mov    $0x1,%esi
mov    %r13d,%edx
mov    %r13,%r12
mov    %r15d,%ecx
mov    $0x1,%r8d
mov    0xc0(%rsp),%r13
mov    %r13d,%r9d
push   0x30(%rsp)
mov    0x10(%rsp),%rax
push   %rax
mov    %rax,%rbx
push   0x58(%rsp)
push   0x30(%rsp)
call   1090 <difeq@plt>
add    $0x18,%rsp
mov    $0x1,%edi
mov    %r13d,%esi
mov    0x30(%rsp),%edx
mov    0x68(%rsp),%ecx
mov    0x2c(%rsp),%r8d
mov    0x64(%rsp),%r9d
push   %rbx
mov    0x20(%rsp),%rbx
push   %rbx
push   %r12
mov    0x60(%rsp),%r12
push   %r12
push   $0x1
push   %r12
push   %r15
call   1050 <red@plt>
add    $0x40,%rsp
mov    $0x1,%edi
mov    %r13d,%esi
mov    0x24(%rsp),%edx
mov    %r15d,%ecx
mov    %r12d,%r8d
mov    %r14d,%r9d
push   0x8(%rsp)
push   %rbx
call   1030 <pinvs@plt>
add    $0x10,%rsp
mov    0x48(%rsp),%r15
mov    %r15d,%edi
mov    %ebp,%esi
mov    %r12d,%edx
mov    $0x1,%ecx
mov    0x50(%rsp),%r13
mov    %r13d,%r8d
mov    %rbx,%r9
call   1060 <bksub@plt>
test   %r15d,%r15d
mov    0x8(%rsp),%rbx
xorps  %xmm4,%xmm4
jle    1650 <solvde+0x4b0>
mov    $0x1,%eax
mov    0x18(%rsp),%r8
mov    0xa0(%rsp),%r9
mov    0x78(%rsp),%r10
mov    0x10(%rsp),%r11
mov    0x70(%rsp),%r12
movaps 0xa8e(%rip),%xmm3
mov    0x98(%rsp),%rbp
jmp    15b9 <solvde+0x419>
nopl   0x0(%rax)
xorps  %xmm0,%xmm0
movss  (%r9,%rax,4),%xmm1
divss  %xmm1,%xmm0
addss  %xmm0,%xmm4
movslq %edx,%rdx
movss  (%rcx,%rdx,4),%xmm0
divss  %xmm1,%xmm0
movss  %xmm0,(%r12,%rax,4)
mov    %edx,(%r10,%rax,4)
add    $0x1,%rax
cmp    0x90(%rsp),%rax
je     165a <solvde+0x4ba>
movslq (%r8,%rax,4),%rcx
mov    (%r11,%rcx,8),%rcx
mov    0x8(%rcx),%rcx
xor    %edx,%edx
test   %r13d,%r13d
jle    1580 <solvde+0x3e0>
xorps  %xmm1,%xmm1
mov    $0x1,%esi
xorps  %xmm0,%xmm0
cmp    $0x2,%r14d
je     1623 <solvde+0x483>
nopl   (%rax)
movss  (%rcx,%rsi,4),%xmm2
andps  %xmm3,%xmm2
ucomiss %xmm1,%xmm2
mov    %esi,%edi
cmovbe %edx,%edi
addss  %xmm2,%xmm0
maxss  %xmm1,%xmm2
movss  0x4(%rcx,%rsi,4),%xmm1
andps  %xmm3,%xmm1
ucomiss %xmm2,%xmm1
lea    0x1(%rsi),%rbx
mov    %ebx,%edx
cmovbe %edi,%edx
addss  %xmm1,%xmm0
maxss  %xmm2,%xmm1
add    $0x2,%rsi
cmp    %rbp,%rbx
jne    15e0 <solvde+0x440>
mov    0x8(%rsp),%rbx
testb  $0x1,0x88(%rsp)
01
je     1583 <solvde+0x3e3>
movss  (%rcx,%rsi,4),%xmm2
andps  %xmm3,%xmm2
ucomiss %xmm1,%xmm2
cmova  %esi,%edx
addss  %xmm2,%xmm0
jmp    1583 <solvde+0x3e3>
nopl   0x0(%rax,%rax,1)
mov    0x18(%rsp),%r8
mov    0x10(%rsp),%r11
divss  0x58(%rsp),%xmm4
ucomiss 0x20(%rsp),%xmm4
movsd  0x9b3(%rip),%xmm0
jbe    167d <solvde+0x4dd>
movss  0x20(%rsp),%xmm0
divss  %xmm4,%xmm0
cvtss2sd %xmm0,%xmm0
movss  %xmm4,0x2c(%rsp)
cvtsd2ss %xmm0,%xmm7
test   %r15d,%r15d
mov    0x30(%rsp),%r9
mov    0x88(%rsp),%r10
mov    0xb8(%rsp),%r12
jle    1870 <solvde+0x6d0>
movaps %xmm7,%xmm0
xorps  0x961(%rip),%xmm0
pshufd $0x0,%xmm0,%xmm1
mov    $0x1,%eax
jmp    16d2 <solvde+0x532>
nopl   0x0(%rax,%rax,1)
add    $0x1,%rax
cmp    0x90(%rsp),%rax
je     1870 <solvde+0x6d0>
test   %r13d,%r13d
jle    16c0 <solvde+0x520>
movslq (%r8,%rax,4),%rcx
mov    (%r11,%rcx,8),%rcx
mov    0x8(%rcx),%rcx
mov    (%r9,%rax,8),%rdx
mov    $0x1,%esi
cmp    $0x8,%r10
jb     17f0 <solvde+0x650>
lea    0x4(%rdx),%rdi
lea    (%rcx,%r14,4),%rbx
cmp    %rbx,%rdi
mov    0x8(%rsp),%rbx
jae    171e <solvde+0x57e>
lea    (%rdx,%r14,4),%rdi
lea    0x4(%rcx),%rbx
cmp    %rdi,%rbx
mov    0x8(%rsp),%rbx
jb     17f0 <solvde+0x650>
cmp    $0x8,%r12
jae    1728 <solvde+0x588>
xor    %esi,%esi
jmp    179e <solvde+0x5fe>
mov    0xc8(%rsp),%rdi
xor    %esi,%esi
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
movups 0x4(%rcx,%rsi,4),%xmm2
movups 0x14(%rcx,%rsi,4),%xmm3
mulps  %xmm1,%xmm2
movups 0x4(%rdx,%rsi,4),%xmm4
addps  %xmm2,%xmm4
movups 0x14(%rdx,%rsi,4),%xmm2
movups 0x24(%rdx,%rsi,4),%xmm5
movups 0x34(%rdx,%rsi,4),%xmm6
mulps  %xmm1,%xmm3
addps  %xmm2,%xmm3
movups %xmm4,0x4(%rdx,%rsi,4)
movups %xmm3,0x14(%rdx,%rsi,4)
movups 0x24(%rcx,%rsi,4),%xmm2
movups 0x34(%rcx,%rsi,4),%xmm3
mulps  %xmm1,%xmm2
addps  %xmm5,%xmm2
mulps  %xmm1,%xmm3
addps  %xmm6,%xmm3
movups %xmm2,0x24(%rdx,%rsi,4)
movups %xmm3,0x34(%rdx,%rsi,4)
add    $0x10,%rsi
add    $0xfffffffffffffffe,%rdi
jne    1740 <solvde+0x5a0>
testb  $0x1,0xe0(%rsp)
01
je     17d7 <solvde+0x637>
shl    $0x2,%rsi
or     $0x4,%rsi
movups (%rcx,%rsi,1),%xmm2
movups 0x10(%rcx,%rsi,1),%xmm3
movups (%rdx,%rsi,1),%xmm4
movups 0x10(%rdx,%rsi,1),%xmm5
mulps  %xmm1,%xmm2
addps  %xmm4,%xmm2
mulps  %xmm1,%xmm3
addps  %xmm5,%xmm3
movups %xmm2,(%rdx,%rsi,1)
movups %xmm3,0x10(%rdx,%rsi,1)
mov    0xd0(%rsp),%rsi
cmp    0xd8(%rsp),%r10
je     16c0 <solvde+0x520>
nopl   (%rax)
mov    %rsi,%rdi
test   $0x1,%r14b
jne    1810 <solvde+0x670>
movss  (%rcx,%rsi,4),%xmm2
mulss  %xmm0,%xmm2
addss  (%rdx,%rsi,4),%xmm2
movss  %xmm2,(%rdx,%rsi,4)
lea    0x1(%rsi),%rdi
not    %rsi
cmp    0xe8(%rsp),%rsi
je     16c0 <solvde+0x520>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movss  (%rcx,%rdi,4),%xmm2
mulss  %xmm0,%xmm2
addss  (%rdx,%rdi,4),%xmm2
movss  %xmm2,(%rdx,%rdi,4)
movss  0x4(%rcx,%rdi,4),%xmm2
mulss  %xmm0,%xmm2
addss  0x4(%rdx,%rdi,4),%xmm2
movss  %xmm2,0x4(%rdx,%rdi,4)
add    $0x2,%rdi
cmp    %rdi,%r14
jne    1830 <solvde+0x690>
jmp    16c0 <solvde+0x520>
nopw   0x0(%rax,%rax,1)
lea    0x7b1(%rip),%rdi
lea    0x7b8(%rip),%rsi
lea    0x7b7(%rip),%rdx
lea    0x7b6(%rip),%rcx
xor    %eax,%eax
movaps %xmm7,0xf0(%rsp)
call   1070 <printf@plt>
movss  0x2c(%rsp),%xmm0
cvtss2sd %xmm0,%xmm0
movaps 0xf0(%rsp),%xmm1
cvtss2sd %xmm1,%xmm1
lea    0x78e(%rip),%rdi
mov    0x80(%rsp),%r12
mov    %r12d,%esi
mov    $0x2,%al
call   1070 <printf@plt>
movss  0x68(%rsp),%xmm0
ucomiss 0x2c(%rsp),%xmm0
ja     1919 <solvde+0x779>
lea    0x1(%r12),%eax
cmp    0x64(%rsp),%r12d
mov    %rax,0x80(%rsp)
mov    0xa8(%rsp),%r8
mov    0x40(%rsp),%rax
jne    1340 <solvde+0x1a0>
lea    0x756(%rip),%rdi
add    $0x108,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
jmp    1040 <nrerror@plt>
mov    0x70(%rsp),%rdi
mov    $0x1,%esi
mov    %r15d,%edx
call   1080 <free_vector@plt>
mov    0x78(%rsp),%rdi
mov    $0x1,%esi
mov    %r15d,%edx
add    $0x108,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
jmp    10b0 <free_ivector@plt>

LLVM IR (-O0):
@.str = [14 x i8] c"\0A%8s %9s %9s\0A\00"
@.str.1 = [6 x i8] c"Iter.\00"
@.str.2 = [6 x i8] c"Error\00"
@.str.3 = [4 x i8] c"FAC\00"
@.str.4 = [19 x i8] c"%6d %12.6f %11.6f\0A\00"
@.str.5 = [30 x i8] c"Too many iterations in solvde\00"
define void @solvde(i32 %0, float %1, float %2, float* %3, i32* %4, i32 %5, i32 %6, i32 %7, float** %8, float*** %9, float** %10) {
%12 = alloca i32
%13 = alloca float
%14 = alloca float
%15 = alloca float*
%16 = alloca i32*
%17 = alloca i32
%18 = alloca i32
%19 = alloca i32
%20 = alloca float**
%21 = alloca float***
%22 = alloca float**
%23 = alloca i32
%24 = alloca i32
%25 = alloca i32
%26 = alloca i32
%27 = alloca i32
%28 = alloca i32
%29 = alloca i32
%30 = alloca i32
%31 = alloca i32
%32 = alloca i32
%33 = alloca i32
%34 = alloca i32
%35 = alloca i32
%36 = alloca i32
%37 = alloca i32
%38 = alloca i32
%39 = alloca i32
%40 = alloca i32
%41 = alloca i32
%42 = alloca i32
%43 = alloca i32
%44 = alloca i32
%45 = alloca i32
%46 = alloca i32
%47 = alloca i32*
%48 = alloca float
%49 = alloca float
%50 = alloca float
%51 = alloca float
%52 = alloca float
%53 = alloca float*
store i32 %0, i32* %12
store float %1, float* %13
store float %2, float* %14
store float* %3, float** %15
store i32* %4, i32** %16
store i32 %5, i32* %17
store i32 %6, i32* %18
store i32 %7, i32* %19
store float** %8, float*** %20
store float*** %9, float**** %21
store float** %10, float*** %22
%54 = load i32, i32* %17
%55 = call i32* @ivector(i32 1, i32 %54)
store i32* %55, i32** %47
%56 = load i32, i32* %17
%57 = call float* @vector(i32 1, i32 %56)
store float* %57, float** %53
store i32 1, i32* %42
%58 = load i32, i32* %19
store i32 %58, i32* %43
%59 = load i32, i32* %17
%60 = load i32, i32* %19
%61 = mul nsw i32 %59, %60
store i32 %61, i32* %46
store i32 1, i32* %29
%62 = load i32, i32* %18
store i32 %62, i32* %30
%63 = load i32, i32* %18
%64 = add nsw i32 %63, 1
store i32 %64, i32* %31
%65 = load i32, i32* %17
store i32 %65, i32* %32
%66 = load i32, i32* %32
%67 = load i32, i32* %29
%68 = add nsw i32 %66, %67
store i32 %68, i32* %33
%69 = load i32, i32* %32
%70 = load i32, i32* %30
%71 = add nsw i32 %69, %70
store i32 %71, i32* %34
%72 = load i32, i32* %32
%73 = load i32, i32* %31
%74 = add nsw i32 %72, %73
store i32 %74, i32* %35
%75 = load i32, i32* %32
%76 = load i32, i32* %32
%77 = add nsw i32 %75, %76
store i32 %77, i32* %36
%78 = load i32, i32* %36
%79 = load i32, i32* %29
%80 = add nsw i32 %78, %79
store i32 %80, i32* %37
store i32 1, i32* %23
%81 = load i32, i32* %17
%82 = load i32, i32* %18
%83 = sub nsw i32 %81, %82
store i32 %83, i32* %24
%84 = load i32, i32* %24
%85 = add nsw i32 %84, 1
store i32 %85, i32* %25
%86 = load i32, i32* %17
store i32 %86, i32* %26
store i32 1, i32* %38
%87 = load i32, i32* %25
store i32 %87, i32* %39
store i32 1, i32* %27
br label %88
88: ; preds = %358, %11
%89 = load i32, i32* %27
%90 = load i32, i32* %12
%91 = icmp sle i32 %89, %90
br i1 %91, label %92, label %361
92: ; preds = %88
%93 = load i32, i32* %42
store i32 %93, i32* %41
%94 = load i32, i32* %41
%95 = load i32, i32* %42
%96 = load i32, i32* %43
%97 = load i32, i32* %37
%98 = load i32, i32* %25
%99 = load i32, i32* %26
%100 = load i32*, i32** %16
%101 = load i32, i32* %17
%102 = load float**, float*** %22
%103 = load float**, float*** %20
%104 = load i32, i32* %25
%105 = load i32, i32* %26
%106 = load i32, i32* %33
%107 = load i32, i32* %37
%108 = load i32, i32* %38
%109 = load i32, i32* %42
%110 = load float***, float**** %21
%111 = load float**, float*** %22
%112 = load i32, i32* %42
%113 = add nsw i32 %112, 1
store i32 %113, i32* %41
br label %114
114: ; preds = %152, %92
%115 = load i32, i32* %41
%116 = load i32, i32* %43
%117 = icmp sle i32 %115, %116
br i1 %117, label %118, label %155
118: ; preds = %114
%119 = load i32, i32* %41
%120 = sub nsw i32 %119, 1
store i32 %120, i32* %45
%121 = load i32, i32* %41
%122 = load i32, i32* %42
%123 = load i32, i32* %43
%124 = load i32, i32* %37
%125 = load i32, i32* %23
%126 = load i32, i32* %26
%127 = load i32*, i32** %16
%128 = load i32, i32* %17
%129 = load float**, float*** %22
%130 = load float**, float*** %20
%131 = load i32, i32* %23
%132 = load i32, i32* %26
%133 = load i32, i32* %29
%134 = load i32, i32* %30
%135 = load i32, i32* %31
%136 = load i32, i32* %32
%137 = load i32, i32* %37
%138 = load i32, i32* %25
%139 = load i32, i32* %38
%140 = load i32, i32* %39
%141 = load i32, i32* %45
%142 = load float***, float**** %21
%143 = load float**, float*** %22
%144 = load i32, i32* %23
%145 = load i32, i32* %26
%146 = load i32, i32* %31
%147 = load i32, i32* %37
%148 = load i32, i32* %38
%149 = load i32, i32* %41
%150 = load float***, float**** %21
%151 = load float**, float*** %22
br label %152
152: ; preds = %118
%153 = load i32, i32* %41
%154 = add nsw i32 %153, 1
store i32 %154, i32* %41
br label %114
155: ; preds = %114
%156 = load i32, i32* %43
%157 = add nsw i32 %156, 1
store i32 %157, i32* %41
%158 = load i32, i32* %41
%159 = load i32, i32* %42
%160 = load i32, i32* %43
%161 = load i32, i32* %37
%162 = load i32, i32* %23
%163 = load i32, i32* %24
%164 = load i32*, i32** %16
%165 = load i32, i32* %17
%166 = load float**, float*** %22
%167 = load float**, float*** %20
%168 = load i32, i32* %23
%169 = load i32, i32* %24
%170 = load i32, i32* %33
%171 = load i32, i32* %34
%172 = load i32, i32* %35
%173 = load i32, i32* %36
%174 = load i32, i32* %37
%175 = load i32, i32* %25
%176 = load i32, i32* %38
%177 = load i32, i32* %39
%178 = load i32, i32* %43
%179 = load float***, float**** %21
%180 = load float**, float*** %22
%181 = load i32, i32* %23
%182 = load i32, i32* %24
%183 = load i32, i32* %35
%184 = load i32, i32* %37
%185 = load i32, i32* %39
%186 = load i32, i32* %43
%187 = add nsw i32 %186, 1
%188 = load float***, float**** %21
%189 = load float**, float*** %22
%190 = load i32, i32* %17
%191 = load i32, i32* %18
%192 = load i32, i32* %39
%193 = load i32, i32* %42
%194 = load i32, i32* %43
%195 = load float***, float**** %21
store float 0.000000e+00, float* %48
store i32 1, i32* %28
br label %196
196: ; preds = %273, %155
%197 = load i32, i32* %28
%198 = load i32, i32* %17
%199 = icmp sle i32 %197, %198
br i1 %199, label %200, label %276
200: ; preds = %196
%201 = load i32*, i32** %16
%202 = load i32, i32* %28
%203 = sext i32 %202 to i64
%204 = getelementptr inbounds i32, i32* %201, i64 %203
%205 = load i32, i32* %204
store i32 %205, i32* %40
store float 0.000000e+00, float* %51
store float 0.000000e+00, float* %49
store i32 0, i32* %44
%206 = load i32, i32* %42
store i32 %206, i32* %41
br label %207
207: ; preds = %234, %200
%208 = load i32, i32* %41
%209 = load i32, i32* %43
%210 = icmp sle i32 %208, %209
br i1 %210, label %211, label %237
211: ; preds = %207
%212 = load float***, float**** %21
%213 = load i32, i32* %40
%214 = sext i32 %213 to i64
%215 = getelementptr inbounds float**, float*** %212, i64 %214
%216 = load float**, float*** %215
%217 = getelementptr inbounds float*, float** %216, i64 1
%218 = load float*, float** %217
%219 = load i32, i32* %41
%220 = sext i32 %219 to i64
%221 = getelementptr inbounds float, float* %218, i64 %220
%222 = load float, float* %221
%223 = call float @fabs(float %222)
store float %223, float* %52
%224 = load float, float* %52
%225 = load float, float* %51
%226 = fcmp ogt float %224, %225
br i1 %226, label %227, label %230
227: ; preds = %211
%228 = load float, float* %52
store float %228, float* %51
%229 = load i32, i32* %41
store i32 %229, i32* %44
br label %230
230: ; preds = %227, %211
%231 = load float, float* %52
%232 = load float, float* %49
%233 = fadd float %232, %231
store float %233, float* %49
br label %234
234: ; preds = %230
%235 = load i32, i32* %41
%236 = add nsw i32 %235, 1
store i32 %236, i32* %41
br label %207
237: ; preds = %207
%238 = load float, float* %49
%239 = load float*, float** %15
%240 = load i32, i32* %28
%241 = sext i32 %240 to i64
%242 = getelementptr inbounds float, float* %239, i64 %241
%243 = load float, float* %242
%244 = fdiv float %238, %243
%245 = load float, float* %48
%246 = fadd float %245, %244
store float %246, float* %48
%247 = load float***, float**** %21
%248 = load i32, i32* %40
%249 = sext i32 %248 to i64
%250 = getelementptr inbounds float**, float*** %247, i64 %249
%251 = load float**, float*** %250
%252 = getelementptr inbounds float*, float** %251, i64 1
%253 = load float*, float** %252
%254 = load i32, i32* %44
%255 = sext i32 %254 to i64
%256 = getelementptr inbounds float, float* %253, i64 %255
%257 = load float, float* %256
%258 = load float*, float** %15
%259 = load i32, i32* %28
%260 = sext i32 %259 to i64
%261 = getelementptr inbounds float, float* %258, i64 %260
%262 = load float, float* %261
%263 = fdiv float %257, %262
%264 = load float*, float** %53
%265 = load i32, i32* %28
%266 = sext i32 %265 to i64
%267 = getelementptr inbounds float, float* %264, i64 %266
store float %263, float* %267
%268 = load i32, i32* %44
%269 = load i32*, i32** %47
%270 = load i32, i32* %28
%271 = sext i32 %270 to i64
%272 = getelementptr inbounds i32, i32* %269, i64 %271
store i32 %268, i32* %272
br label %273
273: ; preds = %237
%274 = load i32, i32* %28
%275 = add nsw i32 %274, 1
store i32 %275, i32* %28
br label %196
276: ; preds = %196
%277 = load i32, i32* %46
%278 = sitofp i32 %277 to float
%279 = load float, float* %48
%280 = fdiv float %279, %278
store float %280, float* %48
%281 = load float, float* %48
%282 = load float, float* %14
%283 = fcmp ogt float %281, %282
br i1 %283, label %284, label %289
284: ; preds = %276
%285 = load float, float* %14
%286 = load float, float* %48
%287 = fdiv float %285, %286
%288 = fpext float %287 to double
br label %290
289: ; preds = %276
br label %290
290: ; preds = %289, %284
%291 = phi double [ %288, %284 ], [ 1.000000e+00, %289 ]
%292 = fptrunc double %291 to float
store float %292, float* %50
store i32 1, i32* %28
br label %293
293: ; preds = %336, %290
%294 = load i32, i32* %28
%295 = load i32, i32* %17
%296 = icmp sle i32 %294, %295
br i1 %296, label %297, label %339
297: ; preds = %293
%298 = load i32*, i32** %16
%299 = load i32, i32* %28
%300 = sext i32 %299 to i64
%301 = getelementptr inbounds i32, i32* %298, i64 %300
%302 = load i32, i32* %301
store i32 %302, i32* %40
%303 = load i32, i32* %42
store i32 %303, i32* %41
br label %304
304: ; preds = %332, %297
%305 = load i32, i32* %41
%306 = load i32, i32* %43
%307 = icmp sle i32 %305, %306
br i1 %307, label %308, label %335
308: ; preds = %304
%309 = load float, float* %50
%310 = load float***, float**** %21
%311 = load i32, i32* %40
%312 = sext i32 %311 to i64
%313 = getelementptr inbounds float**, float*** %310, i64 %312
%314 = load float**, float*** %313
%315 = getelementptr inbounds float*, float** %314, i64 1
%316 = load float*, float** %315
%317 = load i32, i32* %41
%318 = sext i32 %317 to i64
%319 = getelementptr inbounds float, float* %316, i64 %318
%320 = load float, float* %319
%321 = load float**, float*** %20
%322 = load i32, i32* %28
%323 = sext i32 %322 to i64
%324 = getelementptr inbounds float*, float** %321, i64 %323
%325 = load float*, float** %324
%326 = load i32, i32* %41
%327 = sext i32 %326 to i64
%328 = getelementptr inbounds float, float* %325, i64 %327
%329 = load float, float* %328
%330 = fneg float %309
%331 = call float @llvm.fmuladd.f32(float %330, float %320, float %329)
store float %331, float* %328
br label %332
332: ; preds = %308
%333 = load i32, i32* %41
%334 = add nsw i32 %333, 1
store i32 %334, i32* %41
br label %304
335: ; preds = %304
br label %336
336: ; preds = %335
%337 = load i32, i32* %28
%338 = add nsw i32 %337, 1
store i32 %338, i32* %28
br label %293
339: ; preds = %293
%340 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([14 x i8], [14 x i8]* @.str, i64 0, i64 0), i8* getelementptr inbounds ([6 x i8], [6 x i8]* @.str.1, i64 0, i64 0), i8* getelementptr inbounds ([6 x i8], [6 x i8]* @.str.2, i64 0, i64 0), i8* getelementptr inbounds ([4 x i8], [4 x i8]* @.str.3, i64 0, i64 0))
%341 = load i32, i32* %27
%342 = load float, float* %48
%343 = fpext float %342 to double
%344 = load float, float* %50
%345 = fpext float %344 to double
%346 = call i32 (i8*, ...) @printf(i8* getelementptr inbounds ([19 x i8], [19 x i8]* @.str.4, i64 0, i64 0), i32 %341, double %343, double %345)
%347 = load float, float* %48
%348 = load float, float* %13
%349 = fcmp olt float %347, %348
br i1 %349, label %350, label %357
350: ; preds = %339
%351 = load float*, float** %53
%352 = load i32, i32* %17
%353 = call i32 @free_vector(float* %351, i32 1, i32 %352)
%354 = load i32*, i32** %47
%355 = load i32, i32* %17
%356 = call i32 @free_ivector(i32* %354, i32 1, i32 %355)
br label %363
357: ; preds = %339
br label %358
358: ; preds = %357
%359 = load i32, i32* %27
%360 = add nsw i32 %359, 1
store i32 %360, i32* %27
br label %88
361: ; preds = %88
%362 = call i32 @nrerror(i8* getelementptr inbounds ([30 x i8], [30 x i8]* @.str.5, i64 0, i64 0))
br label %363
363: ; preds = %361, %350
ret void
}

LLVM IR (-O2):
@.str = [14 x i8] c"\0A%8s %9s %9s\0A\00"
@.str.1 = [6 x i8] c"Iter.\00"
@.str.2 = [6 x i8] c"Error\00"
@.str.3 = [4 x i8] c"FAC\00"
@.str.4 = [19 x i8] c"%6d %12.6f %11.6f\0A\00"
@.str.5 = [30 x i8] c"Too many iterations in solvde\00"
define void @solvde(i32 %0, float %1, float %2, float* nocapture %3, i32* %4, i32 %5, i32 %6, i32 %7, float** %8, float*** %9, float** %10) {
%12 = tail call i32* @ivector(i32 1, i32 %5)
%13 = tail call float* @vector(i32 1, i32 %5)
%14 = add nsw i32 %6, 1
%15 = add i32 %5, 1
%16 = add nsw i32 %6, %5
%17 = add nsw i32 %14, %5
%18 = shl nsw i32 %5, 1
%19 = or i32 %18, 1
%20 = sub nsw i32 %5, %6
%21 = add nsw i32 %20, 1
%22 = icmp slt i32 %0, 1
br i1 %22, label %270, label %23
23: ; preds = %11
%24 = mul nsw i32 %7, %5
%25 = icmp slt i32 %7, 2
%26 = add i32 %7, 1
%27 = icmp slt i32 %5, 1
%28 = icmp slt i32 %7, 1
%29 = sitofp i32 %24 to float
%30 = zext i32 %15 to i64
%31 = zext i32 %26 to i64
%32 = zext i32 %15 to i64
%33 = zext i32 %26 to i64
%34 = add nsw i64 %31, -1
%35 = add nsw i64 %31, -1
%36 = add nsw i64 %31, -9
%37 = lshr i64 %36, 3
%38 = add nuw nsw i64 %37, 1
%39 = and i64 %35, 1
%40 = icmp eq i32 %26, 2
%41 = and i64 %35, -2
%42 = icmp eq i64 %39, 0
%43 = icmp ult i64 %34, 8
%44 = and i64 %34, -8
%45 = or i64 %44, 1
%46 = and i64 %38, 1
%47 = icmp ult i64 %36, 8
%48 = and i64 %38, 4611686018427387902
%49 = icmp eq i64 %46, 0
%50 = icmp eq i64 %34, %44
%51 = and i64 %31, 1
%52 = icmp eq i64 %51, 0
%53 = sub nsw i64 0, %31
br label %57
54: ; preds = %261
%55 = add nuw i32 %58, 1
%56 = icmp eq i32 %58, %0
br i1 %56, label %270, label %57
57: ; preds = %23, %54
%58 = phi i32 [ 1, %23 ], [ %55, %54 ]
tail call void @difeq(i32 1, i32 1, i32 %7, i32 %19, i32 %21, i32 %5, i32* %4, i32 %5, float** %10, float** %8)
tail call void @pinvs(i32 %21, i32 %5, i32 %15, i32 %19, i32 1, i32 1, float*** %9, float** %10)
br i1 %25, label %64, label %59
59: ; preds = %57, %59
%60 = phi i32 [ %62, %59 ], [ 2, %57 ]
%61 = add nsw i32 %60, -1
tail call void @difeq(i32 %60, i32 1, i32 %7, i32 %19, i32 1, i32 %5, i32* %4, i32 %5, float** %10, float** %8)
tail call void @red(i32 1, i32 %5, i32 1, i32 %6, i32 %14, i32 %5, i32 %19, i32 %21, i32 1, i32 %21, i32 %61, float*** %9, float** %10)
tail call void @pinvs(i32 1, i32 %5, i32 %14, i32 %19, i32 1, i32 %60, float*** %9, float** %10)
%62 = add nuw i32 %60, 1
%63 = icmp eq i32 %60, %7
br i1 %63, label %64, label %59
64: ; preds = %59, %57
tail call void @difeq(i32 %26, i32 1, i32 %7, i32 %19, i32 1, i32 %20, i32* %4, i32 %5, float** %10, float** %8)
tail call void @red(i32 1, i32 %20, i32 %15, i32 %16, i32 %17, i32 %18, i32 %19, i32 %21, i32 1, i32 %21, i32 %7, float*** %9, float** %10)
tail call void @pinvs(i32 1, i32 %20, i32 %17, i32 %19, i32 %21, i32 %26, float*** %9, float** %10)
tail call void @bksub(i32 %5, i32 %6, i32 %21, i32 1, i32 %7, float*** %9)
br i1 %27, label %132, label %65
65: ; preds = %64, %117
%66 = phi i64 [ %130, %117 ], [ 1, %64 ]
%67 = phi float [ %123, %117 ], [ 0.000000e+00, %64 ]
%68 = getelementptr inbounds i32, i32* %4, i64 %66
%69 = load i32, i32* %68
%70 = sext i32 %69 to i64
%71 = getelementptr inbounds float**, float*** %9, i64 %70
%72 = load float**, float*** %71
%73 = getelementptr inbounds float*, float** %72, i64 1
%74 = load float*, float** %73
br i1 %28, label %117, label %75
75: ; preds = %65
br i1 %40, label %102, label %76
76: ; preds = %75, %76
%77 = phi i64 [ %99, %76 ], [ 1, %75 ]
%78 = phi float [ %97, %76 ], [ 0.000000e+00, %75 ]
%79 = phi float [ %98, %76 ], [ 0.000000e+00, %75 ]
%80 = phi i32 [ %96, %76 ], [ 0, %75 ]
%81 = phi i64 [ %100, %76 ], [ 0, %75 ]
%82 = getelementptr inbounds float, float* %74, i64 %77
%83 = load float, float* %82
%84 = tail call float @llvm.fabs.f32(float %83)
%85 = fcmp ogt float %84, %78
%86 = trunc i64 %77 to i32
%87 = select i1 %85, i32 %86, i32 %80
%88 = select i1 %85, float %84, float %78
%89 = fadd float %79, %84
%90 = add nuw nsw i64 %77, 1
%91 = getelementptr inbounds float, float* %74, i64 %90
%92 = load float, float* %91
%93 = tail call float @llvm.fabs.f32(float %92)
%94 = fcmp ogt float %93, %88
%95 = trunc i64 %90 to i32
%96 = select i1 %94, i32 %95, i32 %87
%97 = select i1 %94, float %93, float %88
%98 = fadd float %89, %93
%99 = add nuw nsw i64 %77, 2
%100 = add i64 %81, 2
%101 = icmp eq i64 %100, %41
br i1 %101, label %102, label %76
102: ; preds = %76, %75
%103 = phi i32 [ undef, %75 ], [ %96, %76 ]
%104 = phi float [ undef, %75 ], [ %98, %76 ]
%105 = phi i64 [ 1, %75 ], [ %99, %76 ]
%106 = phi float [ 0.000000e+00, %75 ], [ %97, %76 ]
%107 = phi float [ 0.000000e+00, %75 ], [ %98, %76 ]
%108 = phi i32 [ 0, %75 ], [ %96, %76 ]
br i1 %42, label %117, label %109
109: ; preds = %102
%110 = getelementptr inbounds float, float* %74, i64 %105
%111 = load float, float* %110
%112 = tail call float @llvm.fabs.f32(float %111)
%113 = fcmp ogt float %112, %106
%114 = trunc i64 %105 to i32
%115 = select i1 %113, i32 %114, i32 %108
%116 = fadd float %107, %112
br label %117
117: ; preds = %109, %102, %65
%118 = phi i32 [ 0, %65 ], [ %103, %102 ], [ %115, %109 ]
%119 = phi float [ 0.000000e+00, %65 ], [ %104, %102 ], [ %116, %109 ]
%120 = getelementptr inbounds float, float* %3, i64 %66
%121 = load float, float* %120
%122 = fdiv float %119, %121
%123 = fadd float %67, %122
%124 = sext i32 %118 to i64
%125 = getelementptr inbounds float, float* %74, i64 %124
%126 = load float, float* %125
%127 = fdiv float %126, %121
%128 = getelementptr inbounds float, float* %13, i64 %66
store float %127, float* %128
%129 = getelementptr inbounds i32, i32* %12, i64 %66
store i32 %118, i32* %129
%130 = add nuw nsw i64 %66, 1
%131 = icmp eq i64 %130, %30
br i1 %131, label %132, label %65
132: ; preds = %117, %64
%133 = phi float [ 0.000000e+00, %64 ], [ %123, %117 ]
%134 = fdiv float %133, %29
%135 = fcmp ogt float %134, %2
br i1 %135, label %136, label %139
136: ; preds = %132
%137 = fdiv float %2, %134
%138 = fpext float %137 to double
br label %139
139: ; preds = %132, %136
%140 = phi double [ %138, %136 ], [ 1.000000e+00, %132 ]
%141 = fptrunc double %140 to float
br i1 %27, label %261, label %142
142: ; preds = %139
%143 = fneg float %141
%144 = insertelement <4 x float> poison, float %143, i64 0
%145 = shufflevector <4 x float> %144, <4 x float> poison, <4 x i32> zeroinitializer
%146 = insertelement <4 x float> poison, float %143, i64 0
%147 = shufflevector <4 x float> %146, <4 x float> poison, <4 x i32> zeroinitializer
br label %148
148: ; preds = %142, %258
%149 = phi i64 [ 1, %142 ], [ %259, %258 ]
br i1 %28, label %258, label %150
150: ; preds = %148
%151 = getelementptr inbounds i32, i32* %4, i64 %149
%152 = load i32, i32* %151
%153 = sext i32 %152 to i64
%154 = getelementptr inbounds float**, float*** %9, i64 %153
%155 = load float**, float*** %154
%156 = getelementptr inbounds float*, float** %155, i64 1
%157 = load float*, float** %156
%158 = getelementptr inbounds float*, float** %8, i64 %149
%159 = load float*, float** %158
br i1 %43, label %230, label %160
160: ; preds = %150
%161 = getelementptr float, float* %159, i64 1
%162 = getelementptr float, float* %159, i64 %31
%163 = getelementptr float, float* %157, i64 1
%164 = getelementptr float, float* %157, i64 %31
%165 = icmp ult float* %161, %164
%166 = icmp ult float* %163, %162
%167 = and i1 %165, %166
br i1 %167, label %230, label %168
168: ; preds = %160
br i1 %47, label %209, label %169
169: ; preds = %168, %169
%170 = phi i64 [ %206, %169 ], [ 0, %168 ]
%171 = phi i64 [ %207, %169 ], [ 0, %168 ]
%172 = or i64 %170, 1
%173 = getelementptr inbounds float, float* %157, i64 %172
%174 = bitcast float* %173 to <4 x float>*
%175 = load <4 x float>, <4 x float>* %174
%176 = getelementptr inbounds float, float* %173, i64 4
%177 = bitcast float* %176 to <4 x float>*
%178 = load <4 x float>, <4 x float>* %177
%179 = getelementptr inbounds float, float* %159, i64 %172
%180 = bitcast float* %179 to <4 x float>*
%181 = load <4 x float>, <4 x float>* %180
%182 = getelementptr inbounds float, float* %179, i64 4
%183 = bitcast float* %182 to <4 x float>*
%184 = load <4 x float>, <4 x float>* %183
%185 = call <4 x float> @llvm.fmuladd.v4f32(<4 x float> %145, <4 x float> %175, <4 x float> %181)
%186 = call <4 x float> @llvm.fmuladd.v4f32(<4 x float> %147, <4 x float> %178, <4 x float> %184)
%187 = bitcast float* %179 to <4 x float>*
store <4 x float> %185, <4 x float>* %187
%188 = bitcast float* %182 to <4 x float>*
store <4 x float> %186, <4 x float>* %188
%189 = or i64 %170, 9
%190 = getelementptr inbounds float, float* %157, i64 %189
%191 = bitcast float* %190 to <4 x float>*
%192 = load <4 x float>, <4 x float>* %191
%193 = getelementptr inbounds float, float* %190, i64 4
%194 = bitcast float* %193 to <4 x float>*
%195 = load <4 x float>, <4 x float>* %194
%196 = getelementptr inbounds float, float* %159, i64 %189
%197 = bitcast float* %196 to <4 x float>*
%198 = load <4 x float>, <4 x float>* %197
%199 = getelementptr inbounds float, float* %196, i64 4
%200 = bitcast float* %199 to <4 x float>*
%201 = load <4 x float>, <4 x float>* %200
%202 = call <4 x float> @llvm.fmuladd.v4f32(<4 x float> %145, <4 x float> %192, <4 x float> %198)
%203 = call <4 x float> @llvm.fmuladd.v4f32(<4 x float> %147, <4 x float> %195, <4 x float> %201)
%204 = bitcast float* %196 to <4 x float>*
store <4 x float> %202, <4 x float>* %204
%205 = bitcast float* %199 to <4 x float>*
store <4 x float> %203, <4 x float>* %205
%206 = add nuw i64 %170, 16
%207 = add i64 %171, 2
%208 = icmp eq i64 %207, %48
br i1 %208, label %209, label %169
209: ; preds = %169, %168
%210 = phi i64 [ 0, %168 ], [ %206, %169 ]
br i1 %49, label %229, label %211
211: ; preds = %209
%212 = or i64 %210, 1
%213 = getelementptr inbounds float, float* %157, i64 %212
%214 = bitcast float* %213 to <4 x float>*
%215 = load <4 x float>, <4 x float>* %214
%216 = getelementptr inbounds float, float* %213, i64 4
%217 = bitcast float* %216 to <4 x float>*
%218 = load <4 x float>, <4 x float>* %217
%219 = getelementptr inbounds float, float* %159, i64 %212
%220 = bitcast float* %219 to <4 x float>*
%221 = load <4 x float>, <4 x float>* %220
%222 = getelementptr inbounds float, float* %219, i64 4
%223 = bitcast float* %222 to <4 x float>*
%224 = load <4 x float>, <4 x float>* %223
%225 = call <4 x float> @llvm.fmuladd.v4f32(<4 x float> %145, <4 x float> %215, <4 x float> %221)
%226 = call <4 x float> @llvm.fmuladd.v4f32(<4 x float> %147, <4 x float> %218, <4 x float> %224)
%227 = bitcast float* %219 to <4 x float>*
store <4 x float> %225, <4 x float>* %227
%228 = bitcast float* %222 to <4 x float>*
store <4 x float> %226, <4 x float>* %228
br label %229
229: ; preds = %209, %211
br i1 %50, label %258, label %230
230: ; preds = %160, %150, %229
%231 = phi i64 [ 1, %160 ], [ 1, %150 ], [ %45, %229 ]
%232 = xor i64 %231, -1
br i1 %52, label %233, label %240
233: ; preds = %230
%234 = getelementptr inbounds float, float* %157, i64 %231
%235 = load float, float* %234
%236 = getelementptr inbounds float, float* %159, i64 %231
%237 = load float, float* %236
%238 = tail call float @llvm.fmuladd.f32(float %143, float %235, float %237)
store float %238, float* %236
%239 = add nuw nsw i64 %231, 1
br label %240
240: ; preds = %233, %230
%241 = phi i64 [ %231, %230 ], [ %239, %233 ]
%242 = icmp eq i64 %232, %53
br i1 %242, label %258, label %243
243: ; preds = %240, %243
%244 = phi i64 [ %256, %243 ], [ %241, %240 ]
%245 = getelementptr inbounds float, float* %157, i64 %244
%246 = load float, float* %245
%247 = getelementptr inbounds float, float* %159, i64 %244
%248 = load float, float* %247
%249 = tail call float @llvm.fmuladd.f32(float %143, float %246, float %248)
store float %249, float* %247
%250 = add nuw nsw i64 %244, 1
%251 = getelementptr inbounds float, float* %157, i64 %250
%252 = load float, float* %251
%253 = getelementptr inbounds float, float* %159, i64 %250
%254 = load float, float* %253
%255 = tail call float @llvm.fmuladd.f32(float %143, float %252, float %254)
store float %255, float* %253
%256 = add nuw nsw i64 %244, 2
%257 = icmp eq i64 %256, %33
br i1 %257, label %258, label %243
258: ; preds = %240, %243, %229, %148
%259 = add nuw nsw i64 %149, 1
%260 = icmp eq i64 %259, %32
br i1 %260, label %261, label %148
261: ; preds = %258, %139
%262 = tail call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([14 x i8], [14 x i8]* @.str, i64 0, i64 0), i8* getelementptr inbounds ([6 x i8], [6 x i8]* @.str.1, i64 0, i64 0), i8* getelementptr inbounds ([6 x i8], [6 x i8]* @.str.2, i64 0, i64 0), i8* getelementptr inbounds ([4 x i8], [4 x i8]* @.str.3, i64 0, i64 0))
%263 = fpext float %134 to double
%264 = fpext float %141 to double
%265 = tail call i32 (i8*, ...) @printf(i8* dereferenceable(1) getelementptr inbounds ([19 x i8], [19 x i8]* @.str.4, i64 0, i64 0), i32 %58, double %263, double %264)
%266 = fcmp olt float %134, %1
br i1 %266, label %267, label %54
267: ; preds = %261
%268 = tail call i32 @free_vector(float* %13, i32 1, i32 %5)
%269 = tail call i32 @free_ivector(i32* %12, i32 1, i32 %5)
br label %272
270: ; preds = %54, %11
%271 = tail call i32 @nrerror(i8* getelementptr inbounds ([30 x i8], [30 x i8]* @.str.5, i64 0, i64 0))
br label %272
272: ; preds = %270, %267
ret void
}