void f_coef(float *f, int l, int m, float a, float b) {
  float 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    $0xb0,%rsp
mov    %rdi,-0x8(%rbp)
mov    %esi,-0xc(%rbp)
mov    %edx,-0x10(%rbp)
movss  %xmm0,-0x14(%rbp)
movss  %xmm1,-0x18(%rbp)
movslq -0xc(%rbp),%rax
movss  0xed5(%rip),%xmm0
movss  %xmm0,-0xa0(%rbp,%rax,4)
ff ff
mov    -0xc(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x128(%rbp)
cmpl   $0x0,-0x128(%rbp)
jl     11ac <f_coef+0xac>
movss  -0x14(%rbp),%xmm0
mov    -0x128(%rbp),%eax
add    $0x1,%eax
cvtsi2ss %eax,%xmm1
mulss  %xmm1,%xmm0
mov    -0xc(%rbp),%eax
sub    -0x128(%rbp),%eax
cvtsi2ss %eax,%xmm1
divss  %xmm1,%xmm0
mov    -0x128(%rbp),%eax
add    $0x1,%eax
cltq
mulss  -0xa0(%rbp,%rax,4),%xmm0
ff ff
movslq -0x128(%rbp),%rax
movss  %xmm0,-0xa0(%rbp,%rax,4)
ff ff
mov    -0x128(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x128(%rbp)
jmp    1140 <f_coef+0x40>
movslq -0x10(%rbp),%rax
movss  0xe48(%rip),%xmm0
movss  %xmm0,-0x120(%rbp,%rax,4)
ff ff
mov    -0x10(%rbp),%eax
sub    $0x1,%eax
mov    %eax,-0x12c(%rbp)
cmpl   $0x0,-0x12c(%rbp)
jl     1239 <f_coef+0x139>
movss  -0x18(%rbp),%xmm0
mov    -0x12c(%rbp),%eax
add    $0x1,%eax
cvtsi2ss %eax,%xmm1
mulss  %xmm1,%xmm0
mov    -0x10(%rbp),%eax
sub    -0x12c(%rbp),%eax
cvtsi2ss %eax,%xmm1
divss  %xmm1,%xmm0
mov    -0x12c(%rbp),%eax
add    $0x1,%eax
cltq
mulss  -0x120(%rbp,%rax,4),%xmm0
ff ff
movslq -0x12c(%rbp),%rax
movss  %xmm0,-0x120(%rbp,%rax,4)
ff ff
mov    -0x12c(%rbp),%eax
add    $0xffffffff,%eax
mov    %eax,-0x12c(%rbp)
jmp    11cd <f_coef+0xcd>
movl   $0x0,-0x124(%rbp)
mov    -0x124(%rbp),%eax
mov    -0xc(%rbp),%ecx
add    -0x10(%rbp),%ecx
cmp    %ecx,%eax
jg     127e <f_coef+0x17e>
mov    -0x8(%rbp),%rax
movslq -0x124(%rbp),%rcx
xorps  %xmm0,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x124(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x124(%rbp)
jmp    1243 <f_coef+0x143>
movl   $0x0,-0x128(%rbp)
mov    -0x128(%rbp),%eax
mov    -0xc(%rbp),%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jge    1336 <f_coef+0x236>
movl   $0x0,-0x12c(%rbp)
mov    -0x12c(%rbp),%eax
mov    -0x10(%rbp),%ecx
add    $0x1,%ecx
cmp    %ecx,%eax
jge    131d <f_coef+0x21d>
mov    -0x128(%rbp),%eax
add    -0x12c(%rbp),%eax
mov    %eax,-0x124(%rbp)
movslq -0x128(%rbp),%rax
movss  -0xa0(%rbp,%rax,4),%xmm0
ff ff
movslq -0x12c(%rbp),%rax
movss  -0x120(%rbp,%rax,4),%xmm2
ff ff
mov    -0x8(%rbp),%rax
movslq -0x124(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,(%rax,%rcx,4)
mov    -0x12c(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x12c(%rbp)
jmp    12a6 <f_coef+0x1a6>
jmp    1322 <f_coef+0x222>
mov    -0x128(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x128(%rbp)
jmp    1288 <f_coef+0x188>
add    $0xb0,%rsp
pop    %rbp
ret

x86 assembly (clang, -O2):
<f_coef>:
push   %r15
push   %r14
push   %r12
push   %rbx
sub    $0x108,%rsp
mov    %edx,%r14d
mov    %esi,%r12d
mov    %rdi,%r15
movslq %esi,%rax
movl   $0x3f800000,0x80(%rsp,%rax,4)
test   %eax,%eax
jle    119a <f_coef+0x8a>
mov    %r12d,%esi
cmp    $0x4,%r12d
jae    137a <f_coef+0x26a>
mov    %rsi,%rax
mov    %r12d,%ecx
sub    %eax,%ecx
add    $0x1,%ecx
cs nopw 0x0(%rax,%rax,1)
nopl   0x0(%rax,%rax,1)
xorps  %xmm2,%xmm2
cvtsi2ss %eax,%xmm2
mulss  %xmm0,%xmm2
xorps  %xmm3,%xmm3
cvtsi2ss %ecx,%xmm3
divss  %xmm3,%xmm2
mulss  0x80(%rsp,%rax,4),%xmm2
lea    -0x1(%rax),%rdx
mov    %edx,%esi
movss  %xmm2,0x80(%rsp,%rsi,4)
add    $0x1,%ecx
cmp    $0x1,%rax
mov    %rdx,%rax
ja     1160 <f_coef+0x50>
movslq %r14d,%rax
movl   $0x3f800000,(%rsp,%rax,4)
test   %eax,%eax
jle    11f2 <f_coef+0xe2>
mov    %r14d,%esi
cmp    $0x4,%r14d
jae    14dd <f_coef+0x3cd>
mov    %rsi,%rax
mov    %r14d,%ecx
sub    %eax,%ecx
add    $0x1,%ecx
xorps  %xmm0,%xmm0
cvtsi2ss %eax,%xmm0
mulss  %xmm1,%xmm0
xorps  %xmm2,%xmm2
cvtsi2ss %ecx,%xmm2
divss  %xmm2,%xmm0
mulss  (%rsp,%rax,4),%xmm0
lea    -0x1(%rax),%rdx
mov    %edx,%esi
movss  %xmm0,(%rsp,%rsi,4)
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    0x4(,%rax,4),%rdx
mov    %r15,%rdi
xor    %esi,%esi
call   1030 <memset@plt>
test   %r12d,%r12d
js     136b <f_coef+0x25b>
lea    0x1(%r14),%eax
add    $0x1,%r12d
mov    %eax,%r11d
and    $0xfffffff8,%r11d
lea    -0x8(%r11),%r9
mov    %r9,%r10
shr    $0x3,%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    $0x4,%rdi
add    $0x4,%rbx
cmp    %r12,%rdx
je     136b <f_coef+0x25b>
test   %r14d,%r14d
js     1250 <f_coef+0x140>
movss  0x80(%rsp,%rdx,4),%xmm0
cmp    $0x8,%eax
jae    1280 <f_coef+0x170>
xor    %ecx,%ecx
jmp    1340 <f_coef+0x230>
nop
movaps %xmm0,%xmm1
shufps $0x0,%xmm0,%xmm1
test   %r9,%r9
je     1361 <f_coef+0x251>
mov    %r8,%rcx
xor    %esi,%esi
cs nopw 0x0(%rax,%rax,1)
nop
movaps (%rsp,%rsi,4),%xmm2
mulps  %xmm1,%xmm2
movups -0x30(%rdi,%rsi,4),%xmm3
addps  %xmm2,%xmm3
movups -0x20(%rdi,%rsi,4),%xmm2
movups -0x10(%rdi,%rsi,4),%xmm4
movups (%rdi,%rsi,4),%xmm5
movaps 0x10(%rsp,%rsi,4),%xmm6
mulps  %xmm1,%xmm6
addps  %xmm2,%xmm6
movups %xmm3,-0x30(%rdi,%rsi,4)
movups %xmm6,-0x20(%rdi,%rsi,4)
movaps 0x20(%rsp,%rsi,4),%xmm2
mulps  %xmm1,%xmm2
addps  %xmm4,%xmm2
movaps 0x30(%rsp,%rsi,4),%xmm3
mulps  %xmm1,%xmm3
addps  %xmm5,%xmm3
movups %xmm2,-0x10(%rdi,%rsi,4)
movups %xmm3,(%rdi,%rsi,4)
add    $0x10,%rsi
add    $0xfffffffffffffffe,%rcx
jne    12a0 <f_coef+0x190>
test   $0x1,%r10b
je     1330 <f_coef+0x220>
lea    (%rsi,%rdx,1),%rcx
movups (%r15,%rcx,4),%xmm2
movups 0x10(%r15,%rcx,4),%xmm3
movaps %xmm1,%xmm4
mulps  (%rsp,%rsi,4),%xmm4
addps  %xmm2,%xmm4
movaps 0x10(%rsp,%rsi,4),%xmm2
mulps  %xmm1,%xmm2
addps  %xmm3,%xmm2
movups %xmm4,(%r15,%rcx,4)
movups %xmm2,0x10(%r15,%rcx,4)
mov    %r11,%rcx
cmp    %rax,%r11
je     1250 <f_coef+0x140>
nopl   0x0(%rax)
movss  (%rsp,%rcx,4),%xmm1
mulss  %xmm0,%xmm1
addss  (%rbx,%rcx,4),%xmm1
movss  %xmm1,(%rbx,%rcx,4)
add    $0x1,%rcx
cmp    %rcx,%rax
jne    1340 <f_coef+0x230>
jmp    1250 <f_coef+0x140>
xor    %esi,%esi
test   $0x1,%r10b
jne    1301 <f_coef+0x1f1>
jmp    1330 <f_coef+0x220>
add    $0x108,%rsp
pop    %rbx
pop    %r12
pop    %r14
pop    %r15
ret
lea    -0x1(%rsi),%rcx
lea    -0x1(%r12),%edi
mov    $0x4,%edx
mov    %rcx,%rax
mul    %rdx
seto   %dl
cmp    %ecx,%edi
jb     1146 <f_coef+0x36>
shr    $0x20,%rcx
jne    1146 <f_coef+0x36>
lea    (%rsp,%rdi,4),%rbx
add    $0x80,%rbx
mov    %rbx,%rcx
sub    %rax,%rcx
cmp    %rbx,%rcx
ja     1146 <f_coef+0x36>
test   %dl,%dl
jne    1146 <f_coef+0x36>
lea    (%rsp,%rsi,4),%rbx
add    $0x80,%rbx
mov    %rbx,%rcx
sub    %rax,%rcx
cmp    %rbx,%rcx
ja     1146 <f_coef+0x36>
test   %dl,%dl
jne    1146 <f_coef+0x36>
lea    0x1(%rdi),%rax
sub    %rsi,%rax
lea    (%rsp,%rax,4),%rax
add    $0x80,%rax
lea    (%rsp,%rsi,4),%rcx
add    $0x84,%rcx
cmp    %rcx,%rax
jae    1423 <f_coef+0x313>
lea    (%rsp,%rdi,4),%rax
add    $0x84,%rax
lea    0x84(%rsp),%rcx
cmp    %rax,%rcx
jb     1146 <f_coef+0x36>
mov    %esi,%ecx
and    $0xfffffffc,%ecx
mov    %esi,%eax
and    $0x3,%eax
movd   %r12d,%xmm2
pshufd $0x0,%xmm2,%xmm7
movdqa 0xbc1(%rip),%xmm3
paddd  %xmm7,%xmm3
movaps %xmm0,%xmm4
shufps $0x0,%xmm0,%xmm4
movdqa 0xbbd(%rip),%xmm8
movdqa 0xbc4(%rip),%xmm9
movdqa 0xbcb(%rip),%xmm10
mov    %rsi,%rdx
mov    %rcx,%rdi
nopl   0x0(%rax,%rax,1)
movq   %rdx,%xmm5
pshufd $0x44,%xmm5,%xmm5
movdqa %xmm5,%xmm6
paddq  %xmm8,%xmm6
paddq  %xmm9,%xmm5
movq   %xmm5,%rbx
shufps $0x88,%xmm6,%xmm5
movd   %ebx,%xmm6
movss  %xmm6,%xmm5
cvtdq2ps %xmm3,%xmm6
mulps  %xmm4,%xmm6
movdqa %xmm7,%xmm2
psubd  %xmm5,%xmm2
cvtdq2ps %xmm2,%xmm2
divps  %xmm2,%xmm6
movups 0x74(%rsp,%rdx,4),%xmm2
mov    %ebx,%ebx
shufps $0x1b,%xmm6,%xmm6
mulps  %xmm2,%xmm6
movups %xmm6,0x74(%rsp,%rbx,4)
paddd  %xmm10,%xmm3
add    $0xfffffffffffffffc,%rdx
add    $0xfffffffffffffffc,%rdi
jne    1470 <f_coef+0x360>
cmp    %rsi,%rcx
jne    1149 <f_coef+0x39>
jmp    119a <f_coef+0x8a>
lea    -0x1(%rsi),%rcx
lea    -0x1(%r14),%edi
mov    $0x4,%edx
mov    %rcx,%rax
mul    %rdx
seto   %dl
cmp    %ecx,%edi
jb     11b5 <f_coef+0xa5>
shr    $0x20,%rcx
jne    11b5 <f_coef+0xa5>
lea    (%rsp,%rdi,4),%rbx
mov    %rbx,%rcx
sub    %rax,%rcx
cmp    %rbx,%rcx
ja     11b5 <f_coef+0xa5>
test   %dl,%dl
jne    11b5 <f_coef+0xa5>
lea    (%rsp,%rsi,4),%rbx
mov    %rbx,%rcx
sub    %rax,%rcx
cmp    %rbx,%rcx
ja     11b5 <f_coef+0xa5>
test   %dl,%dl
jne    11b5 <f_coef+0xa5>
lea    0x1(%rdi),%rax
sub    %rsi,%rax
lea    (%rsp,%rax,4),%rax
lea    (%rsp,%rsi,4),%rcx
add    $0x4,%rcx
cmp    %rcx,%rax
jae    1569 <f_coef+0x459>
lea    (%rsp,%rdi,4),%rax
add    $0x4,%rax
lea    0x4(%rsp),%rcx
cmp    %rax,%rcx
jb     11b5 <f_coef+0xa5>
mov    %esi,%ecx
and    $0xfffffffc,%ecx
mov    %esi,%eax
and    $0x3,%eax
movd   %r14d,%xmm0
pshufd $0x0,%xmm0,%xmm5
movdqa 0xa7b(%rip),%xmm2
paddd  %xmm5,%xmm2
movaps %xmm1,%xmm3
shufps $0x0,%xmm1,%xmm3
movdqa 0xa77(%rip),%xmm8
movdqa 0xa7e(%rip),%xmm9
movdqa 0xa86(%rip),%xmm6
mov    %rsi,%rdx
mov    %rcx,%rdi
movq   %rdx,%xmm7
pshufd $0x44,%xmm7,%xmm7
movdqa %xmm7,%xmm4
paddq  %xmm8,%xmm4
paddq  %xmm9,%xmm7
movq   %xmm7,%rbx
shufps $0x88,%xmm4,%xmm7
movd   %ebx,%xmm4
movss  %xmm4,%xmm7
cvtdq2ps %xmm2,%xmm4
mulps  %xmm3,%xmm4
movdqa %xmm5,%xmm0
psubd  %xmm7,%xmm0
cvtdq2ps %xmm0,%xmm0
divps  %xmm0,%xmm4
movups -0xc(%rsp,%rdx,4),%xmm0
mov    %ebx,%ebx
shufps $0x1b,%xmm4,%xmm4
mulps  %xmm0,%xmm4
movups %xmm4,-0xc(%rsp,%rbx,4)
paddd  %xmm6,%xmm2
add    $0xfffffffffffffffc,%rdx
add    $0xfffffffffffffffc,%rdi
jne    15b0 <f_coef+0x4a0>
cmp    %rsi,%rcx
jne    11b8 <f_coef+0xa8>
jmp    11f2 <f_coef+0xe2>

LLVM IR (-O0):
define void @f_coef(float* %0, i32 %1, i32 %2, float %3, float %4) {
%6 = alloca float*
%7 = alloca i32
%8 = alloca i32
%9 = alloca float
%10 = alloca float
%11 = alloca [32 x float]
%12 = alloca [32 x float]
%13 = alloca i32
%14 = alloca i32
%15 = alloca i32
store float* %0, float** %6
store i32 %1, i32* %7
store i32 %2, i32* %8
store float %3, float* %9
store float %4, float* %10
%16 = load i32, i32* %7
%17 = sext i32 %16 to i64
%18 = getelementptr inbounds [32 x float], [32 x float]* %11, i64 0, i64 %17
store float 1.000000e+00, float* %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 float, float* %9
%26 = load i32, i32* %14
%27 = add nsw i32 %26, 1
%28 = sitofp i32 %27 to float
%29 = fmul float %25, %28
%30 = load i32, i32* %7
%31 = load i32, i32* %14
%32 = sub nsw i32 %30, %31
%33 = sitofp i32 %32 to float
%34 = fdiv float %29, %33
%35 = load i32, i32* %14
%36 = add nsw i32 %35, 1
%37 = sext i32 %36 to i64
%38 = getelementptr inbounds [32 x float], [32 x float]* %11, i64 0, i64 %37
%39 = load float, float* %38
%40 = fmul float %34, %39
%41 = load i32, i32* %14
%42 = sext i32 %41 to i64
%43 = getelementptr inbounds [32 x float], [32 x float]* %11, i64 0, i64 %42
store float %40, float* %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 float], [32 x float]* %12, i64 0, i64 %49
store float 1.000000e+00, float* %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 float, float* %10
%58 = load i32, i32* %15
%59 = add nsw i32 %58, 1
%60 = sitofp i32 %59 to float
%61 = fmul float %57, %60
%62 = load i32, i32* %8
%63 = load i32, i32* %15
%64 = sub nsw i32 %62, %63
%65 = sitofp i32 %64 to float
%66 = fdiv float %61, %65
%67 = load i32, i32* %15
%68 = add nsw i32 %67, 1
%69 = sext i32 %68 to i64
%70 = getelementptr inbounds [32 x float], [32 x float]* %12, i64 0, i64 %69
%71 = load float, float* %70
%72 = fmul float %66, %71
%73 = load i32, i32* %15
%74 = sext i32 %73 to i64
%75 = getelementptr inbounds [32 x float], [32 x float]* %12, i64 0, i64 %74
store float %72, float* %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 float*, float** %6
%88 = load i32, i32* %13
%89 = sext i32 %88 to i64
%90 = getelementptr inbounds float, float* %87, i64 %89
store float 0.000000e+00, float* %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 float], [32 x float]* %11, i64 0, i64 %111
%113 = load float, float* %112
%114 = load i32, i32* %15
%115 = sext i32 %114 to i64
%116 = getelementptr inbounds [32 x float], [32 x float]* %12, i64 0, i64 %115
%117 = load float, float* %116
%118 = load float*, float** %6
%119 = load i32, i32* %13
%120 = sext i32 %119 to i64
%121 = getelementptr inbounds float, float* %118, i64 %120
%122 = load float, float* %121
%123 = call float @llvm.fmuladd.f32(float %113, float %117, float %122)
store float %123, float* %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(float* nocapture %0, i32 %1, i32 %2, float %3, float %4) {
%6 = bitcast float* %0 to i8*
%7 = alloca [32 x float]
%8 = alloca [32 x float]
%9 = bitcast [32 x float]* %7 to i8*
%10 = bitcast [32 x float]* %8 to i8*
%11 = sext i32 %1 to i64
%12 = getelementptr inbounds [32 x float], [32 x float]* %7, i64 0, i64 %11
store float 1.000000e+00, float* %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 ult i32 %1, 4
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 float], [32 x float]* %7, i64 0, i64 %24
%26 = call { i64, i1 } @llvm.umul.with.overflow.i64(i64 %18, i64 4)
%27 = extractvalue { i64, i1 } %26, 0
%28 = extractvalue { i64, i1 } %26, 1
%29 = bitcast float* %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 float], [32 x float]* %7, i64 0, i64 %15
%35 = call { i64, i1 } @llvm.umul.with.overflow.i64(i64 %18, i64 4)
%36 = extractvalue { i64, i1 } %35, 0
%37 = extractvalue { i64, i1 } %35, 1
%38 = bitcast float* %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 float], [32 x float]* %7, i64 0, i64 %49
%51 = add nuw nsw i64 %47, 1
%52 = getelementptr [32 x float], [32 x float]* %7, i64 0, i64 %51
%53 = getelementptr inbounds [32 x float], [32 x float]* %7, i64 0, i64 1
%54 = add nuw nsw i64 %15, 1
%55 = getelementptr [32 x float], [32 x float]* %7, i64 0, i64 %54
%56 = icmp ult float* %50, %55
%57 = icmp ult float* %53, %52
%58 = and i1 %56, %57
br i1 %58, label %99, label %59
59: ; preds = %45
%60 = and i64 %15, 4294967292
%61 = and i64 %15, 3
%62 = insertelement <4 x i32> poison, i32 %1, i64 0
%63 = shufflevector <4 x i32> %62, <4 x i32> poison, <4 x i32> zeroinitializer
%64 = add <4 x i32> %63, <i32 0, i32 -1, i32 -2, i32 -3>
%65 = insertelement <4 x float> poison, float %3, i64 0
%66 = shufflevector <4 x float> %65, <4 x float> poison, <4 x i32> zeroinitializer
%67 = insertelement <4 x i32> poison, i32 %1, i64 0
%68 = shufflevector <4 x i32> %67, <4 x i32> poison, <4 x i32> zeroinitializer
br label %69
69: ; preds = %69, %59
%70 = phi i64 [ 0, %59 ], [ %94, %69 ]
%71 = phi <4 x i32> [ %64, %59 ], [ %95, %69 ]
%72 = sub i64 %15, %70
%73 = insertelement <4 x i64> poison, i64 %72, i64 0
%74 = shufflevector <4 x i64> %73, <4 x i64> poison, <4 x i32> zeroinitializer
%75 = add <4 x i64> %74, <i64 -1, i64 -2, i64 -3, i64 -4>
%76 = trunc <4 x i64> %75 to <4 x i32>
%77 = sitofp <4 x i32> %71 to <4 x float>
%78 = fmul <4 x float> %66, %77
%79 = sub nsw <4 x i32> %68, %76
%80 = sitofp <4 x i32> %79 to <4 x float>
%81 = fdiv <4 x float> %78, %80
%82 = getelementptr inbounds [32 x float], [32 x float]* %7, i64 0, i64 %72
%83 = getelementptr inbounds float, float* %82, i64 -3
%84 = bitcast float* %83 to <4 x float>*
%85 = load <4 x float>, <4 x float>* %84
%86 = shufflevector <4 x float> %85, <4 x float> poison, <4 x i32> <i32 3, i32 2, i32 1, i32 0>
%87 = fmul <4 x float> %81, %86
%88 = extractelement <4 x i64> %75, i64 0
%89 = and i64 %88, 4294967295
%90 = getelementptr inbounds [32 x float], [32 x float]* %7, i64 0, i64 %89
%91 = shufflevector <4 x float> %87, <4 x float> poison, <4 x i32> <i32 3, i32 2, i32 1, i32 0>
%92 = getelementptr inbounds float, float* %90, i64 -3
%93 = bitcast float* %92 to <4 x float>*
store <4 x float> %91, <4 x float>* %93
%94 = add nuw i64 %70, 4
%95 = add <4 x i32> %71, <i32 -4, i32 -4, i32 -4, i32 -4>
%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 float
%107 = fmul float %106, %3
%108 = sub nsw i32 %1, %104
%109 = sitofp i32 %108 to float
%110 = fdiv float %107, %109
%111 = getelementptr inbounds [32 x float], [32 x float]* %7, i64 0, i64 %102
%112 = load float, float* %111
%113 = fmul float %110, %112
%114 = and i64 %103, 4294967295
%115 = getelementptr inbounds [32 x float], [32 x float]* %7, i64 0, i64 %114
store float %113, float* %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 float], [32 x float]* %8, i64 0, i64 %118
store float 1.000000e+00, float* %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 ult i32 %2, 4
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 float], [32 x float]* %8, i64 0, i64 %131
%133 = call { i64, i1 } @llvm.umul.with.overflow.i64(i64 %125, i64 4)
%134 = extractvalue { i64, i1 } %133, 0
%135 = extractvalue { i64, i1 } %133, 1
%136 = bitcast float* %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 float], [32 x float]* %8, i64 0, i64 %122
%142 = call { i64, i1 } @llvm.umul.with.overflow.i64(i64 %125, i64 4)
%143 = extractvalue { i64, i1 } %142, 0
%144 = extractvalue { i64, i1 } %142, 1
%145 = bitcast float* %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 float], [32 x float]* %8, i64 0, i64 %156
%158 = add nuw nsw i64 %154, 1
%159 = getelementptr [32 x float], [32 x float]* %8, i64 0, i64 %158
%160 = getelementptr inbounds [32 x float], [32 x float]* %8, i64 0, i64 1
%161 = add nuw nsw i64 %122, 1
%162 = getelementptr [32 x float], [32 x float]* %8, i64 0, i64 %161
%163 = icmp ult float* %157, %162
%164 = icmp ult float* %160, %159
%165 = and i1 %163, %164
br i1 %165, label %206, label %166
166: ; preds = %152
%167 = and i64 %122, 4294967292
%168 = and i64 %122, 3
%169 = insertelement <4 x i32> poison, i32 %2, i64 0
%170 = shufflevector <4 x i32> %169, <4 x i32> poison, <4 x i32> zeroinitializer
%171 = add <4 x i32> %170, <i32 0, i32 -1, i32 -2, i32 -3>
%172 = insertelement <4 x float> poison, float %4, i64 0
%173 = shufflevector <4 x float> %172, <4 x float> poison, <4 x i32> zeroinitializer
%174 = insertelement <4 x i32> poison, i32 %2, i64 0
%175 = shufflevector <4 x i32> %174, <4 x i32> poison, <4 x i32> zeroinitializer
br label %176
176: ; preds = %176, %166
%177 = phi i64 [ 0, %166 ], [ %201, %176 ]
%178 = phi <4 x i32> [ %171, %166 ], [ %202, %176 ]
%179 = sub i64 %122, %177
%180 = insertelement <4 x i64> poison, i64 %179, i64 0
%181 = shufflevector <4 x i64> %180, <4 x i64> poison, <4 x i32> zeroinitializer
%182 = add <4 x i64> %181, <i64 -1, i64 -2, i64 -3, i64 -4>
%183 = trunc <4 x i64> %182 to <4 x i32>
%184 = sitofp <4 x i32> %178 to <4 x float>
%185 = fmul <4 x float> %173, %184
%186 = sub nsw <4 x i32> %175, %183
%187 = sitofp <4 x i32> %186 to <4 x float>
%188 = fdiv <4 x float> %185, %187
%189 = getelementptr inbounds [32 x float], [32 x float]* %8, i64 0, i64 %179
%190 = getelementptr inbounds float, float* %189, i64 -3
%191 = bitcast float* %190 to <4 x float>*
%192 = load <4 x float>, <4 x float>* %191
%193 = shufflevector <4 x float> %192, <4 x float> poison, <4 x i32> <i32 3, i32 2, i32 1, i32 0>
%194 = fmul <4 x float> %188, %193
%195 = extractelement <4 x i64> %182, i64 0
%196 = and i64 %195, 4294967295
%197 = getelementptr inbounds [32 x float], [32 x float]* %8, i64 0, i64 %196
%198 = shufflevector <4 x float> %194, <4 x float> poison, <4 x i32> <i32 3, i32 2, i32 1, i32 0>
%199 = getelementptr inbounds float, float* %197, i64 -3
%200 = bitcast float* %199 to <4 x float>*
store <4 x float> %198, <4 x float>* %200
%201 = add nuw i64 %177, 4
%202 = add <4 x i32> %178, <i32 -4, i32 -4, i32 -4, i32 -4>
%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, 2
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 float
%221 = fmul float %220, %4
%222 = sub nsw i32 %2, %218
%223 = sitofp i32 %222 to float
%224 = fdiv float %221, %223
%225 = getelementptr inbounds [32 x float], [32 x float]* %8, i64 0, i64 %216
%226 = load float, float* %225
%227 = fmul float %224, %226
%228 = and i64 %217, 4294967295
%229 = getelementptr inbounds [32 x float], [32 x float]* %8, i64 0, i64 %228
store float %227, float* %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, 4294967288
%240 = add nsw i64 %239, -8
%241 = lshr exact i64 %240, 3
%242 = add nuw nsw i64 %241, 1
%243 = icmp ult i32 %235, 8
%244 = and i64 %238, 4294967288
%245 = and i64 %242, 1
%246 = icmp eq i64 %240, 0
%247 = and i64 %242, 4611686018427387902
%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 float], [32 x float]* %7, i64 0, i64 %251
%254 = load float, float* %253
br i1 %243, label %322, label %255
255: ; preds = %252
%256 = insertelement <4 x float> poison, float %254, i64 0
%257 = shufflevector <4 x float> %256, <4 x float> poison, <4 x i32> zeroinitializer
%258 = insertelement <4 x float> poison, float %254, i64 0
%259 = shufflevector <4 x float> %258, <4 x float> poison, <4 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 float], [32 x float]* %8, i64 0, i64 %261
%265 = bitcast float* %264 to <4 x float>*
%266 = load <4 x float>, <4 x float>* %265
%267 = getelementptr inbounds float, float* %264, i64 4
%268 = bitcast float* %267 to <4 x float>*
%269 = load <4 x float>, <4 x float>* %268
%270 = getelementptr inbounds float, float* %0, i64 %263
%271 = bitcast float* %270 to <4 x float>*
%272 = load <4 x float>, <4 x float>* %271
%273 = getelementptr inbounds float, float* %270, i64 4
%274 = bitcast float* %273 to <4 x float>*
%275 = load <4 x float>, <4 x float>* %274
%276 = call <4 x float> @llvm.fmuladd.v4f32(<4 x float> %257, <4 x float> %266, <4 x float> %272)
%277 = call <4 x float> @llvm.fmuladd.v4f32(<4 x float> %259, <4 x float> %269, <4 x float> %275)
%278 = bitcast float* %270 to <4 x float>*
store <4 x float> %276, <4 x float>* %278
%279 = bitcast float* %273 to <4 x float>*
store <4 x float> %277, <4 x float>* %279
%280 = or i64 %261, 8
%281 = add nuw nsw i64 %280, %251
%282 = getelementptr inbounds [32 x float], [32 x float]* %8, i64 0, i64 %280
%283 = bitcast float* %282 to <4 x float>*
%284 = load <4 x float>, <4 x float>* %283
%285 = getelementptr inbounds float, float* %282, i64 4
%286 = bitcast float* %285 to <4 x float>*
%287 = load <4 x float>, <4 x float>* %286
%288 = getelementptr inbounds float, float* %0, i64 %281
%289 = bitcast float* %288 to <4 x float>*
%290 = load <4 x float>, <4 x float>* %289
%291 = getelementptr inbounds float, float* %288, i64 4
%292 = bitcast float* %291 to <4 x float>*
%293 = load <4 x float>, <4 x float>* %292
%294 = call <4 x float> @llvm.fmuladd.v4f32(<4 x float> %257, <4 x float> %284, <4 x float> %290)
%295 = call <4 x float> @llvm.fmuladd.v4f32(<4 x float> %259, <4 x float> %287, <4 x float> %293)
%296 = bitcast float* %288 to <4 x float>*
store <4 x float> %294, <4 x float>* %296
%297 = bitcast float* %291 to <4 x float>*
store <4 x float> %295, <4 x float>* %297
%298 = add nuw i64 %261, 16
%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 float], [32 x float]* %8, i64 0, i64 %302
%306 = bitcast float* %305 to <4 x float>*
%307 = load <4 x float>, <4 x float>* %306
%308 = getelementptr inbounds float, float* %305, i64 4
%309 = bitcast float* %308 to <4 x float>*
%310 = load <4 x float>, <4 x float>* %309
%311 = getelementptr inbounds float, float* %0, i64 %304
%312 = bitcast float* %311 to <4 x float>*
%313 = load <4 x float>, <4 x float>* %312
%314 = getelementptr inbounds float, float* %311, i64 4
%315 = bitcast float* %314 to <4 x float>*
%316 = load <4 x float>, <4 x float>* %315
%317 = call <4 x float> @llvm.fmuladd.v4f32(<4 x float> %257, <4 x float> %307, <4 x float> %313)
%318 = call <4 x float> @llvm.fmuladd.v4f32(<4 x float> %259, <4 x float> %310, <4 x float> %316)
%319 = bitcast float* %311 to <4 x float>*
store <4 x float> %317, <4 x float>* %319
%320 = bitcast float* %314 to <4 x float>*
store <4 x float> %318, <4 x float>* %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 float], [32 x float]* %8, i64 0, i64 %325
%328 = load float, float* %327
%329 = getelementptr inbounds float, float* %0, i64 %326
%330 = load float, float* %329
%331 = tail call float @llvm.fmuladd.f32(float %254, float %328, float %330)
store float %331, float* %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
}