void estimate_a(float *alpha, float *beta, int nstates, int length,
                float *gamma_sum, float *xi_sum, float *a) {
  float sum_ab = 0;
  float sum_a;
  int i;
  int j;
  for (i = 0; i < nstates; i++) {
    sum_ab += alpha[((length - 1) * nstates) + i] *
              beta[((length - 1) * nstates) + i];
  }
  for (i = 0; i < nstates; i++) {
    for (j = 0; j < nstates; j++) {
      a[(i * nstates) + j] =
          xi_sum[(i * nstates) + j] /
          (gamma_sum[i] - (alpha[((length - 1) * nstates) + i] *
                           beta[((length - 1) * nstates) + i] / sum_ab));
    }
    sum_a = 0;
    for (j = 0; j < nstates; j++) {
      sum_a += a[(i * nstates) + j];
    }
    for (j = 0; j < nstates; j++) {
      a[(i * nstates) + j] = a[(i * nstates) + j] / sum_a;
    }
  }
}


x86 assembly (clang, -O0):
<estimate_a>:
push   %rbp
mov    %rsp,%rbp
mov    0x10(%rbp),%rax
mov    %rdi,-0x8(%rbp)
mov    %rsi,-0x10(%rbp)
mov    %edx,-0x14(%rbp)
mov    %ecx,-0x18(%rbp)
mov    %r8,-0x20(%rbp)
mov    %r9,-0x28(%rbp)
xorps  %xmm0,%xmm0
movss  %xmm0,-0x2c(%rbp)
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    118b <estimate_a+0x8b>
mov    -0x8(%rbp),%rax
mov    -0x18(%rbp),%ecx
sub    $0x1,%ecx
imul   -0x14(%rbp),%ecx
add    -0x34(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x10(%rbp),%rax
mov    -0x18(%rbp),%ecx
sub    $0x1,%ecx
imul   -0x14(%rbp),%ecx
add    -0x34(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm2
movss  -0x2c(%rbp),%xmm1
mulss  %xmm2,%xmm0
addss  %xmm1,%xmm0
movss  %xmm0,-0x2c(%rbp)
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    112d <estimate_a+0x2d>
movl   $0x0,-0x34(%rbp)
mov    -0x34(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    12e5 <estimate_a+0x1e5>
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1237 <estimate_a+0x137>
mov    -0x28(%rbp),%rax
mov    -0x34(%rbp),%ecx
imul   -0x14(%rbp),%ecx
add    -0x38(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
mov    -0x20(%rbp),%rax
movslq -0x34(%rbp),%rcx
movss  (%rax,%rcx,4),%xmm1
mov    -0x8(%rbp),%rax
mov    -0x18(%rbp),%ecx
sub    $0x1,%ecx
imul   -0x14(%rbp),%ecx
add    -0x34(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm2
mov    -0x10(%rbp),%rax
mov    -0x18(%rbp),%ecx
sub    $0x1,%ecx
imul   -0x14(%rbp),%ecx
add    -0x34(%rbp),%ecx
movslq %ecx,%rcx
mulss  (%rax,%rcx,4),%xmm2
divss  -0x2c(%rbp),%xmm2
subss  %xmm2,%xmm1
divss  %xmm1,%xmm0
mov    0x10(%rbp),%rax
mov    -0x34(%rbp),%ecx
imul   -0x14(%rbp),%ecx
add    -0x38(%rbp),%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    11a5 <estimate_a+0xa5>
xorps  %xmm0,%xmm0
movss  %xmm0,-0x30(%rbp)
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    1280 <estimate_a+0x180>
mov    0x10(%rbp),%rax
mov    -0x34(%rbp),%ecx
imul   -0x14(%rbp),%ecx
add    -0x38(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
addss  -0x30(%rbp),%xmm0
movss  %xmm0,-0x30(%rbp)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    1246 <estimate_a+0x146>
movl   $0x0,-0x38(%rbp)
mov    -0x38(%rbp),%eax
cmp    -0x14(%rbp),%eax
jge    12d2 <estimate_a+0x1d2>
mov    0x10(%rbp),%rax
mov    -0x34(%rbp),%ecx
imul   -0x14(%rbp),%ecx
add    -0x38(%rbp),%ecx
movslq %ecx,%rcx
movss  (%rax,%rcx,4),%xmm0
divss  -0x30(%rbp),%xmm0
mov    0x10(%rbp),%rax
mov    -0x34(%rbp),%ecx
imul   -0x14(%rbp),%ecx
add    -0x38(%rbp),%ecx
movslq %ecx,%rcx
movss  %xmm0,(%rax,%rcx,4)
mov    -0x38(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x38(%rbp)
jmp    1287 <estimate_a+0x187>
jmp    12d7 <estimate_a+0x1d7>
mov    -0x34(%rbp),%eax
add    $0x1,%eax
mov    %eax,-0x34(%rbp)
jmp    1192 <estimate_a+0x92>
pop    %rbp
ret

x86 assembly (clang, -O2):
<estimate_a>:
push   %rbp
push   %r15
push   %r14
push   %r13
push   %r12
push   %rbx
mov    %r9,-0x40(%rsp)
test   %edx,%edx
jle    1509 <estimate_a+0x409>
add    $0xffffffff,%ecx
imul   %edx,%ecx
movslq %ecx,%r14
mov    %edx,%ecx
lea    -0x1(%rcx),%r13
mov    %ecx,%r9d
and    $0x3,%r9d
cmp    $0x3,%r13
jae    113a <estimate_a+0x3a>
xorps  %xmm0,%xmm0
xor    %ebp,%ebp
jmp    11a7 <estimate_a+0xa7>
mov    %ecx,%r11d
and    $0xfffffffc,%r11d
lea    (%rsi,%r14,4),%rax
add    $0xc,%rax
lea    (%rdi,%r14,4),%rbx
add    $0xc,%rbx
xorps  %xmm0,%xmm0
xor    %ebp,%ebp
cs nopw 0x0(%rax,%rax,1)
movss  -0xc(%rbx,%rbp,4),%xmm1
movss  -0x8(%rbx,%rbp,4),%xmm2
mulss  -0xc(%rax,%rbp,4),%xmm1
mulss  -0x8(%rax,%rbp,4),%xmm2
addss  %xmm0,%xmm1
movss  -0x4(%rbx,%rbp,4),%xmm3
mulss  -0x4(%rax,%rbp,4),%xmm3
addss  %xmm1,%xmm2
movss  (%rbx,%rbp,4),%xmm0
mulss  (%rax,%rbp,4),%xmm0
addss  %xmm2,%xmm3
addss  %xmm3,%xmm0
add    $0x4,%rbp
cmp    %rbp,%r11
jne    1160 <estimate_a+0x60>
test   %r9,%r9
je     11d8 <estimate_a+0xd8>
add    %r14,%rbp
lea    (%rsi,%rbp,4),%rax
lea    (%rdi,%rbp,4),%rbp
xor    %ebx,%ebx
nopl   0x0(%rax)
movss  0x0(%rbp,%rbx,4),%xmm1
mulss  (%rax,%rbx,4),%xmm1
addss  %xmm1,%xmm0
add    $0x1,%rbx
cmp    %rbx,%r9
jne    11c0 <estimate_a+0xc0>
test   %edx,%edx
jle    1509 <estimate_a+0x409>
mov    0x38(%rsp),%r11
lea    0x1(%r14),%rax
mov    %rax,-0x18(%rsp)
mov    %ecx,%r10d
and    $0xfffffffc,%r10d
lea    -0x4(%r10),%rax
mov    %rax,-0x20(%rsp)
mov    %rax,%rbx
shr    $0x2,%rbx
add    $0x1,%rbx
mov    %ecx,%eax
and    $0x3,%eax
mov    %rbx,%rbp
mov    %rbx,-0x28(%rsp)
and    $0xfffffffffffffffe,%rbx
mov    %rbx,-0x38(%rsp)
lea    0x0(,%rcx,4),%rbp
mov    %rbp,-0x8(%rsp)
and    $0xfffffffffffffff0,%rbp
mov    %rbp,-0x30(%rsp)
xor    %r15d,%r15d
mov    -0x40(%rsp),%r9
mov    %edx,-0x64(%rsp)
mov    %r14,-0x60(%rsp)
mov    %r13,-0x48(%rsp)
jmp    1268 <estimate_a+0x168>
nopl   (%rax)
add    $0x1,%r15
mov    -0x8(%rsp),%rbx
add    %rbx,%r11
add    %rbx,%r9
cmp    %rcx,%r15
je     1509 <estimate_a+0x409>
lea    (%r15,%r14,1),%r12
cmp    $0x4,%edx
jae    1280 <estimate_a+0x180>
xor    %ebx,%ebx
jmp    13d0 <estimate_a+0x2d0>
nopl   0x0(%rax,%rax,1)
mov    %r15,%rdx
imul   %rcx,%rdx
lea    (%rdx,%rcx,1),%rbx
mov    %rbx,-0x58(%rsp)
mov    -0x18(%rsp),%rbp
lea    (%r15,%rbp,1),%rbx
mov    %rbx,-0x50(%rsp)
mov    0x38(%rsp),%r13
lea    0x0(%r13,%rdx,4),%r14
mov    -0x58(%rsp),%rbp
lea    0x0(%r13,%rbp,4),%r13
lea    (%rdi,%rbx,4),%rbx
mov    %rbx,-0x10(%rsp)
mov    -0x40(%rsp),%rbp
mov    -0x58(%rsp),%rbx
lea    0x0(,%rbx,4),%rbx
add    %rbp,%rbx
cmp    %rbx,%r14
lea    0x4(%r8,%r15,4),%rbx
lea    0x0(%rbp,%rdx,4),%rdx
setb   -0x58(%rsp)
cmp    %r13,%rdx
setb   -0x65(%rsp)
cmp    %rbx,%r14
lea    (%r8,%r15,4),%rdx
setb   %bl
cmp    %r13,%rdx
setb   -0x66(%rsp)
cmp    -0x10(%rsp),%r14
lea    (%rdi,%r12,4),%rdx
setb   %bpl
cmp    %r13,%rdx
mov    -0x50(%rsp),%rdx
lea    (%rsi,%rdx,4),%rdx
setb   -0x50(%rsp)
cmp    %rdx,%r14
lea    (%rsi,%r12,4),%rdx
setb   %r14b
cmp    %r13,%rdx
setb   %r13b
mov    -0x65(%rsp),%dl
test   %dl,-0x58(%rsp)
jne    13bd <estimate_a+0x2bd>
mov    %r13d,%edx
and    -0x66(%rsp),%bl
mov    -0x48(%rsp),%r13
jne    13b0 <estimate_a+0x2b0>
and    -0x50(%rsp),%bpl
jne    13b0 <estimate_a+0x2b0>
mov    $0x0,%ebx
and    %dl,%r14b
mov    -0x64(%rsp),%edx
mov    -0x60(%rsp),%r14
mov    -0x30(%rsp),%rbp
jne    13d0 <estimate_a+0x2d0>
movss  (%r8,%r15,4),%xmm1
movss  (%rdi,%r12,4),%xmm2
mulss  (%rsi,%r12,4),%xmm2
divss  %xmm0,%xmm2
subss  %xmm2,%xmm1
shufps $0x0,%xmm1,%xmm1
xor    %ebx,%ebx
nopw   0x0(%rax,%rax,1)
movups (%r9,%rbx,1),%xmm2
divps  %xmm1,%xmm2
movups %xmm2,(%r11,%rbx,1)
add    $0x10,%rbx
cmp    %rbx,%rbp
jne    1390 <estimate_a+0x290>
mov    %r10,%rbx
cmp    %rcx,%r10
jne    13d0 <estimate_a+0x2d0>
jmp    1403 <estimate_a+0x303>
xor    %ebx,%ebx
mov    -0x64(%rsp),%edx
mov    -0x60(%rsp),%r14
jmp    13d0 <estimate_a+0x2d0>
xor    %ebx,%ebx
mov    -0x64(%rsp),%edx
mov    -0x60(%rsp),%r14
mov    -0x48(%rsp),%r13
nopl   (%rax)
movss  (%r9,%rbx,4),%xmm1
movss  (%r8,%r15,4),%xmm2
movss  (%rdi,%r12,4),%xmm3
mulss  (%rsi,%r12,4),%xmm3
divss  %xmm0,%xmm3
subss  %xmm3,%xmm2
divss  %xmm2,%xmm1
movss  %xmm1,(%r11,%rbx,4)
add    $0x1,%rbx
cmp    %rbx,%rcx
jne    13d0 <estimate_a+0x2d0>
xorps  %xmm1,%xmm1
xor    %ebx,%ebx
cmp    $0x3,%r13
jb     1434 <estimate_a+0x334>
xchg   %ax,%ax
addss  (%r11,%rbx,4),%xmm1
addss  0x4(%r11,%rbx,4),%xmm1
addss  0x8(%r11,%rbx,4),%xmm1
addss  0xc(%r11,%rbx,4),%xmm1
add    $0x4,%rbx
cmp    %rbx,%r10
jne    1410 <estimate_a+0x310>
test   %rax,%rax
je     144e <estimate_a+0x34e>
lea    (%r11,%rbx,4),%rbx
xor    %ebp,%ebp
nop
addss  (%rbx,%rbp,4),%xmm1
add    $0x1,%rbp
cmp    %rbp,%rax
jne    1440 <estimate_a+0x340>
cmp    $0x4,%edx
jae    1460 <estimate_a+0x360>
xor    %ebx,%ebx
jmp    14e0 <estimate_a+0x3e0>
nopw   0x0(%rax,%rax,1)
movaps %xmm1,%xmm2
shufps $0x0,%xmm1,%xmm2
cmpq   $0x0,-0x20(%rsp)
je     14fe <estimate_a+0x3fe>
mov    -0x38(%rsp),%rbp
xor    %ebx,%ebx
nopw   0x0(%rax,%rax,1)
movups (%r11,%rbx,4),%xmm3
movups 0x10(%r11,%rbx,4),%xmm4
divps  %xmm2,%xmm3
movups %xmm3,(%r11,%rbx,4)
divps  %xmm2,%xmm4
movups %xmm4,0x10(%r11,%rbx,4)
add    $0x8,%rbx
add    $0xfffffffffffffffe,%rbp
jne    1480 <estimate_a+0x380>
testb  $0x1,-0x28(%rsp)
je     14c9 <estimate_a+0x3c9>
mov    %r15,%rbp
imul   %rcx,%rbp
add    %rbp,%rbx
mov    0x38(%rsp),%rbp
movups 0x0(%rbp,%rbx,4),%xmm3
divps  %xmm2,%xmm3
movups %xmm3,0x0(%rbp,%rbx,4)
mov    %r10,%rbx
cmp    %rcx,%r10
je     1250 <estimate_a+0x150>
cs nopw 0x0(%rax,%rax,1)
nop
movss  (%r11,%rbx,4),%xmm2
divss  %xmm1,%xmm2
movss  %xmm2,(%r11,%rbx,4)
add    $0x1,%rbx
cmp    %rbx,%rcx
jne    14e0 <estimate_a+0x3e0>
jmp    1250 <estimate_a+0x150>
xor    %ebx,%ebx
testb  $0x1,-0x28(%rsp)
jne    14ad <estimate_a+0x3ad>
jmp    14c9 <estimate_a+0x3c9>
pop    %rbx
pop    %r12
pop    %r13
pop    %r14
pop    %r15
pop    %rbp
ret

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

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