#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 */
double fabs(double);
double sqrt(double);
int jacobi(double a[], double v[], int l, int m, int *nr, double eps) {
  int n, i, j, r, c, k1, k2, s1, s2, s3, s5, s6;
  double wmax, w, a1, a2, a3, t1, ta, si, co;
  if (m < 2 || *nr < 1 || eps <= 0.0)
    return 999;
  n = 0;
  for (i = 0; i < m - 1; i++) {
    k1 = i * l;
    for (j = i + 1; j < m; j++)
      if (a[k1 + j] != a[j * l + i])
        return 999;
  }
  for (i = 0; i < m; i++) {
    k1 = i * l;
    for (j = 0; j < m; j++)
      if (i != j)
        v[k1 + j] = 0.0;
    v[k1 + i] = 1.0;
  }
  while (1) {
    wmax = 0.0;
    for (i = 0; i < m; i++) {
      k1 = i * l;
      for (j = i + 1; j < m; j++) {
        w = fabs(a[k1 + j]);
        if (w > wmax) {
          wmax = w;
          r = i;
          c = j;
        }
      }
    }
    if (wmax <= eps) {
      *nr = n;
      return 0;
    }
    if (n >= *nr) {
      *nr = n;
      return 1;
    }
    k1 = r * l;
    k2 = c * l;
    s1 = k1 + r;
    s2 = k2 + c;
    s3 = k1 + c;
    n++;
    a1 = a[s1];
    a2 = a[s2];
    a3 = a[s3];
    t1 = fabs(a1 - a2);
    ta = 2.0 * a3 / (t1 + sqrt(t1 * t1 + 4.0 * a3 * a3));
    if (a1 < a2)
      ta = -ta;
    co = sqrt(1.0 / (ta * ta + 1.0));
    si = ta * co;
    for (i = 0; i < m; i++) {
      s5 = i * l + r;
      s6 = i * l + c;
      w = v[s5];
      v[s5] = w * co + v[s6] * si;
      v[s6] = -w * si + v[s6] * co;
      if (i != r && i != c) {
        w = a[s5];
        a[s5] = w * co + a[s6] * si;
        a[s6] = -w * si + a[s6] * co;
        a[k1 + i] = a[s5];
        a[k2 + i] = a[s6];
      }
    }
    a[s1] = a1 * co * co + a2 * si * si + 2.0 * a3 * co * si;
    a[s2] = a1 + a2 - a[s1];
    a[s3] = 0.0;
    a[k2 + r] = 0.0;
  }
}


x86 assembly (clang, -O0):
<jacobi>:
push   %rbp
mov    %rsp,%rbp
sub    $0xc0,%rsp
mov    %rdi,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %edx,-0x1c(%rbp)
mov    %ecx,-0x20(%rbp)
mov    %r8,-0x28(%rbp)
movsd  %xmm0,-0x30(%rbp)
cmpl   $0x2,-0x20(%rbp)
jl     1157 <jacobi+0x47>
mov    -0x28(%rbp),%rax
cmpl   $0x1,(%rax)
jl     1157 <jacobi+0x47>
xorps  %xmm0,%xmm0
ucomisd -0x30(%rbp),%xmm0
jb     1163 <jacobi+0x53>
movl   $0x3e7,-0x4(%rbp)
jmp    17a5 <jacobi+0x695>
movl   $0x0,-0x34(%rbp)
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
mov    -0x20(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    120c <jacobi+0xfc>
mov    -0x38(%rbp),%eax
imul   -0x1c(%rbp),%eax
mov    %eax,-0x48(%rbp)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x20(%rbp),%eax
jge    11f9 <jacobi+0xe9>
mov    -0x10(%rbp),%rax
mov    -0x48(%rbp),%ecx
add    -0x3c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x3c(%rbp),%ecx
imul   -0x1c(%rbp),%ecx
add    -0x38(%rbp),%ecx
movslq %ecx,%rcx
ucomisd (%rax,%rcx,8),%xmm0
jne    11da <jacobi+0xca>
jp     11da <jacobi+0xca>
jmp    11e6 <jacobi+0xd6>
movl   $0x3e7,-0x4(%rbp)
jmp    17a5 <jacobi+0x695>
jmp    11eb <jacobi+0xdb>
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    1195 <jacobi+0x85>
jmp    11fe <jacobi+0xee>
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    1171 <jacobi+0x61>
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x20(%rbp),%eax
jge    1298 <jacobi+0x188>
mov    -0x38(%rbp),%eax
imul   -0x1c(%rbp),%eax
mov    %eax,-0x48(%rbp)
movl   $0x0,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x20(%rbp),%eax
jge    1270 <jacobi+0x160>
mov    -0x38(%rbp),%eax
cmp    -0x3c(%rbp),%eax
je     125d <jacobi+0x14d>
mov    -0x18(%rbp),%rax
mov    -0x48(%rbp),%ecx
add    -0x3c(%rbp),%ecx
movslq %ecx,%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    1262 <jacobi+0x152>
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    1230 <jacobi+0x120>
mov    -0x18(%rbp),%rax
mov    -0x48(%rbp),%ecx
add    -0x38(%rbp),%ecx
movslq %ecx,%rcx
movsd  0xd8b(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    1213 <jacobi+0x103>
jmp    129d <jacobi+0x18d>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x68(%rbp)
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x20(%rbp),%eax
jge    1345 <jacobi+0x235>
mov    -0x38(%rbp),%eax
imul   -0x1c(%rbp),%eax
mov    %eax,-0x48(%rbp)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x20(%rbp),%eax
jge    1332 <jacobi+0x222>
mov    -0x10(%rbp),%rax
mov    -0x48(%rbp),%ecx
add    -0x3c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
movaps 0xd30(%rip),%xmm1
pand   %xmm1,%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x70(%rbp),%xmm0
ucomisd -0x68(%rbp),%xmm0
jbe    131f <jacobi+0x20f>
movsd  -0x70(%rbp),%xmm0
movsd  %xmm0,-0x68(%rbp)
mov    -0x38(%rbp),%eax
mov    %eax,-0x40(%rbp)
mov    -0x3c(%rbp),%eax
mov    %eax,-0x44(%rbp)
jmp    1324 <jacobi+0x214>
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    12cb <jacobi+0x1bb>
jmp    1337 <jacobi+0x227>
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    12ac <jacobi+0x19c>
movsd  -0x68(%rbp),%xmm1
movsd  -0x30(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jb     136e <jacobi+0x25e>
mov    -0x34(%rbp),%ecx
mov    -0x28(%rbp),%rax
mov    %ecx,(%rax)
movl   $0x0,-0x4(%rbp)
jmp    17a5 <jacobi+0x695>
mov    -0x34(%rbp),%eax
mov    -0x28(%rbp),%rcx
cmp    (%rcx),%eax
jl     1392 <jacobi+0x282>
mov    -0x34(%rbp),%ecx
mov    -0x28(%rbp),%rax
mov    %ecx,(%rax)
movl   $0x1,-0x4(%rbp)
jmp    17a5 <jacobi+0x695>
mov    -0x40(%rbp),%eax
imul   -0x1c(%rbp),%eax
mov    %eax,-0x48(%rbp)
mov    -0x44(%rbp),%eax
imul   -0x1c(%rbp),%eax
mov    %eax,-0x4c(%rbp)
mov    -0x48(%rbp),%eax
add    -0x40(%rbp),%eax
mov    %eax,-0x50(%rbp)
mov    -0x4c(%rbp),%eax
add    -0x44(%rbp),%eax
mov    %eax,-0x54(%rbp)
mov    -0x48(%rbp),%eax
add    -0x44(%rbp),%eax
mov    %eax,-0x58(%rbp)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x50(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x78(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x54(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x80(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x58(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
movsd  -0x78(%rbp),%xmm0
subsd  -0x80(%rbp),%xmm0
movaps 0xc0c(%rip),%xmm1
pand   %xmm1,%xmm0
movsd  %xmm0,-0x90(%rbp)
ff
movsd  0xbe0(%rip),%xmm0
mulsd  -0x88(%rbp),%xmm0
ff
movsd  %xmm0,-0xb0(%rbp)
ff
movsd  -0x90(%rbp),%xmm0
ff
movsd  %xmm0,-0xb8(%rbp)
ff
movsd  -0x90(%rbp),%xmm0
ff
movsd  -0x90(%rbp),%xmm2
ff
movsd  0xba0(%rip),%xmm1
mulsd  -0x88(%rbp),%xmm1
ff
mulsd  -0x88(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
call   1030 <sqrt@plt>
movsd  -0xb8(%rbp),%xmm1
ff
movaps %xmm0,%xmm2
movsd  -0xb0(%rbp),%xmm0
ff
addsd  %xmm2,%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
movsd  -0x78(%rbp),%xmm1
movsd  -0x80(%rbp),%xmm0
ucomisd %xmm1,%xmm0
jbe    14db <jacobi+0x3cb>
movsd  -0x98(%rbp),%xmm0
ff
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  %xmm0,-0x98(%rbp)
ff
movsd  -0x98(%rbp),%xmm1
ff
movsd  -0x98(%rbp),%xmm0
ff
mulsd  %xmm0,%xmm1
movsd  0xb19(%rip),%xmm0
addsd  %xmm0,%xmm1
movsd  0xb0d(%rip),%xmm0
divsd  %xmm1,%xmm0
call   1030 <sqrt@plt>
movsd  %xmm0,-0xa8(%rbp)
ff
movsd  -0x98(%rbp),%xmm0
ff
mulsd  -0xa8(%rbp),%xmm0
ff
movsd  %xmm0,-0xa0(%rbp)
ff
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x20(%rbp),%eax
jge    16f0 <jacobi+0x5e0>
mov    -0x38(%rbp),%eax
imul   -0x1c(%rbp),%eax
add    -0x40(%rbp),%eax
mov    %eax,-0x5c(%rbp)
mov    -0x38(%rbp),%eax
imul   -0x1c(%rbp),%eax
add    -0x44(%rbp),%eax
mov    %eax,-0x60(%rbp)
mov    -0x18(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x70(%rbp),%xmm0
movsd  -0xa8(%rbp),%xmm2
ff
mov    -0x18(%rbp),%rax
movslq -0x60(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  -0xa0(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x70(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  -0xa0(%rbp),%xmm2
ff
mov    -0x18(%rbp),%rax
movslq -0x60(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  -0xa8(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
movslq -0x60(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x38(%rbp),%eax
cmp    -0x40(%rbp),%eax
je     16dd <jacobi+0x5cd>
mov    -0x38(%rbp),%eax
cmp    -0x44(%rbp),%eax
je     16dd <jacobi+0x5cd>
mov    -0x10(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x70(%rbp),%xmm0
movsd  -0xa8(%rbp),%xmm2
ff
mov    -0x10(%rbp),%rax
movslq -0x60(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  -0xa0(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x70(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  -0xa0(%rbp),%xmm2
ff
mov    -0x10(%rbp),%rax
movslq -0x60(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  -0xa8(%rbp),%xmm1
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x60(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%rax
movslq -0x5c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x48(%rbp),%ecx
add    -0x38(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%rax
movslq -0x60(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x4c(%rbp),%ecx
add    -0x38(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
jmp    16e2 <jacobi+0x5d2>
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    1533 <jacobi+0x423>
movsd  -0x78(%rbp),%xmm1
mulsd  -0xa8(%rbp),%xmm1
ff
movsd  -0xa8(%rbp),%xmm2
ff
movsd  -0x80(%rbp),%xmm0
mulsd  -0xa0(%rbp),%xmm0
ff
mulsd  -0xa0(%rbp),%xmm0
ff
mulsd  %xmm2,%xmm1
addsd  %xmm0,%xmm1
movsd  0x8de(%rip),%xmm0
mulsd  -0x88(%rbp),%xmm0
ff
mulsd  -0xa8(%rbp),%xmm0
ff
movsd  -0xa0(%rbp),%xmm2
ff
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x50(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x78(%rbp),%xmm0
addsd  -0x80(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x50(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x54(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%rax
movslq -0x58(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%rax
mov    -0x4c(%rbp),%ecx
add    -0x40(%rbp),%ecx
movslq %ecx,%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    129d <jacobi+0x18d>
mov    -0x4(%rbp),%eax
add    $0xc0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<jacobi>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0xa8,%rsp
mov    $0x3e7,%eax
cmp    $0x2,%ecx
jge    113d <jacobi+0x2d>
add    $0xa8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
movapd %xmm0,%xmm9
xorpd  %xmm0,%xmm0
ucomisd %xmm9,%xmm0
jae    112b <jacobi+0x1b>
cmpl   $0x0,(%r8)
jle    112b <jacobi+0x1b>
lea    -0x1(%rcx),%r10d
movslq %edx,%rbp
mov    %ecx,%r14d
lea    (%rdi,%rbp,8),%r15
mov    %rdi,%rbx
lea    0x0(,%rbp,8),%r11
mov    %rbp,0x58(%rsp)
lea    0x8(,%rbp,8),%rdi
mov    %rdi,0x28(%rsp)
xor    %r9d,%r9d
mov    %r14,%r12
mov    %rbx,0x20(%rsp)
jmp    11a4 <jacobi+0x94>
nopl   0x0(%rax,%rax,1)
mov    0x28(%rsp),%rdi
add    %rdi,%r15
add    %rdi,%rbx
add    $0xffffffffffffffff,%r12
cmp    %r10,%r9
je     11e4 <jacobi+0xd4>
add    $0x1,%r9
cmp    %r14,%r9
jae    1190 <jacobi+0x80>
mov    $0x1,%edi
mov    %r15,%rbp
cs nopw 0x0(%rax,%rax,1)
nop
movsd  (%rbx,%rdi,8),%xmm0
ucomisd 0x0(%rbp),%xmm0
jne    112b <jacobi+0x1b>
jp     112b <jacobi+0x1b>
add    %r11,%rbp
add    $0x1,%rdi
cmp    %rdi,%r12
jne    11c0 <jacobi+0xb0>
jmp    1190 <jacobi+0x80>
test   %ecx,%ecx
jle    129d <jacobi+0x18d>
mov    %r14d,%r10d
and    $0xfffffffe,%r10d
lea    0x8(%rsi),%r15
xor    %eax,%eax
movabs $0x3ff0000000000000,%r9
jmp    1226 <jacobi+0x116>
cs nopw 0x0(%rax,%rax,1)
nop
movslq %ebp,%rdi
add    %rax,%rdi
mov    %r9,(%rsi,%rdi,8)
add    $0x1,%rax
add    %r11,%r15
cmp    %r14,%rax
je     129d <jacobi+0x18d>
cmp    $0x1,%ecx
jne    1230 <jacobi+0x120>
xor    %edi,%edi
jmp    1274 <jacobi+0x164>
nop
xor    %ebp,%ebp
jmp    1249 <jacobi+0x139>
cs nopw 0x0(%rax,%rax,1)
xchg   %ax,%ax
add    $0x1,%rbp
cmp    %r10,%rbp
je     1270 <jacobi+0x160>
mov    %rbp,%rdi
cmp    %rbp,%rax
je     125a <jacobi+0x14a>
movq   $0x0,-0x8(%r15,%rdi,8)
lea    0x1(%rdi),%rbp
cmp    %rbp,%rax
je     1240 <jacobi+0x130>
movq   $0x0,(%r15,%rdi,8)
jmp    1240 <jacobi+0x130>
nopl   (%rax)
add    $0x2,%rdi
mov    %eax,%ebp
imul   %edx,%ebp
test   $0x1,%r14b
je     1210 <jacobi+0x100>
cmp    %rdi,%rax
je     1210 <jacobi+0x100>
mov    %rax,%rbx
imul   0x58(%rsp),%rbx
add    %rbx,%rdi
movq   $0x0,(%rsi,%rdi,8)
jmp    1210 <jacobi+0x100>
mov    %rsi,0x78(%rsp)
mov    %edx,0x34(%rsp)
movslq %ecx,%r12
mov    %r14,%rax
neg    %rax
mov    %rax,0xa0(%rsp)
lea    0x0(,%r14,8),%r15
xor    %eax,%eax
mov    %rax,0x8(%rsp)
movapd 0xd32(%rip),%xmm5
xorpd  %xmm8,%xmm8
mov    %r8,0x80(%rsp)
movsd  %xmm9,0x38(%rsp)
mov    %rcx,0x18(%rsp)
mov    0x20(%rsp),%r10
jmp    138b <jacobi+0x27b>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    0x8(%rsp),%rax
add    $0x1,%eax
mov    %rax,0x8(%rsp)
movapd %xmm10,%xmm1
mulsd  %xmm0,%xmm1
addsd  %xmm4,%xmm10
mulsd  %xmm6,%xmm4
mulsd  %xmm6,%xmm4
mulsd  %xmm0,%xmm1
addsd  %xmm4,%xmm1
mulsd  %xmm0,%xmm7
mulsd  %xmm6,%xmm7
addsd  %xmm1,%xmm7
mov    0x20(%rsp),%r10
mov    0x90(%rsp),%rax
movsd  %xmm7,(%r10,%rax,8)
subsd  %xmm7,%xmm10
mov    0x98(%rsp),%rax
movsd  %xmm10,(%r10,%rax,8)
mov    0x88(%rsp),%rax
movq   $0x0,(%r10,%rax,8)
mov    0x50(%rsp),%rax
add    %ebx,%eax
cltq
movq   $0x0,(%r10,%rax,8)
mov    0x80(%rsp),%r8
mov    0x18(%rsp),%rcx
mov    0x28(%rsp),%r11
test   %ecx,%ecx
jle    13fd <jacobi+0x2ed>
xorpd  %xmm1,%xmm1
mov    %r10,%rax
xor    %r9d,%r9d
jmp    13b0 <jacobi+0x2a0>
movapd %xmm1,%xmm0
add    %r11,%rax
movapd %xmm0,%xmm1
cmp    %r14,%r9
je     1401 <jacobi+0x2f1>
mov    %r9,%rdx
add    $0x1,%r9
cmp    %r12,%r9
jge    13a0 <jacobi+0x290>
mov    0xa0(%rsp),%rcx
add    %rdx,%rcx
mov    $0x1,%edi
nopl   0x0(%rax)
movsd  (%rax,%rdi,8),%xmm0
andpd  %xmm5,%xmm0
ucomisd %xmm1,%xmm0
lea    (%rdx,%rdi,1),%ebp
cmova  %ebp,%r13d
cmova  %edx,%ebx
maxsd  %xmm1,%xmm0
add    $0x1,%rdi
mov    %rcx,%rsi
movapd %xmm0,%xmm1
add    %rdi,%rsi
jne    13d0 <jacobi+0x2c0>
jmp    13a4 <jacobi+0x294>
xorpd  %xmm0,%xmm0
ucomisd %xmm0,%xmm9
jae    1674 <jacobi+0x564>
mov    0x8(%rsp),%rax
cmp    (%r8),%eax
jge    1678 <jacobi+0x568>
mov    %ebx,%ebp
mov    0x34(%rsp),%eax
imul   %eax,%ebp
mov    %r13d,%ecx
imul   %eax,%ecx
lea    (%rbx,%rbp,1),%eax
mov    %rcx,0x50(%rsp)
add    %r13d,%ecx
lea    0x0(%rbp,%r13,1),%edx
movslq %eax,%rsi
movsd  (%r10,%rsi,8),%xmm10
movslq %ecx,%rax
movsd  (%r10,%rax,8),%xmm4
movslq %edx,%rcx
movsd  (%r10,%rcx,8),%xmm7
movapd %xmm10,%xmm2
subsd  %xmm4,%xmm2
movapd %xmm7,%xmm1
mulsd  0xbb7(%rip),%xmm1
mulsd  %xmm7,%xmm1
movapd %xmm2,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
ucomisd %xmm8,%xmm0
movsd  %xmm10,0x48(%rsp)
mov    %rax,0x98(%rsp)
mov    %rsi,0x90(%rsp)
mov    %rcx,0x88(%rsp)
movsd  %xmm4,0x40(%rsp)
jb     14ab <jacobi+0x39b>
sqrtsd %xmm0,%xmm0
jmp    14e2 <jacobi+0x3d2>
movsd  %xmm7,0x10(%rsp)
movapd %xmm2,0x60(%rsp)
call   1030 <sqrt@plt>
movapd 0x60(%rsp),%xmm2
movsd  0x40(%rsp),%xmm4
movsd  0x10(%rsp),%xmm7
movsd  0x48(%rsp),%xmm10
xorpd  %xmm8,%xmm8
movapd 0xb1e(%rip),%xmm5
andpd  %xmm5,%xmm2
addsd  %xmm7,%xmm7
addsd  %xmm0,%xmm2
movapd %xmm7,%xmm6
divsd  %xmm2,%xmm6
movapd %xmm10,%xmm0
cmpltsd %xmm4,%xmm0
movapd %xmm0,%xmm1
andnpd %xmm6,%xmm1
xorpd  0xb00(%rip),%xmm6
andpd  %xmm0,%xmm6
orpd   %xmm1,%xmm6
movapd %xmm6,%xmm1
mulsd  %xmm6,%xmm1
movsd  0xb00(%rip),%xmm0
addsd  %xmm0,%xmm1
divsd  %xmm1,%xmm0
ucomisd %xmm8,%xmm0
jb     153d <jacobi+0x42d>
sqrtsd %xmm0,%xmm0
jmp    1574 <jacobi+0x464>
movsd  %xmm7,0x10(%rsp)
movapd %xmm6,0x60(%rsp)
call   1030 <sqrt@plt>
movapd 0x60(%rsp),%xmm6
movsd  0x40(%rsp),%xmm4
movsd  0x10(%rsp),%xmm7
movsd  0x48(%rsp),%xmm10
xorpd  %xmm8,%xmm8
movapd 0xa8c(%rip),%xmm5
movsd  0x38(%rsp),%xmm9
mov    0x18(%rsp),%rax
mulsd  %xmm0,%xmm6
test   %eax,%eax
jle    1300 <jacobi+0x1f0>
mov    %r13d,%r10d
movslq %ebp,%r9
movslq 0x50(%rsp),%rcx
movslq %ebx,%rbp
movslq %r13d,%rsi
mov    %ebp,%r11d
shl    $0x3,%r11
mov    0x78(%rsp),%rdi
lea    (%rdi,%rsi,8),%rdx
mov    0x20(%rsp),%rax
lea    (%rax,%rsi,8),%rsi
lea    (%rdi,%rbp,8),%rdi
lea    (%rax,%rbp,8),%rbp
lea    (%rax,%rcx,8),%r8
lea    (%rax,%r9,8),%r9
shl    $0x3,%r10
xor    %eax,%eax
xor    %ecx,%ecx
jmp    15e2 <jacobi+0x4d2>
add    $0x8,%rcx
add    0x58(%rsp),%rax
cmp    %rcx,%r15
je     1300 <jacobi+0x1f0>
movsd  (%rdi,%rax,8),%xmm1
movsd  (%rdx,%rax,8),%xmm2
mulsd  %xmm6,%xmm2
movapd %xmm1,%xmm3
mulsd  %xmm0,%xmm3
addsd  %xmm2,%xmm3
movsd  %xmm3,(%rdi,%rax,8)
movsd  (%rdx,%rax,8),%xmm2
mulsd  %xmm0,%xmm2
mulsd  %xmm6,%xmm1
subsd  %xmm1,%xmm2
movsd  %xmm2,(%rdx,%rax,8)
cmp    %rcx,%r11
je     15d0 <jacobi+0x4c0>
cmp    %rcx,%r10
je     15d0 <jacobi+0x4c0>
movsd  0x0(%rbp,%rax,8),%xmm1
movsd  (%rsi,%rax,8),%xmm2
mulsd  %xmm6,%xmm2
movapd %xmm1,%xmm3
mulsd  %xmm0,%xmm3
addsd  %xmm2,%xmm3
movsd  %xmm3,0x0(%rbp,%rax,8)
movsd  (%rsi,%rax,8),%xmm2
mulsd  %xmm0,%xmm2
mulsd  %xmm6,%xmm1
subsd  %xmm1,%xmm2
movsd  %xmm2,(%rsi,%rax,8)
movsd  0x0(%rbp,%rax,8),%xmm1
movsd  %xmm1,(%r9,%rcx,1)
movsd  (%rsi,%rax,8),%xmm1
movsd  %xmm1,(%r8,%rcx,1)
jmp    15d0 <jacobi+0x4c0>
xor    %eax,%eax
jmp    167d <jacobi+0x56d>
mov    $0x1,%eax
mov    0x8(%rsp),%rcx
mov    %ecx,(%r8)
jmp    112b <jacobi+0x1b>

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

LLVM IR (-O2):
define i32 @jacobi(double* nocapture %0, double* nocapture %1, i32 %2, i32 %3, i32* nocapture %4, double %5) {
%7 = icmp slt i32 %3, 2
br i1 %7, label %225, label %8
8: ; preds = %6
%9 = load i32, i32* %4
%10 = icmp slt i32 %9, 1
%11 = fcmp ole double %5, 0.000000e+00
%12 = or i1 %11, %10
br i1 %12, label %225, label %13
13: ; preds = %8
%14 = add nsw i32 %3, -1
%15 = sext i32 %2 to i64
%16 = zext i32 %3 to i64
%17 = zext i32 %14 to i64
%18 = zext i32 %3 to i64
br label %31
19: ; preds = %37, %31
%20 = add nuw nsw i64 %33, 1
%21 = icmp eq i64 %35, %17
br i1 %21, label %22, label %31
22: ; preds = %19
%23 = icmp sgt i32 %3, 0
%24 = sext i32 %2 to i64
br i1 %23, label %25, label %50
25: ; preds = %22
%26 = zext i32 %3 to i64
%27 = and i64 %18, 1
%28 = icmp eq i32 %3, 1
%29 = and i64 %18, 4294967294
%30 = icmp eq i64 %27, 0
br label %54
31: ; preds = %13, %19
%32 = phi i64 [ 0, %13 ], [ %35, %19 ]
%33 = phi i64 [ 1, %13 ], [ %20, %19 ]
%34 = mul nsw i64 %32, %15
%35 = add nuw nsw i64 %32, 1
%36 = icmp ult i64 %35, %16
br i1 %36, label %39, label %19
37: ; preds = %39
%38 = icmp eq i64 %49, %18
br i1 %38, label %19, label %39
39: ; preds = %31, %37
%40 = phi i64 [ %49, %37 ], [ %33, %31 ]
%41 = add nsw i64 %40, %34
%42 = getelementptr inbounds double, double* %0, i64 %41
%43 = load double, double* %42
%44 = mul nsw i64 %40, %15
%45 = add nsw i64 %44, %32
%46 = getelementptr inbounds double, double* %0, i64 %45
%47 = load double, double* %46
%48 = fcmp une double %43, %47
%49 = add nuw nsw i64 %40, 1
br i1 %48, label %225, label %37
50: ; preds = %83, %22
%51 = sext i32 %3 to i64
%52 = zext i32 %3 to i64
%53 = zext i32 %3 to i64
br label %89
54: ; preds = %83, %25
%55 = phi i64 [ 0, %25 ], [ %87, %83 ]
%56 = trunc i64 %55 to i32
%57 = mul nsw i64 %55, %24
%58 = mul nsw i32 %56, %2
br i1 %28, label %76, label %59
59: ; preds = %54, %72
%60 = phi i64 [ %73, %72 ], [ 0, %54 ]
%61 = phi i64 [ %74, %72 ], [ 0, %54 ]
%62 = icmp eq i64 %55, %60
br i1 %62, label %66, label %63
63: ; preds = %59
%64 = add nsw i64 %60, %57
%65 = getelementptr inbounds double, double* %1, i64 %64
store double 0.000000e+00, double* %65
br label %66
66: ; preds = %59, %63
%67 = or i64 %60, 1
%68 = icmp eq i64 %55, %67
br i1 %68, label %72, label %69
69: ; preds = %66
%70 = add nsw i64 %67, %57
%71 = getelementptr inbounds double, double* %1, i64 %70
store double 0.000000e+00, double* %71
br label %72
72: ; preds = %69, %66
%73 = add nuw nsw i64 %60, 2
%74 = add i64 %61, 2
%75 = icmp eq i64 %74, %29
br i1 %75, label %76, label %59
76: ; preds = %72, %54
%77 = phi i64 [ 0, %54 ], [ %73, %72 ]
%78 = icmp eq i64 %55, %77
%79 = select i1 %30, i1 true, i1 %78
br i1 %79, label %83, label %80
80: ; preds = %76
%81 = add nsw i64 %77, %57
%82 = getelementptr inbounds double, double* %1, i64 %81
store double 0.000000e+00, double* %82
br label %83
83: ; preds = %80, %76
%84 = sext i32 %58 to i64
%85 = add nsw i64 %55, %84
%86 = getelementptr inbounds double, double* %1, i64 %85
store double 1.000000e+00, double* %86
%87 = add nuw nsw i64 %55, 1
%88 = icmp eq i64 %87, %26
br i1 %88, label %50, label %54
89: ; preds = %50, %211
%90 = phi i32 [ %128, %211 ], [ undef, %50 ]
%91 = phi i32 [ %129, %211 ], [ undef, %50 ]
%92 = phi i32 [ %140, %211 ], [ 0, %50 ]
br i1 %23, label %99, label %126
93: ; preds = %110, %99
%94 = phi double [ %104, %99 ], [ %120, %110 ]
%95 = phi i32 [ %103, %99 ], [ %122, %110 ]
%96 = phi i32 [ %102, %99 ], [ %123, %110 ]
%97 = add nuw nsw i64 %101, 1
%98 = icmp eq i64 %106, %52
br i1 %98, label %126, label %99
99: ; preds = %89, %93
%100 = phi i64 [ %106, %93 ], [ 0, %89 ]
%101 = phi i64 [ %97, %93 ], [ 1, %89 ]
%102 = phi i32 [ %96, %93 ], [ %91, %89 ]
%103 = phi i32 [ %95, %93 ], [ %90, %89 ]
%104 = phi double [ %94, %93 ], [ 0.000000e+00, %89 ]
%105 = mul nsw i64 %100, %24
%106 = add nuw nsw i64 %100, 1
%107 = icmp slt i64 %106, %51
br i1 %107, label %108, label %93
108: ; preds = %99
%109 = trunc i64 %100 to i32
br label %110
110: ; preds = %108, %110
%111 = phi i64 [ %101, %108 ], [ %124, %110 ]
%112 = phi i32 [ %102, %108 ], [ %123, %110 ]
%113 = phi i32 [ %103, %108 ], [ %122, %110 ]
%114 = phi double [ %104, %108 ], [ %120, %110 ]
%115 = add nsw i64 %111, %105
%116 = getelementptr inbounds double, double* %0, i64 %115
%117 = load double, double* %116
%118 = tail call double @llvm.fabs.f64(double %117)
%119 = fcmp ogt double %118, %114
%120 = select i1 %119, double %118, double %114
%121 = trunc i64 %111 to i32
%122 = select i1 %119, i32 %121, i32 %113
%123 = select i1 %119, i32 %109, i32 %112
%124 = add nuw nsw i64 %111, 1
%125 = icmp eq i64 %124, %52
br i1 %125, label %93, label %110
126: ; preds = %93, %89
%127 = phi double [ 0.000000e+00, %89 ], [ %94, %93 ]
%128 = phi i32 [ %90, %89 ], [ %95, %93 ]
%129 = phi i32 [ %91, %89 ], [ %96, %93 ]
%130 = fcmp ugt double %127, %5
br i1 %130, label %131, label %223
131: ; preds = %126
%132 = load i32, i32* %4
%133 = icmp slt i32 %92, %132
br i1 %133, label %134, label %223
134: ; preds = %131
%135 = mul nsw i32 %129, %2
%136 = mul nsw i32 %128, %2
%137 = add nsw i32 %135, %129
%138 = add nsw i32 %136, %128
%139 = add nsw i32 %135, %128
%140 = add nuw nsw i32 %92, 1
%141 = sext i32 %137 to i64
%142 = getelementptr inbounds double, double* %0, i64 %141
%143 = load double, double* %142
%144 = sext i32 %138 to i64
%145 = getelementptr inbounds double, double* %0, i64 %144
%146 = load double, double* %145
%147 = sext i32 %139 to i64
%148 = getelementptr inbounds double, double* %0, i64 %147
%149 = load double, double* %148
%150 = fsub double %143, %146
%151 = tail call double @llvm.fabs.f64(double %150)
%152 = fmul double %149, 2.000000e+00
%153 = fmul double %149, 4.000000e+00
%154 = fmul double %149, %153
%155 = tail call double @llvm.fmuladd.f64(double %150, double %150, double %154)
%156 = tail call double @sqrt(double %155)
%157 = fadd double %151, %156
%158 = fdiv double %152, %157
%159 = fcmp olt double %143, %146
%160 = fneg double %158
%161 = select i1 %159, double %160, double %158
%162 = tail call double @llvm.fmuladd.f64(double %161, double %161, double 1.000000e+00)
%163 = fdiv double 1.000000e+00, %162
%164 = tail call double @sqrt(double %163)
%165 = fmul double %164, %161
br i1 %23, label %166, label %211
166: ; preds = %134
%167 = zext i32 %128 to i64
%168 = sext i32 %135 to i64
%169 = sext i32 %136 to i64
%170 = sext i32 %129 to i64
%171 = sext i32 %128 to i64
%172 = zext i32 %129 to i64
br label %173
173: ; preds = %166, %208
%174 = phi i64 [ 0, %166 ], [ %209, %208 ]
%175 = mul nsw i64 %174, %24
%176 = add nsw i64 %175, %170
%177 = add nsw i64 %175, %171
%178 = getelementptr inbounds double, double* %1, i64 %176
%179 = load double, double* %178
%180 = getelementptr inbounds double, double* %1, i64 %177
%181 = load double, double* %180
%182 = fmul double %165, %181
%183 = tail call double @llvm.fmuladd.f64(double %179, double %164, double %182)
store double %183, double* %178
%184 = fneg double %179
%185 = load double, double* %180
%186 = fmul double %164, %185
%187 = tail call double @llvm.fmuladd.f64(double %184, double %165, double %186)
store double %187, double* %180
%188 = icmp eq i64 %174, %172
%189 = icmp eq i64 %174, %167
%190 = select i1 %188, i1 true, i1 %189
br i1 %190, label %208, label %191
191: ; preds = %173
%192 = getelementptr inbounds double, double* %0, i64 %176
%193 = load double, double* %192
%194 = getelementptr inbounds double, double* %0, i64 %177
%195 = load double, double* %194
%196 = fmul double %165, %195
%197 = tail call double @llvm.fmuladd.f64(double %193, double %164, double %196)
store double %197, double* %192
%198 = fneg double %193
%199 = load double, double* %194
%200 = fmul double %164, %199
%201 = tail call double @llvm.fmuladd.f64(double %198, double %165, double %200)
store double %201, double* %194
%202 = load double, double* %192
%203 = add nsw i64 %174, %168
%204 = getelementptr inbounds double, double* %0, i64 %203
store double %202, double* %204
%205 = load double, double* %194
%206 = add nsw i64 %174, %169
%207 = getelementptr inbounds double, double* %0, i64 %206
store double %205, double* %207
br label %208
208: ; preds = %173, %191
%209 = add nuw nsw i64 %174, 1
%210 = icmp eq i64 %209, %53
br i1 %210, label %211, label %173
211: ; preds = %208, %134
%212 = fmul double %143, %164
%213 = fmul double %146, %165
%214 = fmul double %165, %213
%215 = tail call double @llvm.fmuladd.f64(double %212, double %164, double %214)
%216 = fmul double %152, %164
%217 = tail call double @llvm.fmuladd.f64(double %216, double %165, double %215)
store double %217, double* %142
%218 = fadd double %143, %146
%219 = fsub double %218, %217
store double %219, double* %145
store double 0.000000e+00, double* %148
%220 = add nsw i32 %136, %129
%221 = sext i32 %220 to i64
%222 = getelementptr inbounds double, double* %0, i64 %221
store double 0.000000e+00, double* %222
br label %89
223: ; preds = %131, %126
%224 = phi i32 [ 0, %126 ], [ 1, %131 ]
store i32 %92, i32* %4
br label %225
225: ; preds = %39, %223, %6, %8
%226 = phi i32 [ 999, %8 ], [ 999, %6 ], [ %224, %223 ], [ 999, %39 ]
ret i32 %226
}