#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);
int /*<<< orphan*/ fprintf(int /*<<< orphan*/, char *);
double sqrt(int);
int /*<<< orphan*/ stderr;
int jacobi(double (*a)[3], double d[3], double (*v)[3], int *nrot) {
  int j, iq, ip, i;
  double tresh, theta, tau, t, sm, s, h, g, c;
  static double b[3];
  static double z[3];
  for (ip = 0; ip < 3; ip++) {
    for (iq = 0; iq < 3; iq++) {
      v[ip][iq] = 0.0;
    }
    v[ip][ip] = 1.0;
  }
  for (ip = 0; ip < 3; ip++) {
    b[ip] = d[ip] = a[ip][ip];
    z[ip] = 0.0;
  }
  *nrot = 0;
  for (i = 0; i < (50); i++) {
    sm = 0.0;
    for (ip = 0; ip < 3 - 1; ip++) {
      for (iq = ip + 1; iq < 3; iq++) {
        sm += fabs(a[ip][iq]);
      }
    }
    if (sm == 0.0) {
      return (0);
    }
    if (i < 3) {
      tresh = 0.2 * sm / (3 * 3);
    } else {
      tresh = 0.0;
    }
    for (ip = 0; ip < 3 - 1; ip++) {
      for (iq = ip + 1; iq < 3; iq++) {
        g = 100.0 * fabs(a[ip][iq]);
        if (i > 3 && (fabs(d[ip]) + g) == fabs(d[ip]) &&
            (fabs(d[iq]) + g) == fabs(d[iq])) {
          a[ip][iq] = 0.0;
        } else if (fabs(a[ip][iq]) > tresh) {
          h = d[iq] - d[ip];
          if ((fabs(h) + g) == fabs(h)) {
            t = (a[ip][iq]) / h;
          } else {
            theta = 0.5 * h / (a[ip][iq]);
            t = 1.0 / (fabs(theta) + sqrt(1.0 + theta * theta));
            if (theta < 0.0) {
              t = -t;
            }
          }
          c = 1.0 / sqrt(1 + t * t);
          s = t * c;
          tau = s / (1.0 + c);
          h = t * a[ip][iq];
          z[ip] -= h;
          z[iq] += h;
          d[ip] -= h;
          d[iq] += h;
          a[ip][iq] = 0.0;
          for (j = 0; j <= ip - 1; j++) {
            {
              g = a[j][ip];
              h = a[j][iq];
              a[j][ip] = g - s * (h + g * tau);
              a[j][iq] = h + s * (g - h * tau);
            };
          }
          for (j = ip + 1; j <= iq - 1; j++) {
            {
              g = a[ip][j];
              h = a[j][iq];
              a[ip][j] = g - s * (h + g * tau);
              a[j][iq] = h + s * (g - h * tau);
            };
          }
          for (j = iq + 1; j < 3; j++) {
            {
              g = a[ip][j];
              h = a[iq][j];
              a[ip][j] = g - s * (h + g * tau);
              a[iq][j] = h + s * (g - h * tau);
            };
          }
          for (j = 0; j < 3; j++) {
            {
              g = v[j][ip];
              h = v[j][iq];
              v[j][ip] = g - s * (h + g * tau);
              v[j][iq] = h + s * (g - h * tau);
            };
          }
          ++(*nrot);
        }
      }
    }
    for (ip = 0; ip < 3; ip++) {
      b[ip] += z[ip];
      d[ip] = b[ip];
      z[ip] = 0.0;
    }
  }
  fprintf(stderr, "jacobi(): too many iterations\n");
  return (-1);
}


x86 assembly (clang, -O0):
<jacobi>:
push   %rbp
mov    %rsp,%rbp
sub    $0x90,%rsp
mov    %rdi,-0x10(%rbp)
mov    %rsi,-0x18(%rbp)
mov    %rdx,-0x20(%rbp)
mov    %rcx,-0x28(%rbp)
movl   $0x0,-0x34(%rbp)
cmpl   $0x3,-0x34(%rbp)
jge    11b4 <jacobi+0x94>
movl   $0x0,-0x30(%rbp)
cmpl   $0x3,-0x30(%rbp)
jge    1186 <jacobi+0x66>
mov    -0x20(%rbp),%rax
movslq -0x34(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movslq -0x30(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1153 <jacobi+0x33>
mov    -0x20(%rbp),%rax
movslq -0x34(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movslq -0x34(%rbp),%rcx
movsd  0xe77(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1142 <jacobi+0x22>
movl   $0x0,-0x34(%rbp)
cmpl   $0x3,-0x34(%rbp)
jge    121b <jacobi+0xfb>
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movslq -0x34(%rbp),%rcx
lea    0x2e4b(%rip),%rax
movsd  %xmm0,(%rax,%rcx,8)
movslq -0x34(%rbp),%rcx
lea    0x2e5b(%rip),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    11bb <jacobi+0x9b>
mov    -0x28(%rbp),%rax
movl   $0x0,(%rax)
movl   $0x0,-0x38(%rbp)
cmpl   $0x32,-0x38(%rbp)
jge    1b84 <jacobi+0xa64>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x60(%rbp)
movl   $0x0,-0x34(%rbp)
cmpl   $0x2,-0x34(%rbp)
jge    12b0 <jacobi+0x190>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
cmpl   $0x3,-0x30(%rbp)
jge    129d <jacobi+0x17d>
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movaps 0xdaf(%rip),%xmm1
pand   %xmm1,%xmm0
addsd  -0x60(%rbp),%xmm0
movsd  %xmm0,-0x60(%rbp)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1258 <jacobi+0x138>
jmp    12a2 <jacobi+0x182>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1245 <jacobi+0x125>
movsd  -0x60(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    12d4 <jacobi+0x1b4>
jp     12d4 <jacobi+0x1b4>
movl   $0x0,-0x4(%rbp)
jmp    1ba0 <jacobi+0xa80>
cmpl   $0x3,-0x38(%rbp)
jge    1301 <jacobi+0x1e1>
movsd  0xd22(%rip),%xmm0
mulsd  -0x60(%rbp),%xmm0
movsd  0xd0d(%rip),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x40(%rbp)
jmp    1309 <jacobi+0x1e9>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x40(%rbp)
movl   $0x0,-0x34(%rbp)
cmpl   $0x2,-0x34(%rbp)
jge    1af6 <jacobi+0x9d6>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
cmpl   $0x3,-0x30(%rbp)
jge    1ae3 <jacobi+0x9c3>
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movaps 0xce4(%rip),%xmm0
pand   %xmm0,%xmm1
movsd  0xcb8(%rip),%xmm0
mulsd  %xmm1,%xmm0
movsd  %xmm0,-0x78(%rbp)
cmpl   $0x3,-0x38(%rbp)
jle    1415 <jacobi+0x2f5>
mov    -0x18(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movaps 0xcb1(%rip),%xmm1
pand   %xmm1,%xmm0
addsd  -0x78(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movaps 0xc94(%rip),%xmm2
pand   %xmm2,%xmm1
ucomisd %xmm1,%xmm0
jne    1415 <jacobi+0x2f5>
jp     1415 <jacobi+0x2f5>
mov    -0x18(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movaps 0xc6c(%rip),%xmm1
pand   %xmm1,%xmm0
addsd  -0x78(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
movaps 0xc4f(%rip),%xmm2
pand   %xmm2,%xmm1
ucomisd %xmm1,%xmm0
jne    1415 <jacobi+0x2f5>
jp     1415 <jacobi+0x2f5>
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movslq -0x30(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
jmp    1ad0 <jacobi+0x9b0>
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movaps 0xbfc(%rip),%xmm1
pand   %xmm1,%xmm0
ucomisd -0x40(%rbp),%xmm0
jbe    1acb <jacobi+0x9ab>
mov    -0x18(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x34(%rbp),%rcx
subsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x70(%rbp),%xmm0
movaps 0xbc2(%rip),%xmm1
pand   %xmm1,%xmm0
addsd  -0x78(%rbp),%xmm0
movsd  -0x70(%rbp),%xmm1
movaps 0xbad(%rip),%xmm2
pand   %xmm2,%xmm1
ucomisd %xmm1,%xmm0
jne    14be <jacobi+0x39e>
jp     14be <jacobi+0x39e>
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
divsd  -0x70(%rbp),%xmm0
movsd  %xmm0,-0x58(%rbp)
jmp    1574 <jacobi+0x454>
movsd  0xb5a(%rip),%xmm0
mulsd  -0x70(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movslq -0x30(%rbp),%rcx
divsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x48(%rbp)
movsd  -0x48(%rbp),%xmm0
movaps 0xb3c(%rip),%xmm1
pand   %xmm1,%xmm0
movsd  %xmm0,-0x88(%rbp)
ff
movsd  -0x48(%rbp),%xmm0
movsd  -0x48(%rbp),%xmm1
mulsd  %xmm1,%xmm0
movsd  0xb02(%rip),%xmm1
addsd  %xmm1,%xmm0
cvttsd2si %xmm0,%edi
call   1040 <sqrt@plt>
movsd  -0x88(%rbp),%xmm1
ff
addsd  %xmm0,%xmm1
movsd  0xae1(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x58(%rbp)
xorps  %xmm0,%xmm0
ucomisd -0x48(%rbp),%xmm0
jbe    156f <jacobi+0x44f>
movsd  -0x58(%rbp),%xmm0
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
movsd  %xmm0,-0x58(%rbp)
jmp    1574 <jacobi+0x454>
movsd  -0x58(%rbp),%xmm0
movsd  -0x58(%rbp),%xmm1
mulsd  %xmm1,%xmm0
movsd  0xa8e(%rip),%xmm1
addsd  %xmm1,%xmm0
cvttsd2si %xmm0,%edi
call   1040 <sqrt@plt>
movaps %xmm0,%xmm1
movsd  0xa76(%rip),%xmm0
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x80(%rbp)
movsd  -0x58(%rbp),%xmm0
mulsd  -0x80(%rbp),%xmm0
movsd  %xmm0,-0x68(%rbp)
movsd  -0x68(%rbp),%xmm0
movsd  0xa51(%rip),%xmm1
addsd  -0x80(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  %xmm0,-0x50(%rbp)
movsd  -0x58(%rbp),%xmm0
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movslq -0x30(%rbp),%rcx
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x70(%rbp),%xmm1
movslq -0x34(%rbp),%rcx
lea    0x2a59(%rip),%rax
movsd  (%rax,%rcx,8),%xmm0
subsd  %xmm1,%xmm0
lea    0x2a49(%rip),%rax
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x70(%rbp),%xmm0
movslq -0x30(%rbp),%rcx
lea    0x2a34(%rip),%rax
addsd  (%rax,%rcx,8),%xmm0
lea    0x2a28(%rip),%rax
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x70(%rbp),%xmm1
mov    -0x18(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
subsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x70(%rbp),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x30(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movslq -0x30(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
mov    -0x34(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jg     179a <jacobi+0x67a>
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x78(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x78(%rbp),%xmm1
movsd  -0x68(%rbp),%xmm0
movsd  -0x70(%rbp),%xmm3
movsd  -0x78(%rbp),%xmm2
movsd  -0x50(%rbp),%xmm4
mulsd  %xmm4,%xmm2
addsd  %xmm3,%xmm2
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movslq -0x34(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x70(%rbp),%xmm1
movsd  -0x68(%rbp),%xmm0
movsd  -0x78(%rbp),%xmm3
movsd  -0x70(%rbp),%xmm2
movsd  -0x50(%rbp),%xmm4
movq   %xmm2,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm2
mulsd  %xmm4,%xmm2
addsd  %xmm3,%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movslq -0x30(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1691 <jacobi+0x571>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
mov    -0x30(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jg     18ac <jacobi+0x78c>
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x78(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x78(%rbp),%xmm1
movsd  -0x68(%rbp),%xmm0
movsd  -0x70(%rbp),%xmm3
movsd  -0x78(%rbp),%xmm2
movsd  -0x50(%rbp),%xmm4
mulsd  %xmm4,%xmm2
addsd  %xmm3,%xmm2
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movslq -0x2c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x70(%rbp),%xmm1
movsd  -0x68(%rbp),%xmm0
movsd  -0x78(%rbp),%xmm3
movsd  -0x70(%rbp),%xmm2
movsd  -0x50(%rbp),%xmm4
movq   %xmm2,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm2
mulsd  %xmm4,%xmm2
addsd  %xmm3,%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x2c(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movslq -0x30(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    17a3 <jacobi+0x683>
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
cmpl   $0x3,-0x2c(%rbp)
jge    19b7 <jacobi+0x897>
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x78(%rbp)
mov    -0x10(%rbp),%rax
movslq -0x30(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movslq -0x2c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x78(%rbp),%xmm1
movsd  -0x68(%rbp),%xmm0
movsd  -0x70(%rbp),%xmm3
movsd  -0x78(%rbp),%xmm2
movsd  -0x50(%rbp),%xmm4
mulsd  %xmm4,%xmm2
addsd  %xmm3,%xmm2
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x34(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movslq -0x2c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x70(%rbp),%xmm1
movsd  -0x68(%rbp),%xmm0
movsd  -0x78(%rbp),%xmm3
movsd  -0x70(%rbp),%xmm2
movsd  -0x50(%rbp),%xmm4
movq   %xmm2,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm2
mulsd  %xmm4,%xmm2
addsd  %xmm3,%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x10(%rbp),%rax
movslq -0x30(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movslq -0x2c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    18b5 <jacobi+0x795>
movl   $0x0,-0x2c(%rbp)
cmpl   $0x3,-0x2c(%rbp)
jge    1ac0 <jacobi+0x9a0>
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movslq -0x34(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x78(%rbp)
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x70(%rbp)
movsd  -0x78(%rbp),%xmm1
movsd  -0x68(%rbp),%xmm0
movsd  -0x70(%rbp),%xmm3
movsd  -0x78(%rbp),%xmm2
movsd  -0x50(%rbp),%xmm4
mulsd  %xmm4,%xmm2
addsd  %xmm3,%xmm2
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movslq -0x34(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movsd  -0x70(%rbp),%xmm1
movsd  -0x68(%rbp),%xmm0
movsd  -0x78(%rbp),%xmm3
movsd  -0x70(%rbp),%xmm2
movsd  -0x50(%rbp),%xmm4
movq   %xmm2,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm2
mulsd  %xmm4,%xmm2
addsd  %xmm3,%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
movslq -0x2c(%rbp),%rcx
imul   $0x18,%rcx,%rcx
add    %rcx,%rax
movslq -0x30(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    19be <jacobi+0x89e>
mov    -0x28(%rbp),%rax
mov    (%rax),%ecx
add    $0x1,%ecx
mov    %ecx,(%rax)
jmp    1ad0 <jacobi+0x9b0>
jmp    1ad5 <jacobi+0x9b5>
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1323 <jacobi+0x203>
jmp    1ae8 <jacobi+0x9c8>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1310 <jacobi+0x1f0>
movl   $0x0,-0x34(%rbp)
cmpl   $0x3,-0x34(%rbp)
jge    1b71 <jacobi+0xa51>
movslq -0x34(%rbp),%rcx
lea    0x254e(%rip),%rax
movsd  (%rax,%rcx,8),%xmm0
movslq -0x34(%rbp),%rcx
lea    0x251e(%rip),%rax
addsd  (%rax,%rcx,8),%xmm0
lea    0x2512(%rip),%rax
movsd  %xmm0,(%rax,%rcx,8)
movslq -0x34(%rbp),%rcx
lea    0x2502(%rip),%rax
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x34(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
movslq -0x34(%rbp),%rcx
lea    0x2505(%rip),%rax
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1afd <jacobi+0x9dd>
jmp    1b76 <jacobi+0xa56>
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    122c <jacobi+0x10c>
mov    0x2455(%rip),%rax
mov    (%rax),%edi
lea    0x4ac(%rip),%rsi
call   1030 <fprintf@plt>
movl   $0xffffffff,-0x4(%rbp)
mov    -0x4(%rbp),%eax
add    $0x90,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<jacobi>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x68,%rsp
mov    %rdx,%r12
mov    %rdi,%rbx
xorps  %xmm0,%xmm0
movups %xmm0,0x8(%rdx)
movabs $0x3ff0000000000000,%rax
mov    %rax,(%rdx)
movups %xmm0,0x18(%rdx)
movq   $0x0,0x28(%rdx)
mov    %rax,0x20(%rdx)
movups %xmm0,0x30(%rdx)
mov    %rax,0x40(%rdx)
movaps %xmm0,0x2f09(%rip)
movq   $0x0,0x2f0e(%rip)
movsd  (%rdi),%xmm0
movsd  %xmm0,(%rsi)
movsd  %xmm0,0x2ebe(%rip)
movsd  0x20(%rdi),%xmm0
movsd  %xmm0,0x8(%rsi)
movsd  %xmm0,0x2ebc(%rip)
movsd  0x40(%rdi),%xmm0
movsd  %xmm0,0x10(%rsi)
movsd  %xmm0,0x2eba(%rip)
mov    %rcx,0x30(%rsp)
movl   $0x0,(%rcx)
movdqa 0xe47(%rip),%xmm6
movsd  0xe7f(%rip),%xmm4
xor    %ecx,%ecx
mov    %rsi,0x28(%rsp)
jmp    1240 <jacobi+0x120>
nopw   0x0(%rax,%rax,1)
movsd  0x2e98(%rip),%xmm0
addsd  0x2e60(%rip),%xmm0
movsd  %xmm0,0x2e58(%rip)
movsd  %xmm0,(%rsi)
movupd 0x2e84(%rip),%xmm0
movsd  0x2e54(%rip),%xmm1
movhpd 0x2e5c(%rip),%xmm1
addpd  %xmm0,%xmm1
movlpd %xmm1,0x2e40(%rip)
xorpd  %xmm0,%xmm0
movapd %xmm0,0x2e54(%rip)
movhpd %xmm1,0x2e3c(%rip)
movupd %xmm1,0x8(%rsi)
movq   $0x0,0x2e4c(%rip)
add    $0x1,%ecx
cmp    $0x32,%ecx
je     16a4 <jacobi+0x584>
movdqu 0x8(%rbx),%xmm0
pand   %xmm6,%xmm0
pshufd $0xee,%xmm0,%xmm1
addsd  %xmm0,%xmm1
movq   0x28(%rbx),%xmm0
pand   %xmm6,%xmm0
addsd  %xmm1,%xmm0
ucomisd 0xdf1(%rip),%xmm0
jne    126f <jacobi+0x14f>
jnp    16c0 <jacobi+0x5a0>
xorpd  %xmm5,%xmm5
cmp    $0x2,%ecx
ja     128c <jacobi+0x16c>
mulsd  0xdb0(%rip),%xmm0
divsd  0xdb0(%rip),%xmm0
movapd %xmm0,%xmm5
mov    $0x1,%ebp
mov    $0x1,%al
xor    %edx,%edx
mov    %rdx,0x8(%rsp)
xor    %r14d,%r14d
mov    %ecx,0x4(%rsp)
movsd  %xmm5,0x38(%rsp)
jmp    12d5 <jacobi+0x1b5>
nopl   0x0(%rax)
mov    0x20(%rsp),%rbp
add    $0x1,%rbp
mov    $0x1,%r14d
mov    $0x1,%al
mov    %rax,0x8(%rsp)
testb  $0x1,(%rsp)
mov    $0x0,%eax
je     11d0 <jacobi+0xb0>
mov    %eax,(%rsp)
lea    0x1(%r14),%rax
mov    %rax,0x48(%rsp)
lea    (%r14,%r14,2),%rdx
lea    (%rbx,%rdx,8),%r13
add    $0x8,%r13
mov    %rbp,0x20(%rsp)
mov    %rdx,0x40(%rsp)
jmp    144b <jacobi+0x32b>
nopl   0x0(%rax)
lea    0x0(,%rbp,2),%rax
add    %rbp,%rax
movapd %xmm0,%xmm2
xorpd  0xcf9(%rip),%xmm2
movsd  0x8(%r13),%xmm3
movsd  0x10(%rbx,%rax,8),%xmm4
movapd %xmm3,%xmm5
mulsd  %xmm1,%xmm5
addsd  %xmm4,%xmm5
mulsd  %xmm0,%xmm5
movapd %xmm3,%xmm6
subsd  %xmm5,%xmm6
movsd  %xmm6,0x8(%r13)
movapd %xmm4,%xmm5
mulsd  %xmm1,%xmm5
subsd  %xmm5,%xmm3
mulsd  %xmm0,%xmm3
addsd  %xmm4,%xmm3
movsd  %xmm3,0x10(%rbx,%rax,8)
movapd 0xc9d(%rip),%xmm6
movsd  (%r12,%r14,8),%xmm3
movsd  (%r12,%rbp,8),%xmm4
movapd %xmm3,%xmm5
mulsd  %xmm1,%xmm5
addsd  %xmm4,%xmm5
mulsd  %xmm2,%xmm5
addsd  %xmm3,%xmm5
movsd  %xmm5,(%r12,%r14,8)
movapd %xmm4,%xmm5
mulsd  %xmm1,%xmm5
subsd  %xmm5,%xmm3
mulsd  %xmm0,%xmm3
addsd  %xmm4,%xmm3
movsd  %xmm3,(%r12,%rbp,8)
movsd  0x18(%r12,%r14,8),%xmm3
movsd  0x18(%r12,%rbp,8),%xmm4
movapd %xmm3,%xmm5
mulsd  %xmm1,%xmm5
addsd  %xmm4,%xmm5
mulsd  %xmm2,%xmm5
addsd  %xmm3,%xmm5
movsd  %xmm5,0x18(%r12,%r14,8)
movapd %xmm4,%xmm5
mulsd  %xmm1,%xmm5
subsd  %xmm5,%xmm3
mulsd  %xmm0,%xmm3
addsd  %xmm4,%xmm3
movsd  %xmm3,0x18(%r12,%rbp,8)
movsd  0x30(%r12,%r14,8),%xmm3
movsd  0x30(%r12,%rbp,8),%xmm4
movapd %xmm3,%xmm5
mulsd  %xmm1,%xmm5
addsd  %xmm4,%xmm5
mulsd  %xmm2,%xmm5
addsd  %xmm3,%xmm5
movsd  %xmm5,0x30(%r12,%r14,8)
mulsd  %xmm4,%xmm1
subsd  %xmm1,%xmm3
mulsd  %xmm0,%xmm3
addsd  %xmm4,%xmm3
movsd  %xmm3,0x30(%r12,%rbp,8)
mov    0x30(%rsp),%rax
addl   $0x1,(%rax)
movsd  0xc09(%rip),%xmm4
movsd  0x38(%rsp),%xmm5
add    $0x1,%rbp
cmp    $0x3,%rbp
je     12b0 <jacobi+0x190>
lea    (%rbx,%rdx,8),%rax
lea    (%rax,%rbp,8),%r15
movq   (%rax,%rbp,8),%xmm7
movdqa %xmm7,%xmm1
pand   %xmm6,%xmm1
movdqa %xmm1,%xmm0
mulsd  %xmm4,%xmm0
cmp    $0x4,%ecx
jb     14b0 <jacobi+0x390>
movq   (%rsi,%r14,8),%xmm2
pand   %xmm6,%xmm2
movapd %xmm0,%xmm3
addsd  %xmm2,%xmm3
ucomisd %xmm2,%xmm3
jne    14b0 <jacobi+0x390>
jp     14b0 <jacobi+0x390>
movq   (%rsi,%rbp,8),%xmm2
pand   %xmm6,%xmm2
movapd %xmm0,%xmm3
addsd  %xmm2,%xmm3
ucomisd %xmm2,%xmm3
jne    14b0 <jacobi+0x390>
jp     14b0 <jacobi+0x390>
movq   $0x0,(%r15)
jmp    143d <jacobi+0x31d>
nopl   0x0(%rax)
ucomisd %xmm5,%xmm1
jbe    143d <jacobi+0x31d>
movsd  (%rsi,%rbp,8),%xmm2
subsd  (%rsi,%r14,8),%xmm2
movapd %xmm2,%xmm1
andpd  %xmm6,%xmm1
addsd  %xmm1,%xmm0
ucomisd %xmm1,%xmm0
jne    14db <jacobi+0x3bb>
jp     14db <jacobi+0x3bb>
divsd  %xmm2,%xmm7
jmp    1540 <jacobi+0x420>
mulsd  0xb65(%rip),%xmm2
divsd  %xmm7,%xmm2
movapd %xmm2,0x50(%rsp)
movapd %xmm2,%xmm0
andpd  %xmm6,%xmm0
movapd %xmm0,0x10(%rsp)
movapd %xmm2,%xmm0
mulsd  %xmm2,%xmm0
movsd  0xb45(%rip),%xmm1
addsd  %xmm1,%xmm0
cvttsd2si %xmm0,%edi
call   1040 <sqrt@plt>
addsd  0x10(%rsp),%xmm0
movsd  0xb2a(%rip),%xmm7
divsd  %xmm0,%xmm7
xorpd  %xmm0,%xmm0
ucomisd 0x50(%rsp),%xmm0
jbe    1540 <jacobi+0x420>
xorpd  0xad2(%rip),%xmm7
xchg   %ax,%ax
movapd %xmm7,0x10(%rsp)
movapd %xmm7,%xmm0
mulsd  %xmm7,%xmm0
movsd  0xafa(%rip),%xmm1
addsd  %xmm1,%xmm0
cvttsd2si %xmm0,%edi
call   1040 <sqrt@plt>
movsd  0xae5(%rip),%xmm1
movapd %xmm1,%xmm2
divsd  %xmm0,%xmm2
movapd 0x10(%rsp),%xmm3
movapd %xmm3,%xmm0
mulsd  %xmm2,%xmm0
mulsd  (%r15),%xmm3
addsd  %xmm1,%xmm2
lea    0x2adf(%rip),%rax
movsd  (%rax,%r14,8),%xmm1
subsd  %xmm3,%xmm1
movsd  %xmm1,(%rax,%r14,8)
movsd  (%rax,%rbp,8),%xmm1
addsd  %xmm3,%xmm1
movsd  %xmm1,(%rax,%rbp,8)
mov    0x28(%rsp),%rsi
movsd  (%rsi,%r14,8),%xmm1
subsd  %xmm3,%xmm1
movsd  %xmm1,(%rsi,%r14,8)
movapd %xmm0,%xmm1
addsd  (%rsi,%rbp,8),%xmm3
divsd  %xmm2,%xmm1
movsd  %xmm3,(%rsi,%rbp,8)
movq   $0x0,(%r15)
testb  $0x1,0x8(%rsp)
je     1626 <jacobi+0x506>
movsd  (%rbx,%r14,8),%xmm2
movsd  (%rbx,%rbp,8),%xmm3
movapd %xmm2,%xmm4
mulsd  %xmm1,%xmm4
addsd  %xmm3,%xmm4
mulsd  %xmm0,%xmm4
movapd %xmm2,%xmm5
subsd  %xmm4,%xmm5
movsd  %xmm5,(%rbx,%r14,8)
movapd %xmm3,%xmm4
mulsd  %xmm1,%xmm4
subsd  %xmm4,%xmm2
mulsd  %xmm0,%xmm2
addsd  %xmm3,%xmm2
movsd  %xmm2,(%rbx,%rbp,8)
cmp    %rbp,0x48(%rsp)
jae    1680 <jacobi+0x560>
movapd %xmm0,%xmm2
xorpd  0x9d7(%rip),%xmm2
movsd  0x18(%rbx,%rbp,8),%xmm3
nop
movsd  0x0(%r13),%xmm4
movapd %xmm4,%xmm5
mulsd  %xmm1,%xmm5
addsd  %xmm3,%xmm5
mulsd  %xmm2,%xmm5
addsd  %xmm4,%xmm5
movsd  %xmm5,0x0(%r13)
movapd %xmm3,%xmm5
mulsd  %xmm1,%xmm5
subsd  %xmm5,%xmm4
mulsd  %xmm0,%xmm4
addsd  %xmm4,%xmm3
movsd  %xmm3,0x18(%rbx,%rbp,8)
cmp    $0x2,%rbp
ja     1640 <jacobi+0x520>
cmp    $0x2,%rbp
mov    0x4(%rsp),%ecx
mov    0x40(%rsp),%rdx
jb     1300 <jacobi+0x1e0>
movapd %xmm0,%xmm2
xorpd  0x971(%rip),%xmm2
jmp    135b <jacobi+0x23b>
mov    0x2935(%rip),%rax
mov    (%rax),%edi
lea    0x9ac(%rip),%rsi
call   1030 <fprintf@plt>
mov    $0xffffffff,%eax
jmp    16c2 <jacobi+0x5a2>
xor    %eax,%eax
add    $0x68,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
@jacobi.b = internal global [3 x double] zeroinitializer
@jacobi.z = internal global [3 x double] zeroinitializer
@stderr = global i32 0
@.str = [31 x i8] c"jacobi(): too many iterations\0A\00"
define i32 @jacobi([3 x double]* %0, double* %1, [3 x double]* %2, i32* %3) {
%5 = alloca i32
%6 = alloca [3 x double]*
%7 = alloca double*
%8 = alloca [3 x double]*
%9 = alloca i32*
%10 = alloca i32
%11 = alloca i32
%12 = alloca i32
%13 = alloca i32
%14 = alloca double
%15 = alloca double
%16 = alloca double
%17 = alloca double
%18 = alloca double
%19 = alloca double
%20 = alloca double
%21 = alloca double
%22 = alloca double
store [3 x double]* %0, [3 x double]** %6
store double* %1, double** %7
store [3 x double]* %2, [3 x double]** %8
store i32* %3, i32** %9
store i32 0, i32* %12
br label %23
23: ; preds = %49, %4
%24 = load i32, i32* %12
%25 = icmp slt i32 %24, 3
br i1 %25, label %26, label %52
26: ; preds = %23
store i32 0, i32* %11
br label %27
27: ; preds = %38, %26
%28 = load i32, i32* %11
%29 = icmp slt i32 %28, 3
br i1 %29, label %30, label %41
30: ; preds = %27
%31 = load [3 x double]*, [3 x double]** %8
%32 = load i32, i32* %12
%33 = sext i32 %32 to i64
%34 = getelementptr inbounds [3 x double], [3 x double]* %31, i64 %33
%35 = load i32, i32* %11
%36 = sext i32 %35 to i64
%37 = getelementptr inbounds [3 x double], [3 x double]* %34, i64 0, i64 %36
store double 0.000000e+00, double* %37
br label %38
38: ; preds = %30
%39 = load i32, i32* %11
%40 = add nsw i32 %39, 1
store i32 %40, i32* %11
br label %27
41: ; preds = %27
%42 = load [3 x double]*, [3 x double]** %8
%43 = load i32, i32* %12
%44 = sext i32 %43 to i64
%45 = getelementptr inbounds [3 x double], [3 x double]* %42, i64 %44
%46 = load i32, i32* %12
%47 = sext i32 %46 to i64
%48 = getelementptr inbounds [3 x double], [3 x double]* %45, i64 0, i64 %47
store double 1.000000e+00, double* %48
br label %49
49: ; preds = %41
%50 = load i32, i32* %12
%51 = add nsw i32 %50, 1
store i32 %51, i32* %12
br label %23
52: ; preds = %23
store i32 0, i32* %12
br label %53
53: ; preds = %75, %52
%54 = load i32, i32* %12
%55 = icmp slt i32 %54, 3
br i1 %55, label %56, label %78
56: ; preds = %53
%57 = load [3 x double]*, [3 x double]** %6
%58 = load i32, i32* %12
%59 = sext i32 %58 to i64
%60 = getelementptr inbounds [3 x double], [3 x double]* %57, i64 %59
%61 = load i32, i32* %12
%62 = sext i32 %61 to i64
%63 = getelementptr inbounds [3 x double], [3 x double]* %60, i64 0, i64 %62
%64 = load double, double* %63
%65 = load double*, double** %7
%66 = load i32, i32* %12
%67 = sext i32 %66 to i64
%68 = getelementptr inbounds double, double* %65, i64 %67
store double %64, double* %68
%69 = load i32, i32* %12
%70 = sext i32 %69 to i64
%71 = getelementptr inbounds [3 x double], [3 x double]* @jacobi.b, i64 0, i64 %70
store double %64, double* %71
%72 = load i32, i32* %12
%73 = sext i32 %72 to i64
%74 = getelementptr inbounds [3 x double], [3 x double]* @jacobi.z, i64 0, i64 %73
store double 0.000000e+00, double* %74
br label %75
75: ; preds = %56
%76 = load i32, i32* %12
%77 = add nsw i32 %76, 1
store i32 %77, i32* %12
br label %53
78: ; preds = %53
%79 = load i32*, i32** %9
store i32 0, i32* %79
store i32 0, i32* %13
br label %80
80: ; preds = %578, %78
%81 = load i32, i32* %13
%82 = icmp slt i32 %81, 50
br i1 %82, label %83, label %581
83: ; preds = %80
store double 0.000000e+00, double* %18
store i32 0, i32* %12
br label %84
84: ; preds = %109, %83
%85 = load i32, i32* %12
%86 = icmp slt i32 %85, 2
br i1 %86, label %87, label %112
87: ; preds = %84
%88 = load i32, i32* %12
%89 = add nsw i32 %88, 1
store i32 %89, i32* %11
br label %90
90: ; preds = %105, %87
%91 = load i32, i32* %11
%92 = icmp slt i32 %91, 3
br i1 %92, label %93, label %108
93: ; preds = %90
%94 = load [3 x double]*, [3 x double]** %6
%95 = load i32, i32* %12
%96 = sext i32 %95 to i64
%97 = getelementptr inbounds [3 x double], [3 x double]* %94, i64 %96
%98 = load i32, i32* %11
%99 = sext i32 %98 to i64
%100 = getelementptr inbounds [3 x double], [3 x double]* %97, i64 0, i64 %99
%101 = load double, double* %100
%102 = call double @llvm.fabs.f64(double %101)
%103 = load double, double* %18
%104 = fadd double %103, %102
store double %104, double* %18
br label %105
105: ; preds = %93
%106 = load i32, i32* %11
%107 = add nsw i32 %106, 1
store i32 %107, i32* %11
br label %90
108: ; preds = %90
br label %109
109: ; preds = %108
%110 = load i32, i32* %12
%111 = add nsw i32 %110, 1
store i32 %111, i32* %12
br label %84
112: ; preds = %84
%113 = load double, double* %18
%114 = fcmp oeq double %113, 0.000000e+00
br i1 %114, label %115, label %116
115: ; preds = %112
store i32 0, i32* %5
br label %584
116: ; preds = %112
%117 = load i32, i32* %13
%118 = icmp slt i32 %117, 3
br i1 %118, label %119, label %123
119: ; preds = %116
%120 = load double, double* %18
%121 = fmul double 2.000000e-01, %120
%122 = fdiv double %121, 9.000000e+00
store double %122, double* %14
br label %124
123: ; preds = %116
store double 0.000000e+00, double* %14
br label %124
124: ; preds = %123, %119
store i32 0, i32* %12
br label %125
125: ; preds = %546, %124
%126 = load i32, i32* %12
%127 = icmp slt i32 %126, 2
br i1 %127, label %128, label %549
128: ; preds = %125
%129 = load i32, i32* %12
%130 = add nsw i32 %129, 1
store i32 %130, i32* %11
br label %131
131: ; preds = %542, %128
%132 = load i32, i32* %11
%133 = icmp slt i32 %132, 3
br i1 %133, label %134, label %545
134: ; preds = %131
%135 = load [3 x double]*, [3 x double]** %6
%136 = load i32, i32* %12
%137 = sext i32 %136 to i64
%138 = getelementptr inbounds [3 x double], [3 x double]* %135, i64 %137
%139 = load i32, i32* %11
%140 = sext i32 %139 to i64
%141 = getelementptr inbounds [3 x double], [3 x double]* %138, i64 0, i64 %140
%142 = load double, double* %141
%143 = call double @llvm.fabs.f64(double %142)
%144 = fmul double 1.000000e+02, %143
store double %144, double* %21
%145 = load i32, i32* %13
%146 = icmp sgt i32 %145, 3
br i1 %146, label %147, label %187
147: ; preds = %134
%148 = load double*, double** %7
%149 = load i32, i32* %12
%150 = sext i32 %149 to i64
%151 = getelementptr inbounds double, double* %148, i64 %150
%152 = load double, double* %151
%153 = call double @llvm.fabs.f64(double %152)
%154 = load double, double* %21
%155 = fadd double %153, %154
%156 = load double*, double** %7
%157 = load i32, i32* %12
%158 = sext i32 %157 to i64
%159 = getelementptr inbounds double, double* %156, i64 %158
%160 = load double, double* %159
%161 = call double @llvm.fabs.f64(double %160)
%162 = fcmp oeq double %155, %161
br i1 %162, label %163, label %187
163: ; preds = %147
%164 = load double*, double** %7
%165 = load i32, i32* %11
%166 = sext i32 %165 to i64
%167 = getelementptr inbounds double, double* %164, i64 %166
%168 = load double, double* %167
%169 = call double @llvm.fabs.f64(double %168)
%170 = load double, double* %21
%171 = fadd double %169, %170
%172 = load double*, double** %7
%173 = load i32, i32* %11
%174 = sext i32 %173 to i64
%175 = getelementptr inbounds double, double* %172, i64 %174
%176 = load double, double* %175
%177 = call double @llvm.fabs.f64(double %176)
%178 = fcmp oeq double %171, %177
br i1 %178, label %179, label %187
179: ; preds = %163
%180 = load [3 x double]*, [3 x double]** %6
%181 = load i32, i32* %12
%182 = sext i32 %181 to i64
%183 = getelementptr inbounds [3 x double], [3 x double]* %180, i64 %182
%184 = load i32, i32* %11
%185 = sext i32 %184 to i64
%186 = getelementptr inbounds [3 x double], [3 x double]* %183, i64 0, i64 %185
store double 0.000000e+00, double* %186
br label %541
187: ; preds = %163, %147, %134
%188 = load [3 x double]*, [3 x double]** %6
%189 = load i32, i32* %12
%190 = sext i32 %189 to i64
%191 = getelementptr inbounds [3 x double], [3 x double]* %188, i64 %190
%192 = load i32, i32* %11
%193 = sext i32 %192 to i64
%194 = getelementptr inbounds [3 x double], [3 x double]* %191, i64 0, i64 %193
%195 = load double, double* %194
%196 = call double @llvm.fabs.f64(double %195)
%197 = load double, double* %14
%198 = fcmp ogt double %196, %197
br i1 %198, label %199, label %540
199: ; preds = %187
%200 = load double*, double** %7
%201 = load i32, i32* %11
%202 = sext i32 %201 to i64
%203 = getelementptr inbounds double, double* %200, i64 %202
%204 = load double, double* %203
%205 = load double*, double** %7
%206 = load i32, i32* %12
%207 = sext i32 %206 to i64
%208 = getelementptr inbounds double, double* %205, i64 %207
%209 = load double, double* %208
%210 = fsub double %204, %209
store double %210, double* %20
%211 = load double, double* %20
%212 = call double @llvm.fabs.f64(double %211)
%213 = load double, double* %21
%214 = fadd double %212, %213
%215 = load double, double* %20
%216 = call double @llvm.fabs.f64(double %215)
%217 = fcmp oeq double %214, %216
br i1 %217, label %218, label %229
218: ; preds = %199
%219 = load [3 x double]*, [3 x double]** %6
%220 = load i32, i32* %12
%221 = sext i32 %220 to i64
%222 = getelementptr inbounds [3 x double], [3 x double]* %219, i64 %221
%223 = load i32, i32* %11
%224 = sext i32 %223 to i64
%225 = getelementptr inbounds [3 x double], [3 x double]* %222, i64 0, i64 %224
%226 = load double, double* %225
%227 = load double, double* %20
%228 = fdiv double %226, %227
store double %228, double* %17
br label %256
229: ; preds = %199
%230 = load double, double* %20
%231 = fmul double 5.000000e-01, %230
%232 = load [3 x double]*, [3 x double]** %6
%233 = load i32, i32* %12
%234 = sext i32 %233 to i64
%235 = getelementptr inbounds [3 x double], [3 x double]* %232, i64 %234
%236 = load i32, i32* %11
%237 = sext i32 %236 to i64
%238 = getelementptr inbounds [3 x double], [3 x double]* %235, i64 0, i64 %237
%239 = load double, double* %238
%240 = fdiv double %231, %239
store double %240, double* %15
%241 = load double, double* %15
%242 = call double @llvm.fabs.f64(double %241)
%243 = load double, double* %15
%244 = load double, double* %15
%245 = call double @llvm.fmuladd.f64(double %243, double %244, double 1.000000e+00)
%246 = fptosi double %245 to i32
%247 = call double @sqrt(i32 %246)
%248 = fadd double %242, %247
%249 = fdiv double 1.000000e+00, %248
store double %249, double* %17
%250 = load double, double* %15
%251 = fcmp olt double %250, 0.000000e+00
br i1 %251, label %252, label %255
252: ; preds = %229
%253 = load double, double* %17
%254 = fneg double %253
store double %254, double* %17
br label %255
255: ; preds = %252, %229
br label %256
256: ; preds = %255, %218
%257 = load double, double* %17
%258 = load double, double* %17
%259 = call double @llvm.fmuladd.f64(double %257, double %258, double 1.000000e+00)
%260 = fptosi double %259 to i32
%261 = call double @sqrt(i32 %260)
%262 = fdiv double 1.000000e+00, %261
store double %262, double* %22
%263 = load double, double* %17
%264 = load double, double* %22
%265 = fmul double %263, %264
store double %265, double* %19
%266 = load double, double* %19
%267 = load double, double* %22
%268 = fadd double 1.000000e+00, %267
%269 = fdiv double %266, %268
store double %269, double* %16
%270 = load double, double* %17
%271 = load [3 x double]*, [3 x double]** %6
%272 = load i32, i32* %12
%273 = sext i32 %272 to i64
%274 = getelementptr inbounds [3 x double], [3 x double]* %271, i64 %273
%275 = load i32, i32* %11
%276 = sext i32 %275 to i64
%277 = getelementptr inbounds [3 x double], [3 x double]* %274, i64 0, i64 %276
%278 = load double, double* %277
%279 = fmul double %270, %278
store double %279, double* %20
%280 = load double, double* %20
%281 = load i32, i32* %12
%282 = sext i32 %281 to i64
%283 = getelementptr inbounds [3 x double], [3 x double]* @jacobi.z, i64 0, i64 %282
%284 = load double, double* %283
%285 = fsub double %284, %280
store double %285, double* %283
%286 = load double, double* %20
%287 = load i32, i32* %11
%288 = sext i32 %287 to i64
%289 = getelementptr inbounds [3 x double], [3 x double]* @jacobi.z, i64 0, i64 %288
%290 = load double, double* %289
%291 = fadd double %290, %286
store double %291, double* %289
%292 = load double, double* %20
%293 = load double*, double** %7
%294 = load i32, i32* %12
%295 = sext i32 %294 to i64
%296 = getelementptr inbounds double, double* %293, i64 %295
%297 = load double, double* %296
%298 = fsub double %297, %292
store double %298, double* %296
%299 = load double, double* %20
%300 = load double*, double** %7
%301 = load i32, i32* %11
%302 = sext i32 %301 to i64
%303 = getelementptr inbounds double, double* %300, i64 %302
%304 = load double, double* %303
%305 = fadd double %304, %299
store double %305, double* %303
%306 = load [3 x double]*, [3 x double]** %6
%307 = load i32, i32* %12
%308 = sext i32 %307 to i64
%309 = getelementptr inbounds [3 x double], [3 x double]* %306, i64 %308
%310 = load i32, i32* %11
%311 = sext i32 %310 to i64
%312 = getelementptr inbounds [3 x double], [3 x double]* %309, i64 0, i64 %311
store double 0.000000e+00, double* %312
store i32 0, i32* %10
br label %313
313: ; preds = %365, %256
%314 = load i32, i32* %10
%315 = load i32, i32* %12
%316 = sub nsw i32 %315, 1
%317 = icmp sle i32 %314, %316
br i1 %317, label %318, label %368
318: ; preds = %313
%319 = load [3 x double]*, [3 x double]** %6
%320 = load i32, i32* %10
%321 = sext i32 %320 to i64
%322 = getelementptr inbounds [3 x double], [3 x double]* %319, i64 %321
%323 = load i32, i32* %12
%324 = sext i32 %323 to i64
%325 = getelementptr inbounds [3 x double], [3 x double]* %322, i64 0, i64 %324
%326 = load double, double* %325
store double %326, double* %21
%327 = load [3 x double]*, [3 x double]** %6
%328 = load i32, i32* %10
%329 = sext i32 %328 to i64
%330 = getelementptr inbounds [3 x double], [3 x double]* %327, i64 %329
%331 = load i32, i32* %11
%332 = sext i32 %331 to i64
%333 = getelementptr inbounds [3 x double], [3 x double]* %330, i64 0, i64 %332
%334 = load double, double* %333
store double %334, double* %20
%335 = load double, double* %21
%336 = load double, double* %19
%337 = load double, double* %20
%338 = load double, double* %21
%339 = load double, double* %16
%340 = call double @llvm.fmuladd.f64(double %338, double %339, double %337)
%341 = fneg double %336
%342 = call double @llvm.fmuladd.f64(double %341, double %340, double %335)
%343 = load [3 x double]*, [3 x double]** %6
%344 = load i32, i32* %10
%345 = sext i32 %344 to i64
%346 = getelementptr inbounds [3 x double], [3 x double]* %343, i64 %345
%347 = load i32, i32* %12
%348 = sext i32 %347 to i64
%349 = getelementptr inbounds [3 x double], [3 x double]* %346, i64 0, i64 %348
store double %342, double* %349
%350 = load double, double* %20
%351 = load double, double* %19
%352 = load double, double* %21
%353 = load double, double* %20
%354 = load double, double* %16
%355 = fneg double %353
%356 = call double @llvm.fmuladd.f64(double %355, double %354, double %352)
%357 = call double @llvm.fmuladd.f64(double %351, double %356, double %350)
%358 = load [3 x double]*, [3 x double]** %6
%359 = load i32, i32* %10
%360 = sext i32 %359 to i64
%361 = getelementptr inbounds [3 x double], [3 x double]* %358, i64 %360
%362 = load i32, i32* %11
%363 = sext i32 %362 to i64
%364 = getelementptr inbounds [3 x double], [3 x double]* %361, i64 0, i64 %363
store double %357, double* %364
br label %365
365: ; preds = %318
%366 = load i32, i32* %10
%367 = add nsw i32 %366, 1
store i32 %367, i32* %10
br label %313
368: ; preds = %313
%369 = load i32, i32* %12
%370 = add nsw i32 %369, 1
store i32 %370, i32* %10
br label %371
371: ; preds = %423, %368
%372 = load i32, i32* %10
%373 = load i32, i32* %11
%374 = sub nsw i32 %373, 1
%375 = icmp sle i32 %372, %374
br i1 %375, label %376, label %426
376: ; preds = %371
%377 = load [3 x double]*, [3 x double]** %6
%378 = load i32, i32* %12
%379 = sext i32 %378 to i64
%380 = getelementptr inbounds [3 x double], [3 x double]* %377, i64 %379
%381 = load i32, i32* %10
%382 = sext i32 %381 to i64
%383 = getelementptr inbounds [3 x double], [3 x double]* %380, i64 0, i64 %382
%384 = load double, double* %383
store double %384, double* %21
%385 = load [3 x double]*, [3 x double]** %6
%386 = load i32, i32* %10
%387 = sext i32 %386 to i64
%388 = getelementptr inbounds [3 x double], [3 x double]* %385, i64 %387
%389 = load i32, i32* %11
%390 = sext i32 %389 to i64
%391 = getelementptr inbounds [3 x double], [3 x double]* %388, i64 0, i64 %390
%392 = load double, double* %391
store double %392, double* %20
%393 = load double, double* %21
%394 = load double, double* %19
%395 = load double, double* %20
%396 = load double, double* %21
%397 = load double, double* %16
%398 = call double @llvm.fmuladd.f64(double %396, double %397, double %395)
%399 = fneg double %394
%400 = call double @llvm.fmuladd.f64(double %399, double %398, double %393)
%401 = load [3 x double]*, [3 x double]** %6
%402 = load i32, i32* %12
%403 = sext i32 %402 to i64
%404 = getelementptr inbounds [3 x double], [3 x double]* %401, i64 %403
%405 = load i32, i32* %10
%406 = sext i32 %405 to i64
%407 = getelementptr inbounds [3 x double], [3 x double]* %404, i64 0, i64 %406
store double %400, double* %407
%408 = load double, double* %20
%409 = load double, double* %19
%410 = load double, double* %21
%411 = load double, double* %20
%412 = load double, double* %16
%413 = fneg double %411
%414 = call double @llvm.fmuladd.f64(double %413, double %412, double %410)
%415 = call double @llvm.fmuladd.f64(double %409, double %414, double %408)
%416 = load [3 x double]*, [3 x double]** %6
%417 = load i32, i32* %10
%418 = sext i32 %417 to i64
%419 = getelementptr inbounds [3 x double], [3 x double]* %416, i64 %418
%420 = load i32, i32* %11
%421 = sext i32 %420 to i64
%422 = getelementptr inbounds [3 x double], [3 x double]* %419, i64 0, i64 %421
store double %415, double* %422
br label %423
423: ; preds = %376
%424 = load i32, i32* %10
%425 = add nsw i32 %424, 1
store i32 %425, i32* %10
br label %371
426: ; preds = %371
%427 = load i32, i32* %11
%428 = add nsw i32 %427, 1
store i32 %428, i32* %10
br label %429
429: ; preds = %479, %426
%430 = load i32, i32* %10
%431 = icmp slt i32 %430, 3
br i1 %431, label %432, label %482
432: ; preds = %429
%433 = load [3 x double]*, [3 x double]** %6
%434 = load i32, i32* %12
%435 = sext i32 %434 to i64
%436 = getelementptr inbounds [3 x double], [3 x double]* %433, i64 %435
%437 = load i32, i32* %10
%438 = sext i32 %437 to i64
%439 = getelementptr inbounds [3 x double], [3 x double]* %436, i64 0, i64 %438
%440 = load double, double* %439
store double %440, double* %21
%441 = load [3 x double]*, [3 x double]** %6
%442 = load i32, i32* %11
%443 = sext i32 %442 to i64
%444 = getelementptr inbounds [3 x double], [3 x double]* %441, i64 %443
%445 = load i32, i32* %10
%446 = sext i32 %445 to i64
%447 = getelementptr inbounds [3 x double], [3 x double]* %444, i64 0, i64 %446
%448 = load double, double* %447
store double %448, double* %20
%449 = load double, double* %21
%450 = load double, double* %19
%451 = load double, double* %20
%452 = load double, double* %21
%453 = load double, double* %16
%454 = call double @llvm.fmuladd.f64(double %452, double %453, double %451)
%455 = fneg double %450
%456 = call double @llvm.fmuladd.f64(double %455, double %454, double %449)
%457 = load [3 x double]*, [3 x double]** %6
%458 = load i32, i32* %12
%459 = sext i32 %458 to i64
%460 = getelementptr inbounds [3 x double], [3 x double]* %457, i64 %459
%461 = load i32, i32* %10
%462 = sext i32 %461 to i64
%463 = getelementptr inbounds [3 x double], [3 x double]* %460, i64 0, i64 %462
store double %456, double* %463
%464 = load double, double* %20
%465 = load double, double* %19
%466 = load double, double* %21
%467 = load double, double* %20
%468 = load double, double* %16
%469 = fneg double %467
%470 = call double @llvm.fmuladd.f64(double %469, double %468, double %466)
%471 = call double @llvm.fmuladd.f64(double %465, double %470, double %464)
%472 = load [3 x double]*, [3 x double]** %6
%473 = load i32, i32* %11
%474 = sext i32 %473 to i64
%475 = getelementptr inbounds [3 x double], [3 x double]* %472, i64 %474
%476 = load i32, i32* %10
%477 = sext i32 %476 to i64
%478 = getelementptr inbounds [3 x double], [3 x double]* %475, i64 0, i64 %477
store double %471, double* %478
br label %479
479: ; preds = %432
%480 = load i32, i32* %10
%481 = add nsw i32 %480, 1
store i32 %481, i32* %10
br label %429
482: ; preds = %429
store i32 0, i32* %10
br label %483
483: ; preds = %533, %482
%484 = load i32, i32* %10
%485 = icmp slt i32 %484, 3
br i1 %485, label %486, label %536
486: ; preds = %483
%487 = load [3 x double]*, [3 x double]** %8
%488 = load i32, i32* %10
%489 = sext i32 %488 to i64
%490 = getelementptr inbounds [3 x double], [3 x double]* %487, i64 %489
%491 = load i32, i32* %12
%492 = sext i32 %491 to i64
%493 = getelementptr inbounds [3 x double], [3 x double]* %490, i64 0, i64 %492
%494 = load double, double* %493
store double %494, double* %21
%495 = load [3 x double]*, [3 x double]** %8
%496 = load i32, i32* %10
%497 = sext i32 %496 to i64
%498 = getelementptr inbounds [3 x double], [3 x double]* %495, i64 %497
%499 = load i32, i32* %11
%500 = sext i32 %499 to i64
%501 = getelementptr inbounds [3 x double], [3 x double]* %498, i64 0, i64 %500
%502 = load double, double* %501
store double %502, double* %20
%503 = load double, double* %21
%504 = load double, double* %19
%505 = load double, double* %20
%506 = load double, double* %21
%507 = load double, double* %16
%508 = call double @llvm.fmuladd.f64(double %506, double %507, double %505)
%509 = fneg double %504
%510 = call double @llvm.fmuladd.f64(double %509, double %508, double %503)
%511 = load [3 x double]*, [3 x double]** %8
%512 = load i32, i32* %10
%513 = sext i32 %512 to i64
%514 = getelementptr inbounds [3 x double], [3 x double]* %511, i64 %513
%515 = load i32, i32* %12
%516 = sext i32 %515 to i64
%517 = getelementptr inbounds [3 x double], [3 x double]* %514, i64 0, i64 %516
store double %510, double* %517
%518 = load double, double* %20
%519 = load double, double* %19
%520 = load double, double* %21
%521 = load double, double* %20
%522 = load double, double* %16
%523 = fneg double %521
%524 = call double @llvm.fmuladd.f64(double %523, double %522, double %520)
%525 = call double @llvm.fmuladd.f64(double %519, double %524, double %518)
%526 = load [3 x double]*, [3 x double]** %8
%527 = load i32, i32* %10
%528 = sext i32 %527 to i64
%529 = getelementptr inbounds [3 x double], [3 x double]* %526, i64 %528
%530 = load i32, i32* %11
%531 = sext i32 %530 to i64
%532 = getelementptr inbounds [3 x double], [3 x double]* %529, i64 0, i64 %531
store double %525, double* %532
br label %533
533: ; preds = %486
%534 = load i32, i32* %10
%535 = add nsw i32 %534, 1
store i32 %535, i32* %10
br label %483
536: ; preds = %483
%537 = load i32*, i32** %9
%538 = load i32, i32* %537
%539 = add nsw i32 %538, 1
store i32 %539, i32* %537
br label %540
540: ; preds = %536, %187
br label %541
541: ; preds = %540, %179
br label %542
542: ; preds = %541
%543 = load i32, i32* %11
%544 = add nsw i32 %543, 1
store i32 %544, i32* %11
br label %131
545: ; preds = %131
br label %546
546: ; preds = %545
%547 = load i32, i32* %12
%548 = add nsw i32 %547, 1
store i32 %548, i32* %12
br label %125
549: ; preds = %125
store i32 0, i32* %12
br label %550
550: ; preds = %574, %549
%551 = load i32, i32* %12
%552 = icmp slt i32 %551, 3
br i1 %552, label %553, label %577
553: ; preds = %550
%554 = load i32, i32* %12
%555 = sext i32 %554 to i64
%556 = getelementptr inbounds [3 x double], [3 x double]* @jacobi.z, i64 0, i64 %555
%557 = load double, double* %556
%558 = load i32, i32* %12
%559 = sext i32 %558 to i64
%560 = getelementptr inbounds [3 x double], [3 x double]* @jacobi.b, i64 0, i64 %559
%561 = load double, double* %560
%562 = fadd double %561, %557
store double %562, double* %560
%563 = load i32, i32* %12
%564 = sext i32 %563 to i64
%565 = getelementptr inbounds [3 x double], [3 x double]* @jacobi.b, i64 0, i64 %564
%566 = load double, double* %565
%567 = load double*, double** %7
%568 = load i32, i32* %12
%569 = sext i32 %568 to i64
%570 = getelementptr inbounds double, double* %567, i64 %569
store double %566, double* %570
%571 = load i32, i32* %12
%572 = sext i32 %571 to i64
%573 = getelementptr inbounds [3 x double], [3 x double]* @jacobi.z, i64 0, i64 %572
store double 0.000000e+00, double* %573
br label %574
574: ; preds = %553
%575 = load i32, i32* %12
%576 = add nsw i32 %575, 1
store i32 %576, i32* %12
br label %550
577: ; preds = %550
br label %578
578: ; preds = %577
%579 = load i32, i32* %13
%580 = add nsw i32 %579, 1
store i32 %580, i32* %13
br label %80
581: ; preds = %80
%582 = load i32, i32* @stderr
%583 = call i32 @fprintf(i32 %582, i8* getelementptr inbounds ([31 x i8], [31 x i8]* @.str, i64 0, i64 0))
store i32 -1, i32* %5
br label %584
584: ; preds = %581, %115
%585 = load i32, i32* %5
ret i32 %585
}

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