void polydiv(double p[], int m, double q[], int n, double s[], int k,
             double r[], int l) {
  int i, j, mp, mrp;
  for (i = 0; i <= k - 1; i++)
    s[i] = 0.0;
  if (q[n - 1] != 0.0) {
    for (i = k, mp = 1; i >= 1; i--, mp++) {
      s[i - 1] = p[m - mp] / q[n - 1];
      for (mrp = m - mp, j = 1; j <= n - 1; j++) {
        p[mrp - j] = p[mrp - j] - s[i - 1] * q[n - 1 - j];
      }
    }
    for (i = 0; i <= l - 1; i++)
      r[i] = p[i];
  }
}


x86 assembly (clang, -O0):
<polydiv>:
push   %rbp
mov    %rsp,%rbp
mov    0x18(%rbp),%eax
mov    0x10(%rbp),%rax
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %rdx,-0x18(%rbp)
mov    %ecx,-0x1c(%rbp)
mov    %r8,-0x28(%rbp)
mov    %r9d,-0x2c(%rbp)
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
mov    -0x2c(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jg     1157 <polydiv+0x57>
mov    -0x28(%rbp),%rax
movslq -0x30(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    1128 <polydiv+0x28>
mov    -0x18(%rbp),%rax
mov    -0x1c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
xorps  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    1181 <polydiv+0x81>
jp     1181 <polydiv+0x81>
jmp    12c8 <polydiv+0x1c8>
mov    -0x2c(%rbp),%eax
mov    %eax,-0x30(%rbp)
movl   $0x1,-0x38(%rbp)
cmpl   $0x1,-0x30(%rbp)
jl     1283 <polydiv+0x183>
mov    -0x8(%rbp),%rax
mov    -0xc(%rbp),%ecx
sub    -0x38(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x1c(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
divsd  (%rax,%rcx,8),%xmm0
mov    -0x28(%rbp),%rax
mov    -0x30(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0xc(%rbp),%eax
sub    -0x38(%rbp),%eax
mov    %eax,-0x3c(%rbp)
movl   $0x1,-0x34(%rbp)
mov    -0x34(%rbp),%eax
mov    -0x1c(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jg     1267 <polydiv+0x167>
mov    -0x8(%rbp),%rax
mov    -0x3c(%rbp),%ecx
sub    -0x34(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm1
mov    -0x28(%rbp),%rax
mov    -0x30(%rbp),%ecx
sub    $0x1,%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    -0x18(%rbp),%rax
mov    -0x1c(%rbp),%ecx
sub    $0x1,%ecx
sub    -0x34(%rbp),%ecx
movslq %ecx,%rcx
movsd  (%rax,%rcx,8),%xmm2
movq   %xmm0,%rax
movabs $0x8000000000000000,%rcx
xor    %rcx,%rax
movq   %rax,%xmm0
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
mov    -0x8(%rbp),%rax
mov    -0x3c(%rbp),%ecx
sub    -0x34(%rbp),%ecx
movslq %ecx,%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    11de <polydiv+0xde>
jmp    126c <polydiv+0x16c>
mov    -0x30(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x30(%rbp)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    118e <polydiv+0x8e>
movl   $0x0,-0x30(%rbp)
mov    -0x30(%rbp),%eax
mov    0x18(%rbp),%ecx
sub    $0x1,%ecx
cmp    %ecx,%eax
jg     12c3 <polydiv+0x1c3>
mov    -0x8(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm0
mov    0x10(%rbp),%rax
movslq -0x30(%rbp),%rcx
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x30(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x30(%rbp)
jmp    128a <polydiv+0x18a>
jmp    12c8 <polydiv+0x1c8>
pop    %rbp
ret

x86 assembly (clang, -O2):
<polydiv>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
sub    $0x88,%rsp
mov    %ecx,%r15d
mov    %rdx,%r12
mov    %rdi,%rbx
mov    0xc8(%rsp),%edx
test   %r9d,%r9d
jle    14b0 <polydiv+0x3a0>
mov    %r9d,%ebp
mov    %r8,%r13
mov    %esi,0x10(%rsp)
mov    %r9d,%r14d
lea    0x0(,%r14,8),%rdx
mov    %r8,%rdi
xor    %esi,%esi
call   1030 <memset@plt>
movslq %r15d,%rdx
movsd  -0x8(%r12,%rdx,8),%xmm1
xorpd  %xmm0,%xmm0
ucomisd %xmm0,%xmm1
jne    1173 <polydiv+0x63>
jnp    1649 <polydiv+0x539>
add    $0xffffffffffffffff,%rdx
movslq 0x10(%rsp),%rcx
mov    %rbp,0x78(%rsp)
lea    0x1(%rbp),%eax
mov    %rax,0x80(%rsp)
mov    %r15d,%r11d
lea    -0x3(%r11),%rax
lea    0x8(%r13),%rsi
mov    %rsi,0x70(%rsp)
mov    %rcx,0x38(%rsp)
add    $0xfffffffffffffffe,%rcx
mov    %rcx,0x68(%rsp)
mov    %rdx,%rcx
sub    %r11,%rcx
lea    (%r12,%rcx,8),%rcx
add    $0x8,%rcx
mov    %rcx,0x60(%rsp)
mov    %rdx,0x10(%rsp)
lea    (%r12,%rdx,8),%rdi
add    $0xfffffffffffffff8,%rdi
lea    -0x2(%r11),%rdx
cmp    $0x7ffffffe,%eax
setae  %cl
shr    $0x20,%rax
setne  %al
or     %cl,%al
mov    %al,0xf(%rsp)
mov    %rdx,0x28(%rsp)
and    $0xfffffffffffffffc,%rdx
lea    0x2(%rdx),%rax
mov    %rax,0x50(%rsp)
mov    %r11,%rax
mov    %r11,0x30(%rsp)
neg    %r11
mov    %rdx,0x58(%rsp)
neg    %rdx
mov    %rdx,0x48(%rsp)
mov    $0x1,%r8d
movapd 0xde5(%rip),%xmm0
xor    %r9d,%r9d
mov    %r12,%rbp
mov    %r15d,0x24(%rsp)
mov    %r12,0x40(%rsp)
nopl   0x0(%rax,%rax,1)
mov    0x38(%rsp),%r10
sub    %r8,%r10
add    $0xffffffffffffffff,%r14
movsd  (%rbx,%r10,8),%xmm2
divsd  %xmm1,%xmm2
mov    %r14d,%esi
movsd  %xmm2,0x0(%r13,%rsi,8)
cmp    $0x2,%r15d
jl     1480 <polydiv+0x370>
movsd  -0x8(%rbx,%r10,8),%xmm1
mov    0x10(%rsp),%rax
mulsd  -0x8(%rbp,%rax,8),%xmm2
subsd  %xmm2,%xmm1
movsd  %xmm1,-0x8(%rbx,%r10,8)
je     1480 <polydiv+0x370>
cmpq   $0x4,0x28(%rsp)
setb   %al
or     0xf(%rsp),%al
mov    $0x2,%r15d
je     1364 <polydiv+0x254>
mov    0x30(%rsp),%rax
sub    %r15d,%eax
mov    %r15,%r12
test   $0x1,%al
je     12d5 <polydiv+0x1c5>
movsd  0x0(%r13,%rsi,8),%xmm1
movslq %r15d,%rax
mov    %r10,%rcx
sub    %rax,%rcx
mov    0x10(%rsp),%rax
sub    %r15,%rax
mulsd  0x0(%rbp,%rax,8),%xmm1
movsd  (%rbx,%rcx,8),%xmm2
subsd  %xmm1,%xmm2
movsd  %xmm2,(%rbx,%rcx,8)
lea    0x1(%r15),%r12
not    %r15
cmp    %r11,%r15
je     1480 <polydiv+0x370>
mov    %r12,%rax
neg    %rax
shl    $0x20,%r12
nopl   0x0(%rax,%rax,1)
movsd  0x0(%r13,%rsi,8),%xmm1
mov    %r12,%rcx
sar    $0x20,%rcx
mov    %r10,%rdx
sub    %rcx,%rdx
movsd  (%rbx,%rdx,8),%xmm2
mulsd  0x8(%rdi,%rax,8),%xmm1
subsd  %xmm1,%xmm2
movsd  %xmm2,(%rbx,%rdx,8)
movsd  0x0(%r13,%rsi,8),%xmm1
movabs $0x100000000,%rcx
add    %r12,%rcx
sar    $0x20,%rcx
mov    %r10,%rdx
sub    %rcx,%rdx
movsd  (%rbx,%rdx,8),%xmm2
mulsd  (%rdi,%rax,8),%xmm1
subsd  %xmm1,%xmm2
movsd  %xmm2,(%rbx,%rdx,8)
add    $0xfffffffffffffffe,%rax
movabs $0x200000000,%rcx
add    %rcx,%r12
cmp    %rax,%r11
jne    12f0 <polydiv+0x1e0>
jmp    1480 <polydiv+0x370>
mov    %r9d,%eax
not    %eax
add    0x78(%rsp),%eax
mov    0x30(%rsp),%rcx
lea    (%r9,%rcx,1),%r12
mov    0x38(%rsp),%rcx
sub    %r12,%rcx
mov    0x68(%rsp),%r12
sub    %r9,%r12
lea    0x0(,%rax,8),%rdx
add    %r13,%rdx
mov    %rdx,0x18(%rsp)
mov    0x70(%rsp),%rdx
lea    (%rdx,%rax,8),%rax
lea    (%rbx,%rcx,8),%rcx
lea    (%rbx,%r12,8),%r12
cmp    %r12,0x18(%rsp)
setb   0x18(%rsp)
cmp    %rax,%rcx
setb   %dl
cmp    %r12,0x60(%rsp)
setb   %al
cmp    %rdi,%rcx
setb   %cl
test   %dl,0x18(%rsp)
jne    1296 <polydiv+0x186>
and    %cl,%al
jne    1296 <polydiv+0x186>
movsd  0x0(%r13,%rsi,8),%xmm1
unpcklpd %xmm1,%xmm1
xorpd  %xmm0,%xmm1
xor    %eax,%eax
movabs $0x200000000,%rbp
mov    0x48(%rsp),%r15
movabs $0x400000000,%r12
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
mov    %rbp,%rcx
sar    $0x20,%rcx
mov    %r10,%rdx
sub    %rcx,%rdx
movupd -0x18(%rbx,%rdx,8),%xmm2
movupd -0x8(%rbx,%rdx,8),%xmm3
movupd -0x20(%rdi,%rax,8),%xmm4
movupd -0x10(%rdi,%rax,8),%xmm5
movapd %xmm1,%xmm6
shufpd $0x1,%xmm1,%xmm6
mulpd  %xmm6,%xmm5
addpd  %xmm3,%xmm5
movupd %xmm5,-0x8(%rbx,%rdx,8)
mulpd  %xmm4,%xmm6
addpd  %xmm2,%xmm6
movupd %xmm6,-0x18(%rbx,%rdx,8)
add    %r12,%rbp
add    $0xfffffffffffffffc,%rax
cmp    %rax,%r15
jne    1410 <polydiv+0x300>
mov    0x50(%rsp),%r15
mov    0x58(%rsp),%rax
cmp    %rax,0x28(%rsp)
mov    0x40(%rsp),%rbp
jne    1296 <polydiv+0x186>
add    $0x1,%r8
cmp    0x80(%rsp),%r8
je     14a7 <polydiv+0x397>
mov    0x10(%rsp),%rax
movsd  0x0(%rbp,%rax,8),%xmm1
add    $0x1,%r9
mov    0x24(%rsp),%r15d
jmp    1230 <polydiv+0x120>
mov    0xc8(%rsp),%edx
jmp    14ca <polydiv+0x3ba>
movslq %r15d,%rax
movsd  -0x8(%r12,%rax,8),%xmm0
xorpd  %xmm1,%xmm1
ucomisd %xmm1,%xmm0
jne    14ca <polydiv+0x3ba>
jnp    1649 <polydiv+0x539>
test   %edx,%edx
jle    1649 <polydiv+0x539>
mov    0xc0(%rsp),%rax
mov    %edx,%ecx
cmp    $0x4,%edx
jb     14f7 <polydiv+0x3e7>
lea    (%rbx,%rcx,8),%rdx
cmp    %rax,%rdx
jbe    156c <polydiv+0x45c>
lea    (%rax,%rcx,8),%rdx
cmp    %rbx,%rdx
jbe    156c <polydiv+0x45c>
xor    %edx,%edx
mov    %rdx,%rsi
not    %rsi
add    %rcx,%rsi
mov    %rcx,%rdi
and    $0x3,%rdi
je     1524 <polydiv+0x414>
nopl   0x0(%rax,%rax,1)
movsd  (%rbx,%rdx,8),%xmm0
movsd  %xmm0,(%rax,%rdx,8)
add    $0x1,%rdx
add    $0xffffffffffffffff,%rdi
jne    1510 <polydiv+0x400>
cmp    $0x3,%rsi
jb     1649 <polydiv+0x539>
xchg   %ax,%ax
movsd  (%rbx,%rdx,8),%xmm0
movsd  %xmm0,(%rax,%rdx,8)
movsd  0x8(%rbx,%rdx,8),%xmm0
movsd  %xmm0,0x8(%rax,%rdx,8)
movsd  0x10(%rbx,%rdx,8),%xmm0
movsd  %xmm0,0x10(%rax,%rdx,8)
movsd  0x18(%rbx,%rdx,8),%xmm0
movsd  %xmm0,0x18(%rax,%rdx,8)
add    $0x4,%rdx
cmp    %rdx,%rcx
jne    1530 <polydiv+0x420>
jmp    1649 <polydiv+0x539>
mov    %ecx,%edx
and    $0xfffffffc,%edx
lea    -0x4(%rdx),%rsi
mov    %rsi,%rdi
shr    $0x2,%rdi
add    $0x1,%rdi
mov    %edi,%r8d
and    $0x3,%r8d
cmp    $0xc,%rsi
jae    1591 <polydiv+0x481>
xor    %ebp,%ebp
jmp    15fc <polydiv+0x4ec>
and    $0xfffffffffffffffc,%rdi
xor    %ebp,%ebp
nopw   0x0(%rax,%rax,1)
movups (%rbx,%rbp,8),%xmm0
movups 0x10(%rbx,%rbp,8),%xmm1
movups %xmm0,(%rax,%rbp,8)
movups %xmm1,0x10(%rax,%rbp,8)
movups 0x20(%rbx,%rbp,8),%xmm0
movups 0x30(%rbx,%rbp,8),%xmm1
movups %xmm0,0x20(%rax,%rbp,8)
movups %xmm1,0x30(%rax,%rbp,8)
movups 0x40(%rbx,%rbp,8),%xmm0
movups 0x50(%rbx,%rbp,8),%xmm1
movups %xmm0,0x40(%rax,%rbp,8)
movups %xmm1,0x50(%rax,%rbp,8)
movupd 0x60(%rbx,%rbp,8),%xmm0
movupd 0x70(%rbx,%rbp,8),%xmm1
movupd %xmm0,0x60(%rax,%rbp,8)
movupd %xmm1,0x70(%rax,%rbp,8)
add    $0x10,%rbp
add    $0xfffffffffffffffc,%rdi
jne    15a0 <polydiv+0x490>
test   %r8,%r8
je     1640 <polydiv+0x530>
lea    (%rax,%rbp,8),%rdi
add    $0x10,%rdi
lea    (%rbx,%rbp,8),%rbp
add    $0x10,%rbp
shl    $0x5,%r8
xor    %esi,%esi
nopw   0x0(%rax,%rax,1)
movupd -0x10(%rbp,%rsi,1),%xmm0
movupd 0x0(%rbp,%rsi,1),%xmm1
movupd %xmm0,-0x10(%rdi,%rsi,1)
movupd %xmm1,(%rdi,%rsi,1)
add    $0x20,%rsi
cmp    %rsi,%r8
jne    1620 <polydiv+0x510>
cmp    %rcx,%rdx
jne    14f9 <polydiv+0x3e9>
add    $0x88,%rsp
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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