#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*/ fprintf(int /*<<< orphan*/, char *, int, int);
int /*<<< orphan*/ stderr;
void ILU(int n, const double *a, double *lu) {
  int k;
  for (k = 0; k < n * n; k++) {
    lu[k] = a[k];
  }
  for (k = 0; k < (n - 1); k++) {
    if (lu[k * n + k] == 0.0) {
      fprintf(stderr, "ILU : breakdown a[%d,%d] = 0\n", k, k);
    }
    double d = 1.0 / lu[k * n + k];
    int i;
    for (i = (k + 1); i < n; i++) {
      double e = d * lu[i * n + k];
      lu[i * n + k] = e;
      int j;
      for (j = (k + 1); j < n; j++) {
        lu[i * n + j] -= e * lu[k * n + j];
      }
    }
  }
}


x86 assembly (clang, -O0):
<ILU>:
push   %rbp
mov    %rsp,%rbp
sub    $0x40,%rsp
mov    %edi,-0x4(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %rdx,-0x18(%rbp)
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
mov    -0x4(%rbp),%ecx
imul   -0x4(%rbp),%ecx
cmp    %ecx,%eax
jge    1164 <ILU+0x54>
mov    -0x10(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
movslq -0x1c(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    112a <ILU+0x1a>
movl   $0x0,-0x1c(%rbp)
mov    -0x1c(%rbp),%eax
mov    -0x4(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jge    12d1 <ILU+0x1c1>
mov    -0x18(%rbp),%rax
mov    -0x1c(%rbp),%ecx
imul   -0x4(%rbp),%ecx
add    -0x1c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    11c0 <ILU+0xb0>
jp     11c0 <ILU+0xb0>
mov    0x2e34(%rip),%rax
mov    (%rax),%edi
mov    -0x1c(%rbp),%edx
mov    -0x1c(%rbp),%ecx
lea    0xe4d(%rip),%rsi
call   1030 <fprintf@plt>
mov    -0x18(%rbp),%rax
mov    -0x1c(%rbp),%ecx
imul   -0x4(%rbp),%ecx
add    -0x1c(%rbp),%ecx
movslq %ecx,%rcx
movsd  0xe27(%rip),%xmm0
divsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x28(%rbp)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
mov    -0x2c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    12be <ILU+0x1ae>
movsd  -0x28(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0x4(%rbp),%ecx
add    -0x1c(%rbp),%ecx
movslq %ecx,%rcx
mulsd  (%rax,%rcx,8),%xmm0
movsd  %xmm0,-0x38(%rbp)
movsd  -0x38(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0x4(%rbp),%ecx
add    -0x1c(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
mov    -0x3c(%rbp),%eax
cmp    -0x4(%rbp),%eax
jge    12ab <ILU+0x19b>
movsd  -0x38(%rbp),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x1c(%rbp),%ecx
imul   -0x4(%rbp),%ecx
add    -0x3c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
mov    -0x18(%rbp),%rax
mov    -0x2c(%rbp),%ecx
imul   -0x4(%rbp),%ecx
add    -0x3c(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
movq   %xmm0,%rdx
movabs $0x8000000000000000,%rsi
xor    %rsi,%rdx
movq   %rdx,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x3c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x3c(%rbp)
jmp    123c <ILU+0x12c>
jmp    12b0 <ILU+0x1a0>
mov    -0x2c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x2c(%rbp)
jmp    11ec <ILU+0xdc>
jmp    12c3 <ILU+0x1b3>
mov    -0x1c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x1c(%rbp)
jmp    116b <ILU+0x5b>
add    $0x40,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<ILU>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0xa8,%rsp
mov    %rdx,%r12
mov    %rdi,%rax
mov    %rdi,0x10(%rsp)
mov    %edi,%ecx
imul   %edi,%ecx
test   %ecx,%ecx
je     12a9 <ILU+0x199>
mov    %ecx,%eax
cmp    $0x4,%ecx
jb     1156 <ILU+0x46>
lea    (%rsi,%rax,8),%rcx
cmp    %r12,%rcx
jbe    11d0 <ILU+0xc0>
lea    (%r12,%rax,8),%rcx
cmp    %rsi,%rcx
jbe    11d0 <ILU+0xc0>
xor    %ecx,%ecx
mov    %rcx,%rdx
not    %rdx
add    %rax,%rdx
mov    %rax,%rdi
and    $0x3,%rdi
je     1185 <ILU+0x75>
nopw   0x0(%rax,%rax,1)
movsd  (%rsi,%rcx,8),%xmm0
movsd  %xmm0,(%r12,%rcx,8)
add    $0x1,%rcx
add    $0xffffffffffffffff,%rdi
jne    1170 <ILU+0x60>
cmp    $0x3,%rdx
jb     12a9 <ILU+0x199>
nop
movsd  (%rsi,%rcx,8),%xmm0
movsd  %xmm0,(%r12,%rcx,8)
movsd  0x8(%rsi,%rcx,8),%xmm0
movsd  %xmm0,0x8(%r12,%rcx,8)
movsd  0x10(%rsi,%rcx,8),%xmm0
movsd  %xmm0,0x10(%r12,%rcx,8)
movsd  0x18(%rsi,%rcx,8),%xmm0
movsd  %xmm0,0x18(%r12,%rcx,8)
add    $0x4,%rcx
cmp    %rcx,%rax
jne    1190 <ILU+0x80>
jmp    12a9 <ILU+0x199>
mov    %eax,%ecx
and    $0xfffffffc,%ecx
lea    -0x4(%rcx),%rbp
mov    %rbp,%rdi
shr    $0x2,%rdi
add    $0x1,%rdi
mov    %edi,%edx
and    $0x3,%edx
cmp    $0xc,%rbp
jae    11f3 <ILU+0xe3>
xor    %ebp,%ebp
jmp    1264 <ILU+0x154>
and    $0xfffffffffffffffc,%rdi
xor    %ebp,%ebp
nopl   0x0(%rax)
movups (%rsi,%rbp,8),%xmm0
movups 0x10(%rsi,%rbp,8),%xmm1
movups %xmm0,(%r12,%rbp,8)
movups %xmm1,0x10(%r12,%rbp,8)
movups 0x20(%rsi,%rbp,8),%xmm0
movups 0x30(%rsi,%rbp,8),%xmm1
movups %xmm0,0x20(%r12,%rbp,8)
movups %xmm1,0x30(%r12,%rbp,8)
movups 0x40(%rsi,%rbp,8),%xmm0
movups 0x50(%rsi,%rbp,8),%xmm1
movups %xmm0,0x40(%r12,%rbp,8)
movups %xmm1,0x50(%r12,%rbp,8)
movupd 0x60(%rsi,%rbp,8),%xmm0
movupd 0x70(%rsi,%rbp,8),%xmm1
movupd %xmm0,0x60(%r12,%rbp,8)
movupd %xmm1,0x70(%r12,%rbp,8)
add    $0x10,%rbp
add    $0xfffffffffffffffc,%rdi
jne    1200 <ILU+0xf0>
test   %rdx,%rdx
je     12a0 <ILU+0x190>
lea    (%r12,%rbp,8),%rdi
add    $0x10,%rdi
lea    (%rsi,%rbp,8),%rbp
add    $0x10,%rbp
shl    $0x5,%rdx
xor    %ebx,%ebx
nop
movupd -0x10(%rbp,%rbx,1),%xmm0
movupd 0x0(%rbp,%rbx,1),%xmm1
movupd %xmm0,-0x10(%rdi,%rbx,1)
movupd %xmm1,(%rdi,%rbx,1)
add    $0x20,%rbx
cmp    %rbx,%rdx
jne    1280 <ILU+0x170>
cmp    %rax,%rcx
jne    1158 <ILU+0x48>
cmpl   $0x2,0x10(%rsp)
jl     15eb <ILU+0x4db>
mov    0x10(%rsp),%rax
lea    -0x1(%rax),%ecx
mov    %rcx,0x50(%rsp)
mov    %eax,%r14d
lea    0x1(%r14),%rax
mov    %rax,0x88(%rsp)
lea    0x0(,%r14,8),%rax
mov    %rax,0x38(%rsp)
lea    0x8(,%r14,8),%rax
mov    %rax,0x30(%rsp)
lea    (%r14,%r14,1),%rax
mov    %rax,0x48(%rsp)
lea    0x8(%r12),%rax
mov    %rax,0x40(%rsp)
mov    %r14,%rax
neg    %rax
mov    %rax,0xa0(%rsp)
lea    -0x1(%r14),%rax
mov    %rax,0x28(%rsp)
lea    0x18(%r12),%r15
lea    0x18(%r12,%r14,8),%rax
mov    %rax,0x20(%rsp)
lea    0x8(%r12,%r14,8),%rax
mov    %rax,0x18(%rsp)
mov    $0x1,%eax
mov    %rax,0x8(%rsp)
xor    %ebp,%ebp
movapd 0xccf(%rip),%xmm7
mov    %r12,%rbx
jmp    1386 <ILU+0x276>
cs nopw 0x0(%rax,%rax,1)
addq   $0x1,0x8(%rsp)
addq   $0xffffffffffffffff,0x28(%rsp)
mov    0x30(%rsp),%rax
add    %rax,%r15
add    %rax,0x20(%rsp)
mov    0x38(%rsp),%rax
add    %rax,0x18(%rsp)
add    %rax,%rbx
mov    0x58(%rsp),%rbp
cmp    0x50(%rsp),%rbp
je     15eb <ILU+0x4db>
mov    %ebp,%eax
imul   0x10(%rsp),%eax
movslq %eax,%r13
add    %rbp,%r13
movsd  (%r12,%r13,8),%xmm1
ucomisd 0xc67(%rip),%xmm1
jne    13cc <ILU+0x2bc>
jp     13cc <ILU+0x2bc>
mov    0x2c34(%rip),%rax
mov    (%rax),%edi
lea    0xc6b(%rip),%rsi
mov    %ebp,%edx
mov    %ebp,%ecx
call   1030 <fprintf@plt>
movapd 0xc4a(%rip),%xmm7
movsd  (%r12,%r13,8),%xmm1
lea    0x1(%rbp),%rax
mov    %rax,0x58(%rsp)
cmp    %r14,%rax
jae    1350 <ILU+0x240>
mov    %rbp,%rdx
imul   %r14,%rdx
mov    0x30(%rsp),%rax
imul   %rbp,%rax
lea    (%rdx,%r14,1),%r8
mov    0x28(%rsp),%rsi
and    $0xfffffffffffffffc,%rsi
mov    %rbp,%rdi
not    %rdi
add    %r14,%rdi
mov    0x48(%rsp),%rcx
mov    %rdx,0x80(%rsp)
add    %rdx,%rcx
mov    %rcx,0x78(%rsp)
mov    0x40(%rsp),%rcx
mov    %rax,0x98(%rsp)
add    %rcx,%rax
mov    %rax,0x70(%rsp)
lea    (%r12,%r8,8),%rax
mov    %rax,0x90(%rsp)
movsd  0xbbe(%rip),%xmm0
divsd  %xmm1,%xmm0
mov    %rdi,%rax
and    $0xfffffffffffffffc,%rax
mov    0x8(%rsp),%r10
mov    %rax,0x68(%rsp)
add    %r10,%rax
mov    %rax,0x60(%rsp)
mov    0x18(%rsp),%r8
mov    0x20(%rsp),%r11
xor    %r9d,%r9d
jmp    148c <ILU+0x37c>
xchg   %ax,%ax
add    $0x1,%r10
add    $0x1,%r9
mov    0x38(%rsp),%rax
add    %rax,%r11
add    %rax,%r8
cmp    %r14,%r10
je     1350 <ILU+0x240>
mov    %r10,%r13
imul   %r14,%r13
mov    %r13,%rax
add    %rbp,%rax
movsd  (%r12,%rax,8),%xmm1
mulsd  %xmm0,%xmm1
movsd  %xmm1,(%r12,%rax,8)
xorpd  %xmm7,%xmm1
mov    0x8(%rsp),%rcx
cmp    $0x4,%rdi
jb     1560 <ILU+0x450>
mov    %r9,%rax
imul   %r14,%rax
mov    0x88(%rsp),%rcx
add    %rax,%rcx
lea    (%r12,%rcx,8),%rcx
add    0x98(%rsp),%rcx
cmp    0x90(%rsp),%rcx
jae    14f9 <ILU+0x3e9>
add    0x78(%rsp),%rax
lea    (%r12,%rax,8),%rax
mov    0x8(%rsp),%rcx
cmp    %rax,0x70(%rsp)
jb     1560 <ILU+0x450>
movapd %xmm1,%xmm2
unpcklpd %xmm1,%xmm2
xor    %eax,%eax
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movupd -0x10(%r15,%rax,8),%xmm3
movupd (%r15,%rax,8),%xmm4
movupd -0x10(%r11,%rax,8),%xmm5
movupd (%r11,%rax,8),%xmm6
mulpd  %xmm2,%xmm3
addpd  %xmm5,%xmm3
mulpd  %xmm2,%xmm4
addpd  %xmm6,%xmm4
movupd %xmm3,-0x10(%r11,%rax,8)
movupd %xmm4,(%r11,%rax,8)
add    $0x4,%rax
cmp    %rax,%rsi
jne    1510 <ILU+0x400>
mov    0x60(%rsp),%rcx
cmp    0x68(%rsp),%rdi
je     1470 <ILU+0x360>
mov    %rbp,%rdx
mov    %r14d,%ebp
sub    %ecx,%ebp
mov    %rcx,%rax
test   $0x1,%bpl
je     1599 <ILU+0x489>
mov    0x80(%rsp),%rax
add    %rcx,%rax
add    %rcx,%r13
movsd  (%r12,%rax,8),%xmm2
mulsd  %xmm1,%xmm2
addsd  (%r12,%r13,8),%xmm2
movsd  %xmm2,(%r12,%r13,8)
lea    0x1(%rcx),%rax
not    %rcx
cmp    0xa0(%rsp),%rcx
mov    %rdx,%rbp
je     1470 <ILU+0x360>
nopl   (%rax)
movsd  (%rbx,%rax,8),%xmm2
mulsd  %xmm1,%xmm2
addsd  -0x8(%r8,%rax,8),%xmm2
movsd  %xmm2,-0x8(%r8,%rax,8)
movsd  0x8(%rbx,%rax,8),%xmm2
mulsd  %xmm1,%xmm2
addsd  (%r8,%rax,8),%xmm2
movsd  %xmm2,(%r8,%rax,8)
add    $0x2,%rax
cmp    %rax,%r14
jne    15b0 <ILU+0x4a0>
jmp    1470 <ILU+0x360>
add    $0xa8,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

LLVM IR (-O0):
@stderr = global i32 0
@.str = [30 x i8] c"ILU : breakdown a[%d,%d] = 0\0A\00"
define void @ILU(i32 %0, double* %1, double* %2) {
%4 = alloca i32
%5 = alloca double*
%6 = alloca double*
%7 = alloca i32
%8 = alloca double
%9 = alloca i32
%10 = alloca double
%11 = alloca i32
store i32 %0, i32* %4
store double* %1, double** %5
store double* %2, double** %6
store i32 0, i32* %7
br label %12
12: ; preds = %28, %3
%13 = load i32, i32* %7
%14 = load i32, i32* %4
%15 = load i32, i32* %4
%16 = mul nsw i32 %14, %15
%17 = icmp slt i32 %13, %16
br i1 %17, label %18, label %31
18: ; preds = %12
%19 = load double*, double** %5
%20 = load i32, i32* %7
%21 = sext i32 %20 to i64
%22 = getelementptr inbounds double, double* %19, i64 %21
%23 = load double, double* %22
%24 = load double*, double** %6
%25 = load i32, i32* %7
%26 = sext i32 %25 to i64
%27 = getelementptr inbounds double, double* %24, i64 %26
store double %23, double* %27
br label %28
28: ; preds = %18
%29 = load i32, i32* %7
%30 = add nsw i32 %29, 1
store i32 %30, i32* %7
br label %12
31: ; preds = %12
store i32 0, i32* %7
br label %32
32: ; preds = %127, %31
%33 = load i32, i32* %7
%34 = load i32, i32* %4
%35 = sub nsw i32 %34, 1
%36 = icmp slt i32 %33, %35
br i1 %36, label %37, label %130
37: ; preds = %32
%38 = load double*, double** %6
%39 = load i32, i32* %7
%40 = load i32, i32* %4
%41 = mul nsw i32 %39, %40
%42 = load i32, i32* %7
%43 = add nsw i32 %41, %42
%44 = sext i32 %43 to i64
%45 = getelementptr inbounds double, double* %38, i64 %44
%46 = load double, double* %45
%47 = fcmp oeq double %46, 0.000000e+00
br i1 %47, label %48, label %53
48: ; preds = %37
%49 = load i32, i32* @stderr
%50 = load i32, i32* %7
%51 = load i32, i32* %7
%52 = call i32 @fprintf(i32 %49, i8* getelementptr inbounds ([30 x i8], [30 x i8]* @.str, i64 0, i64 0), i32 %50, i32 %51)
br label %53
53: ; preds = %48, %37
%54 = load double*, double** %6
%55 = load i32, i32* %7
%56 = load i32, i32* %4
%57 = mul nsw i32 %55, %56
%58 = load i32, i32* %7
%59 = add nsw i32 %57, %58
%60 = sext i32 %59 to i64
%61 = getelementptr inbounds double, double* %54, i64 %60
%62 = load double, double* %61
%63 = fdiv double 1.000000e+00, %62
store double %63, double* %8
%64 = load i32, i32* %7
%65 = add nsw i32 %64, 1
store i32 %65, i32* %9
br label %66
66: ; preds = %123, %53
%67 = load i32, i32* %9
%68 = load i32, i32* %4
%69 = icmp slt i32 %67, %68
br i1 %69, label %70, label %126
70: ; preds = %66
%71 = load double, double* %8
%72 = load double*, double** %6
%73 = load i32, i32* %9
%74 = load i32, i32* %4
%75 = mul nsw i32 %73, %74
%76 = load i32, i32* %7
%77 = add nsw i32 %75, %76
%78 = sext i32 %77 to i64
%79 = getelementptr inbounds double, double* %72, i64 %78
%80 = load double, double* %79
%81 = fmul double %71, %80
store double %81, double* %10
%82 = load double, double* %10
%83 = load double*, double** %6
%84 = load i32, i32* %9
%85 = load i32, i32* %4
%86 = mul nsw i32 %84, %85
%87 = load i32, i32* %7
%88 = add nsw i32 %86, %87
%89 = sext i32 %88 to i64
%90 = getelementptr inbounds double, double* %83, i64 %89
store double %82, double* %90
%91 = load i32, i32* %7
%92 = add nsw i32 %91, 1
store i32 %92, i32* %11
br label %93
93: ; preds = %119, %70
%94 = load i32, i32* %11
%95 = load i32, i32* %4
%96 = icmp slt i32 %94, %95
br i1 %96, label %97, label %122
97: ; preds = %93
%98 = load double, double* %10
%99 = load double*, double** %6
%100 = load i32, i32* %7
%101 = load i32, i32* %4
%102 = mul nsw i32 %100, %101
%103 = load i32, i32* %11
%104 = add nsw i32 %102, %103
%105 = sext i32 %104 to i64
%106 = getelementptr inbounds double, double* %99, i64 %105
%107 = load double, double* %106
%108 = load double*, double** %6
%109 = load i32, i32* %9
%110 = load i32, i32* %4
%111 = mul nsw i32 %109, %110
%112 = load i32, i32* %11
%113 = add nsw i32 %111, %112
%114 = sext i32 %113 to i64
%115 = getelementptr inbounds double, double* %108, i64 %114
%116 = load double, double* %115
%117 = fneg double %98
%118 = call double @llvm.fmuladd.f64(double %117, double %107, double %116)
store double %118, double* %115
br label %119
119: ; preds = %97
%120 = load i32, i32* %11
%121 = add nsw i32 %120, 1
store i32 %121, i32* %11
br label %93
122: ; preds = %93
br label %123
123: ; preds = %122
%124 = load i32, i32* %9
%125 = add nsw i32 %124, 1
store i32 %125, i32* %9
br label %66
126: ; preds = %66
br label %127
127: ; preds = %126
%128 = load i32, i32* %7
%129 = add nsw i32 %128, 1
store i32 %129, i32* %7
br label %32
130: ; preds = %32
ret void
}

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