#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 *calloc(int, int);
int /*<<< orphan*/ exit(int);
int /*<<< orphan*/ free(double *);
int /*<<< orphan*/ printf(char *);
double sqrt(double);
double conv_damping(int nBasis, double alpha, double *PA, double *PB) {
  static double *prevPA = NULL;
  static double *prevPB = NULL;
  static int nIter = 0;
  int i, j;
  double avgdP;
  if (nBasis == 0 && alpha == 0.0 && PA == NULL && PB == NULL) {
    nIter = 0;
    if (prevPA != NULL) {
      free(prevPA);
      prevPA = NULL;
    }
    if (prevPB != NULL) {
      free(prevPB);
      prevPB = NULL;
    }
    return 0.0;
  }
  if (alpha < 0.0 || alpha > 1.0) {
    printf("conv_damping : invalid range of alpha\n");
    exit(-1);
  }
  if (nIter == 0) {
    if (prevPA != NULL || prevPB != NULL) {
      printf("conv_damping : detect prior initialized pointers\n");
      exit(-1);
    }
    prevPA = calloc(nBasis * nBasis, sizeof(double));
    prevPB = calloc(nBasis * nBasis, sizeof(double));
    if (prevPA == NULL || prevPB == NULL) {
      printf("conv_damping : cannot allocate memory\n");
      exit(-1);
    }
    for (i = 0; i < nBasis; i++)
      for (j = 0; j < nBasis; j++) {
        prevPA[i * nBasis + j] = PA[i * nBasis + j];
        prevPB[i * nBasis + j] = PB[i * nBasis + j];
      }
    nIter++;
    return 0.0;
  }
  avgdP = 0.0;
  for (i = 0; i < nBasis; i++)
    for (j = 0; j < nBasis; j++) {
      avgdP += (PA[i * nBasis + j] - prevPA[i * nBasis + j]) *
               (PA[i * nBasis + j] - prevPA[i * nBasis + j]);
      avgdP += (PB[i * nBasis + j] - prevPB[i * nBasis + j]) *
               (PB[i * nBasis + j] - prevPB[i * nBasis + j]);
    }
  avgdP = sqrt(avgdP / nBasis / nBasis / 2.0);
  for (i = 0; i < nBasis; i++)
    for (j = 0; j < nBasis; j++) {
      PA[i * nBasis + j] =
          prevPA[i * nBasis + j] +
          alpha * (PA[i * nBasis + j] - prevPA[i * nBasis + j]);
      PB[i * nBasis + j] =
          prevPB[i * nBasis + j] +
          alpha * (PB[i * nBasis + j] - prevPB[i * nBasis + j]);
    }
  for (i = 0; i < nBasis; i++)
    for (j = 0; j < nBasis; j++) {
      prevPA[i * nBasis + j] = PA[i * nBasis + j];
      prevPB[i * nBasis + j] = PB[i * nBasis + j];
    }
  nIter++;
  return avgdP;
}


x86 assembly (clang, -O0):
<conv_damping>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
mov    %edi,-0xc(%rbp)
movsd  %xmm0,-0x18(%rbp)
mov    %rsi,-0x20(%rbp)
mov    %rdx,-0x28(%rbp)
cmpl   $0x0,-0xc(%rbp)
jne    1201 <conv_damping+0xb1>
movsd  -0x18(%rbp),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1201 <conv_damping+0xb1>
jp     1201 <conv_damping+0xb1>
cmpq   $0x0,-0x20(%rbp)
jne    1201 <conv_damping+0xb1>
cmpq   $0x0,-0x28(%rbp)
jne    1201 <conv_damping+0xb1>
movl   $0x0,0x2eb6(%rip)
cmpq   $0x0,0x2e9e(%rip)
je     11cf <conv_damping+0x7f>
mov    0x2e91(%rip),%rdi
call   1030 <free@plt>
movq   $0x0,0x2e81(%rip)
cmpq   $0x0,0x2e81(%rip)
je     11f4 <conv_damping+0xa4>
mov    0x2e74(%rip),%rdi
call   1030 <free@plt>
movq   $0x0,0x2e64(%rip)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    16d5 <conv_damping+0x585>
xorps  %xmm0,%xmm0
ucomisd -0x18(%rbp),%xmm0
ja     1226 <conv_damping+0xd6>
movsd  -0x18(%rbp),%xmm0
movsd  0xde4(%rip),%xmm1
ucomisd %xmm1,%xmm0
jbe    123c <conv_damping+0xec>
lea    0xde3(%rip),%rdi
call   1040 <printf@plt>
mov    $0xffffffff,%edi
call   1070 <exit@plt>
cmpl   $0x0,0x2e1d(%rip)
jne    139e <conv_damping+0x24e>
cmpq   $0x0,0x2dff(%rip)
jne    1265 <conv_damping+0x115>
cmpq   $0x0,0x2df9(%rip)
je     127b <conv_damping+0x12b>
lea    0xdcb(%rip),%rdi
call   1040 <printf@plt>
mov    $0xffffffff,%edi
call   1070 <exit@plt>
mov    -0xc(%rbp),%edi
imul   -0xc(%rbp),%edi
mov    $0x8,%esi
call   1050 <calloc@plt>
mov    %rax,0x2dbd(%rip)
mov    -0xc(%rbp),%edi
imul   -0xc(%rbp),%edi
mov    $0x8,%esi
call   1050 <calloc@plt>
mov    %rax,0x2dad(%rip)
cmpq   $0x0,0x2d9d(%rip)
je     12c7 <conv_damping+0x177>
cmpq   $0x0,0x2d97(%rip)
jne    12dd <conv_damping+0x18d>
lea    0xd9b(%rip),%rdi
call   1040 <printf@plt>
mov    $0xffffffff,%edi
call   1070 <exit@plt>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1382 <conv_damping+0x232>
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    136f <conv_damping+0x21f>
mov    -0x20(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x2d30(%rip),%rax
mov    -0x2c(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x2d09(%rip),%rax
mov    -0x2c(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    12f7 <conv_damping+0x1a7>
jmp    1374 <conv_damping+0x224>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    12e4 <conv_damping+0x194>
mov    0x2cd8(%rip),%eax
add    $0x1,%eax
mov    %eax,0x2ccf(%rip)
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x8(%rbp)
jmp    16d5 <conv_damping+0x585>
xorps  %xmm0,%xmm0
movsd  %xmm0,-0x38(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    14cd <conv_damping+0x37d>
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    14ba <conv_damping+0x36a>
mov    -0x20(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x2c67(%rip),%rax
mov    -0x2c(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    0x2c38(%rip),%rax
mov    -0x2c(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm2
movsd  -0x38(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x28(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x2bff(%rip),%rax
mov    -0x2c(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    0x2bd0(%rip),%rax
mov    -0x2c(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm2
movsd  -0x38(%rbp),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,-0x38(%rbp)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    13c0 <conv_damping+0x270>
jmp    14bf <conv_damping+0x36f>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    13ad <conv_damping+0x25d>
movsd  -0x38(%rbp),%xmm0
cvtsi2sdl -0xc(%rbp),%xmm1
divsd  %xmm1,%xmm0
cvtsi2sdl -0xc(%rbp),%xmm1
divsd  %xmm1,%xmm0
movsd  0xb1c(%rip),%xmm1
divsd  %xmm1,%xmm0
call   1060 <sqrt@plt>
movsd  %xmm0,-0x38(%rbp)
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1617 <conv_damping+0x4c7>
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    1604 <conv_damping+0x4b4>
mov    0x2b29(%rip),%rax
mov    -0x2c(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x18(%rbp),%xmm0
mov    -0x20(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    0x2af5(%rip),%rax
mov    -0x2c(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x20(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    0x2ac6(%rip),%rax
mov    -0x2c(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
movsd  -0x18(%rbp),%xmm0
mov    -0x28(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    0x2a92(%rip),%rax
mov    -0x2c(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
subsd  (%rax,%rcx,8),%xmm2
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x28(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1514 <conv_damping+0x3c4>
jmp    1609 <conv_damping+0x4b9>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    1501 <conv_damping+0x3b1>
movl   $0x0,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    16bc <conv_damping+0x56c>
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    16a9 <conv_damping+0x559>
mov    -0x20(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x29f6(%rip),%rax
mov    -0x2c(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x28(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x29cf(%rip),%rax
mov    -0x2c(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x30(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1631 <conv_damping+0x4e1>
jmp    16ae <conv_damping+0x55e>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    161e <conv_damping+0x4ce>
mov    0x299e(%rip),%eax
add    $0x1,%eax
mov    %eax,0x2995(%rip)
movsd  -0x38(%rbp),%xmm0
movsd  %xmm0,-0x8(%rbp)
movsd  -0x8(%rbp),%xmm0
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<conv_damping>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0xb8,%rsp
mov    %rdx,%r14
movapd %xmm0,%xmm5
mov    %edi,%r12d
test   %edi,%edi
jne    11e8 <conv_damping+0x88>
xorpd  %xmm0,%xmm0
ucomisd %xmm0,%xmm5
jne    11e8 <conv_damping+0x88>
jp     11e8 <conv_damping+0x88>
test   %rsi,%rsi
jne    11e8 <conv_damping+0x88>
test   %r14,%r14
jne    11e8 <conv_damping+0x88>
movl   $0x0,0x2ec9(%rip)
mov    0x2eb2(%rip),%rdi
test   %rdi,%rdi
je     11bb <conv_damping+0x5b>
call   1030 <free@plt>
movq   $0x0,0x2e9d(%rip)
mov    0x2e9e(%rip),%rdi
xorpd  %xmm3,%xmm3
test   %rdi,%rdi
je     1b3f <conv_damping+0x9df>
call   1030 <free@plt>
xorpd  %xmm3,%xmm3
movq   $0x0,0x2e7d(%rip)
jmp    1b3f <conv_damping+0x9df>
xorpd  %xmm0,%xmm0
ucomisd %xmm5,%xmm0
ja     1b55 <conv_damping+0x9f5>
ucomisd 0xe02(%rip),%xmm5
ja     1b55 <conv_damping+0x9f5>
cmpl   $0x0,0x2e5d(%rip)
mov    %r14,0x20(%rsp)
mov    %rsi,0x8(%rsp)
mov    %r12,0x18(%rsp)
je     136d <conv_damping+0x20d>
xorpd  %xmm1,%xmm1
mov    %r12d,%r13d
xorpd  %xmm0,%xmm0
test   %r12d,%r12d
jle    1336 <conv_damping+0x1d6>
mov    0x2e1d(%rip),%r8
mov    0x2e1e(%rip),%r9
mov    %r13d,%r15d
and    $0xfffffffe,%r15d
lea    0x8(%r9),%rdx
lea    0x0(,%r13,8),%r10
lea    0x8(%r14),%rbp
lea    0x8(%r8),%rax
lea    0x8(%rsi),%rcx
xorpd  %xmm0,%xmm0
xor    %r11d,%r11d
jmp    1289 <conv_damping+0x129>
nopw   0x0(%rax,%rax,1)
add    $0x1,%r11
add    %r10,%rdx
add    %r10,%rbp
add    %r10,%rax
add    %r10,%rcx
cmp    %r13,%r11
je     1336 <conv_damping+0x1d6>
xor    %edi,%edi
cmp    $0x1,%r12d
je     12f6 <conv_damping+0x196>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movsd  -0x8(%rcx,%rdi,8),%xmm2
movsd  (%rcx,%rdi,8),%xmm3
subsd  -0x8(%rax,%rdi,8),%xmm2
mulsd  %xmm2,%xmm2
addsd  %xmm0,%xmm2
movsd  -0x8(%rbp,%rdi,8),%xmm4
movsd  0x0(%rbp,%rdi,8),%xmm0
subsd  -0x8(%rdx,%rdi,8),%xmm4
mulsd  %xmm4,%xmm4
addsd  %xmm2,%xmm4
subsd  (%rax,%rdi,8),%xmm3
mulsd  %xmm3,%xmm3
addsd  %xmm4,%xmm3
subsd  (%rdx,%rdi,8),%xmm0
mulsd  %xmm0,%xmm0
addsd  %xmm3,%xmm0
add    $0x2,%rdi
cmp    %rdi,%r15
jne    12a0 <conv_damping+0x140>
test   $0x1,%r13b
je     1270 <conv_damping+0x110>
mov    %r11,%rbx
imul   %r13,%rbx
add    %rbx,%rdi
movsd  (%rsi,%rdi,8),%xmm2
subsd  (%r8,%rdi,8),%xmm2
mulsd  %xmm2,%xmm2
addsd  %xmm0,%xmm2
movsd  (%r14,%rdi,8),%xmm0
subsd  (%r9,%rdi,8),%xmm0
mulsd  %xmm0,%xmm0
addsd  %xmm2,%xmm0
jmp    1270 <conv_damping+0x110>
xorps  %xmm2,%xmm2
cvtsi2sd %r12d,%xmm2
divsd  %xmm2,%xmm0
divsd  %xmm2,%xmm0
mulsd  0xcba(%rip),%xmm0
ucomisd %xmm1,%xmm0
jb     1527 <conv_damping+0x3c7>
xorps  %xmm3,%xmm3
sqrtsd %xmm0,%xmm3
test   %r12d,%r12d
jg     1545 <conv_damping+0x3e5>
jmp    1b38 <conv_damping+0x9d8>
mov    0x2ce4(%rip),%rax
or     0x2ce5(%rip),%rax
jne    1b5e <conv_damping+0x9fe>
mov    %r12d,%ebx
imul   %r12d,%ebx
mov    %ebx,%edi
mov    $0x8,%esi
call   1050 <calloc@plt>
mov    %rax,%r15
mov    %rax,0x2cba(%rip)
mov    %ebx,%edi
mov    $0x8,%esi
call   1050 <calloc@plt>
mov    %rax,0x2caf(%rip)
mov    %r15,0x40(%rsp)
test   %r15,%r15
je     1b67 <conv_damping+0xa07>
mov    %rax,%rbp
test   %rax,%rax
je     1b67 <conv_damping+0xa07>
test   %r12d,%r12d
jle    1514 <conv_damping+0x3b4>
mov    %r12,%rax
mov    %eax,%ecx
lea    0x0(,%rcx,8),%rbx
xor    %r12d,%r12d
cmp    $0x1,%eax
je     14d2 <conv_damping+0x372>
mov    0x18(%rsp),%rax
add    %eax,%eax
mov    %eax,0x38(%rsp)
mov    %rcx,0x58(%rsp)
and    $0xfffffffe,%ecx
neg    %rcx
mov    %rcx,0x50(%rsp)
xor    %r12d,%r12d
mov    %rbx,0x28(%rsp)
xor    %r13d,%r13d
mov    0x28(%rsp),%r15
nopl   0x0(%rax,%rax,1)
movslq %r13d,%r13
mov    0x8(%rsp),%r14
lea    (%r14,%r13,8),%rsi
mov    0x40(%rsp),%rax
lea    (%rax,%r13,8),%rdi
mov    0x20(%rsp),%rax
lea    (%rax,%r13,8),%rax
mov    %rax,0x10(%rsp)
lea    0x0(,%r13,8),%rbx
add    %rbp,%rbx
mov    %r15,%rdx
call   1060 <memcpy@plt>
mov    %rbx,%rdi
mov    0x10(%rsp),%rsi
mov    0x20(%rsp),%rbx
mov    %r15,%rdx
call   1060 <memcpy@plt>
mov    0x18(%rsp),%rax
add    %r13d,%eax
cltq
lea    (%r14,%rax,8),%rsi
mov    %rbx,%r14
mov    0x40(%rsp),%rcx
lea    (%rcx,%rax,8),%rdi
lea    (%rbx,%rax,8),%rcx
mov    %rcx,0x10(%rsp)
lea    0x0(,%rax,8),%rbx
add    %rbp,%rbx
mov    %r15,%rdx
call   1060 <memcpy@plt>
mov    %rbx,%rdi
mov    0x10(%rsp),%rsi
mov    %r15,%rdx
call   1060 <memcpy@plt>
add    0x38(%rsp),%r13d
add    $0xfffffffffffffffe,%r12
cmp    %r12,0x50(%rsp)
jne    1420 <conv_damping+0x2c0>
neg    %r12d
mov    0x28(%rsp),%rbx
mov    0x58(%rsp),%rcx
test   $0x1,%cl
mov    0x18(%rsp),%rax
je     1514 <conv_damping+0x3b4>
imul   %eax,%r12d
movslq %r12d,%rax
mov    0x8(%rsp),%rcx
lea    (%rcx,%rax,8),%rsi
mov    0x40(%rsp),%rcx
lea    (%rcx,%rax,8),%rdi
lea    (%r14,%rax,8),%r14
lea    0x0(%rbp,%rax,8),%rbp
mov    %rbx,%rdx
call   1060 <memcpy@plt>
mov    %rbp,%rdi
mov    %r14,%rsi
mov    %rbx,%rdx
call   1060 <memcpy@plt>
movl   $0x1,0x2b4a(%rip)
xorpd  %xmm3,%xmm3
jmp    1b3f <conv_damping+0x9df>
movapd %xmm5,0x40(%rsp)
call   1070 <sqrt@plt>
movapd 0x40(%rsp),%xmm5
movapd %xmm0,%xmm3
test   %r12d,%r12d
jle    1b38 <conv_damping+0x9d8>
mov    0x2b0c(%rip),%rax
mov    0x2b0d(%rip),%rcx
mov    %r13d,%edx
and    $0xfffffffe,%edx
mov    %rdx,0x68(%rsp)
movapd %xmm5,%xmm0
unpcklpd %xmm5,%xmm0
lea    0x0(,%r13,8),%rdi
mov    %rdi,%r11
and    $0xfffffffffffffff0,%r11
xor    %ebx,%ebx
mov    %rax,0x58(%rsp)
mov    0x8(%rsp),%r8
mov    %r8,%rdx
mov    %rcx,0x70(%rsp)
mov    0x20(%rsp),%rsi
mov    0x18(%rsp),%r9
mov    %rdi,0x30(%rsp)
jmp    15b9 <conv_damping+0x459>
nopw   0x0(%rax,%rax,1)
add    $0x1,%rbx
add    %rdi,%rsi
add    %rdi,%rcx
add    %rdi,%rdx
add    %rdi,%rax
cmp    %r13,%rbx
je     1760 <conv_damping+0x600>
cmp    $0x2,%r9d
jae    1620 <conv_damping+0x4c0>
xor    %ebp,%ebp
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
movsd  (%rax,%rbp,8),%xmm1
movsd  (%rdx,%rbp,8),%xmm2
subsd  %xmm1,%xmm2
mulsd  %xmm5,%xmm2
addsd  %xmm1,%xmm2
movsd  %xmm2,(%rdx,%rbp,8)
movsd  (%rcx,%rbp,8),%xmm1
movsd  (%rsi,%rbp,8),%xmm2
subsd  %xmm1,%xmm2
mulsd  %xmm5,%xmm2
addsd  %xmm1,%xmm2
movsd  %xmm2,(%rsi,%rbp,8)
add    $0x1,%rbp
cmp    %rbp,%r13
jne    15d0 <conv_damping+0x470>
jmp    15a0 <conv_damping+0x440>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    %rbx,0x40(%rsp)
imul   %r13,%rbx
lea    (%rbx,%r13,1),%rbp
lea    (%r8,%rbx,8),%rdi
lea    (%r8,%rbp,8),%r8
mov    0x20(%rsp),%r9
lea    (%r9,%rbx,8),%r15
lea    (%r9,%rbp,8),%r10
mov    0x58(%rsp),%r12
lea    (%r12,%rbx,8),%r14
lea    (%r12,%rbp,8),%r12
mov    0x70(%rsp),%r9
lea    (%r9,%rbp,8),%rbp
cmp    %r10,%rdi
setb   0x10(%rsp)
cmp    %r8,%r15
setb   0x28(%rsp)
cmp    %r12,%rdi
setb   0x38(%rsp)
cmp    %r8,%r14
setb   0x50(%rsp)
cmp    %rbp,%rdi
lea    (%r9,%rbx,8),%rbx
setb   %dil
cmp    %r8,%rbx
setb   %r8b
cmp    %r12,%r15
setb   %r12b
cmp    %r10,%r14
setb   %r9b
cmp    %rbp,%r15
setb   %r15b
cmp    %r10,%rbx
setb   %r10b
mov    0x28(%rsp),%bl
test   %bl,0x10(%rsp)
jne    1745 <conv_damping+0x5e5>
mov    0x50(%rsp),%bl
and    %bl,0x38(%rsp)
jne    1745 <conv_damping+0x5e5>
and    %r8b,%dil
jne    1745 <conv_damping+0x5e5>
and    %r9b,%r12b
jne    1745 <conv_damping+0x5e5>
mov    $0x0,%ebp
and    %r10b,%r15b
mov    0x8(%rsp),%r8
mov    0x18(%rsp),%r9
mov    0x30(%rsp),%rdi
mov    0x40(%rsp),%rbx
jne    15d0 <conv_damping+0x470>
xor    %ebp,%ebp
movupd (%rax,%rbp,1),%xmm1
movupd (%rdx,%rbp,1),%xmm2
subpd  %xmm1,%xmm2
mulpd  %xmm0,%xmm2
addpd  %xmm1,%xmm2
movupd %xmm2,(%rdx,%rbp,1)
movupd (%rcx,%rbp,1),%xmm1
movupd (%rsi,%rbp,1),%xmm2
subpd  %xmm1,%xmm2
mulpd  %xmm0,%xmm2
addpd  %xmm1,%xmm2
movupd %xmm2,(%rsi,%rbp,1)
add    $0x10,%rbp
cmp    %rbp,%r11
jne    16f0 <conv_damping+0x590>
mov    0x68(%rsp),%r10
mov    %r10,%rbp
cmp    %r13,%r10
jne    15d0 <conv_damping+0x470>
jmp    15a0 <conv_damping+0x440>
xor    %ebp,%ebp
mov    0x8(%rsp),%r8
mov    0x18(%rsp),%r9
mov    0x30(%rsp),%rdi
mov    0x40(%rsp),%rbx
jmp    15d0 <conv_damping+0x470>
test   %r9d,%r9d
jle    1b38 <conv_damping+0x9d8>
mov    0x28e8(%rip),%rdi
mov    0x28e9(%rip),%r9
mov    %r13d,%eax
and    $0xfffffffc,%eax
mov    %rax,0x98(%rsp)
add    $0xfffffffffffffffc,%rax
mov    %rax,0x90(%rsp)
mov    %rax,%rcx
shr    $0x2,%rcx
add    $0x1,%rcx
mov    %rcx,%rax
mov    %rcx,0x88(%rsp)
and    $0xfffffffffffffffe,%rcx
mov    %rcx,0x80(%rsp)
mov    %r13,%rax
neg    %rax
mov    %rax,0xa0(%rsp)
lea    0x30(%r9),%r12
lea    0x0(,%r13,8),%r11
mov    0x20(%rsp),%rsi
lea    0x30(%rsi),%rdx
lea    0x30(%rdi),%rbp
mov    0x8(%rsp),%rcx
lea    0x30(%rcx),%rax
add    $0x8,%rcx
mov    %rdi,0x28(%rsp)
add    $0x8,%rdi
lea    0x8(%rsi),%r8
mov    %r9,0x38(%rsp)
add    $0x8,%r9
xor    %r10d,%r10d
mov    %r11,0x78(%rsp)
jmp    1835 <conv_damping+0x6d5>
nopw   0x0(%rax,%rax,1)
add    $0x1,%r10
add    %r11,%r12
add    %r11,%rdx
add    %r11,%rbp
add    %r11,%rax
add    %r11,%rcx
add    %r11,%rdi
add    %r11,%r8
add    %r11,%r9
cmp    %r13,%r10
je     1b38 <conv_damping+0x9d8>
mov    %r10,%rsi
imul   %r13,%rsi
mov    %rsi,0x40(%rsp)
cmpl   $0x4,0x18(%rsp)
jae    18f0 <conv_damping+0x790>
xor    %r15d,%r15d
mov    %r15,%r14
test   $0x1,%r13b
je     188f <conv_damping+0x72f>
mov    0x40(%rsp),%rbx
add    %r15,%rbx
mov    0x8(%rsp),%rsi
movsd  (%rsi,%rbx,8),%xmm0
mov    0x28(%rsp),%rsi
movsd  %xmm0,(%rsi,%rbx,8)
mov    0x20(%rsp),%rsi
movsd  (%rsi,%rbx,8),%xmm0
mov    0x38(%rsp),%rsi
movsd  %xmm0,(%rsi,%rbx,8)
mov    %r15,%r14
or     $0x1,%r14
not    %r15
cmp    0xa0(%rsp),%r15
je     1810 <conv_damping+0x6b0>
movsd  -0x8(%rcx,%r14,8),%xmm0
movsd  %xmm0,-0x8(%rdi,%r14,8)
movsd  -0x8(%r8,%r14,8),%xmm0
movsd  %xmm0,-0x8(%r9,%r14,8)
movsd  (%rcx,%r14,8),%xmm0
movsd  %xmm0,(%rdi,%r14,8)
movsd  (%r8,%r14,8),%xmm0
movsd  %xmm0,(%r9,%r14,8)
add    $0x2,%r14
cmp    %r14,%r13
jne    18a0 <conv_damping+0x740>
jmp    1810 <conv_damping+0x6b0>
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax)
mov    %r10,0x50(%rsp)
mov    %r10,%rsi
imul   %r13,%rsi
lea    (%rsi,%r13,1),%r11
mov    0x28(%rsp),%r10
lea    (%r10,%rsi,8),%r15
lea    (%r10,%r11,8),%rbx
mov    %rbx,0x60(%rsp)
mov    0x38(%rsp),%r10
lea    (%r10,%rsi,8),%rbx
mov    %rbx,0xa8(%rsp)
lea    (%r10,%r11,8),%r14
mov    %r14,0x30(%rsp)
mov    0x8(%rsp),%r10
lea    (%r10,%rsi,8),%rbx
mov    %rbx,0x68(%rsp)
lea    (%r10,%r11,8),%rbx
mov    %rbx,0x10(%rsp)
mov    0x20(%rsp),%rbx
lea    (%rbx,%r11,8),%rbx
mov    %rbx,0xb0(%rsp)
cmp    %r14,%r15
setb   0x58(%rsp)
mov    0xa8(%rsp),%r11
mov    0x60(%rsp),%r14
cmp    %r14,%r11
setb   0x70(%rsp)
cmp    0x10(%rsp),%r15
setb   0x7(%rsp)
mov    0x68(%rsp),%r10
cmp    %r14,%r10
setb   0x6(%rsp)
cmp    %rbx,%r15
mov    0x20(%rsp),%rbx
lea    (%rbx,%rsi,8),%r15
setb   0x5(%rsp)
cmp    %r14,%r15
setb   0x60(%rsp)
cmp    0x10(%rsp),%r11
setb   %r14b
mov    0x30(%rsp),%rsi
cmp    %rsi,%r10
setb   0x10(%rsp)
cmp    0xb0(%rsp),%r11
setb   %r10b
cmp    %rsi,%r15
setb   %r15b
mov    0x70(%rsp),%bl
test   %bl,0x58(%rsp)
jne    1b26 <conv_damping+0x9c6>
mov    %r10d,%esi
mov    0x6(%rsp),%bl
and    %bl,0x7(%rsp)
mov    0x78(%rsp),%r11
mov    0x50(%rsp),%r10
jne    184c <conv_damping+0x6ec>
mov    0x60(%rsp),%bl
and    %bl,0x5(%rsp)
jne    184c <conv_damping+0x6ec>
and    0x10(%rsp),%r14b
jne    184c <conv_damping+0x6ec>
mov    %r15d,%ebx
mov    $0x0,%r15d
and    %bl,%sil
jne    184f <conv_damping+0x6ef>
cmpq   $0x0,0x90(%rsp)
je     1b17 <conv_damping+0x9b7>
mov    0x80(%rsp),%rsi
xor    %r14d,%r14d
nopl   0x0(%rax,%rax,1)
movups -0x30(%rax,%r14,8),%xmm0
movups -0x20(%rax,%r14,8),%xmm1
movups %xmm0,-0x30(%rbp,%r14,8)
movups %xmm1,-0x20(%rbp,%r14,8)
movups -0x30(%rdx,%r14,8),%xmm0
movups -0x20(%rdx,%r14,8),%xmm1
movups %xmm0,-0x30(%r12,%r14,8)
movups %xmm1,-0x20(%r12,%r14,8)
movups -0x10(%rax,%r14,8),%xmm0
movups (%rax,%r14,8),%xmm1
movups %xmm0,-0x10(%rbp,%r14,8)
movups %xmm1,0x0(%rbp,%r14,8)
movupd -0x10(%rdx,%r14,8),%xmm0
movupd (%rdx,%r14,8),%xmm1
movupd %xmm0,-0x10(%r12,%r14,8)
movupd %xmm1,(%r12,%r14,8)
add    $0x8,%r14
add    $0xfffffffffffffffe,%rsi
jne    1a40 <conv_damping+0x8e0>
testb  $0x1,0x88(%rsp)
01
je     1afe <conv_damping+0x99e>
add    0x40(%rsp),%r14
mov    0x8(%rsp),%rsi
movups (%rsi,%r14,8),%xmm0
movups 0x10(%rsi,%r14,8),%xmm1
mov    0x28(%rsp),%rsi
movups %xmm0,(%rsi,%r14,8)
movups %xmm1,0x10(%rsi,%r14,8)
mov    0x20(%rsp),%rsi
movupd (%rsi,%r14,8),%xmm0
movupd 0x10(%rsi,%r14,8),%xmm1
mov    0x38(%rsp),%rsi
movupd %xmm0,(%rsi,%r14,8)
movupd %xmm1,0x10(%rsi,%r14,8)
mov    0x98(%rsp),%rsi
mov    %rsi,%r15
cmp    %r13,%rsi
je     1810 <conv_damping+0x6b0>
jmp    184f <conv_damping+0x6ef>
xor    %r14d,%r14d
testb  $0x1,0x88(%rsp)
01
jne    1ab5 <conv_damping+0x955>
jmp    1afe <conv_damping+0x99e>
xor    %r15d,%r15d
mov    0x78(%rsp),%r11
mov    0x50(%rsp),%r10
jmp    184f <conv_damping+0x6ef>
addl   $0x1,0x2529(%rip)
movapd %xmm3,%xmm0
add    $0xb8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret
lea    0x50b(%rip),%rdi
jmp    1b6e <conv_damping+0xa0e>
lea    0x4d1(%rip),%rdi
jmp    1b6e <conv_damping+0xa0e>
lea    0x4a2(%rip),%rdi
call   1040 <puts@plt>
mov    $0xffffffff,%edi
call   1080 <exit@plt>

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

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