#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 */
int /*<<< orphan*/ memcpy(double *, double const *, int);
int /*<<< orphan*/ memset(double *, int /*<<< orphan*/, int);
int /*<<< orphan*/ myinvWarm(double *, double *, double *, int *, int *, int);
void DenmanBeaversSqrtm(const double *A, int n, int nSteps, double *solution,
                        char *workingBuffer) {
  int i, j;
  memcpy(solution, A, n * n * sizeof(double));
  double *Y = (double *)workingBuffer;
  double *Xk = Y + n * n;
  double *Yk = Xk + n * n;
  double *tmp = Yk + n * n;
  double *bPtr = tmp + n * n;
  int *permPtr = (int *)(bPtr + n * n);
  int *perPtr2 = permPtr + n * n;
  memset(Y, 0, n * n * sizeof(double));
  for (i = 0; i < n; i++)
    Y[i * n + i] = 1.0;
  for (j = 0; j < nSteps; j++) {
    memcpy(Xk, solution, n * n * sizeof(double));
    memcpy(Yk, Y, n * n * sizeof(double));
    myinvWarm(Yk, tmp, bPtr, permPtr, perPtr2, n);
    for (i = 0; i < n * n; i++)
      solution[i] = (Xk[i] + tmp[i]) * 0.5;
    myinvWarm(Xk, tmp, bPtr, permPtr, perPtr2, n);
    for (i = 0; i < n * n; i++)
      Y[i] = (Yk[i] + tmp[i]) * 0.5;
  }
}


x86 assembly (clang, -O0):
<DenmanBeaversSqrtm>:
push   %rbp
mov    %rsp,%rbp
sub    $0x60,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
mov    %rcx,-0x18(%rbp)
mov    %r8,-0x20(%rbp)
mov    -0x18(%rbp),%rdi
mov    -0x8(%rbp),%rsi
mov    -0xc(%rbp),%eax
imul   -0xc(%rbp),%eax
cltq
shl    $0x3,%rax
mov    %eax,%edx
call   1050 <memcpy@plt>
mov    -0x20(%rbp),%rax
mov    %rax,-0x30(%rbp)
mov    -0x30(%rbp),%rax
mov    -0xc(%rbp),%ecx
imul   -0xc(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x38(%rbp)
mov    -0x38(%rbp),%rax
mov    -0xc(%rbp),%ecx
imul   -0xc(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x40(%rbp)
mov    -0x40(%rbp),%rax
mov    -0xc(%rbp),%ecx
imul   -0xc(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x48(%rbp)
mov    -0x48(%rbp),%rax
mov    -0xc(%rbp),%ecx
imul   -0xc(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x50(%rbp)
mov    -0x50(%rbp),%rax
mov    -0xc(%rbp),%ecx
imul   -0xc(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x3,%rcx
add    %rcx,%rax
mov    %rax,-0x58(%rbp)
mov    -0x58(%rbp),%rax
mov    -0xc(%rbp),%ecx
imul   -0xc(%rbp),%ecx
movslq %ecx,%rcx
shl    $0x2,%rcx
add    %rcx,%rax
mov    %rax,-0x60(%rbp)
mov    -0x30(%rbp),%rdi
mov    -0xc(%rbp),%eax
imul   -0xc(%rbp),%eax
cltq
shl    $0x3,%rax
mov    %eax,%edx
xor    %esi,%esi
call   1040 <memset@plt>
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
cmp    -0xc(%rbp),%eax
jge    125d <DenmanBeaversSqrtm+0x12d>
mov    -0x30(%rbp),%rax
mov    -0x24(%rbp),%ecx
imul   -0xc(%rbp),%ecx
add    -0x24(%rbp),%ecx
movslq %ecx,%rcx
movsd  0xdbe(%rip),%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1225 <DenmanBeaversSqrtm+0xf5>
movl   $0x0,-0x28(%rbp)
mov    -0x28(%rbp),%eax
cmp    -0x10(%rbp),%eax
jge    13a9 <DenmanBeaversSqrtm+0x279>
mov    -0x38(%rbp),%rdi
mov    -0x18(%rbp),%rsi
mov    -0xc(%rbp),%eax
imul   -0xc(%rbp),%eax
cltq
shl    $0x3,%rax
mov    %eax,%edx
call   1050 <memcpy@plt>
mov    -0x40(%rbp),%rdi
mov    -0x30(%rbp),%rsi
mov    -0xc(%rbp),%eax
imul   -0xc(%rbp),%eax
cltq
shl    $0x3,%rax
mov    %eax,%edx
call   1050 <memcpy@plt>
mov    -0x40(%rbp),%rdi
mov    -0x48(%rbp),%rsi
mov    -0x50(%rbp),%rdx
mov    -0x58(%rbp),%rcx
mov    -0x60(%rbp),%r8
mov    -0xc(%rbp),%r9d
call   1030 <myinvWarm@plt>
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
mov    -0xc(%rbp),%ecx
imul   -0xc(%rbp),%ecx
cmp    %ecx,%eax
jge    131f <DenmanBeaversSqrtm+0x1ef>
mov    -0x38(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x48(%rbp),%rax
movslq -0x24(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  0xd00(%rip),%xmm1
mulsd  %xmm1,%xmm0
mov    -0x18(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    12cc <DenmanBeaversSqrtm+0x19c>
mov    -0x38(%rbp),%rdi
mov    -0x48(%rbp),%rsi
mov    -0x50(%rbp),%rdx
mov    -0x58(%rbp),%rcx
mov    -0x60(%rbp),%r8
mov    -0xc(%rbp),%r9d
call   1030 <myinvWarm@plt>
movl   $0x0,-0x24(%rbp)
mov    -0x24(%rbp),%eax
mov    -0xc(%rbp),%ecx
imul   -0xc(%rbp),%ecx
cmp    %ecx,%eax
jge    1396 <DenmanBeaversSqrtm+0x266>
mov    -0x40(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x48(%rbp),%rax
movslq -0x24(%rbp),%rcx
addsd  (%rax,%rcx,8),%xmm0
movsd  0xc89(%rip),%xmm1
mulsd  %xmm1,%xmm0
mov    -0x30(%rbp),%rax
movslq -0x24(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x24(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x24(%rbp)
jmp    1343 <DenmanBeaversSqrtm+0x213>
jmp    139b <DenmanBeaversSqrtm+0x26b>
mov    -0x28(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x28(%rbp)
jmp    1264 <DenmanBeaversSqrtm+0x134>
add    $0x60,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<DenmanBeaversSqrtm>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0xa8,%rsp
mov    %r8,%r14
mov    %rcx,%r13
mov    %edx,0x14(%rsp)
mov    %esi,%ebp
mov    %rdi,%rsi
mov    %ebp,%r15d
imul   %r15d,%r15d
lea    0x0(,%r15,8),%ebx
mov    %rcx,%rdi
mov    %ebx,%edx
call   1050 <memcpy@plt>
mov    %r14,%rdi
xor    %esi,%esi
mov    %ebx,0x4c(%rsp)
mov    %ebx,%edx
call   1040 <memset@plt>
test   %ebp,%ebp
mov    %rbp,0x40(%rsp)
mov    %r15,0x78(%rsp)
jle    1276 <DenmanBeaversSqrtm+0x146>
mov    %ebp,%r11d
lea    -0x1(%r11),%rax
mov    %r11,%rcx
and    $0x3,%r11d
xor    %edx,%edx
movabs $0x3ff0000000000000,%rbx
cmp    $0x3,%rax
jb     1255 <DenmanBeaversSqrtm+0x125>
lea    0x0(%rbp,%rbp,1),%edi
lea    0x0(,%rbp,4),%eax
mov    %eax,0x38(%rsp)
lea    0x0(%rbp,%rbp,2),%esi
and    $0xfffffffc,%ecx
neg    %rcx
mov    %rcx,0x8(%rsp)
mov    $0x2,%r15d
mov    $0x1,%r12d
mov    %r14,%r10
mov    $0x3,%r14d
xor    %r9d,%r9d
mov    %ebp,%eax
xor    %ecx,%ecx
nopl   0x0(%rax)
movslq %edi,%rdi
cltq
movslq %edx,%rdx
lea    (%rcx,%rdx,1),%r8
mov    %rbx,(%r10,%r8,8)
lea    (%r12,%rax,1),%rbp
mov    %rbx,(%r10,%rbp,8)
lea    (%r15,%rdi,1),%rbp
movslq %esi,%rsi
mov    %rbx,(%r10,%rbp,8)
lea    (%r14,%rsi,1),%rbp
mov    %rbx,(%r10,%rbp,8)
add    $0x4,%rcx
mov    0x38(%rsp),%ebp
add    %ebp,%edx
add    $0x4,%r15
add    %ebp,%edi
add    $0x4,%r12
add    %ebp,%eax
add    $0x4,%r14
add    %ebp,%esi
add    $0xfffffffffffffffc,%r9
cmp    %r9,0x8(%rsp)
jne    11f0 <DenmanBeaversSqrtm+0xc0>
neg    %r9
mov    %r9,%rdx
mov    0x40(%rsp),%rbp
mov    0x78(%rsp),%r15
mov    %r10,%r14
test   %r11,%r11
je     1276 <DenmanBeaversSqrtm+0x146>
mov    %ebp,%eax
imul   %edx,%eax
nop
cltq
lea    (%rdx,%rax,1),%rcx
mov    %rbx,(%r14,%rcx,8)
add    $0x1,%rdx
add    %ebp,%eax
add    $0xffffffffffffffff,%r11
jne    1260 <DenmanBeaversSqrtm+0x130>
cmpl   $0x0,0x14(%rsp)
jle    179c <DenmanBeaversSqrtm+0x66c>
lea    (%r14,%r15,8),%rbp
lea    0x0(,%r15,8),%rax
add    %rbp,%rax
mov    %rax,0x80(%rsp)
lea    (%rax,%r15,8),%rax
mov    %rax,0x8(%rsp)
lea    (%rax,%r15,8),%rax
mov    %rax,0x30(%rsp)
lea    (%rax,%r15,8),%r11
lea    0x0(%r13,%r15,8),%r10
lea    0x0(,%r15,8),%rcx
lea    (%rcx,%rcx,2),%rdi
lea    (%r14,%rdi,1),%r8
mov    %r15,%rdx
shl    $0x5,%rdx
add    %r14,%rdx
mov    %r15,%rbx
shl    $0x4,%rbx
lea    (%r14,%rbx,1),%r9
mov    %r15d,%ecx
and    $0x7ffffffc,%ecx
mov    %rcx,0x68(%rsp)
add    $0xfffffffffffffffc,%rcx
mov    %rcx,0x60(%rsp)
shr    $0x2,%rcx
add    $0x1,%rcx
cmp    %r13,%rdx
seta   %dl
cmp    %r10,%r8
setb   %al
and    %dl,%al
cmp    %r13,%r9
seta   %r8b
lea    (%r11,%r15,4),%rsi
mov    %rsi,0x28(%rsp)
cmp    %r10,%rbp
setb   %dl
and    %r8b,%dl
or     %al,%dl
mov    %dl,0x7(%rsp)
mov    %rcx,%rax
mov    %rcx,0x58(%rsp)
and    $0xfffffffffffffffe,%rcx
mov    %rcx,0x50(%rsp)
mov    %r15,%rax
neg    %rax
mov    %rax,0x90(%rsp)
lea    (%r14,%rdi,1),%r12
add    $0x30,%r12
lea    (%r14,%rbx,1),%rax
add    $0x30,%rax
mov    %rax,0x88(%rsp)
xor    %esi,%esi
mov    %r14,0x20(%rsp)
mov    %rdi,0xa0(%rsp)
mov    %rbx,0x98(%rsp)
mov    %r11,0x70(%rsp)
mov    %r13,0x18(%rsp)
jmp    13bd <DenmanBeaversSqrtm+0x28d>
xchg   %ax,%ax
mov    %rbp,%rdi
mov    0x8(%rsp),%rsi
mov    %r15,%rdx
mov    %r13,%r8
mov    0x70(%rsp),%r13
mov    %r13,%rcx
mov    %r14d,%r9d
call   1030 <myinvWarm@plt>
mov    0x20(%rsp),%r14
mov    %r13,%r11
mov    0x38(%rsp),%rsi
mov    0x18(%rsp),%r13
add    $0x1,%esi
cmp    0x14(%rsp),%esi
je     179c <DenmanBeaversSqrtm+0x66c>
mov    %rsi,0x38(%rsp)
mov    %rbp,%rdi
mov    %r13,%rsi
mov    %r11,%r13
mov    0x4c(%rsp),%r15d
mov    %r15d,%edx
call   1050 <memcpy@plt>
mov    0x80(%rsp),%rbx
mov    %rbx,%rdi
mov    %r14,%rsi
mov    %r15d,%edx
call   1050 <memcpy@plt>
mov    %rbx,%rdi
mov    0x8(%rsp),%rsi
mov    0x30(%rsp),%r15
mov    %r15,%rdx
mov    %r13,%rcx
mov    0x28(%rsp),%r13
mov    %r13,%r8
mov    0x78(%rsp),%rbx
mov    0x40(%rsp),%r14
mov    %r14d,%r9d
call   1030 <myinvWarm@plt>
test   %rbx,%rbx
je     1380 <DenmanBeaversSqrtm+0x250>
mov    0x70(%rsp),%rsi
cmp    $0x4,%ebx
setb   %al
or     0x7(%rsp),%al
mov    %rbx,%r15
je     1460 <DenmanBeaversSqrtm+0x330>
xor    %edi,%edi
mov    0x20(%rsp),%r14
mov    0x30(%rsp),%rdx
mov    %rsi,%rcx
mov    0x28(%rsp),%r8
movsd  0xbbc(%rip),%xmm3
mov    0x18(%rsp),%r13
jmp    1567 <DenmanBeaversSqrtm+0x437>
xchg   %ax,%ax
cmpq   $0x0,0x60(%rsp)
je     1751 <DenmanBeaversSqrtm+0x621>
mov    0x50(%rsp),%rdi
xor    %eax,%eax
mov    0x20(%rsp),%r14
mov    0x30(%rsp),%rdx
mov    %rsi,%rcx
mov    0x28(%rsp),%r8
movsd  0xb83(%rip),%xmm3
movapd 0xb6b(%rip),%xmm4
mov    0x18(%rsp),%r13
nopw   0x0(%rax,%rax,1)
movupd 0x0(%rbp,%rax,8),%xmm0
movupd 0x10(%rbp,%rax,8),%xmm1
movupd -0x30(%r12,%rax,8),%xmm2
addpd  %xmm0,%xmm2
movupd -0x20(%r12,%rax,8),%xmm0
addpd  %xmm1,%xmm0
mulpd  %xmm4,%xmm2
mulpd  %xmm4,%xmm0
movupd %xmm2,0x0(%r13,%rax,8)
movupd %xmm0,0x10(%r13,%rax,8)
movupd 0x20(%rbp,%rax,8),%xmm0
movupd 0x30(%rbp,%rax,8),%xmm1
movupd -0x10(%r12,%rax,8),%xmm2
addpd  %xmm0,%xmm2
movupd (%r12,%rax,8),%xmm0
addpd  %xmm1,%xmm0
mulpd  %xmm4,%xmm2
mulpd  %xmm4,%xmm0
movupd %xmm2,0x20(%r13,%rax,8)
movupd %xmm0,0x30(%r13,%rax,8)
add    $0x8,%rax
add    $0xfffffffffffffffe,%rdi
jne    14a0 <DenmanBeaversSqrtm+0x370>
testb  $0x1,0x58(%rsp)
je     155a <DenmanBeaversSqrtm+0x42a>
movupd 0x0(%rbp,%rax,8),%xmm0
movupd 0x10(%rbp,%rax,8),%xmm1
mov    0x8(%rsp),%rsi
movupd (%rsi,%rax,8),%xmm2
addpd  %xmm0,%xmm2
movupd 0x10(%rsi,%rax,8),%xmm0
addpd  %xmm1,%xmm0
mulpd  %xmm4,%xmm2
mulpd  %xmm4,%xmm0
movupd %xmm2,0x0(%r13,%rax,8)
movupd %xmm0,0x10(%r13,%rax,8)
mov    0x68(%rsp),%rax
mov    %rax,%rdi
cmp    %r15,%rax
je     15d9 <DenmanBeaversSqrtm+0x4a9>
mov    %rdi,%rax
test   $0x1,%r15b
je     1592 <DenmanBeaversSqrtm+0x462>
movsd  0x0(%rbp,%rdi,8),%xmm0
mov    0x8(%rsp),%rax
addsd  (%rax,%rdi,8),%xmm0
mulsd  %xmm3,%xmm0
movsd  %xmm0,0x0(%r13,%rdi,8)
mov    %rdi,%rax
or     $0x1,%rax
not    %rdi
cmp    0x90(%rsp),%rdi
je     15d9 <DenmanBeaversSqrtm+0x4a9>
nop
movsd  0x0(%rbp,%rax,8),%xmm0
addsd  -0x30(%r12,%rax,8),%xmm0
mulsd  %xmm3,%xmm0
movsd  %xmm0,0x0(%r13,%rax,8)
movsd  0x8(%rbp,%rax,8),%xmm0
addsd  -0x28(%r12,%rax,8),%xmm0
mulsd  %xmm3,%xmm0
movsd  %xmm0,0x8(%r13,%rax,8)
add    $0x2,%rax
cmp    %rax,%r15
jne    15a0 <DenmanBeaversSqrtm+0x470>
mov    %rbp,%rdi
mov    0x8(%rsp),%rsi
mov    0x40(%rsp),%r9
mov    %rcx,%rbx
call   1030 <myinvWarm@plt>
mov    %rbx,%r11
movapd 0xa07(%rip),%xmm4
movsd  0xa0f(%rip),%xmm3
test   %r15,%r15
mov    0xa0(%rsp),%rdi
mov    0x98(%rsp),%rbx
mov    0x38(%rsp),%rsi
je     13b0 <DenmanBeaversSqrtm+0x280>
cmp    $0x4,%r15d
jae    162c <DenmanBeaversSqrtm+0x4fc>
xor    %edx,%edx
jmp    171f <DenmanBeaversSqrtm+0x5ef>
cmpq   $0x0,0x60(%rsp)
je     178a <DenmanBeaversSqrtm+0x65a>
mov    0x50(%rsp),%rcx
xor    %eax,%eax
mov    0x88(%rsp),%rdx
nopw   0x0(%rax,%rax,1)
movupd -0x30(%rdx,%rax,8),%xmm0
movupd -0x20(%rdx,%rax,8),%xmm1
movupd -0x30(%r12,%rax,8),%xmm2
addpd  %xmm0,%xmm2
movupd -0x20(%r12,%rax,8),%xmm0
addpd  %xmm1,%xmm0
mulpd  %xmm4,%xmm2
mulpd  %xmm4,%xmm0
movupd %xmm2,(%r14,%rax,8)
movupd %xmm0,0x10(%r14,%rax,8)
movupd -0x10(%rdx,%rax,8),%xmm0
movupd (%rdx,%rax,8),%xmm1
movupd -0x10(%r12,%rax,8),%xmm2
addpd  %xmm0,%xmm2
movupd (%r12,%rax,8),%xmm0
addpd  %xmm1,%xmm0
mulpd  %xmm4,%xmm2
mulpd  %xmm4,%xmm0
movupd %xmm2,0x20(%r14,%rax,8)
movupd %xmm0,0x30(%r14,%rax,8)
add    $0x8,%rax
add    $0xfffffffffffffffe,%rcx
jne    1650 <DenmanBeaversSqrtm+0x520>
testb  $0x1,0x58(%rsp)
je     170e <DenmanBeaversSqrtm+0x5de>
mov    0x80(%rsp),%rcx
movupd (%rcx,%rax,8),%xmm0
movupd 0x10(%rcx,%rax,8),%xmm1
mov    0x8(%rsp),%rcx
movupd (%rcx,%rax,8),%xmm2
addpd  %xmm0,%xmm2
movupd 0x10(%rcx,%rax,8),%xmm0
addpd  %xmm1,%xmm0
mulpd  %xmm4,%xmm2
mulpd  %xmm4,%xmm0
movupd %xmm2,(%r14,%rax,8)
movupd %xmm0,0x10(%r14,%rax,8)
mov    0x68(%rsp),%rax
mov    %rax,%rdx
cmp    %r15,%rax
je     13b0 <DenmanBeaversSqrtm+0x280>
lea    (%r14,%rdx,8),%rax
mov    %r15,%rcx
sub    %rdx,%rcx
nopl   0x0(%rax)
movsd  (%rax,%rbx,1),%xmm0
addsd  (%rax,%rdi,1),%xmm0
mulsd  %xmm3,%xmm0
movsd  %xmm0,(%rax)
add    $0x8,%rax
add    $0xffffffffffffffff,%rcx
jne    1730 <DenmanBeaversSqrtm+0x600>
jmp    13b0 <DenmanBeaversSqrtm+0x280>
xor    %eax,%eax
mov    0x20(%rsp),%r14
mov    0x30(%rsp),%rdx
mov    %rsi,%rcx
mov    0x28(%rsp),%r8
movsd  0x8a3(%rip),%xmm3
movapd 0x88b(%rip),%xmm4
mov    0x18(%rsp),%r13
testb  $0x1,0x58(%rsp)
jne    1520 <DenmanBeaversSqrtm+0x3f0>
jmp    155a <DenmanBeaversSqrtm+0x42a>
xor    %eax,%eax
testb  $0x1,0x58(%rsp)
jne    16ce <DenmanBeaversSqrtm+0x59e>
jmp    170e <DenmanBeaversSqrtm+0x5de>
add    $0xa8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
define void @DenmanBeaversSqrtm(double* %0, i32 %1, i32 %2, double* %3, i8* %4) {
%6 = alloca double*
%7 = alloca i32
%8 = alloca i32
%9 = alloca double*
%10 = alloca i8*
%11 = alloca i32
%12 = alloca i32
%13 = alloca double*
%14 = alloca double*
%15 = alloca double*
%16 = alloca double*
%17 = alloca double*
%18 = alloca i32*
%19 = alloca i32*
store double* %0, double** %6
store i32 %1, i32* %7
store i32 %2, i32* %8
store double* %3, double** %9
store i8* %4, i8** %10
%20 = load double*, double** %9
%21 = load double*, double** %6
%22 = load i32, i32* %7
%23 = load i32, i32* %7
%24 = mul nsw i32 %22, %23
%25 = sext i32 %24 to i64
%26 = mul i64 %25, 8
%27 = trunc i64 %26 to i32
%28 = call i32 @memcpy(double* %20, double* %21, i32 %27)
%29 = load i8*, i8** %10
%30 = bitcast i8* %29 to double*
store double* %30, double** %13
%31 = load double*, double** %13
%32 = load i32, i32* %7
%33 = load i32, i32* %7
%34 = mul nsw i32 %32, %33
%35 = sext i32 %34 to i64
%36 = getelementptr inbounds double, double* %31, i64 %35
store double* %36, double** %14
%37 = load double*, double** %14
%38 = load i32, i32* %7
%39 = load i32, i32* %7
%40 = mul nsw i32 %38, %39
%41 = sext i32 %40 to i64
%42 = getelementptr inbounds double, double* %37, i64 %41
store double* %42, double** %15
%43 = load double*, double** %15
%44 = load i32, i32* %7
%45 = load i32, i32* %7
%46 = mul nsw i32 %44, %45
%47 = sext i32 %46 to i64
%48 = getelementptr inbounds double, double* %43, i64 %47
store double* %48, double** %16
%49 = load double*, double** %16
%50 = load i32, i32* %7
%51 = load i32, i32* %7
%52 = mul nsw i32 %50, %51
%53 = sext i32 %52 to i64
%54 = getelementptr inbounds double, double* %49, i64 %53
store double* %54, double** %17
%55 = load double*, double** %17
%56 = load i32, i32* %7
%57 = load i32, i32* %7
%58 = mul nsw i32 %56, %57
%59 = sext i32 %58 to i64
%60 = getelementptr inbounds double, double* %55, i64 %59
%61 = bitcast double* %60 to i32*
store i32* %61, i32** %18
%62 = load i32*, i32** %18
%63 = load i32, i32* %7
%64 = load i32, i32* %7
%65 = mul nsw i32 %63, %64
%66 = sext i32 %65 to i64
%67 = getelementptr inbounds i32, i32* %62, i64 %66
store i32* %67, i32** %19
%68 = load double*, double** %13
%69 = load i32, i32* %7
%70 = load i32, i32* %7
%71 = mul nsw i32 %69, %70
%72 = sext i32 %71 to i64
%73 = mul i64 %72, 8
%74 = trunc i64 %73 to i32
%75 = call i32 @memset(double* %68, i32 0, i32 %74)
store i32 0, i32* %11
br label %76
76: ; preds = %89, %5
%77 = load i32, i32* %11
%78 = load i32, i32* %7
%79 = icmp slt i32 %77, %78
br i1 %79, label %80, label %92
80: ; preds = %76
%81 = load double*, double** %13
%82 = load i32, i32* %11
%83 = load i32, i32* %7
%84 = mul nsw i32 %82, %83
%85 = load i32, i32* %11
%86 = add nsw i32 %84, %85
%87 = sext i32 %86 to i64
%88 = getelementptr inbounds double, double* %81, i64 %87
store double 1.000000e+00, double* %88
br label %89
89: ; preds = %80
%90 = load i32, i32* %11
%91 = add nsw i32 %90, 1
store i32 %91, i32* %11
br label %76
92: ; preds = %76
store i32 0, i32* %12
br label %93
93: ; preds = %184, %92
%94 = load i32, i32* %12
%95 = load i32, i32* %8
%96 = icmp slt i32 %94, %95
br i1 %96, label %97, label %187
97: ; preds = %93
%98 = load double*, double** %14
%99 = load double*, double** %9
%100 = load i32, i32* %7
%101 = load i32, i32* %7
%102 = mul nsw i32 %100, %101
%103 = sext i32 %102 to i64
%104 = mul i64 %103, 8
%105 = trunc i64 %104 to i32
%106 = call i32 @memcpy(double* %98, double* %99, i32 %105)
%107 = load double*, double** %15
%108 = load double*, double** %13
%109 = load i32, i32* %7
%110 = load i32, i32* %7
%111 = mul nsw i32 %109, %110
%112 = sext i32 %111 to i64
%113 = mul i64 %112, 8
%114 = trunc i64 %113 to i32
%115 = call i32 @memcpy(double* %107, double* %108, i32 %114)
%116 = load double*, double** %15
%117 = load double*, double** %16
%118 = load double*, double** %17
%119 = load i32*, i32** %18
%120 = load i32*, i32** %19
%121 = load i32, i32* %7
%122 = call i32 @myinvWarm(double* %116, double* %117, double* %118, i32* %119, i32* %120, i32 %121)
store i32 0, i32* %11
br label %123
123: ; preds = %146, %97
%124 = load i32, i32* %11
%125 = load i32, i32* %7
%126 = load i32, i32* %7
%127 = mul nsw i32 %125, %126
%128 = icmp slt i32 %124, %127
br i1 %128, label %129, label %149
129: ; preds = %123
%130 = load double*, double** %14
%131 = load i32, i32* %11
%132 = sext i32 %131 to i64
%133 = getelementptr inbounds double, double* %130, i64 %132
%134 = load double, double* %133
%135 = load double*, double** %16
%136 = load i32, i32* %11
%137 = sext i32 %136 to i64
%138 = getelementptr inbounds double, double* %135, i64 %137
%139 = load double, double* %138
%140 = fadd double %134, %139
%141 = fmul double %140, 5.000000e-01
%142 = load double*, double** %9
%143 = load i32, i32* %11
%144 = sext i32 %143 to i64
%145 = getelementptr inbounds double, double* %142, i64 %144
store double %141, double* %145
br label %146
146: ; preds = %129
%147 = load i32, i32* %11
%148 = add nsw i32 %147, 1
store i32 %148, i32* %11
br label %123
149: ; preds = %123
%150 = load double*, double** %14
%151 = load double*, double** %16
%152 = load double*, double** %17
%153 = load i32*, i32** %18
%154 = load i32*, i32** %19
%155 = load i32, i32* %7
%156 = call i32 @myinvWarm(double* %150, double* %151, double* %152, i32* %153, i32* %154, i32 %155)
store i32 0, i32* %11
br label %157
157: ; preds = %180, %149
%158 = load i32, i32* %11
%159 = load i32, i32* %7
%160 = load i32, i32* %7
%161 = mul nsw i32 %159, %160
%162 = icmp slt i32 %158, %161
br i1 %162, label %163, label %183
163: ; preds = %157
%164 = load double*, double** %15
%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 = load double*, double** %16
%170 = load i32, i32* %11
%171 = sext i32 %170 to i64
%172 = getelementptr inbounds double, double* %169, i64 %171
%173 = load double, double* %172
%174 = fadd double %168, %173
%175 = fmul double %174, 5.000000e-01
%176 = load double*, double** %13
%177 = load i32, i32* %11
%178 = sext i32 %177 to i64
%179 = getelementptr inbounds double, double* %176, i64 %178
store double %175, double* %179
br label %180
180: ; preds = %163
%181 = load i32, i32* %11
%182 = add nsw i32 %181, 1
store i32 %182, i32* %11
br label %157
183: ; preds = %157
br label %184
184: ; preds = %183
%185 = load i32, i32* %12
%186 = add nsw i32 %185, 1
store i32 %186, i32* %12
br label %93
187: ; preds = %93
ret void
}

LLVM IR (-O2):
define void @DenmanBeaversSqrtm(double* %0, i32 %1, i32 %2, double* %3, i8* %4) {
%6 = bitcast double* %3 to i8*
%7 = mul nsw i32 %1, %1
%8 = zext i32 %7 to i64
%9 = shl i32 %7, 3
%10 = tail call i32 @memcpy(double* %3, double* %0, i32 %9)
%11 = bitcast i8* %4 to double*
%12 = getelementptr inbounds double, double* %11, i64 %8
%13 = getelementptr inbounds double, double* %12, i64 %8
%14 = getelementptr inbounds double, double* %13, i64 %8
%15 = getelementptr inbounds double, double* %14, i64 %8
%16 = getelementptr inbounds double, double* %15, i64 %8
%17 = bitcast double* %16 to i32*
%18 = getelementptr inbounds i32, i32* %17, i64 %8
%19 = tail call i32 @memset(double* %11, i32 0, i32 %9)
%20 = icmp sgt i32 %1, 0
br i1 %20, label %21, label %42
21: ; preds = %5
%22 = zext i32 %1 to i64
%23 = add nsw i64 %22, -1
%24 = and i64 %22, 3
%25 = icmp ult i64 %23, 3
br i1 %25, label %28, label %26
26: ; preds = %21
%27 = and i64 %22, 4294967292
br label %84
28: ; preds = %84, %21
%29 = phi i64 [ 0, %21 ], [ %110, %84 ]
%30 = icmp eq i64 %24, 0
br i1 %30, label %42, label %31
31: ; preds = %28, %31
%32 = phi i64 [ %39, %31 ], [ %29, %28 ]
%33 = phi i64 [ %40, %31 ], [ 0, %28 ]
%34 = trunc i64 %32 to i32
%35 = mul nsw i32 %34, %1
%36 = sext i32 %35 to i64
%37 = add nsw i64 %32, %36
%38 = getelementptr inbounds double, double* %11, i64 %37
store double 1.000000e+00, double* %38
%39 = add nuw nsw i64 %32, 1
%40 = add i64 %33, 1
%41 = icmp eq i64 %40, %24
br i1 %41, label %42, label %31
42: ; preds = %28, %31, %5
%43 = icmp sgt i32 %2, 0
br i1 %43, label %44, label %319
44: ; preds = %42
%45 = icmp eq i32 %7, 0
%46 = getelementptr double, double* %3, i64 %8
%47 = bitcast double* %46 to i8*
%48 = mul nuw nsw i64 %8, 24
%49 = getelementptr i8, i8* %4, i64 %48
%50 = shl nuw nsw i64 %8, 5
%51 = getelementptr i8, i8* %4, i64 %50
%52 = shl nuw nsw i64 %8, 3
%53 = getelementptr i8, i8* %4, i64 %52
%54 = shl nuw nsw i64 %8, 4
%55 = getelementptr i8, i8* %4, i64 %54
%56 = and i64 %8, 2147483644
%57 = add nsw i64 %56, -4
%58 = lshr exact i64 %57, 2
%59 = add nuw nsw i64 %58, 1
%60 = icmp ult i32 %7, 4
%61 = icmp ugt i8* %51, %6
%62 = icmp ult i8* %49, %47
%63 = and i1 %61, %62
%64 = icmp ugt i8* %55, %6
%65 = icmp ult i8* %53, %47
%66 = and i1 %64, %65
%67 = or i1 %63, %66
%68 = and i64 %8, 2147483644
%69 = and i64 %59, 1
%70 = icmp eq i64 %57, 0
%71 = and i64 %59, 9223372036854775806
%72 = icmp eq i64 %69, 0
%73 = icmp eq i64 %68, %8
%74 = and i64 %8, 1
%75 = icmp eq i64 %74, 0
%76 = sub nsw i64 0, %8
%77 = icmp ult i32 %7, 4
%78 = and i64 %8, 2147483644
%79 = and i64 %59, 1
%80 = icmp eq i64 %57, 0
%81 = and i64 %59, 9223372036854775806
%82 = icmp eq i64 %79, 0
%83 = icmp eq i64 %78, %8
br label %113
84: ; preds = %84, %26
%85 = phi i64 [ 0, %26 ], [ %110, %84 ]
%86 = phi i64 [ 0, %26 ], [ %111, %84 ]
%87 = trunc i64 %85 to i32
%88 = mul nsw i32 %87, %1
%89 = sext i32 %88 to i64
%90 = add nsw i64 %85, %89
%91 = getelementptr inbounds double, double* %11, i64 %90
store double 1.000000e+00, double* %91
%92 = or i64 %85, 1
%93 = trunc i64 %92 to i32
%94 = mul nsw i32 %93, %1
%95 = sext i32 %94 to i64
%96 = add nsw i64 %92, %95
%97 = getelementptr inbounds double, double* %11, i64 %96
store double 1.000000e+00, double* %97
%98 = or i64 %85, 2
%99 = trunc i64 %98 to i32
%100 = mul nsw i32 %99, %1
%101 = sext i32 %100 to i64
%102 = add nsw i64 %98, %101
%103 = getelementptr inbounds double, double* %11, i64 %102
store double 1.000000e+00, double* %103
%104 = or i64 %85, 3
%105 = trunc i64 %104 to i32
%106 = mul nsw i32 %105, %1
%107 = sext i32 %106 to i64
%108 = add nsw i64 %104, %107
%109 = getelementptr inbounds double, double* %11, i64 %108
store double 1.000000e+00, double* %109
%110 = add nuw nsw i64 %85, 4
%111 = add i64 %86, 4
%112 = icmp eq i64 %111, %27
br i1 %112, label %28, label %84
113: ; preds = %44, %316
%114 = phi i32 [ 0, %44 ], [ %317, %316 ]
%115 = tail call i32 @memcpy(double* %12, double* %3, i32 %9)
%116 = tail call i32 @memcpy(double* %13, double* %11, i32 %9)
%117 = tail call i32 @myinvWarm(double* %13, double* %14, double* %15, i32* %17, i32* %18, i32 %1)
br i1 %45, label %207, label %118
118: ; preds = %113
%119 = select i1 %60, i1 true, i1 %67
br i1 %119, label %192, label %120
120: ; preds = %118
br i1 %70, label %168, label %121
121: ; preds = %120, %121
%122 = phi i64 [ %165, %121 ], [ 0, %120 ]
%123 = phi i64 [ %166, %121 ], [ 0, %120 ]
%124 = getelementptr inbounds double, double* %12, i64 %122
%125 = bitcast double* %124 to <2 x double>*
%126 = load <2 x double>, <2 x double>* %125
%127 = getelementptr inbounds double, double* %124, i64 2
%128 = bitcast double* %127 to <2 x double>*
%129 = load <2 x double>, <2 x double>* %128
%130 = getelementptr inbounds double, double* %14, i64 %122
%131 = bitcast double* %130 to <2 x double>*
%132 = load <2 x double>, <2 x double>* %131
%133 = getelementptr inbounds double, double* %130, i64 2
%134 = bitcast double* %133 to <2 x double>*
%135 = load <2 x double>, <2 x double>* %134
%136 = fadd <2 x double> %126, %132
%137 = fadd <2 x double> %129, %135
%138 = fmul <2 x double> %136, <double 5.000000e-01, double 5.000000e-01>
%139 = fmul <2 x double> %137, <double 5.000000e-01, double 5.000000e-01>
%140 = getelementptr inbounds double, double* %3, i64 %122
%141 = bitcast double* %140 to <2 x double>*
store <2 x double> %138, <2 x double>* %141
%142 = getelementptr inbounds double, double* %140, i64 2
%143 = bitcast double* %142 to <2 x double>*
store <2 x double> %139, <2 x double>* %143
%144 = or i64 %122, 4
%145 = getelementptr inbounds double, double* %12, i64 %144
%146 = bitcast double* %145 to <2 x double>*
%147 = load <2 x double>, <2 x double>* %146
%148 = getelementptr inbounds double, double* %145, i64 2
%149 = bitcast double* %148 to <2 x double>*
%150 = load <2 x double>, <2 x double>* %149
%151 = getelementptr inbounds double, double* %14, i64 %144
%152 = bitcast double* %151 to <2 x double>*
%153 = load <2 x double>, <2 x double>* %152
%154 = getelementptr inbounds double, double* %151, i64 2
%155 = bitcast double* %154 to <2 x double>*
%156 = load <2 x double>, <2 x double>* %155
%157 = fadd <2 x double> %147, %153
%158 = fadd <2 x double> %150, %156
%159 = fmul <2 x double> %157, <double 5.000000e-01, double 5.000000e-01>
%160 = fmul <2 x double> %158, <double 5.000000e-01, double 5.000000e-01>
%161 = getelementptr inbounds double, double* %3, i64 %144
%162 = bitcast double* %161 to <2 x double>*
store <2 x double> %159, <2 x double>* %162
%163 = getelementptr inbounds double, double* %161, i64 2
%164 = bitcast double* %163 to <2 x double>*
store <2 x double> %160, <2 x double>* %164
%165 = add nuw i64 %122, 8
%166 = add i64 %123, 2
%167 = icmp eq i64 %166, %71
br i1 %167, label %168, label %121
168: ; preds = %121, %120
%169 = phi i64 [ 0, %120 ], [ %165, %121 ]
br i1 %72, label %191, label %170
170: ; preds = %168
%171 = getelementptr inbounds double, double* %12, i64 %169
%172 = bitcast double* %171 to <2 x double>*
%173 = load <2 x double>, <2 x double>* %172
%174 = getelementptr inbounds double, double* %171, i64 2
%175 = bitcast double* %174 to <2 x double>*
%176 = load <2 x double>, <2 x double>* %175
%177 = getelementptr inbounds double, double* %14, i64 %169
%178 = bitcast double* %177 to <2 x double>*
%179 = load <2 x double>, <2 x double>* %178
%180 = getelementptr inbounds double, double* %177, i64 2
%181 = bitcast double* %180 to <2 x double>*
%182 = load <2 x double>, <2 x double>* %181
%183 = fadd <2 x double> %173, %179
%184 = fadd <2 x double> %176, %182
%185 = fmul <2 x double> %183, <double 5.000000e-01, double 5.000000e-01>
%186 = fmul <2 x double> %184, <double 5.000000e-01, double 5.000000e-01>
%187 = getelementptr inbounds double, double* %3, i64 %169
%188 = bitcast double* %187 to <2 x double>*
store <2 x double> %185, <2 x double>* %188
%189 = getelementptr inbounds double, double* %187, i64 2
%190 = bitcast double* %189 to <2 x double>*
store <2 x double> %186, <2 x double>* %190
br label %191
191: ; preds = %168, %170
br i1 %73, label %228, label %192
192: ; preds = %118, %191
%193 = phi i64 [ 0, %118 ], [ %68, %191 ]
%194 = xor i64 %193, -1
br i1 %75, label %204, label %195
195: ; preds = %192
%196 = getelementptr inbounds double, double* %12, i64 %193
%197 = load double, double* %196
%198 = getelementptr inbounds double, double* %14, i64 %193
%199 = load double, double* %198
%200 = fadd double %197, %199
%201 = fmul double %200, 5.000000e-01
%202 = getelementptr inbounds double, double* %3, i64 %193
store double %201, double* %202
%203 = or i64 %193, 1
br label %204
204: ; preds = %195, %192
%205 = phi i64 [ %193, %192 ], [ %203, %195 ]
%206 = icmp eq i64 %194, %76
br i1 %206, label %228, label %209
207: ; preds = %113
%208 = tail call i32 @myinvWarm(double* %12, double* %14, double* %15, i32* %17, i32* %18, i32 %1)
br label %316
209: ; preds = %204, %209
%210 = phi i64 [ %226, %209 ], [ %205, %204 ]
%211 = getelementptr inbounds double, double* %12, i64 %210
%212 = load double, double* %211
%213 = getelementptr inbounds double, double* %14, i64 %210
%214 = load double, double* %213
%215 = fadd double %212, %214
%216 = fmul double %215, 5.000000e-01
%217 = getelementptr inbounds double, double* %3, i64 %210
store double %216, double* %217
%218 = add nuw nsw i64 %210, 1
%219 = getelementptr inbounds double, double* %12, i64 %218
%220 = load double, double* %219
%221 = getelementptr inbounds double, double* %14, i64 %218
%222 = load double, double* %221
%223 = fadd double %220, %222
%224 = fmul double %223, 5.000000e-01
%225 = getelementptr inbounds double, double* %3, i64 %218
store double %224, double* %225
%226 = add nuw nsw i64 %210, 2
%227 = icmp eq i64 %226, %8
br i1 %227, label %228, label %209
228: ; preds = %204, %209, %191
%229 = tail call i32 @myinvWarm(double* %12, double* %14, double* %15, i32* %17, i32* %18, i32 %1)
br i1 %45, label %316, label %230
230: ; preds = %228
br i1 %77, label %303, label %231
231: ; preds = %230
br i1 %80, label %279, label %232
232: ; preds = %231, %232
%233 = phi i64 [ %276, %232 ], [ 0, %231 ]
%234 = phi i64 [ %277, %232 ], [ 0, %231 ]
%235 = getelementptr inbounds double, double* %13, i64 %233
%236 = bitcast double* %235 to <2 x double>*
%237 = load <2 x double>, <2 x double>* %236
%238 = getelementptr inbounds double, double* %235, i64 2
%239 = bitcast double* %238 to <2 x double>*
%240 = load <2 x double>, <2 x double>* %239
%241 = getelementptr inbounds double, double* %14, i64 %233
%242 = bitcast double* %241 to <2 x double>*
%243 = load <2 x double>, <2 x double>* %242
%244 = getelementptr inbounds double, double* %241, i64 2
%245 = bitcast double* %244 to <2 x double>*
%246 = load <2 x double>, <2 x double>* %245
%247 = fadd <2 x double> %237, %243
%248 = fadd <2 x double> %240, %246
%249 = fmul <2 x double> %247, <double 5.000000e-01, double 5.000000e-01>
%250 = fmul <2 x double> %248, <double 5.000000e-01, double 5.000000e-01>
%251 = getelementptr inbounds double, double* %11, i64 %233
%252 = bitcast double* %251 to <2 x double>*
store <2 x double> %249, <2 x double>* %252
%253 = getelementptr inbounds double, double* %251, i64 2
%254 = bitcast double* %253 to <2 x double>*
store <2 x double> %250, <2 x double>* %254
%255 = or i64 %233, 4
%256 = getelementptr inbounds double, double* %13, i64 %255
%257 = bitcast double* %256 to <2 x double>*
%258 = load <2 x double>, <2 x double>* %257
%259 = getelementptr inbounds double, double* %256, i64 2
%260 = bitcast double* %259 to <2 x double>*
%261 = load <2 x double>, <2 x double>* %260
%262 = getelementptr inbounds double, double* %14, i64 %255
%263 = bitcast double* %262 to <2 x double>*
%264 = load <2 x double>, <2 x double>* %263
%265 = getelementptr inbounds double, double* %262, i64 2
%266 = bitcast double* %265 to <2 x double>*
%267 = load <2 x double>, <2 x double>* %266
%268 = fadd <2 x double> %258, %264
%269 = fadd <2 x double> %261, %267
%270 = fmul <2 x double> %268, <double 5.000000e-01, double 5.000000e-01>
%271 = fmul <2 x double> %269, <double 5.000000e-01, double 5.000000e-01>
%272 = getelementptr inbounds double, double* %11, i64 %255
%273 = bitcast double* %272 to <2 x double>*
store <2 x double> %270, <2 x double>* %273
%274 = getelementptr inbounds double, double* %272, i64 2
%275 = bitcast double* %274 to <2 x double>*
store <2 x double> %271, <2 x double>* %275
%276 = add nuw i64 %233, 8
%277 = add i64 %234, 2
%278 = icmp eq i64 %277, %81
br i1 %278, label %279, label %232
279: ; preds = %232, %231
%280 = phi i64 [ 0, %231 ], [ %276, %232 ]
br i1 %82, label %302, label %281
281: ; preds = %279
%282 = getelementptr inbounds double, double* %13, i64 %280
%283 = bitcast double* %282 to <2 x double>*
%284 = load <2 x double>, <2 x double>* %283
%285 = getelementptr inbounds double, double* %282, i64 2
%286 = bitcast double* %285 to <2 x double>*
%287 = load <2 x double>, <2 x double>* %286
%288 = getelementptr inbounds double, double* %14, i64 %280
%289 = bitcast double* %288 to <2 x double>*
%290 = load <2 x double>, <2 x double>* %289
%291 = getelementptr inbounds double, double* %288, i64 2
%292 = bitcast double* %291 to <2 x double>*
%293 = load <2 x double>, <2 x double>* %292
%294 = fadd <2 x double> %284, %290
%295 = fadd <2 x double> %287, %293
%296 = fmul <2 x double> %294, <double 5.000000e-01, double 5.000000e-01>
%297 = fmul <2 x double> %295, <double 5.000000e-01, double 5.000000e-01>
%298 = getelementptr inbounds double, double* %11, i64 %280
%299 = bitcast double* %298 to <2 x double>*
store <2 x double> %296, <2 x double>* %299
%300 = getelementptr inbounds double, double* %298, i64 2
%301 = bitcast double* %300 to <2 x double>*
store <2 x double> %297, <2 x double>* %301
br label %302
302: ; preds = %279, %281
br i1 %83, label %316, label %303
303: ; preds = %230, %302
%304 = phi i64 [ 0, %230 ], [ %78, %302 ]
br label %305
305: ; preds = %303, %305
%306 = phi i64 [ %314, %305 ], [ %304, %303 ]
%307 = getelementptr inbounds double, double* %13, i64 %306
%308 = load double, double* %307
%309 = getelementptr inbounds double, double* %14, i64 %306
%310 = load double, double* %309
%311 = fadd double %308, %310
%312 = fmul double %311, 5.000000e-01
%313 = getelementptr inbounds double, double* %11, i64 %306
store double %312, double* %313
%314 = add nuw nsw i64 %306, 1
%315 = icmp eq i64 %314, %8
br i1 %315, label %316, label %305
316: ; preds = %305, %302, %207, %228
%317 = add nuw nsw i32 %114, 1
%318 = icmp eq i32 %317, %2
br i1 %318, label %319, label %113
319: ; preds = %316, %42
ret void
}