void f_coef(double *f, int l, int m, double a, double b) {
  double Ca[32], Cb[32];
  int i, j, k;
  Ca[l] = 1.0;
  for (j = l - 1; j >= 0; j--)
    Ca[j] = a * (j + 1) / (l - j) * Ca[j + 1];
  Cb[m] = 1.0;
  for (k = m - 1; k >= 0; k--)
    Cb[k] = b * (k + 1) / (m - k) * Cb[k + 1];
  for (i = 0; i <= (l + m); i++)
    f[i] = 0.0;
  for (j = 0; j < l + 1; j++)
    for (k = 0; k < m + 1; k++) {
      i = j + k;
      f[i] += Ca[j] * Cb[k];
    }
  return;
}


x86 assembly (clang, -O0):
<f_coef>:
push   %rbp
mov    %rsp,%rbp
sub    $0x1b0,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
movsd  %xmm0,-0x18(%rbp)
movsd  %xmm1,-0x20(%rbp)
movslq -0xc(%rbp),%rax
movsd  0xed5(%rip),%xmm0
movsd  %xmm0,-0x120(%rbp,%rax,8)
ff ff
mov    -0xc(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x228(%rbp)
cmpl   $0x0,-0x228(%rbp)
jl     11ac <f_coef+0xac>
movsd  -0x18(%rbp),%xmm0
mov    -0x228(%rbp),%eax
add    $0x1,%eax
cvtsi2sd %eax,%xmm1
mulsd  %xmm1,%xmm0
mov    -0xc(%rbp),%eax
sub    -0x228(%rbp),%eax
cvtsi2sd %eax,%xmm1
divsd  %xmm1,%xmm0
mov    -0x228(%rbp),%eax
add    $0x1,%eax
cltq
mulsd  -0x120(%rbp,%rax,8),%xmm0
ff ff
movslq -0x228(%rbp),%rax
movsd  %xmm0,-0x120(%rbp,%rax,8)
ff ff
mov    -0x228(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x228(%rbp)
jmp    1140 <f_coef+0x40>
movslq -0x10(%rbp),%rax
movsd  0xe48(%rip),%xmm0
movsd  %xmm0,-0x220(%rbp,%rax,8)
ff ff
mov    -0x10(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x22c(%rbp)
cmpl   $0x0,-0x22c(%rbp)
jl     1239 <f_coef+0x139>
movsd  -0x20(%rbp),%xmm0
mov    -0x22c(%rbp),%eax
add    $0x1,%eax
cvtsi2sd %eax,%xmm1
mulsd  %xmm1,%xmm0
mov    -0x10(%rbp),%eax
sub    -0x22c(%rbp),%eax
cvtsi2sd %eax,%xmm1
divsd  %xmm1,%xmm0
mov    -0x22c(%rbp),%eax
add    $0x1,%eax
cltq
mulsd  -0x220(%rbp,%rax,8),%xmm0
ff ff
movslq -0x22c(%rbp),%rax
movsd  %xmm0,-0x220(%rbp,%rax,8)
ff ff
mov    -0x22c(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x22c(%rbp)
jmp    11cd <f_coef+0xcd>
movl   $0x0,-0x224(%rbp)
mov    -0x224(%rbp),%eax
mov    -0xc(%rbp),%ecx
add    -0x10(%rbp),%ecx
cmp    %ecx,%eax
jg     127e <f_coef+0x17e>
mov    -0x8(%rbp),%rax
movslq -0x224(%rbp),%rcx
xorps  %xmm0,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x224(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x224(%rbp)
jmp    1243 <f_coef+0x143>
movl   $0x0,-0x228(%rbp)
mov    -0x228(%rbp),%eax
mov    -0xc(%rbp),%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jge    1336 <f_coef+0x236>
movl   $0x0,-0x22c(%rbp)
mov    -0x22c(%rbp),%eax
mov    -0x10(%rbp),%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jge    131d <f_coef+0x21d>
mov    -0x228(%rbp),%eax
add    -0x22c(%rbp),%eax
mov    %eax,-0x224(%rbp)
movslq -0x228(%rbp),%rax
movsd  -0x120(%rbp,%rax,8),%xmm0
ff ff
movslq -0x22c(%rbp),%rax
movsd  -0x220(%rbp,%rax,8),%xmm2
ff ff
mov    -0x8(%rbp),%rax
movslq -0x224(%rbp),%rcx
movsd  (%rax,%rcx,8),%xmm1
mulsd  %xmm2,%xmm0
addsd  %xmm1,%xmm0
movsd  %xmm0,(%rax,%rcx,8)
mov    -0x22c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x22c(%rbp)
jmp    12a6 <f_coef+0x1a6>
jmp    1322 <f_coef+0x222>
mov    -0x228(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x228(%rbp)
jmp    1288 <f_coef+0x188>
add    $0x1b0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<f_coef>:
push   %r15
push   %r14
push   %r12
push   %rbx
sub    $0x208,%rsp
mov    %edx,%r14d
mov    %esi,%r12d
mov    %rdi,%r15
movslq %esi,%rax
movabs $0x3ff0000000000000,%r8
mov    %r8,0x100(%rsp,%rax,8)
test   %eax,%eax
jle    119a <f_coef+0x8a>
mov    %r12d,%edi
cmp    $0x1,%r12d
jne    139a <f_coef+0x28a>
mov    %rdi,%rax
mov    %r12d,%ecx
sub    %eax,%ecx
add    $0x1,%ecx
nopl   0x0(%rax,%rax,1)
xorps  %xmm2,%xmm2
cvtsi2sd %eax,%xmm2
mulsd  %xmm0,%xmm2
xorps  %xmm3,%xmm3
cvtsi2sd %ecx,%xmm3
divsd  %xmm3,%xmm2
mulsd  0x100(%rsp,%rax,8),%xmm2
lea    -0x1(%rax),%rdx
mov    %edx,%esi
movsd  %xmm2,0x100(%rsp,%rsi,8)
add    $0x1,%ecx
cmp    $0x1,%rax
mov    %rdx,%rax
ja     1160 <f_coef+0x50>
movslq %r14d,%rax
mov    %r8,(%rsp,%rax,8)
test   %eax,%eax
jle    11f2 <f_coef+0xe2>
mov    %r14d,%esi
cmp    $0x1,%r14d
jne    1503 <f_coef+0x3f3>
mov    %rsi,%rax
mov    %r14d,%ecx
sub    %eax,%ecx
add    $0x1,%ecx
nopl   (%rax)
xorps  %xmm0,%xmm0
cvtsi2sd %eax,%xmm0
mulsd  %xmm1,%xmm0
xorps  %xmm2,%xmm2
cvtsi2sd %ecx,%xmm2
divsd  %xmm2,%xmm0
mulsd  (%rsp,%rax,8),%xmm0
lea    -0x1(%rax),%rdx
mov    %edx,%esi
movsd  %xmm0,(%rsp,%rsi,8)
add    $0x1,%ecx
cmp    $0x1,%rax
mov    %rdx,%rax
ja     11c0 <f_coef+0xb0>
mov    %r14d,%eax
add    %r12d,%eax
js     120e <f_coef+0xfe>
mov    %eax,%eax
lea    0x8(,%rax,8),%rdx
mov    %r15,%rdi
xor    %esi,%esi
call   1030 <memset@plt>
test   %r12d,%r12d
js     138b <f_coef+0x27b>
lea    0x1(%r14),%eax
add    $0x1,%r12d
mov    %eax,%r11d
and    $0xfffffffc,%r11d
lea    -0x4(%r11),%r9
mov    %r9,%r10
shr    $0x2,%r10
add    $0x1,%r10
mov    %r10,%r8
and    $0xfffffffffffffffe,%r8
lea    0x30(%r15),%rdi
xor    %edx,%edx
mov    %r15,%rbx
jmp    1265 <f_coef+0x155>
nopw   0x0(%rax,%rax,1)
add    $0x1,%rdx
add    $0x8,%rdi
add    $0x8,%rbx
cmp    %r12,%rdx
je     138b <f_coef+0x27b>
test   %r14d,%r14d
js     1250 <f_coef+0x140>
movsd  0x100(%rsp,%rdx,8),%xmm0
cmp    $0x4,%eax
jae    1280 <f_coef+0x170>
xor    %ecx,%ecx
jmp    1360 <f_coef+0x250>
nop
movapd %xmm0,%xmm1
unpcklpd %xmm0,%xmm1
test   %r9,%r9
je     1381 <f_coef+0x271>
mov    %r8,%rcx
xor    %esi,%esi
cs nopw 0x0(%rax,%rax,1)
movapd (%rsp,%rsi,8),%xmm2
mulpd  %xmm1,%xmm2
movupd -0x30(%rdi,%rsi,8),%xmm3
addpd  %xmm2,%xmm3
movupd -0x20(%rdi,%rsi,8),%xmm2
movupd -0x10(%rdi,%rsi,8),%xmm4
movupd (%rdi,%rsi,8),%xmm5
movapd 0x10(%rsp,%rsi,8),%xmm6
mulpd  %xmm1,%xmm6
addpd  %xmm2,%xmm6
movupd %xmm3,-0x30(%rdi,%rsi,8)
movupd %xmm6,-0x20(%rdi,%rsi,8)
movapd 0x20(%rsp,%rsi,8),%xmm2
mulpd  %xmm1,%xmm2
addpd  %xmm4,%xmm2
movapd 0x30(%rsp,%rsi,8),%xmm3
mulpd  %xmm1,%xmm3
addpd  %xmm5,%xmm3
movupd %xmm2,-0x10(%rdi,%rsi,8)
movupd %xmm3,(%rdi,%rsi,8)
add    $0x8,%rsi
add    $0xfffffffffffffffe,%rcx
jne    12a0 <f_coef+0x190>
test   $0x1,%r10b
je     134a <f_coef+0x23a>
lea    (%rsi,%rdx,1),%rcx
movupd (%r15,%rcx,8),%xmm2
movupd 0x10(%r15,%rcx,8),%xmm3
movapd (%rsp,%rsi,8),%xmm4
mulpd  %xmm1,%xmm4
mulpd  0x10(%rsp,%rsi,8),%xmm1
addpd  %xmm2,%xmm4
addpd  %xmm3,%xmm1
movupd %xmm4,(%r15,%rcx,8)
movupd %xmm1,0x10(%r15,%rcx,8)
mov    %r11,%rcx
cmp    %rax,%r11
je     1250 <f_coef+0x140>
cs nopw 0x0(%rax,%rax,1)
movsd  (%rsp,%rcx,8),%xmm1
mulsd  %xmm0,%xmm1
addsd  (%rbx,%rcx,8),%xmm1
movsd  %xmm1,(%rbx,%rcx,8)
add    $0x1,%rcx
cmp    %rcx,%rax
jne    1360 <f_coef+0x250>
jmp    1250 <f_coef+0x140>
xor    %esi,%esi
test   $0x1,%r10b
jne    1315 <f_coef+0x205>
jmp    134a <f_coef+0x23a>
add    $0x208,%rsp
pop    %rbx
pop    %r12
pop    %r14
pop    %r15
ret
lea    -0x1(%rdi),%rcx
lea    -0x1(%r12),%ebx
mov    $0x8,%edx
mov    %rcx,%rax
mul    %rdx
seto   %dl
cmp    %ecx,%ebx
jb     114d <f_coef+0x3d>
shr    $0x20,%rcx
jne    114d <f_coef+0x3d>
lea    (%rsp,%rbx,8),%rsi
add    $0x100,%rsi
mov    %rsi,%rcx
sub    %rax,%rcx
cmp    %rsi,%rcx
ja     114d <f_coef+0x3d>
test   %dl,%dl
jne    114d <f_coef+0x3d>
lea    (%rsp,%rdi,8),%rcx
add    $0x100,%rcx
mov    %rcx,%rsi
sub    %rax,%rsi
cmp    %rcx,%rsi
ja     114d <f_coef+0x3d>
test   %dl,%dl
jne    114d <f_coef+0x3d>
lea    0x1(%rbx),%rax
sub    %rdi,%rax
lea    (%rsp,%rax,8),%rax
add    $0x100,%rax
lea    (%rsp,%rdi,8),%rcx
add    $0x108,%rcx
cmp    %rcx,%rax
jae    1443 <f_coef+0x333>
lea    (%rsp,%rbx,8),%rax
add    $0x108,%rax
lea    0x108(%rsp),%rcx
cmp    %rax,%rcx
jb     114d <f_coef+0x3d>
mov    %edi,%ebx
and    $0xfffffffe,%ebx
mov    %edi,%eax
and    $0x1,%eax
movd   %r12d,%xmm2
pshufd $0x50,%xmm2,%xmm6
movdqa 0xba1(%rip),%xmm3
paddd  %xmm6,%xmm3
movapd %xmm0,%xmm4
unpcklpd %xmm0,%xmm4
movdqa 0xb9c(%rip),%xmm8
movdqa 0xba3(%rip),%xmm9
mov    %rdi,%rdx
mov    %rbx,%rcx
cs nopw 0x0(%rax,%rax,1)
nopl   (%rax)
movq   %rdx,%xmm7
pshufd $0x44,%xmm7,%xmm7
paddq  %xmm8,%xmm7
pshufd $0xfa,%xmm7,%xmm5
movq   %xmm7,%rsi
movd   %esi,%xmm7
movss  %xmm7,%xmm5
cvtdq2pd %xmm3,%xmm7
mulpd  %xmm4,%xmm7
movdqa %xmm6,%xmm2
psubd  %xmm5,%xmm2
cvtdq2pd %xmm2,%xmm2
divpd  %xmm2,%xmm7
movupd 0xf8(%rsp,%rdx,8),%xmm2
mov    %esi,%esi
shufpd $0x1,%xmm7,%xmm7
mulpd  %xmm2,%xmm7
movupd %xmm7,0xf8(%rsp,%rsi,8)
paddd  %xmm9,%xmm3
add    $0xfffffffffffffffe,%rdx
add    $0xfffffffffffffffe,%rcx
jne    1490 <f_coef+0x380>
cmp    %rdi,%rbx
jne    1150 <f_coef+0x40>
jmp    119a <f_coef+0x8a>
lea    -0x1(%rsi),%rcx
lea    -0x1(%r14),%edi
mov    $0x8,%edx
mov    %rcx,%rax
mul    %rdx
seto   %dl
cmp    %ecx,%edi
jb     11b2 <f_coef+0xa2>
shr    $0x20,%rcx
jne    11b2 <f_coef+0xa2>
lea    (%rsp,%rdi,8),%rbx
mov    %rbx,%rcx
sub    %rax,%rcx
cmp    %rbx,%rcx
ja     11b2 <f_coef+0xa2>
test   %dl,%dl
jne    11b2 <f_coef+0xa2>
lea    (%rsp,%rsi,8),%rbx
mov    %rbx,%rcx
sub    %rax,%rcx
cmp    %rbx,%rcx
ja     11b2 <f_coef+0xa2>
test   %dl,%dl
jne    11b2 <f_coef+0xa2>
lea    0x1(%rdi),%rax
sub    %rsi,%rax
lea    (%rsp,%rax,8),%rax
lea    (%rsp,%rsi,8),%rcx
add    $0x8,%rcx
cmp    %rcx,%rax
jae    158f <f_coef+0x47f>
lea    (%rsp,%rdi,8),%rax
add    $0x8,%rax
lea    0x8(%rsp),%rcx
cmp    %rax,%rcx
jb     11b2 <f_coef+0xa2>
mov    %esi,%ecx
and    $0xfffffffe,%ecx
mov    %esi,%eax
and    $0x1,%eax
movd   %r14d,%xmm0
pshufd $0x50,%xmm0,%xmm4
movdqa 0xa55(%rip),%xmm2
paddd  %xmm4,%xmm2
movapd %xmm1,%xmm3
unpcklpd %xmm1,%xmm3
movdqa 0xa50(%rip),%xmm8
movdqa 0xa58(%rip),%xmm5
mov    %rsi,%rdx
mov    %rcx,%rdi
xchg   %ax,%ax
movq   %rdx,%xmm6
pshufd $0x44,%xmm6,%xmm6
paddq  %xmm8,%xmm6
pshufd $0xfa,%xmm6,%xmm7
movq   %xmm6,%rbx
movd   %ebx,%xmm6
movss  %xmm6,%xmm7
cvtdq2pd %xmm2,%xmm6
mulpd  %xmm3,%xmm6
movdqa %xmm4,%xmm0
psubd  %xmm7,%xmm0
cvtdq2pd %xmm0,%xmm0
divpd  %xmm0,%xmm6
movupd -0x8(%rsp,%rdx,8),%xmm0
mov    %ebx,%ebx
shufpd $0x1,%xmm6,%xmm6
mulpd  %xmm0,%xmm6
movupd %xmm6,-0x8(%rsp,%rbx,8)
paddd  %xmm5,%xmm2
add    $0xfffffffffffffffe,%rdx
add    $0xfffffffffffffffe,%rdi
jne    15d0 <f_coef+0x4c0>
cmp    %rsi,%rcx
jne    11b5 <f_coef+0xa5>
jmp    11f2 <f_coef+0xe2>

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

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